!=======================================================================
 $MODEL_GRIDS

   EXPNME   = '{EXPNME}',

   RUNTYPE  = '{RUNTYPE}',    ! Type of run: MAKESFC, INITIAL, HISTORY,
                            !              MAKEVFILE, or MEMORY

   TIMEUNIT = 'h',          ! 'h','m','s' - Time units of TIMMAX, TIMSTR

   TIMMAX   = {TIMMAX},           ! Final time of simulation

!  Start of simulation or ISAN processing

   IMONTH1  ={IMONTH1},      ! Month
   IDATE1   ={IDATE1},     ! Day
   IYEAR1   ={IYEAR1},      ! Year
   ITIME1   ={ITIME1},         ! GMT of model TIME = 0.

! Grid specifications

   NGRIDS   = 1,            ! Number of grids to run
 
   NNXP     = {NNXP},         ! Number of x gridpoints
   NNYP     = {NNYP},      ! Number of y gridpoints
   NNZP     = {NNZP}, 47,         ! Number of z gridpoints
   NZG      = 7,	      ! Number of soil layers
   NZS      = 1,	        ! Maximum number of snow layers
 
   NXTNEST  = 0,1,2,3,          ! Grid number which is the next coarser grid
   
! Coarse grid specifications

   IF_ADAP  = 0,

   IHTRAN   = 1,	        ! 0-Cartesian, 1-Polar stereo
   DELTAX   = {DELTAX},
   DELTAY   = {DELTAY},       ! X and Y grid spacing

   DELTAZ   = 80.,         ! Z grid spacing (set to 0. to use ZZ)
   DZRAT    = 1.08,       ! Vertical grid stretch ratio
   DZMAX    = 850.,        ! Maximum delta Z for vertical stretch

   ZZ       = 0.0,          ! Vertical levels if DELTAZ = 0
     20.0,      46.0,      80.0,     120.0,     165.0,
    220.0,     290.0,     380.0,     480.0,     590.0,
    720.0,     870.0,    1030.0,    1200.0,    1380.0,
   1595.0,    1850.0,    2120.0,    2410.0,    2715.0,
   3030.0,    3400.0,    3840.0,    4380.0,    5020.0,
   5800.0,    6730.0,    7700.0,    8700.0,    9700.0,
  10700., 11700., 12700., 13700., 14700., 15700., 16700.,
  17700., 18700., 19700.,
 
   DTLONG   = {DTLONG}, 	   ! Coarse grid long timestep
   NACOUST  = 4,           ! Small timestep ratio
   IDELTAT  = 0,            ! =0 - constant timesteps
                            ! >0 - initial computation <0 - variable

                            ! Nest ratios between this grid
                            !  and the next coarser grid.
   NSTRATX =   1, 3, 2, 4, 4,		  ! x-direction
   NSTRATY =   1, 3, 2, 4, 4,		  ! y-direction
   NNDTRAT =   1, 4, 2, 4, 4,		  ! time

   NESTZ1   = 0,                          ! Contort coarser grids if negative
   NSTRATZ1 = 2,2,2,1,            !
   NESTZ2   = 0,                          ! Contort coarser grids if negative
   NSTRATZ2 = 3,3,2,1,            !


   POLELAT  = {POLELAT},         ! Latitude of pole point
   POLELON  = {POLELON},         ! Longitude of pole point

   CENTLAT  = {CENTLAT}, -23.0, -23.5, -20.9,
   CENTLON  = {CENTLON}, -46.65, -46.5, -47.5,

 
                            ! Grid point on the next coarser
                            !  nest where the lower southwest
                            !  corner of this nest will start.
                            !  If NINEST or NJNEST = 0, use CENTLAT/LON
   NINEST   = 1, 0, 0, 0,         ! i-point
   NJNEST   = 1, 0, 0, 0,         ! j-point
   NKNEST   = 1, 1, 1, 1,         ! k-point
 
   NNSTTOP  = 1, 1, 1, 1,	   ! Flag (0-no or 1-yes) if this
   NNSTBOT  = 1, 1, 1, 1,	   ! Nest goes the top or bottom of the
                                        !  coarsest nest.
 $END
 
 $CCATT_INFO

  CCATT = 1, ! 1-CCATT environmental model activated 0-off

   CHEMISTRY = 4,  ! -1- only atmospheric model
                   !  0- only tracer transport with atmos. model
		   !  1- activate the loss/production by kinetic chemical
                   !    and photochemical reactions, with solver 1
		   !  2- activate the loss/production by kinetic chemical
                   !     and photochemical reactions, with solver 2
		   !  3- Rosenbrock 2nd order - dynamic timestep
		   !  4- Rosenbrock "RODAS3" 3rd order - dynamic timestep

   SPLIT_METHOD = 'SYMMETRIC',! Splitting operator method
                              !'SYMMETRIC' = sequential symmetric (D->Q->D)
			      !'SEQUENTIAL'= sequential at right (D->Q )
			      !'PARALLEL' = original method (D->)
                              !                             (Q->)
   CHEM_TIMESTEP = {CHEM_TIMESTEP}, ! chemistry timestep integration (seconds)
                         ! should be a multiple of dtlong
                         
   CHEMISTRY_AQ = 0,  !  4dda for chemistry species with mocage data
                      !  1 - ON, 0 - OFF aqueous chemistry

   CHEM_ASSIM = {CHEM_ASSIM}, ! 1- activate 4dda with mocage data, 0 - off
   AER_ASSIM = {AER_ASSIM}, ! 1 - Activate 4dda aer assim   
   
   SRCMAPFN = {SRCMAPFN}, ! Source Map file prefix ou "NONE" !caminho prepchem

   RECYCLE_TRACERS = {RECYCLE_TRACERS}, !1,

   DEF_PROC_SRC = 'LAST_SOURCES', ! what to do if src files are not available 
                        ! (STOP, LAST_SOURCES)
   DIUR_CYCLE = 1,1,1,1,

   NA_EXTRA2D = 0,   !Number of extras 2d arrays
   NA_EXTRA3D = 0,   !Number of extras 3d arrays

   PLUMERISE  =1, !0-Desativated, 1-Activated Plume Rise process
   
   PRFRQ =7200.,

   VOLCANOES = 0, !0 -Desactivated, 1-Activated volcanoes emissiosn/processes

   AEROSOL = 1,    ! 1 -  aerosol model: on
                   ! 0 -  aerosol model: off
                   
   AER_TIMESTEP = {DTLONG}, ! aerosol time step integration (sec)
                        ! should be an multiple of dtlong (max 4)		
                ! 0 - off, 1 - on 

 $END

 $TEB_SPM_INFO

  TEB_SPM  = 0, ! 1-TEB activated 0-off

 $END

 $MODEL_FILE_INFO

