!namelist

 $MODEL_GRIDS

 !  Simulation title (64 chars)
 
    EXPNME   = 'RAMS Simulations',
 
 !  Type of model run to perform
 
    RUNTYPE  = 'INITIAL',         
 !   RUNTYPE  = 'MAKESFC',                       
 !   RUNTYPE  = 'MAKEVFILE',
 !   RUNTYPE  = 'MAKEHFILE',
 !   RUNTYPE  = 'HISTORY',
 !   RUNTYPE  = 'ERROR',
 
    TIMEUNIT = 'h',          ! 'h','m','s' - Time units of TIMMAX, TIMSTR
 
    TIMMAX   = 36,           ! Final time of simulation
 
 !  Start of simulation or ISAN processing
 
    IMONTH1  = 10,           ! Month
    IDATE1   = 01,           ! Day
    IYEAR1   = 2019,         ! Year
    ITIME1   = 0000,         ! GMT of model TIME = 0.
 
 ! Grid specifications
 
    NGRIDS   = 1,            ! Number of grids to run
  
    NNXP     = 96,30,        ! Number of x gridpoints
    NNYP     = 96,50,        ! Number of y gridpoints
    NNZP     = 200,68,40,        ! Number of z gridpoints
    NZG      = 5,           ! Number of soil layers
    NZS      = 1,            ! Maximum number of snow layers
  
    NXTNEST  = 0,1,          ! Grid number which is the next coarser grid
    INESTING = 0,            ! 0-One way nesting, 1-Two way nesting
 
    IPRNTSTMT = 1,           ! Output extra print statements: 0=no,1=yes
 
 ! Coarse grid specifications
 
    IHTRAN   = 0,              ! 0-Cartesian, 1-Polar stereo
    DELTAX   = 62.5,          ! X-Y grid spacing (meters)
    DELTAZ   = 6.25,           ! Z grid spacing (set to 0. to use ZZ)
    DZRAT    = 1.,           ! Vertical grid stretch ratio
    DZMAX    = 1000.,          ! Maximum delta Z for vertical stretch
  
    ZZ       = 3212.875, 3219.125, 3225.375, 3231.625, 3237.875, 3244.125,
        3250.375, 3256.625, 3262.875, 3269.125, 3275.375, 3281.625,
        3287.875, 3294.125, 3300.375, 3306.625, 3312.875, 3319.125,
        3325.375, 3331.625, 3337.875, 3344.125, 3350.375, 3356.625,
        3362.875, 3369.125, 3375.375, 3381.625, 3387.875, 3394.125,
        3400.375, 3406.625, 3412.875, 3419.125, 3425.375, 3431.625,
        3437.875, 3444.125, 3450.375, 3456.625, 3462.875, 3469.125,
        3475.375, 3481.625, 3487.875, 3494.125, 3500.375, 3506.625,
        3512.875, 3519.125, 3525.375, 3531.625, 3537.875, 3544.125,
        3550.375, 3556.625, 3562.875, 3569.125, 3575.375, 3581.625,
        3587.875, 3594.125, 3600.375, 3606.625, 3612.875, 3619.125,
        3625.375, 3631.625, 3637.875, 3644.125, 3650.375, 3656.625,
        3662.875, 3669.125, 3675.375, 3681.625, 3687.875, 3694.125,
        3700.375, 3706.625, 3712.875, 3719.125, 3725.375, 3731.625,
        3737.875, 3744.125, 3750.375, 3756.625, 3762.875, 3769.125,
        3775.375, 3781.625, 3787.875, 3794.125, 3800.375, 3806.625,
        3812.875, 3819.125, 3825.375, 3831.625, 3837.875, 3844.125,
        3850.375, 3856.625, 3862.875, 3869.125, 3875.375, 3881.625,
        3887.875, 3894.125, 3900.375, 3906.625, 3912.875, 3919.125,
        3925.375, 3931.625, 3937.875, 3944.125, 3950.375, 3956.625,
        3962.875, 3969.125, 3975.375, 3981.625, 3987.875, 3994.125,
        4000.375, 4006.625, 4012.875, 4019.125, 4025.375, 4031.625,
        4037.875, 4044.125, 4050.375, 4056.625, 4062.875, 4069.125,
        4075.375, 4081.625, 4087.875, 4094.125, 4100.375, 4106.625,
        4112.875, 4119.125, 4125.375, 4131.625, 4137.875, 4144.125,
        4150.375, 4156.625, 4162.875, 4169.125, 4175.375, 4181.625,
        4187.875, 4194.125, 4200.375, 4206.625, 4212.875, 4219.125,
        4225.375, 4231.625, 4237.875, 4244.125, 4250.375, 4256.625,
        4262.875, 4269.125, 4275.375, 4281.625, 4287.875, 4294.125,
        4300.375, 4306.625, 4312.875, 4319.125, 4325.375, 4331.625,
        4337.875, 4344.125, 4350.375, 4356.625, 4362.875, 4369.125,
        4375.375, 4381.625, 4387.875, 4394.125, 4400.375, 4406.625,
        4412.875, 4419.125, 4425.375, 4431.625, 4437.875, 4444.125,
        4450.375, 4456.625, 4462.875, 4469.125, 4475.375, 4481.625,
        4487.875, 4494.125, 4500.375, 4506.625, 4512.875, 4519.125,
        4525.375, 4531.625, 4537.875, 4544.125, 4550.375, 4556.625,
        4562.875, 4569.125, 4575.375, 4581.625, 4587.875, 4594.125,
        4600.375, 4606.625, 4612.875, 4619.125, 4625.375, 4631.625,
        4637.875, 4644.125, 4650.375, 4656.625, 4662.875, 4669.125,
        4675.375, 4681.625, 4687.875, 4694.125, 4700.375, 4706.625,
        4712.875, 4719.125, 4725.375, 4731.625, 4737.875, 4744.125,
        4750.375, 4756.625, 4762.875, 4769.125, 4775.375, 4781.625,
        4787.875, 4794.125, 4800.375, 4806.625,
  
    DTLONG   = 1.,          ! Coarse grid long timestep
    NACOUST  = 5,           ! 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,2,5,4,         ! x-y-direction
    NNDTRAT  = 1,2,3,3,         ! Time
 
    NESTZ   = 0,                     ! Contort coarser grids if negative
    NSTRATZ = 1,1,3,3,3,3,2,2,2,2,1, !
                             
    POLELAT  = 85.0,         ! If polar, latitude/longitude of pole point
    POLELON  = 0.0,        ! If lambert, lat/lon of grid origin (x=y=0.)
 
    ! Grid point on the next coarser nest where the lower southwest
    !  corner of this nest will start.
    !  If NINEST or NJNEST = 0, use CENTLAT/CENTLON
    
    NINEST   = 1,0,0,0,           ! i-point
    NJNEST   = 1,0,0,0,           ! j-point
    NKNEST   = 1,1,1,1,           ! k-point
 
    CENTLAT =   85.0,37.0, ! Center lat/lon of grids, may or
    CENTLON =  0.0,-97.5, ! may not be same as pole 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
  
  $MODEL_FILE_INFO
 
 ! Variable initialization input
 
    INITIAL  = 1,             ! Initial fields - 1=horiz.homogeneous, 
                              !                  2=init from varfile
                              !                  3=init from history
 
    INORAINTIME = 0,     ! Seconds until rain is allowed
 
    ITEMPNUDGE = 0,           ! Run temperature nudging to initial sounding
                              ! 0 = off, 1 = on
 
    ITNTS = 1800, !3600,             ! Temperature Nudging time scale (seconds) 
 !  ------ Analysis nudging parameters -----------------------------------
    NUD_TYPE = 0,             ! =0 - no analysis nudging
                              ! =1 - nudge from varfiles
                              
    VARFPFX  = './var', ! Varfile prefix
    VWAIT1   = 0.,            ! wait between each VFILE check (s)
    VWAITTOT = 0.,            ! total wait before giving up on a VFILE (s)
    
    NUDLAT   = 5,             ! Number of points in lateral bnd region
    TNUDLAT  = 900.,          ! Nudging time scale(s) at lateral boundary
    TNUDCENT = 0.,            ! Nudging time scale(s) in center of domain
    TNUDTOP  = 00.,           ! Nudging time scale (s) at top of domain
    ZNUDTOP  = 15000.,        ! Nudging at top of domain above height(m)
    
    WT_NUDGE_G =  1., 0.8, 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
    
 !-------------------------------------------------------------------------
 
 !----------- Condensate nudging ------------------------------------------
 ! Only nudge total water where condensate exists. Even if turned on, this
 ! will not occur using traditional varfiles. Requires use of History-Varfiles
 ! made from a MAKEHFILE run where condensate is ingested from History files.
 
    NUD_COND   = 0,
    TCOND_BEG  = 0.,              ! Model time start of cond nudging (sec)
    TCOND_END  = 21600.,          ! Model time end of cond nudging (sec)
    T_NUDGE_RC = 3600.,           ! Cond nudging timescale for r_total
    WT_NUDGEC  =  1.,0.8,0.7,0.5, ! Relative nudging weights for active grids
                                  ! =0., turns off nudging for that grid
 !-------------------------------------------------------------------------
 
 !--------- Observation Data Assimilation (ODA) ---------------------------------
    IF_ODA = 0,                            ! Flag to turn on oda 
    ODA_UPAPREF = './zobs/dp-r',           ! File prefix for upper air obs
    ODA_SFCPREF = './zobs/dt-s',           ! File prefix for surface obs
    
    FRQODA=300.,                           ! Frequency of obs analysis
    TODABEG=0., 
    TODAEND=99999999.,         ! Model time start and end of oda (sec)
 
    TNUDODA= 900.,                         ! Nudging timescale for each grid
    WT_ODA_GRID =  1., 0.8, 0.7, 0.5,      ! Relative nudging weights for active grids
                                           ! =0., turns off nudging for that grid
 
    WT_ODA_UV = 1.,         ! ODA nudging weight for u and v
    WT_ODA_TH = 1.,         ! ODA nudging weight for theta
    WT_ODA_PI = 1.,         ! ODA nudging weight for pi
    WT_ODA_RT = 1.,         ! ODA nudging weight for r_tot
       
   ! Following are radii that affect the "smoothness" of the analyzed fields
   !     The SFCE and UPAE are the radii where the affect falls off to e**(-2)
   !     The SFC0 and UPA0 are the radii where the affect falls off to 0
   !     Values are grid dependent.
    
    RODA_SFCE = 50000.,100.,100.,100.,
    RODA_SFC0 = 100000.,100000.,100000.,100000.,
    RODA_UPAE = 100000.,200.,200.,200.,
    RODA_UPA0 = 200000.,2000.,2000.,2000.,
    
    RODA_HGT = 3000.,3000.,3000.,3000.,  ! Height at which transition from SFC radii 
                                         !   to UPA radii occurs
    
    RODA_ZFAC = 100.,100.,100.,100.,     ! Vertical factor related to dx/dz
    
    ! - Time interpolate limit (TIL)- if the future-past obs time 
    !    is > this limit, do not use to interpolate
    !
    ! - Time extrapolate limit (TEL)- if past/future obs is greater than TIL,
    !    but less than TEL, use the obs
    
    ODA_SFC_TIL=21600.,
    ODA_SFC_TEL=900.,
    ODA_UPA_TIL=43200.,
    ODA_UPA_TEL=21600.,
 !------------------------------------------------------------------------   
 
 ! History start (if RUNTYPE='HISTORY')
 
    !HFILIN   = './output_20170512/nosedim-A-2017-05-12-070000-head.txt',
    HFILIN   = './output/control/control-A-2019-07-03-000000-head.txt',
                              ! Input state file name

 ! Analysis file input for assimilation: 1=yes, 0=no
 
    IPAST_SFC = 0,            ! Initialize LEAF fields from history
 
 !------------------------------------------------------------------------
    ICLOBBER = 1,             ! 0=stop if files exist, 1=overwrite files
 !------------------------------------------------------------------------
 
 !  Analysis file output
  
    IOUTPUT  =  1,             ! 0-no files, 1-write files
 AFILEPREF = './output/20/control',
                              ! File prefix for all analysis files
                              !   state,mean,lite,both
                              
