!namelist

 $MODEL_GRIDS

!  Simulation title (64 chars)

   EXPNME   = 'RAMS Simulations',

!  Type of model run to perform

  RUNTYPE  = 'INITIAL',                                                             
!   RUNTYPE  = 'INITIAL',                       
!   RUNTYPE  = 'MAKEVFILE',
!   RUNTYPE  = 'MAKEHFILE',
!   RUNTYPE  = 'HISTORY',
!   RUNTYPE  = 'ERROR',

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

   TIMMAX   = 9000.,           ! Final time of simulation

!  Start of simulation or ISAN processing

   IMONTH1  = 01,           ! Month
   IDATE1   = 26,           ! Day
   IYEAR1   = 1999,         ! Year
   ITIME1   = 1200,         ! GMT of model TIME = 0.

! Grid specifications

   NGRIDS   = 1,            ! Number of grids to run
 
   NNXP     = 400,        ! Number of x gridpoints
   NNYP     = 400,        ! Number of y gridpoints
   NNZP     = 82,        ! Number of z gridpoints
   NZG      = 8,           ! Number of soil layers
   NZS      = 1,            ! Maximum number of snow layers
 
   NXTNEST  = 0,1,          ! Grid number which is the next coarser grid

   IPRNTSTMT = 1,           ! Output extra print statements: 0=no,1=yes

! Coarse grid specifications

   IHTRAN   = 0,              ! 0-Cartesian, 1-Polar stereo
   DELTAX   = 250.,          ! X-Y grid spacing (meters)
   DELTAZ   = 50.,           ! Z grid spacing (set to 0. to use ZZ)
   DZRAT    = 1.05,           ! Vertical grid stretch ratio
   DZMAX    = 450.,          ! 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   = 2.,          ! Coarse grid long timestep
   NACOUST  = 5,           ! Small timestep ratio

   ! 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  = 10.0,         ! If polar, latitude/longitude of pole point
   POLELON  = -20.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 =   10.0,10.0, ! Center lat/lon of grids, may or
   CENTLON =  -20.0,-20.0, ! 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

!  ------ 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   = '/nobackup/rstorer/controls/aug17/out--A-1999-01-26-140000-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=overwite files
!------------------------------------------------------------------------

!  Analysis file output
 
   IOUTPUT  =  1,             ! 0-no files, 1-write files
   AFILEPREF  = '/nobackup/rstorer/convperts/aug17-control/out-',  
                              ! 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=18,

   LITE_VARS= 'GLAT','GLON','TOPT','PATCH_AREA',
   'UC','VC','WC','RAP', 'RCP', 'RDP', 'RGP', 'RHP', 'RPP', 'RRP',
   'RSP', 'RV', 'THETA', 'PI', 
!-----------------------------------------------------------------------
   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 = '/nobackup/rstorer/convperts/aug17-control/toph', ! Path and prefix for topo files. 
   SFCFILES = '/nobackup/rstorer/convperts/aug17-control/sfch', ! Path and prefix for surface files. 
   SSTFPFX  = '/nobackup/rstorer/convperts/aug17-control/ssth', ! Path and prefix for sst files
   NDVIFPFX = '/nobackup/rstorer/convperts/aug17-control/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.storm',
                             ! Filename to read or write

 $END
 
 $MODEL_OPTIONS

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

   IBND     = 1,             ! Lateral boundary condition flags
   JBND     = 1,             ! 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     = 6,             ! Rayleigh friction - number of points from the top
   DISTIM   = 60.,           !                   - dissipation time scale
 
! Radiation parameters
 
   ISWRTYP  = 3,             ! Shortwave radiation type
   ILWRTYP  = 3,             ! Longwave radiation type
                             !    0-none, 2-Mahrer/Pielke, 1-Chen, 3-Harrington
   RADFRQ   = 600.,          ! 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 and 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    = 1,             ! 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
   SNUDCENT = 0.,            ! Soil moisture interior nudging time scale (sec)

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

   NVGCON   = 0,             ! 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  = 0.0,            ! 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   = .00015,           ! Constant roughness if for all domain
   ALBEDO   = .05,            ! Constant albedo if not running soil model
   SEATMP   = 299.85,          ! 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   = 5.,5.,5.,5.,3.5,2.,.5,-1.,-1.5,-1.8,-2.,
                             ! Initial soil temperature offset
                             !  from lowest atmospheric level
!--------------------------------------------------------------------------
 
!--------- Eddy diffusion coefficient parameters -----------------------
 
   IDIFFK  = 2,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)
                             !  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   = 1.0,1.0,1.5,1.5,      ! Ratio of minimum horizontal eddy
                             !    viscosity coefficientto typical value
                             !    from deformation K