! Variable initialization input

   INITIAL  = 2,             ! Initial fields - 1=horiz.homogeneous, 

!  ------ Analysis nudging parameters -----------------------------------
   NUD_TYPE = 2,             ! =1 - nudge from history files(1-way nest)
                             ! =2 - nudge from varfiles
                             ! =0 - no analysis nudging
                             !                  2=variable
!
   VARFPFX  = '/lustre_xc50/luiz_flavio/oper/{DATE}/dataout/IVAR/OPQUE',      ! Varfile initialization file prefix
   VWAIT1   = 0.,            ! wait between each VFILE check (s)
   VWAITTOT = 0.,            ! total wait befor giving up on a VFILE (s)

   NUD_HFILE = './H/a-H-2001-07-21-000000-head.txt', 
            ! Header file name for history nudging files (only prefix is used)


   NUDLAT   = 20,             ! Number of points in lateral bnd region
   TNUDLAT  = 3600.,          ! Nudging time scale(s) at lateral boundary
   TNUDCENT = 0.,        ! Nudging time scale(s) in center of domain
   TNUDTOP  = 3600.,        ! Nudging time scale (s) at top of domain
   ZNUDTOP  = 15000.,        ! Nudging at top of domain above height(m)

   WT_NUDGE_GRID =  1., 0.7, 0.7, 0.5, ! Relative nudging weights for active grids
                                       ! =0., turns off nudging for that grid

                             ! These weights will multiply the base timescales
                             !   to determine full nudging weight. 
                             !   (Timescales)/(WT_NUDGE_*) 
                             !    must be larger than DTLONG
   WT_NUDGE_UV = 1.,         ! Anal nudging weight for u and v
   WT_NUDGE_TH = 1.,         ! Anal nudging weight for theta
   WT_NUDGE_PI = 1.,         ! Anal nudging weight for pi
   WT_NUDGE_RT = 1.,         ! Anal nudging weight for r_tot
   
!-------------------------------------------------------------------------
! History file input

   TIMSTR   = 0.,           ! Time of history start (see TIMEUNIT)
   HFILIN   = '/scratchout/oper/tempo/BRAMS/ams_20km/dataout/2019020800/HIS/hOPQUE-H-2003-06-15-000000.vfm',
                             ! Input history file name

! Analysis file input for assimilation (currently LEAF variables)

   IPASTIN   = 0,            ! Initialize various fields from analysis file?
                             !  1=yes, 0=no
   PASTFN    = {PASTFN}, ! usar a linha abaixo para o primeiro dia somente

                             ! Input analysis file name

 