!-----------------------------------------------------------------------
   FRQSTATE   = 300.,600.,3600.,1800., ! "state" file write frequency
   !  can be different for each grid
   !  works best if grids are multiples of each other
FRQST_KEEP   =    0.,      ! If > 0., frequency at which "state" files will be 
   !  kept. Others will be written, but deleted on
   !  the subsequent write. 
   !  Intended to be a mutiple of and >= FRQSTATE
!-----------------------------------------------------------------------
FRQLITE  = 0.,            ! Frequency for "lite" files
  !   = 0 : no lite files
NLITE_VARS=0,

LITE_VARS= 'GLAT','GLON','TOPT','PATCH_AREA', !Keep these always
'UP','VP','WP', !This would be useful for hi-frequency trajectories

!-----------------------------------------------------------------------
AVGTIM   = 0.,            ! Averaging time for analysis variables
  !   must be abs(AVGTIM) <= FRQANL
  !   > 0 : averaging is centered at FRQANL
  !   < 0 : averaging ends at FRQANL
  !   = 0 : no averaged files
FRQMEAN  = 0.,            ! Frequency for "averaged" files
FRQBOTH  = 0.,            ! Frequency for Both ("averaged"/"lite") files

! Input topography variables

TOPFILES = './sfc/toph', ! Path and prefix for topo files. 
SFCFILES = './sfc/sfch', ! Path and prefix for surface files. 
SSTFPFX  = './sfc/ssth', ! Path and prefix for sst files
NDVIFPFX = './sfc/ndh',  ! Path and prefix for ndni files

