!namelist

 $MODEL_GRIDS

   ! Simulation title (64 chars)

   EXPNME   = 'LP Supercell 300m',

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

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

   TIMMAX   = 10800.,           ! Final time of simulation

!  Start of simulation or ISAN processing

   IMONTH1  = 04,           ! Month
   IDATE1   = 26,           ! Day
   IYEAR1   = 2000,         ! Year
   ITIME1   = 2100,         ! GMT of model TIME = 0.

! Grid specifications

   NGRIDS   = 1,            ! Number of grids to run
 
   NNXP     = 950,100,      ! Number of x gridpoints
   NNYP     = 750,100,      ! Number of y gridpoints
   NNZP     = 92,35,        ! Number of z gridpoints
   NZG      = 11,           ! Number of soil layers
   NZS      = 1,            ! Maximum number of snow layers
 
   NXTNEST  = 0,1,          ! Grid number which is the next coarser grid

! Coarse grid specifications

   IF_ADAP  = 0,
   
   IHTRAN   = 0,              ! 0-Cartesian, 1-Polar stereo, 2-Lambert-conformal
   DELTAX   = 300.,
   DELTAY   = 300.,           ! X and Y grid spacing

   DELTAZ   = 25.,            ! Z grid spacing (set to 0. to use ZZ)
   DZRAT    = 1.1,            ! Vertical grid stretch ratio
   DZMAX    = 300.,           ! 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   = 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-direction
   NSTRATY  = 1,2,5,4,         ! y-direction
   NNDTRAT  = 1,2,3,3,         ! Time

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

                            
   POLELAT  = 37.0,         ! If polar, latitude/longitude of pole point
   POLELON  = -97.5,        ! If lambert, lat/lon of grid origin (x=y=0.)

   STDLAT1  = 30.,          ! If polar, unused
   STDLAT2  = 60.,          ! If lambert, standard latitudes of projection
                            ! stdlat1 must be <= stdlat2 (not working yet)
   
   ! 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 =   37.0,37.0, ! Center lat/lon of grids, may or
   CENTLON =  -97.5,-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.

   GRIDU    = 0.,0.,0.,0.,       ! u-component for moving grids
   GRIDV    = 0.,0.,0.,0.,       ! v-component for moving grids
                                 !  (not working again!)
 
 $END
 
 $MODEL_FILE_INFO

! Variable initialization input

   INITIAL  = 1,             ! Initial fields - 1=horiz.homogeneous, 
                             !                  2=init from varfile
                             !                  3=init from HFILIN



!  ------ Analysis nudging parameters -----------------------------------
   NUD_TYPE = 0,             ! =1 - nudge from history files(1-way nest)
                             ! =2 - nudge from varfiles
                             ! =0 - no analysis nudging
                             
   VARFPFX  = '/home/smsaleeb/data.storm/var', ! Varfile prefix
   VWAIT1   = 0.,            ! wait between each VFILE check (s)
   VWAITTOT = 0.,            ! total wait before giving up on a VFILE (s)
   
   NUD_HFILE = './a-H-2001-07-21-000000-head.txt', 
            ! Header file name for history nudging files (only prefix is used)
                             
   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_GRID =  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 ------------------------------------------
   NUD_COND    = 0,          ! Only nudge total water where condensate
                             !   exists (from previous history files)
   COND_HFILE = './zhist/a-H-2001-07-21-000000-head.txt', 
            ! Header file name for cond nudging history files (only prefix is used)
   TCOND_BEG=0., 
   TCOND_END=21600.,                    ! Model time start and end of cond nudging (sec)
   T_NUDGE_RC = 3600.,                  ! Cond nudging timescale for r_total
   WT_NUDGEC_GRID =  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_UPAPREFIX = './zobs/dp-r',          ! File prefix for upper air obs
   ODA_SFCPREFIX = './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_ZFACT = 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.,
!------------------------------------------------------------------------   

!----- Cumulus inversion tendency input -----
   IF_CUINV = 0,
   CU_PREFIX = './zmicro/t5-C-',
   
   TNUDCU=900.,
   WT_CU_GRID=1., 1., 1., 1.,
      
   TCU_BEG=0., 
   TCU_END=7200.,
   CU_TEL=3600.,
   CU_TIL=21600.,