! History/analysis file output
 
   IPOS     =  2,                ! 0-no files, 2-grads files, 1-hdf5(not available yet) 
   IOUTPUT  =  0,                ! 0-no files, 1-save ASCII, 2-save binary
   HFILOUT  = {HFILOUT},      ! History file prefix
   AFILOUT  = {AFILOUT},       ! Analysis file prefix
   ICLOBBER = 1,                 ! 0=stop if files exist, 1=overwite files
   IHISTDEL = 1,                 ! 0=keep all hist files, 1=delete previous
   FRQHIS   = 2160000.,           ! History file frequency
   FRQANL   = {FRQANL},            ! Analysis file frequency


!-----------------------------------------------------------------------
   KWRITE   = 0,             ! 1-write,0-don't write scalar K's to anal.

! Printed output controls
 
   FRQPRT   = 21600.,        ! Printout frequency
   INITFLD  = 0,             ! Initial field print flag 0=no prnt,1=prnt
 
! Input topography variables

   TOPFILES = './{DATE}/dataout/SFC/top_OQ3g',  ! File path and prefix for topo files.
   SFCFILES = './{DATE}/dataout/SFC/sfc_OQ3g',   ! File path and prefix for surface
   SSTFPFX  = './{DATE}/dataout/SFC/sst_OQ3g',   ! Path and prefix for sst files
   NDVIFPFX = './{DATE}/dataout/SFC/ndv_OQ3g',  ! Path and prefix for ndvi files

   ITOPTFLG = 1,1,1,1,0,	    ! 2 - Fill data in "rsurf"
   ISSTFLG  = 1,1,1,1,0,	    ! 0 - Interpolate from coarser grid
   IVEGTFLG = 1,1,1,1,0,	    ! 1 - Read from standard Lat/Lon data file
   ISOILFLG = 1,1,1,1,0,            ! soil files not yet available: avoid isoilflg=1
   NDVIFLG  = 1,1,1,1,0,          !
                                    !avoid isoilflg=1
   
   NOFILFLG = 2,2,2,2,       ! 2 - Fill data in "rsurf"
                             ! 0 - Interpolate from coarser grid   

   IUPDNDVI  = 1,             ! 0 - No update of NDVI values during run
   IUPDSST  = 1,             ! 0 - No update of SST values during run
                             ! 1 - Update SST values during run
 
                             ! The following only apply for IxxxxFLG=1


   ITOPTFN  = './datafix/TOPO1KM/EL',
                                      
   ISSTFN   = './datafix/SST_WEEK/W',

   IVEGTFN  = './datafix/GL_OGE_INPE/OGE',
             
   ISOILFN  = './datafix/GL_FAO_INPE/FAO',

   NDVIFN  = './datafix/NDVI-MODIS_GRADS/N',	     
             
! Topography scheme

   ITOPSFLG = 0,0,0,1,       ! 0 = Average Orography
                             ! 1 = Silhouette Orography
                             ! 2 = Envelope Orography
                             ! 3 = Reflected Envelope Orography
   TOPTENH  = 0.,0.,0.,0.,   ! For ITOPSFLG=1, Weighting of topo
                             !   silhouette averaging 
                             ! For ITOPSFLG=2 or 3, Reflected Envelope
                             !   and Envelope Orography enhancement
                             !   factor
   TOPTWVL  = -4.,-2.,2.,2.,   ! Topo wavelength cutoff in filter 
                                
! Surface Roughness scheme

   IZ0FLG   = 0,0,0,         ! 0 = Based off vege, bare soil and water surface
                             ! 1 = Based off subgrid scale topography
   Z0MAX    = 5.,5.,5.,      ! Max zo for IZ0FLG=1
   Z0FACT   = 0.005,         ! Subgrid scale orograhic roughness factor
   
! Microphysics collection tables

   MKCOLTAB = 0,             ! Make table: 0 = no, 1 = yes
   COLTABFN = './tables/micro/ct2.0',
                             ! Filename to read or write
   MAPAOTFILE = './tables/rad_carma/infMapAOT.vfm',
   JULESIN = './jules.in',

 $END
 
 $MODEL_OPTIONS

