!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 = 18., ! Final time of simulation ! Start of simulation or ISAN processing IMONTH1 = 04, ! Month IDATE1 = 26, ! Day IYEAR1 = 2008, ! Year ITIME1 = 1800, ! GMT of model TIME = 0. ! Grid specifications NGRIDS = 1, ! Number of grids to run NNXP = 256,100, ! Number of x gridpoints NNYP = 256,100, ! Number of y gridpoints NNZP = 150,35, ! Number of z gridpoints NZG = 3, ! 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 = 50., ! X-Y grid spacing (meters) DELTAZ = 10., ! Z grid spacing (set to 0. to use ZZ) DZRAT = 1.00, ! Vertical grid stretch ratio DZMAX = 10., ! 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-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 = 71.32, ! If polar, latitude/longitude of pole point POLELON = -156.61, ! 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 = 71.32,37.0, ! Center lat/lon of grids, may or CENTLON = -156.61,-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 ! ------ 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 = './out/z.atex.ccn0400.sst298-A-1999-02-11-131500-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 = './output_100_101_10/a', ! File prefix for all analysis files ! state,mean,lite,both !------------------------------------------------------------------------ FRQSTATE = 900.,10800.,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 = './surface/toph', ! Path and prefix for topo files. SFCFILES = './surface/sfch', ! Path and prefix for surface files. SSTFPFX = './surface/ssth', ! Path and prefix for sst files NDVIFPFX = './surface/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, ! 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,1,1,1, ! 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.isdac', ! 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 = 30, ! 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 = 3., ! 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 ! 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 = .00015, ! Constant roughness if for all domain ALBEDO = .5, ! Constant albedo if not running soil model SEATMP = 267., ! 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 = -.3, -.25,-.20,-.16,-.12,-.09,-.06,-.03, ! Soil grid levels SLMSTR = 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, ! Initial soil moisture STGOFF = 5., 5., 3.5, 2., .5, -1.,-1.5,-1.8, ! Initial soil temperature offset ! from lowest atmospheric level !-------------------------------------------------------------------------- DIVLS = 0.0e-6, !--------- Eddy diffusion coefficient parameters ----------------------- 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,1.,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) IBUBGRD = 1, ! Grid Number On Which To Initialize Bubble IBDXIA = 25, ! DX start grid point IBDXIZ = 30, ! DX end grid point IBDYJA = 25, ! DY start grid point IBDYJZ = 30, ! DY end grid point IBDZK1 = 2, ! Lower k level IBDZK2 = 60, ! 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 = 21, ! 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 !----- 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 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 = 0, !------------------- IRAIN = 5, ! 1 - diagnostic concen. IPRIS = 0, ! 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 = 10e-6, ! Characteristic diameter, # concentration DPARM = .1e6, ! or y-intercept RPARM = 1e-3, ! PPARM = 0., ! Specify values for the given SPARM = 1e3, ! 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., ! Aerosol parameters ! See mic_init.f90 for setting profiles ! ------------------------------------------ IAEROSOL = 1, ! CCN,GCCN: 0=off, 1=userset ISALT = 0, ! Sea Salt: 0=off, 1=userset, 2=saltmodel IDUST = 0, ! Dust: 0=off, 1=userset 2=dustmodel ICCNLEV = 0, ! Aerosol treatment level ! 0 = no source or sink ! 1 = nucleation + deposition sink ! 2 = same as 1 + evaporation regeneration IFORCECCN = 0, ! Force aerosol values ! 0 = no forcing ! 1 = BL decrease, FT constant ! 2 = FT decrease, BL constant ! 3 = Decrease in both BL and FT BLH = 825., !Definition of boundary layer height (m) !for FORCECCN FCCNTS = 0., !Decay timescale for FORCECCN. 0 = no decay. FCCNSTART = 21600., !Time after which to start decay !Before FCCNSTART, aerosol forced to CCN_MAX IIFN = 2, ! 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 = 0, ! Print aerosol init profiles (0=no,1=yes) IAEROHIST = 0, ! Re-initialize aerosols on hist restart:0,1 CIN_MAX = .01, ! Ice Nuclei (#/mg) CCN_MAX = 100, ! CCN (#/cc - RAMS code must convert to #/kg) GCCN_MAX = 0.e-5, ! Giant-CCN (#/mg) DUST1_MAX = 0., ! Small Dust Mode (#/mg) DUST2_MAX = 0., ! Large Dust Mode (#/mg) SALTF_MAX = 0., ! Film Sea Salt (#/mg) SALTJ_MAX = 0., ! Jet Sea Salt (#/mg) SALTS_MAX = 0.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.50,0.95,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,1.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 = 0, ! Tracking aerosol solubility for regeneration ITRKDUST = 0, ! Tracking dust sub-species within hydrometeors ITRKDUSTIFN = 0, ! Tracing dust that served as ice nuclei !Mechanism denial flags ISEDTEST = 1, IEVAPTEST = 0, IRADTEST = 1, ISATADJ = 0, ICLDSED = 1, $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., PS = 1020,10,20,30,40,50,60,70,80, 90,100,110,120,130,140,150,160,170, 180,190,200,210,220,230,240,250,260, 270,280,290,300,310,320,330,340,350, 360,370,380,390,400,410,420,430,440, 450,460,470,480,490,500,510,520,530, 540,550,560,570,580,590,600,610,620, 630,640,650,660,670,680,690,700,710, 720,730,740,750,760,770,780,790,800, 810,820,830,840,850,860,870,880,890, 900,910,920,930,940,950,960,970,980, 990,1000,1010,1020,1030,1040,1050,1060,1070, 1080,1090,1100,1110,1120,1130,1140,1150,1160, 1170,1180,1190,1200,1210,1220,1230,1240,1250, 1260,1270,1280,1290,1300,1310,1320,1330,1340, 1350,1360,1370,1380,1390,1400,1410,1420,1430, 1440,1450,1460,1470,1480,1490,1500,1510,1520, TS = 263.4,263.44,263.48,263.52,263.56,263.6,263.64,263.68,263.72, 263.76,263.8,263.84,263.88,263.92,263.96,264,264.04,264.08, 264.12,264.16,264.2,264.24,264.28,264.32,264.36,264.4,264.44, 264.48,264.52,264.56,264.6,264.64,264.68,264.72,264.76,264.8, 264.84,264.88,264.92,264.96,265,265,265,265,265, 265,265,265,265,265,265,265,265,265, 265,265,265,265,265,265,265,265,265, 265,265,265,265,265,265,265,265,265, 265,265,265,265,265,265,265,265,265, 265,265,267.62,268.25,268.63,268.91,269.13,269.33,269.5, 269.65,269.79,269.92,270.04,270.15,270.26,270.36,270.45,270.54, 270.63,270.71,270.79,270.86,270.94,271.01,271.08,271.14,271.21, 271.27,271.33,271.39,271.45,271.51,271.56,271.62,271.67,271.72, 271.77,271.82,271.87,271.92,271.97,272.01,272.06,272.1,272.14, 272.19,272.23,272.27,272.31,272.35,272.39,272.43,272.47,272.51, 272.55,272.58,272.62,272.66,272.69,272.73,272.76,272.8,272.83, 272.87,272.9,272.93,272.96,273,273.03,273.06,273.09,273.12, RTS = 1.8,1.7925,1.785,1.7775,1.77,1.7625,1.755,1.7475,1.74, 1.7325,1.725,1.7175,1.71,1.7025,1.695,1.6875,1.68,1.6725, 1.665,1.6575,1.65,1.6425,1.635,1.6275,1.62,1.6125,1.605, 1.5975,1.59,1.5825,1.575,1.5675,1.56,1.5525,1.545,1.5375, 1.53,1.5225,1.515,1.5075,1.5,1.5,1.5,1.5,1.5, 1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5, 1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5, 1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5, 1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5, 1.5,1.5,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, 1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2,1.2, US = -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, -7,-7,-7,-7,-7,-7,-7,-7,-7, VS = -2,-1.97,-1.94,-1.91,-1.88,-1.85,-1.82,-1.79,-1.76, -1.73,-1.7,-1.67,-1.64,-1.61,-1.58,-1.55,-1.52,-1.49, -1.46,-1.43,-1.4,-1.37,-1.34,-1.31,-1.28,-1.25,-1.22, -1.19,-1.16,-1.13,-1.1,-1.07,-1.04,-1.01,-0.98,-0.95, -0.92,-0.89,-0.86,-0.83,-0.8,-0.77,-0.74,-0.71,-0.68, -0.65,-0.62,-0.59,-0.56,-0.53,-0.5,-0.47,-0.44,-0.41, -0.38,-0.35,-0.32,-0.29,-0.26,-0.23,-0.2,-0.17,-0.14, -0.11,-0.08,-0.05,-0.02,0.01,0.04,0.07,0.1,0.13, 0.16,0.19,0.22,0.25,0.28,0.31,0.34,0.37,0.4, 0.43,0.46,0.49,0.52,0.55,0.58,0.61,0.64,0.67, 0.7,0.73,0.76,0.79,0.82,0.85,0.88,0.91,0.94, 0.97,1,1.03,1.06,1.09,1.12,1.15,1.18,1.21, 1.24,1.27,1.3,1.33,1.36,1.39,1.42,1.45,1.48, 1.51,1.54,1.57,1.6,1.63,1.66,1.69,1.72,1.75, 1.78,1.81,1.84,1.87,1.9,1.93,1.96,1.99,2.02, 2.05,2.08,2.11,2.14,2.17,2.2,2.23,2.26,2.29, VS = -4.25, -4.25, -4.25, -4.25, -4.25, -4.25, -4.25, -4.25 $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