!------Bubble Simulation--------------------------------------------------
   IBUBBLE = 0,              ! Bubble initialization: 0 = off    
                             ! 1=square bubble, 2=gaussian bubble
                             ! 3=random perturbations (user adjust ruser.f90)
   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 = 2,              ! 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 = 300,          ! I center grid point
   ICJCENT = 200,          ! J center grid point
   CXRAD   = 7000.,        ! I radius width (meters)
   CYRAD   = 21000.,        ! 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 = -1.2e-5,   ! 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   = 24,            ! 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 = 2,             ! 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   = 5,              !-------------------
   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.,2.,2.,2.,2.,2.,2.,4.,

!HUCM parameters
   NDTCOLL = 3,              !Collision-coalescence called once
                             !every NDTCOLL time steps

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

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

   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)
   IIFN_FORMULA = 2,         ! DeMott IN Formulation if IIFN = 2 or 3
                             ! 1 = DeMott(2010) (composite aerosol types)
                             ! 2 = DeMott(2014) (dust formula)

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

   IAERODEP  = 1,            ! 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.01,        ! Ice Nuclei (#/mg)
   CCN_MAX    = 1000.,       ! CCN (#/mg)
   GCCN_MAX   = 10.,         ! Giant-CCN (#/mg)
   DUST1_MAX  = 100.,        ! Small Dust Mode (#/mg)
   DUST2_MAX  = 10.,         ! Large Dust Mode (#/mg)
   SALTF_MAX  = 100.,        ! Film Sea Salt (#/mg)
   SALTJ_MAX  = 1.,          ! Jet Sea Salt (#/mg)
   SALTS_MAX  = 1.e-3,       ! Spume Sea Salt (#/mg)

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

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

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

   ! Aerosol tracking flags: 0=off, 1=on
   ITRKEPSILON = 1,        ! 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   = 0,             ! 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  = 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       = 88.2,
   Sounding:  aug17  control  
   PS = 1000.0, 990.0, 980.0, 970.0, 960.0, 950.0, 940.0, 
   930.0, 920.0, 910.0, 900.0, 890.0, 880.0, 
   870.0, 860.0, 850.0, 840.0, 830.0, 820.0, 
   810.0, 800.0, 790.0, 780.0, 770.0, 760.0, 
   750.0, 740.0, 730.0, 720.0, 710.0, 700.0, 
   690.0, 680.0, 670.0, 660.0, 650.0, 640.0, 
   630.0, 620.0, 610.0, 600.0, 590.0, 580.0, 
   570.0, 560.0, 550.0, 540.0, 530.0, 520.0, 
   510.0, 500.0, 490.0, 480.0, 470.0, 460.0, 
   450.0, 440.0, 430.0, 420.0, 410.0, 400.0, 
   390.0, 380.0, 370.0, 360.0, 350.0, 340.0, 
   330.0, 320.0, 310.0, 300.0, 290.0, 280.0, 
   270.0, 260.0, 250.0, 240.0, 230.0, 220.0, 
   210.0, 200.0, 190.0, 180.0, 170.0, 160.0, 
   150.0, 140.0, 130.0, 120.0, 110.0, 100.0, 
   90.0, 80.0, 70.0, 60.0, 50.0, 40.0, 
   30.0, 20.0, 10.0, 
 
   TS = 25.75, 25.18, 24.61, 24.05, 23.48, 22.91, 22.34, 
   21.77, 21.2, 20.63, 20.07, 19.5, 18.93, 
   18.36, 17.79, 17.22, 16.65, 16.08, 15.52, 
   14.95, 14.38, 13.81, 13.24, 12.67, 12.1, 
   11.54, 10.97, 10.4, 9.83, 9.24, 8.43, 
   7.62, 6.8, 5.99, 5.17, 4.36, 3.54, 
   2.73, 1.92, 1.1, 0.29, -0.53, -1.34, 
   -2.16, -2.97, -3.83, -4.14, -4.76, -5.59, 
   -6.54, -7.48, -8.43, -9.37, -10.32, -11.34, 
   -12.4, -13.51, -14.75, -15.99, -17.23, -18.46, 
   -19.58, -20.44, -21.84, -23.48, -25.12, -26.77, 
   -28.41, -30.1, -31.87, -33.64, -35.41, -37.21, 
   -39.72, -42.23, -44.74, -47.25, -49.76, -52.27, 
   -54.79, -57.3, -59.81, -62.32, -64.83, -67.34, 
   -69.6, -71.77, -72.56, -71.01, -73.49, -75.7, 
   -76.64, -73.95, -71.18, -68.39, -65.6, -62.81, 
   -60.02, -57.23, -54.44, 
 
   RTS = 17.9067, 17.4381, 16.981, 16.5353, 16.1006, 15.6768, 15.2636, 
   14.8607, 14.468, 14.0942, 13.7209, 13.357, 13.0024, 
   12.6568, 12.3201, 11.992, 11.6724, 11.3686, 11.0652, 
   10.7697, 10.4819, 10.2016, 9.9288, 9.6632, 9.4046, 
   9.1592, 8.9141, 8.6756, 8.4436, 8.212, 7.9533, 
   7.697, 7.4487, 7.2082, 6.9754, 6.75, 6.5317, 
   6.3205, 6.116, 5.9182, 5.7269, 5.5459, 5.3668, 
   5.1937, 5.0263, 4.835, 4.1817, 3.8871, 3.7653, 
   3.5312, 3.3131, 3.105, 2.9089, 2.7265, 2.3783, 
   1.9739, 1.6865, 1.5871, 1.4947, 1.4075, 1.3252, 
   1.1696, 0.8986, 0.7629, 0.678, 0.6017, 0.5328, 
   0.4717, 0.3944, 0.2935, 0.2164, 0.158, 0.1148, 
   0.0936, 0.076, 0.0615, 0.0496, 0.0398, 0.0318, 
   0.0253, 0.0201, 0.0159, 0.0125, 0.0098, 0.0076, 
   0.006, 0.0048, 0.0046, 0.0063, 0.0048, 0.0038, 
   0.0002, 0.0002, 0.0002, 0.0003, 0.0003, 0.0004, 
   0.0004, 0.0005, 0.0005, 
 
   US = -8.6219, -8.7553, -8.8984, -9.0322, -9.1666, -9.3103, -9.4445, 
   -9.5885, -9.7232, -9.868, -10.0029, -10.1476, -10.2828, 
   -10.4277, -10.5637, -10.6993, -10.8447, -10.9805, -11.1261, 
   -11.2622, -11.4085, -11.5448, -11.6908, -11.8272, -11.9734, 
   -12.1106, -12.2569, -12.3937, -12.5305, -12.6489, -12.4302, 
   -12.2091, -11.9866, -11.7521, -11.526, -11.2978, -11.0684, 
   -10.8272, -10.5946, -10.3602, -10.1247, -9.8777, -9.6395, 
   -9.3999, -9.1593, -8.9732, -8.9998, -8.9394, -8.6966, 
   -8.3826, -8.0421, -7.6934, -7.3214, -6.9368, -5.9943, 
   -4.6993, -3.7354, -3.5361, -3.2472, -2.879, -2.4272, 
   -1.6127, -0.1689, 0.3819, 0.5389, 0.7085, 0.8909, 
   1.0834, 1.1009, 0.5896, -0.1612, -1.1497, -2.2697, 
   -1.1046, 0.1695, 1.5463, 3.0174, 4.5817, 6.2287, 
   7.9501, 9.7367, 11.58, 13.469, 15.3984, 17.3463, 
   12.2192, 5.8028, 0.1613, -2.8856, -0.6362, 6.1031, 
   -2.6016, -3.3528, -4.5091, -4.6423, -4.7754, -4.9086, 
   -5.0417, -5.1749, -5.308, 
 
   VS = -2.9419, -2.9533, -2.967, -2.9766, -2.9837, -2.9946, -3.0014, 
   -3.0104, -3.0154, -3.0207, -3.0239, -3.029, -3.0303, 
   -3.0335, -3.0311, -3.0296, -3.0299, -3.0266, -3.025, 
   -3.0198, -3.0143, -3.0072, -3.0017, -2.9928, -2.9853, 
   -2.9724, -2.963, -2.9503, -2.9367, -2.9132, -2.7466, 
   -2.5862, -2.4278, -2.2737, -2.1237, -1.9799, -1.8383, 
   -1.7013, -1.5682, -1.4413, -1.3168, -1.1971, -1.0812, 
   -0.9714, -0.8642, -0.6936, -0.0597, 1.0422, 2.3172, 
   2.7431, 3.1306, 3.4948, 3.8243, 4.1237, 4.4909, 
   4.6911, 4.8226, 5.4659, 6.0995, 6.7301, 7.3391, 
   7.8358, 7.9982, 8.2512, 8.6132, 8.9821, 9.3376, 
   9.6896, 10.6533, 12.7864, 14.8991, 16.9511, 18.8942, 
   19.8192, 20.6593, 21.4142, 22.0848, 22.6411, 23.0845, 
   23.4067, 23.6112, 23.6691, 23.585, 23.3617, 22.9776, 
   22.0622, 19.2441, 16.4992, 13.9956, 12.3036, 11.4783, 
   3.9143, 1.6345, 0.4342, -2.4944, -5.423, -8.3515, 
   -11.2801, -14.2087, -17.1373, 
 
	 
 $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