! Numerical schemes
   
   !** this flag is experimental, must be used with caution
   !** only flags 0 and 1 are suitable for parallel runs.
   DYNCORE_FLAG = 2, !0- leapfrog Robert-Asselin filter, 
                     !1- leapfrog Robert-Asselin_Williams filter
		     !2- Runge-Kutta 3rd order
		     !3- Adams-Bashforth-Moulton 3rd order
  
   PD_OR_MNT_CONSTRAINT = 1,
 
   ADVMNT = 1,  ! 0 = Forward 2nd order (non-monotonic)
                ! 1 = Walcek monotonic advection
                ! 2 = hybrid: forward 2nd order for thetail,microphysics, TKE and
		!             Walcek scheme for tracers 

   GHOSTZONELENGTH = 3,     ! Size of ghostzone, only applicable for 
                            ! ADVMNT > 0
 
   ICORFLG  = 1,             ! Coriolis flag/2D v-component  -  0=off, 1=on
   VVELDAMP = 0,
   IEXEV    = 2,             ! Exner tendency
                             ! 1- simplest form
                             ! 2- complete, mass conservative, solution 
			     !    (works only with sigma-z coordinate)
   
   IMASSFLX = 0,             ! Output advective and convective mass flux for
                             ! STILT Lagrangian Particle Dispersion Modelling
			     ! 0 = off, 1 = on

   IBND     = 1,             ! Lateral boundary condition flags
   JBND     = 1,             ! 1-Klemp/Wilhelmson, 2-Klemp/Lilly, 3-Orlanski
                             !  4-cyclic
   CPHAS    = 20.,           ! Phase speed if IBND or JBND = 1
   LSFLG    = 0,             ! Large-scale gradient flag for variables other than
                             !  normal velocity:
                             !  0 = zero gradient inflow and outflow
                             !  1 = zero gradient inflow, radiative b.c. outflow
                             !  2 = constant inflow, radiative b.c. outflow
                             !  3 = constant inflow and outflow
   NFPT     = 0,             ! Rayleigh friction - number of points from the top
   DISTIM   = 400.,          !	          - dissipation time scale
 
! Radiation parameters
 
!AAF-20160706   ISWRTYP  = 6,             ! Shortwave radiation type
!AAF-20160706   ILWRTYP  = 6,             ! Longwave radiation type
   ISWRTYP  = 6,             ! Shortwave radiation type
   ILWRTYP  = 6,             ! Longwave radiation type
                             !    0-none, 2-Mahrer/Pielke, 1-Chen, 3-Harrington
			     !    4- CARMA 
! Parameters needed by CARMA radiation scheme
   RADDATFN  = './tables/rad_carma/rad_param.data',

!   RADFRQ   = 1200.,          ! Freq. of radiation tendency update (s)
   RADFRQ   = 3600.,          ! Freq. of radiation tendency update (s)
   LONRAD   = 1,	      ! Longitudinal variation of shortwave
                              !    (0-no, 1-yes)
! Cumulus parameterization parameters
 
   NNQPARM  = 6,! Convective param. flag 
                       ! 0- off,
                       ! 1- Tremback formulation
		       ! 2- Grell-Deveny scheme
		       ! 3- Grell-3d formulation
                       ! 4- Grell-Deveny scheme as implemented in FIM/NOAA model                      
		       ! 5- Grell-Freitas scheme
   
   CLOSURE_TYPE = 'PB', ! Closure type (for Grell Param.):
                        ! EN: ensemble (all closures)
                        ! GR: Grell
                        ! LO: low level omega
                        ! MC: moisture convergence
                        ! SC: like Fritsch Chappel or Kain Fritsch
                        ! AS: Arakawa-Schubert

   G3D_SPREAD	= 0, ! Lateral spreading of enviromental subsidence (0=off, 1=on)
                     ! (only for NNQPARM  = 3 and grid spacing < 10km)

   NNSHCU   = 2,   ! Shallow Cumulus Param. (0-off, 1-Souza scheme, 2-Grell-Deveny scheme)
   CONFRQ   = 600.,      ! Frequency of conv param. updates (s)
   SHCUFRQ  = 600.,      ! Frequency of Shallow param. updates (s)
   WCLDBS   = .0005,     ! Vertical motion needed at cloud base for
                         !    to trigger convection
 