ITOPTFLG = 2,0,0,0,           ! 2 - Fill data in "leaf3_init"
ISSTFLG  = 2,2,0,0,           ! 0 - Interpolate from coarser grid
IVEGTFLG = 2,2,0,0,           ! 1 - Read from standard Lat/Lon data file
ISOILFLG = 2,2,0,0,           ! 
NDVIFLG  = 2,2,0,0,           ! 

IUPDNDVI  = 0,             ! 0 - No update of NDVI values during run
IUPDSST  = 0,              ! 0 - No update of SST values during run
   ! 1 - Update values during run

  ! The following only apply for IxxxxFLG=1
ITOPTFN  = '/home/smsaleeb/sfctypehdf5/DEM30s-h5/EL',
'/home/smsaleeb/sfctypehdf5/DEM30s-h5/EL',
'/home/smsaleeb/sfctypehdf5/DEM30s-h5/EL',
'/home/smsaleeb/sfctypehdf5/DEM30s-h5/EL',
ISSTFN   = '/home/smsaleeb/sfctypehdf5/sst-h5/S',
'/home/smsaleeb/sfctypehdf5/sst-h5/S',
'/home/smsaleeb/sfctypehdf5/sst-h5/S',
'/home/smsaleeb/sfctypehdf5/sst-h5/S',
IVEGTFN  = '/home/smsaleeb/sfctypehdf5/ogedata-h5/GE',
'/home/smsaleeb/sfctypehdf5/ogedata-h5/GE',
'/home/smsaleeb/sfctypehdf5/ogedata-h5/GE',
'/home/smsaleeb/sfctypehdf5/ogedata-h5/GE',
ISOILFN  = '/home/smsaleeb/sfctypehdf5/FAOdata-h5/FAO',
'/home/smsaleeb/sfctypehdf5/FAOdata-h5/FAO',
'/home/smsaleeb/sfctypehdf5/FAOdata-h5/FAO',
'/home/smsaleeb/sfctypehdf5/FAOdata-h5/FAO',
NDVIFN   = '/home/smsaleeb/sfctypehdf5/ndvi-h5/',           
'/home/smsaleeb/sfctypehdf5/ndvi-h5/',
'/home/smsaleeb/sfctypehdf5/ndvi-h5/',
'/home/smsaleeb/sfctypehdf5/ndvi-h5/',    