!------------------------------------------------------------------------   

! History start (if RUNTYPE='HISTORY')

   HFILIN   = '/glade/scratch/ldgrant/z.CTL/a-A-1991-04-26-220000-head.txt',
                             ! Input state file name
 
! Analysis file input for assimilation (currently LEAF variables)

   IPASTIN   = 0,            ! Initialize various fields from analysis file?
                             !  1=yes, 0=no
   PASTFN    = '/home/smsaleeb/data.storm/a-A-2007-02-11-060000-head.txt',
                             ! Input analysis file name

!------------------------------------------------------------------------
   ICLOBBER = 1,             ! 0=stop if files exist, 1=overwite files
!------------------------------------------------------------------------

!  Analysis file output
 
   IOUTPUT  =  1,             ! 0-no files, 1-write files
   AFILEPREF  = '/glade/scratch/ldgrant/z.CTL/a',  
                              ! File prefix for all analysis files
                              !   state,mean,lite,both
                              
!------------------------------------------------------------------------
   FRQSTATE   = 300.,10800.,3600.,1800., ! "state" file write frequency
                              !  can be different for each grid
                              !  works best if grids are multiples of each other
   FRQSTATE_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
   XLITE    = '/0:0/',       ! nums>0 are absolute grid indexes (not yet)
   YLITE    = '/0:0/',       ! nums<0 count in from the domain edges (not yet)
   ZLITE    = '/0:0/',       ! nums=0 are domain edges (not yet)
   
   NLITE_VARS=5,
   LITE_VARS= 'GLAT','GLON','TOPT','TOPTA','PATCH_AREA',
!        'NUCCLDR','NUCCLDC','NUCICER','NUCICEC','INUCHOMR',
!        'INUCHOMC','INUCCONTR','INUCCONTC','INUCIFNR','INUCIFNC',
!        'INUCHAZR','INUCHAZC','VAPLIQ','VAPICE','VAPCLD',
!        'VAPRAIN','VAPPRIS','VAPSNOW','VAPAGGR','VAPGRAU',
!        'VAPHAIL','VAPDRIZ','MELTICE','MELTPRIS','MELTSNOW',
!        'MELTAGGR','MELTGRAU','MELTHAIL','CLD2RAIN','RIMECLD',
!        'RIMECLDSNOW','RIMECLDAGGR','RIMECLDGRAU','RIMECLDHAIL','RAIN2ICE',
!        'RAIN2PR','RAIN2SN','RAIN2AG','RAIN2GR','RAIN2HA',
!        'RAIN2HA_XTRA','ICE2RAIN','AGGREGATE','AGGRSELFPRIS','AGGRSELFSNOW',
!        'AGGRPRISSNOW',
!-----------------------------------------------------------------------
   
   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

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

   TOPFILES = '/glade/scratch/ldgrant/z.CTL/toph', ! Path and prefix for topo files. 
   SFCFILES = '/glade/scratch/ldgrant/z.CTL/sfch', ! Path and prefix for surface files. 
   SSTFPFX  = '/glade/scratch/ldgrant/z.CTL/ssth', ! Path and prefix for sst files
   NDVIFPFX = '/glade/scratch/ldgrant/z.CTL/ndh',  ! Path and prefix for ndni files

   ITOPTFLG = 2,2,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,           ! 
   
   NOFILFLG = 2,2,2,2,           ! 2 - Fill data in "leaf3_init"
                                 ! 0 - Interpolate from coarser grid   
 
   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

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

 $END
 
 $MODEL_OPTIONS

   NADDSC   = 0,             ! Number of additional scalar species