! Surface layer and soil parameterization
 
   NPATCH   = 4,             ! Number of patches per grid cell (min=2)
                       
   NVEGPAT  = 3,             ! Number of patches per grid cell to be 
                             !    filled from
                             ! vegetation files 
                             !    (min of 1, max of NPATCH-1)

   ISFCL    = 5,             ! Surface layer/soil/veg model
                             !  0-specified surface layer gradients
                             !  1-soil/vegetation model  LEAF
                             !  2-soil/vegetation model LEAF & hidrological model
                             !  3-Reserved
                             !  4-Reserved
                             !  5-JULES

   NVGCON   = 7,             ! Vegetation type (see below)
   
   !   1 -- Crop/mixed farming         2 -- Short grass  
   !   3 -- Evergreen needleleaf tree  4 -- Deciduous needleleaf tree  
   !   5 -- Deciduous broadleaf tree   6 -- Evergreen broadleaf tree
   !   7 -- Tall grass                 8 -- Desert
   !   9 -- Tundra                    10 -- Irrigated crop
   !  11 -- Semi-desert               12 -- Ice cap/glacier
   !  13 -- Bog or marsh              14 -- Inland water
   !  15 -- Ocean                     16 -- Evergreen shrub
   !  17 -- Deciduous shrub           18 -- Mixed woodland

   PCTLCON  = 1.,            ! Constant land % if for all domain
   NSLCON   = 6,             ! Constant soil type if for all domain

   !   1 -- sand             2 -- loamy sand      3 -- sandy loam
   !   4 -- silt loam        5 -- loam            6 -- sandy clay loam
   !   7 -- silty clay loam  8 -- clay loam       9 -- sandy clay
   !  10 -- silty clay      11 -- clay           12 -- peat

   ZROUGH   = .05,           ! Constant roughness if for all domain
   ALBEDO   = .2,            ! Constant albedo if not running soil model
   SEATMP   = 295.,          ! Constant water surface temperature

   DTHCON   = 0.,            ! Constant sfc layer  temp grad for no soil
   DRTCON   = 0.,            ! Constant sfc layer moist grad for no soil
   
   SOIL_MOIST = 'n', ! n => Homogeneous (standard)
                     ! i,h,a => Heterogenous Soil Moist.Init.
                     !          with file.
                     !          i => INITIAL
                     !          h => HISTORY
                     !          a => INITIAL or HISTORY

   SOIL_MOIST_FAIL = 'l', ! In case of Heterogenous Soil Moist.
                          ! file was not found, what to do?
                          ! s => STOP the program
                          ! h => Initialize with Homogenous Soil Moist.
                          ! l => Looking for 5 days old files,
                          ! and if not found again then stop.

   USDATA_IN  = '/lustre_xc50/ioper/data/external/{DATE}00/dataout/Umid_Solo/GL_SM_GPNR.', !gl_sm.gpnr.',  ! Soil Moisture File prefix
   USMODEL_IN = './{DATE}/dataout/UMID/GL_SM_GPNR.', !gl_sm.gpnr.', ! File prefix with pre-calculated data by BRAMS
   
   SLZ     =  -12.25, -7.25, -4.250, -2.25, -1.0, -0.35, -0.10,      ! soil grid levels

   SLMSTR  =   0.45,0.45,0.30,0.25,0.22, 0.22,0.22, ! initial soil moisture
                          
   STGOFF =    .0,  .0,  .0,  .0,   .0, .0, .0,  ! Initial soil temperature offset
                                                      ! from lowest atmospheric level
!--------------------------------------------------------------------------
 
!--------- Eddy diffusion coefficient parameters -----------------------
 
   IDIFFK  = 1,2,2,2,        ! K flag:
                             !  1 - Horiz deform/Vert Mellor-Yamada
                             !  2 - Anisotropic deformormation
                             !      (horiz & vert differ)
                             !  3 - Isotropic deformation 
                             !      (horiz and vert same)
                             !  4 - Deardorff TKE (horiz and vert same)
			     !  7 - Nakanishi TKE 2.5-level model
			     
   IHORGRAD = 1,             !  1 - horiz grad frm decomposed sigma grad
                             !  2 - true horizontal gradient. 
                             !      Non-conserving, but allows small DZ
   CSX     = .32,.32,        ! Deformation horiz. K's coefficient '
   CSZ     = .35,.35,        ! Deformation vert. K's coefficient '
   XKHKM   =  3.,3.,          ! Ratio of horiz K_h to K_m for deformation
   ZKHKM   =  3.,3.,          ! Ratio of vert K_h to K_m for deformation
   AKMIN   =  1.,1.,          ! Ratio of minimum horizontal eddy
                             ! viscosity coefficientto typical value
                             ! from deformation K
!------------------------------------------------------------------------------
 ! Microphysics

   MCPHYS_TYPE = 2,          ! microphysics formulation
                             ! 0 = original brams
			                    ! 1 = 2M rams6.0
   LEVEL   = 3,              ! Moisture complexity level
   IRIME   = 1,              !   0 - Bulk riming, 1 - bin riming
   IPLAWS  = 0,              !   0 - Original 4.3 Power Laws
                             !   1 - New Power Laws (R.Carver,Mitchell 1996)
 
   ICLOUD  = 5,              ! Microphysics flags
   IDRIZ   = 0,              !  
   IRAIN   = 2,              !-------------------
   IPRIS   = 5,              !  1 - diagnostic concen.
   ISNOW   = 2,              !  2 - specified mean diameter
   IAGGR   = 2,              !  3 - specified y-intercept
   IGRAUP  = 2,              !  4 - specified concentration
   IHAIL   = 2,              !  5 - prognostic concentration
                             
   CPARM   = .3e9,           !  Microphysics parameters
   RPARM   = 1.e-3,          !-------------------------
   PPARM   = 1.e-6,             !  Characteristic diameter, # concentration
   SPARM   = 1.e-3,          !    or y-intercept
   APARM   = 1.e-3,
   GPARM   = 1.e-3,
   HPARM   = 3.e-3,
   DPARM   = 1.e-5,          ! (5 = 3D homogeneous 6 = vertical profile)

   GNU     = 2.,2.,2.,2.,2.,2.,2.,2., ! Gamma shape parms for