! Topography scheme

ITOPSFLG = 3,3,3,3,           ! 0 = Average Orography
      ! 1 = Silhouette Orography
      ! 2 = Envelope Orography
      ! 3 = Reflected Envelope Orography

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

IZ0FLG   = 0,0,0,0,           ! 0 = Use veg, bare soil and water surface
      ! 1 = Also use subgrid scale topography
Z0MAX    = 5.,5.,5.,5.,       ! Max zo for IZ0FLG=1
Z0FACT   = 0.005,             ! Subgrid scale orograhic roughness factor

! Microphysics collection tables (Do for LEVEL=3 or 4)

MKCOLTAB = 1,             ! Make table: 0 = no, 1 = yes
COLTABFN = './ct.fthigh2k',
  ! Filename to read or write

$END

$MODEL_OPTIONS

! Numerical schemes

ICORFLG  = 0,             ! Coriolis flag/2D v-component  -  0=off, 1=on

IBND     = 2,             ! Lateral boundary condition flags
JBND     = 2,             ! 1-Klemp/Wilhelmson, 2-cyclic

CPHAS    = 30.,           ! 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     = 16,             ! Rayleigh friction - number of points from the top
DISTIM   = 60.,           !                   - dissipation time scale

! Radiation parameters

ISWRTYP  = 4,             ! Shortwave radiation type
ILWRTYP  = 4,             ! Longwave radiation type
  !    0-none, 2-Mahrer/Pielke, 1-Chen, 3-Harrington