! Case Study or Large Scale Forcing Parameters

   LSTYPE   = 0,             ! Large scale forcing or Case Study flag
                             ! 0 = Off
                             ! 1 = Shallow Cumulus (GCSS RICO - Precipitating)
                             ! 2 = Shallow Cimulus (BOMEX - Nonprecipitating)
                             ! 3 = Keep mean state fixed after specified spinup.
                             !     User defined in rtimi.f90 [lsforcing()]
                             ! 4 = User defined in rtimi.f90 [lsforcing()]
   ISFCFORCING = 0,          ! Force the surface using parameterized fluxes
                             ! and shut off LEAF3 calculations
                             ! 0 = Off
                             ! 1 = On
   IOVERRIDE = 0,            ! Override RAMSIN with data in rtimi.f90?
                             ! This option is only for RUNTYPE=INITIAL
                             ! 1 = yes; 0 = no


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

   IBND     = 1,             ! Lateral boundary condition flags
   JBND     = 1,             ! 1-Klemp/Wilhelmson, 2-Klemp/Lilly, 3-Orlanski
                             !  4-cyclic
   CPHAS    = 30.,           ! Phase speed if IBND or JBND = 1
   LSFLG    = 1,             ! 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  = 0,             ! Shortwave radiation type
   ILWRTYP  = 0,             ! 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    = 0,             ! Surface layer/soil/veg model
                             !  0-specified surface layer gradients
                             !  1-soil/vegetation model

   NVGCON   = 3,             ! 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.,            ! 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   = .00,           ! Constant roughness if for all domain
   ALBEDO   = .2,            ! Constant albedo if not running soil model
   SEATMP   = 300.,          ! 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   = 3.,3.,3.,3.,2.5,1.5,1.,1.,.5,.4,.3,
                             ! Initial soil temperature offset
                             !  from lowest atmospheric level
!--------------------------------------------------------------------------

!----- Urban canopy parameterization --------------------------------

   IF_URBAN_CANOPY =0,

!-------------------------------------------------------------------------
 
!--------- 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 = 2,             !  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   = 2.,1.,1.5,1.5,     ! Ratio of minimum horizontal eddy
                             !    viscosity coefficientto typical value
                             !    from deformation K