!-----------------------------------------------------------------------------

  $END
 
  $MODEL_SOUND

!-----------------------------------
! Sounding specification
!-----------------------------------
 
! Flags for how sounding is specified
 
   IPSFLG   = 1,             ! Specifies what is in PS array
                             !  0-pressure(mb) 1-heights(m) 
                             !  PS(1)=sfc press(mb)
 
   ITSFLG   = 0,             ! Specifies what is in TS array
                             !  0-temp(C) 1-temp(K) 2-pot. temp(K)

   IRTSFLG  = 3,             ! Specifies what is in RTS array
                             !  0-dew pnt.(C) 1-dew pnt.(K) 
                             !  2-mix rat(g/kg) 
                             !  3-relative humidity in %, 
                             !  4-dew pnt depression(K)
 
   IUSFLG   = 0,             ! Specifies what is in US and VS arrays
                             ! 0-u,v component(m/s) 
                             ! 1-umoms-direction, vmoms-speed
 
   HS       = 0.,

   PS       = 1010.,1000.,2000.,3000.,4000.,6000.,8000.,11000.,15000.,20000.,
              25000.,

   TS       = 25.,18.5,12.,4.5,-11.,-24.,-37.,-56.5,-56.5,-56.5,-56.5,

   RTS      = 70.,70.,70.,70.,20.,20.,20.,20.,10.,10.,10.,

   US       = 10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,10.,

   VS       = 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,
 
  $END
 
 
  $MODEL_PRINT

!-----------------------------------
! Specifies the fields to be printed during the simulation
!-----------------------------------

   NPLT     = 0,            ! Number of fields printed at each time
                            !   for various cross-sections (limit of 50)
                            
   IPLFLD   = 'UP','VP','WP','PP','THP','RT',
                            ! Field names - see table below
                            
!   PLFMT(6)  = '0PF7.2',   ! Format spec. if default is unacceptable

   IXSCTN    = 3,3,3,3,3,3,3,3,
                            ! Cross-section type (1=XZ, 2=YZ, 3=XY)

   ISBVAL    = 2,2,2,2,2,2,2,2,
                            ! Grid-point slab value for third direction
 
  ! The following variables can also be set in the  namelist:  IAA,
  ! IAB, JOA, JOB, NAAVG, NOAVG, PLTIT, PLCONLO, PLCONHI, and PLCONIN.
 
!     'UP'    - UP(M/S)     'RC'   - RC(G/KG)    'PCPT' - TOTPRE
!     'VP'    - VP(M/S)     'RR'   - RR(G/KG)    'TKE'  - TKE
!     'WP'    - WP(CM/S)    'RP'   - RP(G/KG)    'HSCL' - HL(M)
!     'PP'   - PRS(MB)      'RA'   - RA(G/KG)    'VSCL' - VL(M)
!     'THP'  - THP(K)
!     'THETA'- THETA(K)     'RL'   - RL(G/KG)    'TG'   - TG (K)
!     'THVP' - THV(K)      'RI'   - RI(G/KG)    'SLM'  - SLM (PCT)
!     'TV'   - TV(K)        'RCOND'- RD(G/KG)    'CONPR'- CON RATE
!     'RT'   - RT(G/KG)     'CP'   - NPRIS       'CONP' - CON PCP
!     'RV'   - RV(G/KG)     'RTP'  - RT'(G/KG)   'CONH' - CON HEAT
!                                                'CONM' - CON MOIS
!     'THIL' - Theta-il (K) 'TEMP' - temperature (K)
!     'TVP'  - Tv' (K)      'THV'  - Theta-v     (K)
!     'RELHUM'-relative humidity (%)           'SPEED'- wind speed (m/s)
!     'FTHRD'- radiative flux convergence (??)
!     'MICRO'- GASPRC
!     'Z0'   - Z0 (M)       'ZI'   - ZI (M)      'ZMAT' - ZMAT (M)
!     'USTARL'-USTARL(M/S)  'USTARW'-USTARW(M/S) 'TSTARL'-TSTARL (K)
!     'TSTARW'-TSTARW(K)    'RSTARL'-RSTARL(G/G) 'RSTARW'-RSTARW(G/G)
!     'UW'   - UW  (M*M/S*S)                     'VW'   - VW (M*M/S*S)
!     'WFZ'  - WFZ (M*M/S*S)                     'TFZ'  - TFZ (K*M/S)
!     'QFZ'  - QFZ (G*M/G*S)                     'RLONG'- RLONG
!     'RSHORT'-RSHORT

  $END

  $ISAN_CONTROL
 