RADFRQ   = 1.,          ! Freq. of radiation tendency update (s)
LONRAD   = 0,             ! Longitudinal variation of shortwave  
  !    (0-no, 1-yes)

! Cumulus parameterization parameters

NNQPARM  = 0,0,0,0,       ! Convective param flag (0-off, 1-Kuo, 2-Kain-Fritsch) 
CONFRQ   = 1200.,         ! Frequency of conv param. updates (s)
WCLDBS   = .01,           ! Vertical motion needed at cloud base for
  !    to trigger convection (Kuo)

! Surface layer DTHCONand soil parameterization

NPATCH   = 2,             ! Number of patches per grid cell (min=2)

NVEGPAT  = 1,             ! Number of patches per grid cell to be 
  !    filled from
  ! vegetation files 
  !    (min of 1, max of NPATCH-1)

ISFCL    = 0,             ! Surface layer/soil/veg model
  !  0-specified surface layer gradients
  !  1-soil/vegetation model

! CO2 vertical profile for SiB submodel initialization. [ppm]
CO2_INIT = 360., 360., 360., 355., 355., 355., 350., 350., 340.,

ISOILDAT = 0,             ! Soil moisture and temperature initialization
  ! 0 = homogeneous default, 1 = varfile ingest
ISNOWDAT = 0,             ! Snow depth and mass initialization
  ! 0 = homogeneous zero default, 1 = varfile ingest

NVGCON   = 2,             ! Vegetation type (see below)

!  0  Ocean                       ! 11  Tundra
!  1  Lakes, rivers, streams      ! 12  Evergreen shrub                      
!  2  Ice cap/glacier             ! 13  Deciduous shrub
!  3  Desert, bare soil           ! 14  Mixed woodland
!  4  Evergreen needleleaf tree   ! 15  Crop/mixed farming, C3 grassland
!  5  Deciduous needleleaf tree   ! 16  Irrigated crop
!  6  Deciduous broadleaf tree    ! 17  Bog or marsh
!  7  Evergreen broadleaf tree    ! 18  Wooded grassland 
!  8  Short grass                 ! 19  Urban and built up
!  9  Tall grass                  ! 20  Wetland evergreen broadleaf tree
! 10  Semi-desert                 ! 21  Very urban

PCTLCON  = 1.0,            ! Constant land % if for all domain
NSLCON   = 11,             ! 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   = .6,            ! Constant albedo if not running soil model
SEATMP   = 275.,          ! Constant water surface temperature

DTHCON   = 0.,            ! Constant sfc layer  temp grad for no soil
DRTCON   = 0.,            ! Constant sfc layer moist grad for no soil

SLZ      = -.50,-.40,-.30,-.25,-.20,-.16,-.12,-.09,-.06,-.03,-.01,
  ! Soil grid levels

SLMSTR   = 0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35,
  ! Initial soil moisture

STGOFF   = -2.,-2.,-2.,-2.,
  ! Initial soil temperature offset
  !  from lowest atmospheric level
!--------------------------------------------------------------------------

!--------- Eddy diffusion coefficient parameters -----------------------

DIVLS   = 5.0e-6,         ! Value from paper = 5e-6

IDIFFK  = 4,4,4,4,        ! 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)
  !  5 - Silvia's TKE-l scheme(horiz and vert same)
  !  6 - Silvia's TKE-eps scheme(horiz and vert same)
IHORGRAD = 1,             !  1 - horiz grad frm decomposed sigma grad
  !  2 - true horizontal gradient. 
  !      Non-conserving, but allows small DZ
CSX     = .2,.2,.2,.2,          ! Deformation horiz. K's coefficient
CSZ     = .2,.2,.2,.2,          ! Deformation vert. K's coefficient
XKHKM   = 3.,3.,3.,3.,          ! Ratio of horiz K_h to K_m for deformation
ZKHKM   = 3.,3.,3.,3.,          ! Ratio of vert K_h to K_m for deformation
AKMIN   = 0.1,2.,1.5,1.5,     ! Ratio of minimum horizontal eddy
  !    viscosity coefficientto typical value
  !    from deformation K