!------Bubble Simulation--------------------------------------------------
   IBUBBLE = 1,              ! Bubble initialization: 0 = off    
                             ! 1=square bubble, 2=gaussian bubble
                             ! 3=random perturbations (user adjust ruser.f90)
			     ! 4=random pert within square
   IBDXIA  = 100,             ! DX start grid point
   IBDXIZ  = 133,            ! DX end grid point
   IBDYJA  = 200,            ! DY start grid point
   IBDYJZ  = 233,            ! DY end grid point
   IBDZK1  = 5,              ! Lower k level
   IBDZK2  = 29,             ! Upper k level
   BTHP    = 2.0,            ! Theta-il perturbation(K) (THP + BTHP)
   BRTP    = 0.0,            ! Mixing ratio perturbation (RTP * (1 + BRTP)

!------MCS Simulation using CONVERGENCE ----------------------------------
   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
   ICICENT = 50,				! I center grid point
   ICJCENT = 25,				! J center grid point
   ICIRAD  = 25,				! I width in grid points
   ICJRAD  = 25,				! J width in grid points
   ICVERT  = 2,				! Convergence vertical structure
                           ! 1=Linear decrease w/height and max at sfc
                           ! 2=Elevated 1/2 vertical wavelength
   ICKMAX  = 12,				! Max K for ICVERT=1 (ignored for ICVERT=2)
   CZRAD   = 2000.,			! Vertical radius for ICVERT=2 (ignored for ICVERT=1)
   ICKCENT = 15,				! K center for ICVERT=2 (ignored for ICVERT=1)
   CAMP    = 3.,				! Amplitude constant
   CTMAX   = 2400.,			! Max time in seconds for convergence zone

!----- Airplane Simulator ------------- ----------------------------------
   IPLANE     = 0,          ! Airplane initialization: 0 = off
                            ! 1=Single sampler

   PFILEPREF  = '/scratch/rseigel/ccn100/FSIM',      ! File prefix for the single output file

   IPNUM      = 3,             ! Number of airplanes to simulate

! The remaining IPLANE variables must have IPNUM values
   IPHROUTE   = 1,1,1,         ! Horizontal flight route/pattern of airplane
                               ! 1 = East to west
                               ! 2 = North to south
   IPVROUTE   = 2,3,4,         ! Vertical flight route/pattern of airplane
                               ! 1 = Constant height (given by IPKSTART)
                               ! 2 = Cloud base (see plane_sim.f90 for defs)
                               ! 3 = Cloud top (see plane_sim.f90 for defs)
                               ! 4 = Cloud center (see plane_sim.f90 for defs)

   PSPEED     = 50.,50.,50.,   ! Airplane ground-relative speed (m/s)
                               ! If set to a negative value, IPROUTE will
                                 fly in opposite direction
   IPISTART   = 250,250,250,      ! i starting location of airplane
   IPJSTART   = 250,250,250,      ! j starting location of airplane
   IPKSTART   = 20,20,20,      ! k starting location of airplane

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

   LEVEL   = 3,              ! Moisture complexity level

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

   IRIME   = 1,              !   0 - Bulk riming, 1 - bin riming 
   IPLAWS  = 1,              !   0 - Original 4.3 Power Laws
                             !   1 - New Power Laws (R.Carver,Mitchell 1996)
 
   ICLOUD  = 6,              ! Microphysics flags
   IDRIZ   = 5,              !-------------------
   IRAIN   = 5,              !  1 - diagnostic concen.
   IPRIS   = 6,              !  2 - specified mean diameter
   ISNOW   = 5,              !  3 - specified y-intercept
   IAGGR   = 5,              !  4 - specified concentration
   IGRAUP  = 5,              !  5 - prognostic concentration
   IHAIL   = 5,

   ICCNLEV = 0,              ! Aerosol treatment level
                             !-------------------
                             ! 0 = no source or sink
                             ! 1 = nucleation + deposition sink
                             ! 2 = same as 1 + evaporation regeneration
                             ! 3 = same as 2 + dust mass tracking and
                             !     surface mass accumulation

   IIFN = 2,                 ! Ice Nuclei treatment level
                             !-------------------
                             ! 0 = no ice nuclei
                             ! 1 = Meyers formula (pparm)
                             ! 2 = DeMott(2010) procedure

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

! Hydrometeor parameters     ! Microphysics parameters
                             ! -----------------------
   CPARM   = .3e9,           ! Characteristic diameter, # concentration
   DPARM   = .1e6,           !   or y-intercept
   RPARM   = 1e-3,           !
   PPARM   = 1e5,            ! If ICLOUD = 5+, Set CCN,GCCN aerosol
   SPARM   = 1e-3,           !   concentration in mic_init.f90
   APARM   = 1e-3,           ! If IPRIS  = 5+, PPARM = Max # IFN (#/kg)
   GPARM   = 1e-3,           ! For ICLOUD or IPRIS >= 5
   HPARM   = 3e-3,           !   (5 = 3D homogeneous 6 = vertical profile)

! Aerosol initial maximum concentrations (#/cm3)
   CCN_MAX    = 100.,       ! CCN
   GCCN_MAX   = 1.e-5,       ! Giant-CCN
   DUST1_MAX  = 100.,        ! Small Dust Mode
   DUST2_MAX  = 10.,         ! Large Dust Mode
   SALTF_MAX  = 100.,        ! Film Sea Salt
   SALTJ_MAX  = 1.,          ! Jet Sea Salt
   SALTS_MAX  = 1.e-3,       ! Spume Sea Salt

! Salt and Dust parameters   ! See mic_init.f90 for 3D dust concentration
                             ! ------------------------------------------
   ISALT   = 0,              ! Sea Salt: 0=off, 1=userset, 2=saltmodel
   IDUST   = 1,              ! Dust: 0=off, 1=userset 2=dustmodel
                             ! If IDUST=2 values below are over-ridden

   GNU     = 2.,2.,2.,2.,2.,2.,2.,2., ! Gamma shape parms for
                                      ! cld rain pris snow aggr graup hail
!-----------------------------------------------------------------------------

 $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.,

! Sounding: WK dry low levels (sounding info from RAMSIN.WK.dryLow
  PS= 1000.0,    250,    500,    750,   1000,  
        1250,   1500,   1750,   2000,   2250,  
        2500,   2750,   3000,   3250,   3500,
        3750,   4000,   4250,   4500,   4750,
        5000,   5250,   5500,   5750,   6000,
        6250,   6500,   6750,   7000,   7250,
        7500,   7750,   8000,   8250,   8500,
        8750,   9000,   9250,   9500,   9750,
       10000,  10250,  10500,  10750,  11000,
       11250,  11500,  11750,  12000,  12250,
       12500,  12750,  13000,  13250,  13500,
       13750,  14000,  14250,  14500,  14750,
       15000,  15250,  15500,  15750,  16000,
       16250,  16500,  16750,  17000,  17250,
       17500,  17750,  18000,  18250,  18500,
       18750,  19000,  19250,  19500,  19750,
       20000,  20250,  20500,  20750,  21000,
       21250,  21500,  21750,  22000,  22250,
       22500,  22750,  23000,  23250,  23500,
       23750,  24000,  25000,  26000,  27000,

! theta for original WK
!  TS= 300.0000, 300.3725, 300.8257, 301.3556, 301.9347,
!      302.5526, 303.2029, 303.8813, 304.5847, 305.3105,

! theta well-mixed in lowest 7 levels (goes with 11g/kg moisture)
!  TS= 303.2029, 303.2029, 303.2029, 303.2029, 303.2029,
!      303.2029, 303.2029, 303.8813, 304.5847, 305.3105,

! theta well-mixed in lowest 6 levels (goes with 12g/kg moisture)
!  TS= 302.5526, 302.5526, 302.5526, 302.5526, 302.5526,
!      302.5526, 303.2029, 303.8813, 304.5847, 305.3105,

! theta well-mixed in lowest 5 levels (goes with 13g/kg moisture)
  TS= 301.9347, 301.9347, 301.9347, 301.9347, 301.9347,
      302.5526, 303.2029, 303.8813, 304.5847, 305.3105,
      306.0570, 306.8224, 307.6055, 308.4052, 309.2204,
      310.0504, 310.8943, 311.7516, 312.6216, 313.5037,
      314.3976, 315.3027, 316.2187, 317.1451, 318.0817,
      319.0281, 319.9841, 320.9492, 321.9234, 322.9063,
      323.8977, 324.8974, 325.9052, 326.9210, 327.9444,
      328.9754, 330.0138, 331.0594, 332.1122, 333.1719,
      334.2384, 335.3116, 336.3913, 337.4775, 338.5701,
      339.6688, 340.7737, 341.8846, 343.6910, 346.8917,
      350.8164, 354.7856, 358.7997, 362.8591, 366.9645,
      371.1164, 375.3152, 379.5616, 383.8559, 388.1989,
      392.5910, 397.0328, 401.5248, 406.0677, 410.6620,
      415.3082, 420.0070, 424.7590, 429.5647, 434.4249,
      439.3400, 444.3107, 449.3376, 454.4214, 459.5628,
      464.7623, 470.0206, 475.3385, 480.7165, 486.1553,
      491.6557, 497.2183, 502.8438, 508.5169, 513.3480,
      518.7270, 524.1070, 529.4870, 534.8670, 540.2460,
      545.6260, 551.0060, 556.3860, 561.7650, 567.1450,
      572.5250, 577.9050, 600.0000, 622.0000, 644.0000,


! well-mixed moisture sounding (drier low levels)
!  RTS= 11.0000,  11.0000,  11.0000,   11.0000,   13.8353,  
!       13.2112,  12.1701,  11.0391,    9.9877,    9.0267,   

! HIGHEST CAPE SOUNDING (original WK? -- first line only)
!  RTS= 16.0000,  16.0000,  16.0000,   16.0000,   13.8353,  
!       13.2112,  12.1701,  11.0391,    9.9877,    9.0267,   

! 11 g/kg moisture well-mixed
!  RTS= 11.0000,  11.0000,  11.0000,   11.0000,   11.0000,  
!       11.0000,  11.0000,  11.0000,    9.9877,    9.0267,   

! 12 g/kg moisture well-mixed
!  RTS= 12.0000,  12.0000,  12.0000,   12.0000,   12.0000,  
!       12.0000,  12.0000,  11.0391,    9.9877,    9.0267,  

! 13 g/kg moisture well-mixed, dry layer ~800-500mb
!  RTS= 13.0000,  13.0000,  13.0000,   13.0000,   13.0000,  
!       13.0000,  12.1701,  11.0391,    8.8,    6.2,   
!        4.7,   3.7,   3.0,    2.5,    2.2,   
!        2.15,   2.1,   2.0,    1.95, 1.93,   
!        1.91, 1.9, 1.87, 1.76,    1.6236,   

! 13 g/kg moisture well-mixed
  RTS= 13.0000,  13.0000,  13.0000,   13.0000,   13.0000,  
       13.0000,  12.1701,  11.0391,    9.9877,    9.0267,   
        8.1485,   7.3462,   6.6135,    5.9450,    5.3355,   
        4.7804,   4.2755,   3.8167,    3.4005,    3.0234,   
        2.6824,   2.3746,   2.0971,    1.8476,    1.6236,   
        1.4229,   1.2437,   1.0838,    0.9417,    0.8156,   
        0.7041,   0.6058,   0.5193,    0.4435,    0.3773,   
        0.3196,   0.2696,   0.2264,    0.1892,    0.1573,   
        0.1301,   0.1070,   0.0875,    0.0710,    0.0573,   
        0.0459,   0.0364,   0.0254,    0.0138,    0.0092,   
        0.0095,   0.0099,   0.0103,    0.0107,    0.0112,   
        0.0116,   0.0121,   0.0126,    0.0131,    0.0136,   
        0.0141,   0.0147,   0.0153,    0.0159,    0.0166,   
        0.0172,   0.0179,   0.0186,    0.0194,    0.0202,   
        0.0210,   0.0218,   0.0227,    0.0236,    0.0245,   
        0.0255,   0.0265,   0.0276,    0.0287,    0.0299,   
        0.0311,   0.0323,   0.0336,    0.0350,    0.0360,    
        0.0370,   0.0380,   0.0390,    0.0400,    0.0410,    
        0.0420,   0.0430,   0.0450,    0.0460,    0.0470,    
        0.0480,   0.0490,   0.0500,    0.0500,    0.0500,

! origonal sounding "u" (changed sfc value to be 0)
  US=  0.0000,  0.19590,   0.7790,   1.7347,   3.0396,  
       4.6615,  6.56060,   8.6900,  10.9973,  13.4258, 
      15.9155,  18.4052,  20.8337,  23.1410,  25.2704, 
      27.1695,  28.7914,  30.0963,  31.0520,  31.6350,  
      32.0310,  32.2310,  32.4310,  32.6310,  32.8310,  
      33.0310,  33.2310,  33.4310,  33.6310,  33.8310,
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310, 
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,   
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,  
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,  
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,  
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310, 
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,  
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,  
      33.8310,  33.8310,  33.8310,  33.8310,  33.8310,  
      33.8310,  33.8310,  33.8310,  33.8310,  25.0000,
      20.0000,  15.0000,  10.0000,   5.0000,   0.0000,
       0.0000,   0.0000,   0.0000,   0.0000,   0.0000,
       0.0000,   0.0000,   0.0000,   0.0000,   0.0000,

! weaker hodograph: v halved in lowest 4.3 km (18 levs), reduced lev 19
!  VS=  0.0000,   1.2597,   2.5082,   3.6255,   4.6549,
!       5.6540,   6.4759,   7.0808,   7.6365,   7.8196,
!       8.0155,   7.8196,   7.6365,   7.0808,   6.4759,
!       5.6540,   4.1549,   3.6255,   3.0182,   2.4897,

! origonal sounding "v" (changed sfc value to be 0)
  VS=  0.0000,   2.4897,   4.9182,   7.2255,   9.3549,
      11.2540,  12.8759,  14.1808,  15.1365,  15.7196,
      15.9155,  15.7196,  15.1365,  14.1808,  12.8759,
      11.2540,   9.3549,   7.2255,   4.9182,   2.4897,
       0.5000,   1.0000,   1.5000,   2.0000,   2.5000,
       3.0000,   3.5000,   4.0000,   4.5000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   5.0000,
       5.0000,   5.0000,   5.0000,   5.0000,   4.0000,
       3.0000,   2.0000,   1.0000,   0.0000,   0.0000,
       0.0000,   0.0000,   0.0000,   0.0000,   0.0000,
       0.0000,   0.0000,   0.0000,   0.0000,   0.0000,

 $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   = 'RC','RC','TOTPRE',
                            ! Field names - see table below
                            
!   PLFMT(1)  = '0PF7.3',   ! Format spec. if default is unacceptable

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

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

! File names and dispose flags 

   VARPFX    = '/home/smsaleeb/data.07Feb11/rams6/varf/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