!-----------------------------------
! Isentropic control 
!-----------------------------------
 
   ISZSTAGE = 1,            ! Main switches for isentropic-sigz
   IVRSTAGE = 1,            !    "varfile" processing
    
   ISAN_INC = 0600,         ! ISAN processing increment (hhmm)
                            !    range controlled by TIMMAX, 
                            !    IYEAR1,...,ITIME1

   GUESS1ST = 'PRESS',      ! Type of first guess input- 'PRESS', 'RAMS'

   I1ST_FLG = 1,            ! What to do if first guess file should be
                            !  used but does not exist.
                            !  1=I know it may not be there, 
                            !        skip this data time
                            !  2=I screwed up, stop the run
                            !  3=interpolate first guess file from 
                            !    nearest surrounding times, stop if unable
                            !     (not yet available)

   IUPA_FLG = 3,            ! UPA-upper air, SFC-surface
   ISFC_FLG = 3,            ! What to do if other data files should be
                            !  uesed, but does not exist.
                            !  1 = I know it may not be there, 
                            !      skip this data time
                            !  2 = I screwed up, stop the run
                            !  3 = Try to continue processing anyway

! Input data file prefixes

   IAPR    = '/lustre_xc50/valter_oliveira/lufla/GFS/DPCHEM/dp-chem-relacs', ! Input press level dataset
   IARAWI  = '', ! Archived rawindsonde file name
   IASRFCE = '', ! Archived surface obs file name

! File names and dispose flags 

   VARPFX    = '/lustre_xc50/luiz_flavio/oper/{DATE}/dataout/IVAR/OPQUE',  ! isan file names prefix
   IOFLGISZ  = 0,           ! Isen-sigz file flag: 0 = no write, 1 = write
   IOFLGVAR  = 1,           ! Var file flag: 0 = no write, 1 = write

  $END
 
 
  $ISAN_ISENTROPIC
 
!-----------------------------------
! Isentropic and sigma-z processing 
!-----------------------------------

!-----------------------------------
! Specify isentropic levels
!-----------------------------------

   NISN     = 43,           ! Number of isentropic levels
   LEVTH    = 280,282,284,286,288,290,292,294,296,298,300,303,306,309,312,
              315,318,321,324,327,330,335,340,345,350,355,360,380,400,420,
              440,460,480,500,520,540,570,600,630,670,700,750,800,
              
!-----------------------------------
! Analyzed grid information:
!-----------------------------------

   NIGRIDS  = 1,            ! Number of RAMS grids to analyze

   TOPSIGZ  = 20000.,       ! Sigma-z coordinates to about this height

   HYBBOT   = 4000.,        ! Bottom (m) of blended sigma-z/isentropic 
                            !    layer in varfiles
   HYBTOP   = 6000.,        ! Top (m) of blended sigma-z/isentropic layr

   SFCINF   = 1000.,        ! Vert influence of sfc observation analysis

   SIGZWT   = 1.,           ! Weight for sigma-z data in varfile: 
                            !  0.= no sigz data, 
                            !  1.=full weight from surface to HYBBOT

   NFEEDVAR = 1,            ! 1=feed back nested grid varfile, 0=dont