!------Bubble Simulation--------------------------------------------------
IBUBBLE = 3,              ! Bubble initialization: 0 = off    
  ! 1=square bubble, 2=gaussian bubble
  ! 3=random perturbations (user adjust ruser.f90)

IBUBBLEALL = 1,           ! 0 = perturbations linearly decrease
  ! 1 = perturbations throughout domain

IBUBGRD = 1,              ! Grid Number On Which To Initialize Bubble
IBDXIA  = 5,              ! DX start grid point
IBDXIZ  = 13,             ! DX end grid point
IBDYJA  = 5,              ! DY start grid point
IBDYJZ  = 13,             ! DY end grid point
IBDZK1  = 2,              ! Lower k level
IBDZK2  = 16,             ! Upper k level
BTHP    = 3.0,            ! Theta-il perturbation(K) (THP + BTHP)
BRTP    = 0.2,            ! Mixing ratio perturbation (RTP * (1 + BRTP)

!------Convergence Forcing -----------------------------------------------
ICONV = 0,              ! Convergence Zone initialization: 0 = off
! 1=Gaussian in x and y, U and V
! 2=Gaussian in x and y, U only
! 3=Gaussian in x and y, V only
! 4=Gaussian only in x, U only
! 5=Gaussian only in y, V only
ICONGR  = 1,            ! Grid for convergence
ICICENT = 100,          ! I center grid point
ICJCENT = 100,          ! J center grid point
CXRAD   = 7000.,        ! I radius width (meters)
CYRAD   = 7000.,        ! J radius width (meters)
ICVERT  = 1,            ! Convergence vertical structure
! 1=Linear decrease w/height and max at sfc
! 2=Elevated 1/2 vertical wavelength
ICKMAX  = 14,           ! Max K for ICVERT=1 (ignored for ICVERT=2)
CZRAD   = 2000.,        ! Vertical radius for ICVERT=2 (ignored for ICVERT=1)
ICKCENT = 17,           ! K center for ICVERT=2 (ignored for ICVERT=1)
CDIVMAX = -6.1224e-6,   ! Divergence amplitude (s^-1) (negative = convergence)
CTAU    = 1.,           ! timescale (s) over which conv. zone strength
!   increases to amplitude CDIVMAX
CTMAX   = 3600.,        ! Max time in seconds for convergence zone
! If ctmax < 0, abs(ctmax) vertical vel thresh past 
!   which to terminate conv zone (m/s)

!----- Non-standard lateral boundary forcing -----------------------------

IAEROLBC   = 0,0,0,       ! Aerosol profile forcing: 0=no, 1=yes
ICO2LBC    = 0,0,0,       ! CO2 profile forcing: 0=no, 1=yes
BCTAU      = 1.,1.,1.,    ! (seconds) Inverse time weight

!----- Extra Tracers and Scalars -----------------------------------------

ITRACER   = 0,            ! Integer number of tracers. Initialize in init_tracer
ITRACHIST = 0,            ! Initialize tracers on history start and/or restart

!----- Microphysics ------------------------------------------------------

LEVEL = 3,                ! Moisture complexity level
  ! (0=dry, 1=vapor only, 2=cloud only,
  !  3=full micro, 4=HUCM-SBM bin micro)

ICHECKMIC = 1,            ! Check Micro Negative & Nans: 0=Off,1=on

IMBUDGET = 3,             ! Micro budgets: 0=Off,1=partial,2=all,3=all+dust

IRIME   = 1,              ! Riming: 0-Bulk, 1-HybridBin 
IPLAWS  = 2,              ! FallSpeed: 0-Original4.3 
  !            1-New (R.Carver,Mitchell96)
  !            2-New With Banded Plaws for sedimentation
ISEDIM  = 1,              ! Sedimentation: 0-RelDisp, 1-True Hybrid-Bin

ICLOUD  = 5,              ! Microphysics flags
IDRIZ   = 0,              !-------------------
IRAIN   = 5,              !  1 - diagnostic concen.
IPRIS   = 5,              !  2 - specified mean diameter
ISNOW   = 5,              !  3 - specified y-intercept
IAGGR   = 5,              !  4 - specified concentration
IGRAUP  = 5,              !  5 - prognostic concentration
IHAIL   = 5,

! Hydrometeor parameters     ! Microphysics parameters
  ! -----------------------
CPARM   = .3e9,           ! Characteristic diameter, # concentration
DPARM   = .1e6,           !   or y-intercept
RPARM   = 1e-3,           !
PPARM   = 0.,             ! Specify values for the given
SPARM   = 1e-3,           !   microphysics flag where appropriate
APARM   = 1e-3,           !
GPARM   = 1e-3,           ! Only used for "microphysics flags" 1-4
HPARM   = 3e-3,           !

!Gamma shape pararmeters     ! cld rain pris snow aggr graup hail driz
GNU     = 4.,4.,4.,4.,4.,4.,4.,4.,

! Aerosol parameters         ! See mic_init.f90 for setting profiles
  ! ------------------------------------------
IAEROSOL = 1,             ! CCN,GCCN: 0=off, 1=userset
ISALT    = 3,             ! Sea Salt: 0=off, 1=userset, 2=saltmodel, 3=separate BL/FT concentratios
IDUST    = 3,             ! Dust: 0=off, 1=userset 2=dustmodel, 3=separate BL/FT concentratios

ICCNLEV = 2,              ! Aerosol treatment level
  ! 0 = no source or sink
  ! 1 = nucleation + deposition sink
  ! 2 = same as 1 + evaporation regeneration

IFORCECCN = 9,            ! Force aerosol values
  ! 0 = no forcing
  ! 1 = BL decrease, FT constant
  ! 2 = FT decrease, BL constant
  ! 3 = Decrease in both BL and FT

BLH = 750.,               !Definition of boundary layer height (m) 
  !for FORCECCN

FCCNTS = 2297.,          !Decay timescale for FORCECCN. 0 = no decay.

FCCNSTART = 7200.,        !Time after which to start decay
  !Before FCCNSTART, aerosol forced to CCN_MAX
IREGENDUST = 1,
IIFN = 3,                 ! Ice Nuclei treatment level
  ! 0 = no ice nuclei
  ! 1 = Meyers formula (CIN_MAX)
  ! 2 = DeMott(2010) (CIN_MAX)
  ! 3 = DeMott(2010) (aerosols D>0.5 microns)

IAERORAD  = 0,            ! Aerosol-Radiation Model: 0 = Off, 1 = On

IAERODEP  = 0,            ! Aerosol wet and dry deposition
  ! (track scavenged aerosols if iccnlev>=2)
  ! 0 = Off, 1 = On

  ! Aerosol initial maximum concentrations
  ! userset option. Default vertical profiles:
  ! MAXNUMBER*exp(-zt(k)/7000.)
  ! ------------------------------------------
IAEROPRNT  = 1,           ! Print aerosol init profiles (0=no,1=yes)
IAEROHIST  = 0,           ! Re-initialize aerosols on hist restart:0,1 
CIN_MAX    = 0.,        ! Ice Nuclei (#/mg)
CCN_MAX    = 0.,       ! CCN (#/mg) normally 170
GCCN_MAX   = 0.,          ! Giant-CCN (#/mg)
DUST1_MAX  = 5.,        ! Small Dust Mode (#/mg)
DUST2_MAX  = 0.,         ! Large Dust Mode (#/mg)
SALTF_MAX  = 20.,        ! Film Sea Salt (#/mg)
SALTJ_MAX  = 0.,          ! Jet Sea Salt (#/mg)
SALTS_MAX  = 0.,       ! Spume Sea Salt (#/mg)

DUST1_MAX_BL = 20.0,
DUST1_MAX_FT = 20.0,
SALTF_MAX_BL = 20.,
SALTF_MAX_FT = 20.,

! Aerosol 9 categories soluble chemistry type: 1=(NH4)2S04, 2=NaCl
IAERO_CHEM   = 2,2,2,2,2,2,2,2,2,

! Aerosol 9 categories solubility fraction: (0->1)
AERO_EPSILON = 0.90,0.90,0.05,0.05,1.00,1.00,1.00,1.00,0.05, 

! Aerosol 9 categories initial median radius (meters)
AERO_MEDRAD  = 0.127e-6,3.00e-6,          !CCN,GCCN
0.50e-6,0.50e-6,          !Dust1,Dust2
0.10e-6,1.00e-6,6.00e-6,  !SeaSalt: film,jet,spume
0.10e-6,0.50e-6,          !Regenerated 1 & 2 default

! Aerosol tracking flags: 0=off, 1=on
ITRKEPSILON = 0,        ! Tracking aerosol solubility for regeneration
ITRKDUST    = 1,        ! Tracking dust sub-species within hydrometeors
ITRKDUSTIFN = 1,        ! Tracing dust that served as ice nuclei

$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   = 2,             ! Specifies what is in TS array
  !  0-temp(C) 1-temp(K) 2-pot. temp(K)

IRTSFLG  = 2,             ! 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.0,

PS= 0.,50.00,100.00,150.00,200.00,250.00,300.00,350.00,400.00,
450.00,500.00,550.00,600.00,650.00,700.00,750.00,800.00,825.00,
850.00,875.00,900.00,925.00,950.00,1000.00,1050.00,1100.00,1150.00,
1200.00,1250.00,1300.00,1350.00,1400.00,1450.00,1500.00,1550.00,
1600.00,1650.00,1700.00,1750.00,1800.00,1850.00,1900.00,2000.00,
2100.00,2200.00,2300.00,2400.00,2500.00,

! original sounding temps
TS= 267.39,267.44,268.02,268.31,268.59,268.68,268.74,268.71,268.78,
268.86,268.92,268.99,269.08,269.25,269.35,269.52,269.57,269.67,269.73,
269.79,270.93,273.27,274.37,275.45,275.96,276.53,277.42,278.23,278.99,
279.62,280.30,281.06,281.87,282.42,283.10,283.67,284.35,284.86,285.07,
285.13,285.14,285.77,286.90,287.95,288.69,289.37,289.81,289.95,

! original sounding dewpoints
RTS= 2.73,2.68,2.71,2.72,2.53,2.50,2.46,2.43,2.42,2.38,2.34,2.31,2.29,
2.31,2.29,2.32,2.25,2.28,2.27,2.25,2.78,2.26,2.27,2.14,2.54,2.23,
1.92,1.95,1.88,1.71,1.57,1.32,1.22,1.19,1.12,1.23,1.04,1.24,1.14,
1.78,1.79,1.43,1.23,7.32,9.78,1.20,1.35,1.89,

! original sounding "u"
US=-1.70,-2.03,-1.93,-1.49,-1.23,-1.28,-1.38,-1.47,-1.65,-1.95,-1.91,
-1.31,-1.08,-1.63,-1.90,-1.72,-1.42,-1.31,-1.28,-1.24,-1.15,0.96,
0.78,0.62,0.89,-1.20,-1.02,0.42,0.14,0.31,0.34,0.05,0.18,0.10,0.05,
0.08,0.46,0.60,0.72,0.94,1.45,2.28,2.70,3.28,2.96,3.84,3.45,3.18, 

! original sounding "v"
VS= -0.77,-0.76,-0.54,-0.60,-0.81,-0.63,-0.39,-0.42,-0.59,-0.53,-0.09,
0.67,0.91,0.53,0.38,0.79,1.25,0.70,-0.46,-1.71,-3.01,-3.99,-4.21,-3.65,-3.47,-4.12,
-5.08,-6.09,-6.32,-5.69,-5.39,-5.89,-6.07,-6.14,-6.58,-7.24,-7.28,
-6.70,-6.91,-7.46,-7.59,-7.19,-6.85,-6.84,-7.46,-8.38,-8.03,-8.11,

$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

I1ST_FLG = 2,            ! What to do if first guess file (IAPR)
 !  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

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
IDATAIN = 0,             ! 1 = NARR, 0 = All Others for now
IAPR    = './dp-p', ! Input press level dataset
IARAWI  = ' ', ! Archived rawindsonde file name
IASRFCE = ' ', ! Archived surface obs file name

! Header file name for history nudging files (only prefix is used)
VAR_HFILE = './a-A-',

! File names and dispose flags 

VARPFX    = './var',  ! 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  = 4,            ! 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=don't

!-----------------------------------
! 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's to be excluded
 !  Prefix with 'r' for rawindsonde, 
 !              's' for surface

USED_FILE = 'none',      ! Filename prefix to output stations actually
 !   used in the analysis

IOBSWIN = 1800,          ! Window (seconds) around analysis time. Obs
 !   outside this window will not be used.

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.


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,.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   = 900.,900.,900.,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  = 300.,300.,300.,300.,  ! Wavelength for surface objective analysis

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

$END