!-----------------------------------
! Observation number limits:
!-----------------------------------

   MAXSTA   = 150,          ! maximum number of rawindsondes 
                            !    (archived + special)
   MAXSFC   = 1000,         ! maximum number of surface observations

 
   NOTSTA   = 0,            ! Number of stations to be excluded
   NOTID    = 'r76458',     ! Station ID to be excluded
                            !  Prefix with 'r' for rawindsonde, 
                            !              's' for surface
   IOBSWIN = 1800,


   STASEP   = .1,           ! Minimum sfc station separation in degrees.
                            !    Any surface obs within this distance
                            !    of another obs will be thrown out
                            !    unless it has less missing data, 
                            !    in which case the other obs will be
                            !    thrown out.ramspost_50


   IGRIDFL  = 3,            ! Grid flag=0 if no grid point, only obs
                            !           1 if all grid point data and obs 
                            !           2 if partial grid point and obs
                            !           3 if only grid data
                            !           4 all data... fast


   GRIDWT   = .01,.01,      ! Relative weight for the gridded press data
                            !   compared to the observational data in 
                            !   the objective analysis


   GOBSEP   = 5.,           ! Grid-observation separation (degrees)
   GOBRAD   = 5.,           ! Grid-obs proximity radius (degrees)

   WVLNTH   = 1200.,900.,   ! Used in S. Barnes objective analysis.
                            !   Wavelength in km to be retained to the
                            !   RESPON % from the data to the upper air 
                            !   grids.
   SWVLNTH  = 750.,300.,    ! Wavelength for surface objective analysis

   RESPON   = .90,.9,       ! Percentage of amplitude to be retained.

   !-----------------------------------
   ! Makevfile information IC = Initial Conditions
   !-----------------------------------

   ICFILETYPE=4, !0- ISAN using dprep, 1- do ISAN using NCEP GFS GRIB2,
                   !2- ISAN using netCDF

   ICPREFIX = "./{DATE}/datain/GRADS/IC",  !IC file prefix

   WIND_U_VARNAME = "UGRD",    !Name of U wind variable inside grib2 descriptions (m/s)

   WIND_V_VARNAME = "VGRD",    !Name of V wind variable inside grib2 descriptions (m/s)

   TEMPERATURE_VARNAME = "TMP", !Name of temperature variable inside grib2 descriptions (k)

   GEO_VARNAME = "HGT", ! Name of geopotential variable inside grib2 descriptions (m)

   UR_VARNAME = "RH", ! Name of humidity variable inside grib2 descriptions (%)

   INITIAL_LATITUDE = -70., !Initial latitude for domain of model (-90 to 90)

   FINAL_LATITUDE  = 29., !Final latitude for domain of model(-90 to 90)

   INITIAL_LONGITUDE = 250., !Initial longitude for domain of model (0 to 360)

   FINAL_LONGITUDE = 358., !Final longitude for domain of model (0 to 360)

   Z_MAX_LEVEL = 23, !Maximum number of levels to be processed

   !          U     V      T      Z       R
   DLIMIT = -500.,-500.,   0.,  -1000.,   0., !Inferior valid values for each var

   !          U     V      T      Z       R
   ULIMIT =  500., 500., 500., 500000., 100., !Superior valid values for each var

   !              U  V  T  Z  R
   SCALE_FACTOR = 1.,1.,1.,1.,0.01, !Factor to multiplie each var

   CCGRADSWRITE = 1, !Is to write a grads with CI?

   ICGRADSPREFIX = "./{DATE}/dataout/IC/icGrads",  !Prefix for grads file with IC


  $END

  $POST
    NVP = 55,
    VP ='topo',         
	'land',        
	'sea_press',    
	'precip',    
	'acccon',       
	'tveg',       
	't2mj',         
	'le',           
	'h',            
	'rshort',       
	'rlong',       
	'albedt',       
	'rlongup',      
	'cape',         
	'cine',         
	'zi',           
	'sst',          
	'pwt',          
	'slp_metar',    
	'td2mj',        
	'u10mj',        
	'v10mj',        
	'smoist',      
	'tempc',        
	'ue_avg',      
	've_avg',      
	'rh',          
	'geo',         
	'cloud',       
	'omeg',        
	'ndvi', 
	'lai'       ,
	'rv'        ,	       
!fim variaveis scamtec
         'sfc_press' ,
         'albedt'    ,
         'tempk'     ,
         'cucldp'    ,
         'curtdp'    ,
         'cuthdp'    ,
         'aot500'    ,
         'aot550'    ,
         'PMINT'     ,
         'CO'        ,
         'PM25'      ,  
	 'PM25wd'    ,       
         'NO'        ,               
         'NO2'       ,       
         'O3'        ,
         'CO_SRC',
         'NO_SRC', 
 'R_BBURN2' ,
 'R_BBURN3' ,
 'R_URBAN2' ,
 'R_URBAN3' ,
 'R_MARIN1' ,

    GPREFIX = './{DATE}/dataout/POST/chem',
    GPREFIX = './{DATE}/dataout/POST/chem',
    CSVFILE = './variables.csv',
    ANL2GRA = 'ONE',

    PROJ='YES',
!   MEAN_TYPE ='VMP', 
    MEAN_TYPE ='BAV', 
    LATI = -90.,-90.,    -90., -90.,
    LATF = +90.,+90.,    +90., +90., 
    LONI = -180.,-180., -180.,-180.,
    LONF =  180.,180.,   180., 180.,
    
    
    ZLEVMAX = 32,32,32,32,
    
    IPRESSLEV = 0,
    INPLEVS = 19,
    IPLEVS = 1000, 925, 900, 850, 800, 750, 700, 650, 600, 550, 500, 450, 400, 350, 300, 250, 200, 150, 100,
    ASCII_DATA ='NO',
    SITE_LAT =-10.76,
    SITE_LON =-62.36,
 $END

 $DIGITALFILTER
 	       applyDigitalFilter	= .false.,
         digitalFilterTimeWindow = 43200.,

 $END
 $METEOGRAM
     applyMeteogram = .false.,
     meteogramFreq  = 1800.,
!    meteogramMap   = '/lustre_xc50/io/models/BRAMS/ams_20km/fix/tables/meteogram/citiesBr.txt',
     meteogramMap   = '/lustre_xc50/luiz_flavio/bin/tables/meteogram/citiesBr.txt',
!    meteogramDir   = '/lustre_xc50/io/data/BRAMS/ams_20km/2020061200G/dataout/HST/HST2020061200'
     meteogramDir   = './HST/HST2020061200'
 $END

