!==========================================================================================! !==========================================================================================! ! RAMSIN. ! ! ! ! This is the file that contains the various namelists that define how BRAMS is to be ! ! run. ! ! ! ! Nota: Todos os comentarios foram retirados deste RAMSIN para evitar problemas na ! ! operacao. Para uma versao mais bem explicada do RAMSIN, procure o arquivo no ! ! diretorio BRAMS406/BRAMS/run/RAMSIN. ! !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls most of the grid and domain parameters. ! !------------------------------------------------------------------------------------------! $MODEL_GRIDS EXPNME = 'Previsao operacional eb00', RUNTYPE = 'INITIAL', ! Type of run: MAKESFC, INITIAL, HISTORY, ! HISTORY, or MAKEVFILE LOAD_BAL = 0, ! Dynamic load balance flag: 1=yes, 0=no ! Start of simulation or ISAN processing IMONTHA = 07, IDATEA = 12, IYEARA = 2016, ITIMEA = 0000, ! End of simulation or ISAN processing IMONTHZ = 07, IDATEZ = 15, IYEARZ = 2016, ITIMEZ = 0000, ! Grid specifications NGRIDS = 1, ! Number of grids to run NNXP = 89, ! Number of x gridpoints NNYP = 89, ! Number of y gridpoints NNZP = 50, ! Number of z gridpoints NZG = 12, ! 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, !Vertical coordinate: 0- sigma_z, 1- shaved ETA IHTRAN = 1, ! 0-Cartesian, 1-Polar stereo DELTAX = 25000., DELTAY = 25000., ! X and Y grid spacing DELTAZ = 120., ! Z grid spacing (set to 0. to use ZZ) DZRAT = 1.05, ! Vertical grid stretch ratio DZMAX = 1000., ! Maximum delta Z for vertical stretc ! Vertical levels if DELTAZ = 0 ZZ = 0.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.0, 11700.0, 12700.0, 13700.0, 14700.0, 15700.0, 16700.0, 17700.0, 18700.0, 19700.0, DTLONG = 30., ! 25 - Coarse grid long timestep NACOUST = 4, ! 3 - Small timestep ratio IDELTAT = 0, ! =0 - constant timesteps ! >0 - initial computation <0 - variable ! Nest ratios between this grid ! and the next coarser grid. NSTRATX = 1,3, ! X direction NSTRATY = 1,3, ! Y direction NNDTRAT = 1,3, ! Time NESTZ1 = 0, ! Contort coarser grids if negative NSTRATZ1 = 3,3,3,3,2,2,2,2,1, ! NESTZ2 = 0, ! Contort coarser grids if negative NSTRATZ2 = 3,3,3,2,2,2,2,1, ! POLELAT = -25.00, ! Latitude of pole point POLELON = -49.00, ! Longitude of pole point CENTLAT = -25.00, ! Latitude of grid centre CENTLON = -49.00, ! Longitude of grid centre ! Grid point on the next coarser ! nest where the lower southwest ! corner of this nest will start. ! If NINEST or NJNEST = 0, use CENTLAT/LON NINEST = 1,0, ! i-point NJNEST = 1,0, ! j-point NKNEST = 1,1, ! k-point NNSTTOP = 1,1, ! Flag (0-no or 1-yes) if this NNSTBOT = 1,1, ! Nest goes the top or bottom of the ! coarsest nest. GRIDU = 0.,0., ! u-component for moving grids GRIDV = 0.,0., ! v-component for moving grids ! (not working again!) $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls parameters relative to the Coupled Aerosol Tracer and ! ! Transport sub-model (CATT). ! !------------------------------------------------------------------------------------------! $CATT_INFO CATT = 0, FIREMAPFN = '/gfs/dk15/fazenda/data/Queima_source', RECYCLE_TRACERS = 0, PLUMERISE = 0, PRFRQ = 3600., $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls parameters relative to the TEB urban model. ! !------------------------------------------------------------------------------------------! $TEB_SPM_INFO TEB_SPM = 0, $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls which initial and boundary atmospheric conditions will be ! ! used by BRAMS, and how they will be used. ! !------------------------------------------------------------------------------------------! $MODEL_FILE_INFO INITIAL = 2, ! Initial fields - 1=horiz.homogeneous, ! 2=init from varfile ! 3=init from history file ! ------ Analysis nudging parameters ----------------------------------- NUD_TYPE = 2, ! =1 - nudge from history files(1-way nest) ! =2 - nudge from varfiles ! =0 - no analysis nudging VARFPFX = '/p1-sonny3/marcos/opera_406/model_qe00/isean/qe00', VWAIT1 = 0., ! wait between each VFILE check (s) VWAITTOT = 0., ! total wait before giving up on a VFILE (s) NUD_HFILE = '/p1-sonny3/marcos/opera_406/model_eb00/histo/h_eb00', ! Prefix for history nudging files NUDLAT = 8 , ! Number of points in lateral bnd region TNUDLAT = 900., ! Nudging time scale(s) at lateral boundary TNUDCENT = 86400., ! Nudging time scale(s) in center of domain TNUDTOP = 10800., ! Nudging time scale (s) at top of domain ZNUDTOP = 18000., ! Nudging at top of domain above height(m) WT_NUDGE_GRID = 1., 0.80, ! 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 (this is awful!) weight for u and v WT_NUDGE_TH = 1., ! Anal nudging weight for theta WT_NUDGE_PI = 4., ! Anal nudging weight for pi WT_NUDGE_RT = 1., ! Anal nudging weight for r_tot WT_NUDGE_CO2 = 1., ! Anal nudging weight for CO2 ! Lower boundary nudging. Nudging intensity at the inner ! part of the domain will decrease quadratically below ! ZB_NUDGE_??. In case ZB_NUDGE_?? = 0., no quadratic ! attenuation will occur. ZB_NUDGE_UV = 0., ! Lower bnd for u and v ZB_NUDGE_TH = 0., ! Lower bnd for theta ZB_NUDGE_PI = 0., ! Lower bnd for pi ZB_NUDGE_RT = 0., ! Lower bnd for r_tot ZB_NUDGE_CO2 = 0., ! Lower bnd for CO2 !------------------------------------------------------------------------- !----------- Condensate nudging ------------------------------------------ NUD_COND = 0, ! Only nudge total water where condensate ! exists (from previous history files) COND_HFILE = './histo/opera_eb00-H-2001-07-21-000000-head.txt', ! Header file name for cond nudging history files (only prefix is used) TCOND_BEG = 0., ! Model time start of cond nudging (sec) TCOND_END = 21600., ! Model time end of cond nudging (sec) WT_NUDGE_RC = 0.01, ! Cond nudging timescale for r_total ZB_NUDGE_RC = 0., ! Lower bnd for r_cond WT_NUDGEC_GRID = 1., 0.8, ! 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 = '/n/Moorcroft_Lab/Lab/data/brams_data/racci/radio/ir', ! File prefix for upper air obs ODA_SFCPREFIX = '/n/Moorcroft_Lab/Lab/data/brams_data/racci/surf/quali/eq', ! File prefix for surface obs FRQODA=21600., ! Frequency of obs analysis TODABEG=0., ! Model time start of oda (sec) TODAEND=99999999., ! Model time end of oda (sec) TNUDODA= 3600., ! Nudging timescale for each grid WT_ODA_GRID = 1., 0.75, ! 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 = 72000., 36000., RODA_SFC0 = 200000., 216000., RODA_UPAE = 150000., 75000., RODA_UPA0 = 450000., 225000., RODA_HGT = 3000.,3000., ! Height at which transition from SFC radii ! to UPA radii occurs RODA_ZFACT = 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=1200., ODA_UPA_TIL=43200., ODA_UPA_TEL=3600., !------------------------------------------------------------------------ !----- Cumulus inversion tendency input ----- IF_CUINV = 0, CU_PREFIX = './t5-C-', TNUDCU=900., WT_CU_GRID=1., 1., .5, TCU_BEG=0., TCU_END=7200., CU_TEL=3600., CU_TIL=21600., !------------------------------------------------------------------------ ! History time input IMONTHH = 10, IDATEH = 25, IYEARH = 2008, ITIMEH = 1200, HFILIN = '/p1-sonny3/marcos/opera_406/model_eb00/histo/h_eb00', ! Input history file name ! Analysis file input for assimilation (currently LEAF variables) IPASTIN = 0, ! Initialize various fields from analysis file? ! 1=yes, 0=no PASTFN = 'a-A-2000-01-09-000000-head.txt', ! Input analysis file name ! History/analysis file output IOUTPUT = 2, ! 0-no files, 1-save ASCII, 2-save binary HFILOUT = '/p1-sonny3/marcos/opera_406/model_qe00/histo/qe00', AFILOUT = '/p1-sonny3/marcos/opera_406/model_qe00/analy/qe00', ICLOBBER = 1, ! 0=stop if files exist, 1=overwite files IHISTDEL = 0, ! 0=keep all hist files, 1=delete previous FRQHIS = 43200., ! History file frequency FRQANL = 3600., ! Analysis file frequency FRQLITE = 0., ! Analysis freq. for "lite" variables ! = 0 : no lite files XLITE = '/0:0/', ! nums>0 are absolute grid indices YLITE = '/0:0/', ! nums<0 count in from the domain edges ZLITE = '/0:0/', ! nums=0 are domain edges NLITE_VARS = 5, LITE_VARS = 'GLAT', 'GLON', 'UP', 'VP', 'WP', !----------------------------------------------------------------------- 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., ! Analysis freq. for "averaged" variables FRQBOTH = 0., ! Analysis freq. for Both "averaged" and ! "lite" variables KWRITE = 0, ! 1-write,0-don't write scalar K's to anal. ! Printed output controls FRQPRT = 21600., ! Printout frequency INITFLD = 1, ! Initial field print flag 0=no prnt,1=prnt ! Input topography variables TOPFILES = '/p1-sonny3/marcos/opera_406/model_qe00/surfa/topo_qe00', SFCFILES = '/p1-sonny3/marcos/opera_406/model_qe00/surfa/surfa_qe00', SSTFPFX = '/p1-sonny3/marcos/opera_406/model_qe00/surfa/sst_qe00', NDVIFPFX = '/p1-sonny3/marcos/opera_406/model_qe00/surfa/ndvi_qe00', ITOPTFLG = 1,1,1, ! 2 - Fill data in "leaf3_init" ISSTFLG = 1,1,1, ! 0 - Interpolate from coarser grid IVEGTFLG = 1,1,1, ! 1 - Read from standard Lat/Lon data file ISOILFLG = 1,1,1, ! NDVIFLG = 1,1,1, ! NOFILFLG = 2,0,0, ! 2 - Fill data in "leaf3_init" ! 0 - Interpolate from coarser grid IUPDNDVI = 0, ! 0 - No update of NDVI values during run IUPDSST = 1, ! 0 - No update of SST values during run ! 1 - Update values during run IUSELAI = 1, ! 0 - Don't substitute LAI data for the NDVI data ! 1 - The NDVI data files are actually LAI data ! The following only apply for IxxxxFLG=1 ITOPTFN = '/p1-sonny2/ramsop/dados_brams/topo/EL', ISSTFN = '/p1-sonny2/ramsop/dados_brams/sst_week/RAMS/W', IVEGTFN = '/p1-sonny3/marcos/dados_brams/veg_MCD12Q1/2012/MCD12Q1LU1_2012_', ISOILFN = '/p1-sonny3/marcos/dados_brams/FAO/FAO_', NDVIFN = '/p1-sonny3/marcos/dados_brams/lai-modis/LAI', ! Topography scheme ITOPSFLG = 0,0, ! 0 = Average Orography ! 1 = Silhouette Orography ! 2 = Envelope Orography ! 3 = Reflected Envelope Orography TOPTENH = 0.,0., ! For ITOPSFLG=1, Weighting of topo ! silhouette averaging ! For ITOPSFLG=2 or 3, Reflected Envelope ! and Envelope Orography enhancement factor TOPTWVL = 3.,2.,2.,4., ! Topo wavelength cutoff in filter ! Surface Roughness scheme IZ0FLG = 0,0, ! 0 = Use veg, bare soil and water surface ! 1 = Also use subgrid scale topography Z0MAX = 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 = '/p1-sonny3/marcos/opera_406/model_qe00/micro/coltab_qe00', ! Filename to read or write $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls the variables used by the cumulus parametrisation. ! !------------------------------------------------------------------------------------------! $CUPARM_OPTIONS !---------------------------------------------------------------------------------------! ! NNQPARM(ngrids) = Convective parameterization flag: 0 - off, 1 - on. ! !---------------------------------------------------------------------------------------! NNQPARM = 1,1, !---------------------------------------------------------------------------------------! ! NCLOUDS: How many clouds do you want to use? ! ! PS1: If you want to run deep convection only, use nclouds=1 ! ! PS2: If you want to run shallow convection only for a certain grid ! ! (or all grids...) set up nclouds to 2 and ndeepest to 0 ! ! PS3: If NCLOUDS >= 3, then the middle ones will necessarily use Grell's ! ! parameterization. ! !---------------------------------------------------------------------------------------! NCLOUDS = 2, !---------------------------------------------------------------------------------------! ! Both define the cumulus scheme to be used for the deepest and shallowest clouds. ! ! NDEEPEST(ngrids) : 0 - turn it off; 1- run Kuo ; 2- run Grell; 3- run old Grell ! ! NSHALLOWEST(ngrids) : 0 - turn it off; 1- run Souza; 2- run Grell 3- run old Grell ! !---------------------------------------------------------------------------------------! NDEEPEST = 2,2, NSHALLOWEST = 2,2, !----- Kuo configuration (will be ignored by Souza and Grell's methods) ----------------! WCLDBS = 0.0, !---------------------------------------------------------------------------------------! ! CONFRQ(nclouds): how often should each cloud size be updated? ! ! CPTIME(nclouds): when I should start computing this cloud. ! PS1: This is NOT for each grid, but for each cloud!!! ! ! PS2: The first value is always assigned for the deepest cloud, even ! ! if you are not running it. ! !---------------------------------------------------------------------------------------! CONFRQ = 1200., CPTIME = 1200., !---------------------------------------------------------------------------------------! !---------------------------------------------------------------------------------------! ! IUPMETHOD -- How to find the level in which updrafts originate: ! ! 1. Maximum moist static energy; ! ! 2. PBL top (if IDIFFK=7, uses the pblhgt from turbulence) ! ! 3. Most turbulent layer ! ! 4. Max w+sigw within the PBL. ! ! 5. Start from the bottom, move upwards (5 is available only if ! ! CAP_MAXS < 0) ! ! ILOUREC -- Should the search for level of origin of updrafts (LOU) be recursive? ! ! 0. No ! ! 1. Yes. ! !---------------------------------------------------------------------------------------! IUPMETHOD = 1, ILOUREC = 0, !---------------------------------------------------------------------------------------! !---------------------------------------------------------------------------------------! ! The following variables have (nclouds) dimension. ! ! ! ! RADIUS - Cloud mean radius [m]. Important: the cloud won't have downdrafts nor ! ! precipitation if its radius is less than min_down_radius ! ! (grell_coms.f90). Default min_down_radius is 900. ! ! DEPTH_MIN - Minimum depth that the cloud must have [m]. ! !---------------------------------------------------------------------------------------! RADIUS = 12000., 600., DEPTH_MIN = 4000., 100., !---------------------------------------------------------------------------------------! !---------------------------------------------------------------------------------------! ! CAP_MAXS controls whether convection can happen or not. Values can be negative or ! ! positive (but never 0). It is a scalar. ! ! Positive: Depth of the inverse capping in hPa. The larger the value is, the less ! ! likely is to develop convection. ! ! Negative: The minimum probability [%] of the air parcel leaving from the level of ! ! origin of the updrafts to reach the level of free convection (if its radius ! ! is greater than min_down_radius) or the minimum depth (if its radius is ! ! less than min_down_radius). The probability is found by using the mean ! ! velocity and sigma-w and assuming normal distribution, so this option works ! ! only when IDIFFK is 1 or 7. ! !---------------------------------------------------------------------------------------! CAP_MAXS = -20., !---------------------------------------------------------------------------------------! !---------------------------------------------------------------------------------------! ! The following variables are scalars. ! !---------------------------------------------------------------------------------------! CLD2PREC = -0.30, ! Ratio of conversion of condensates to precipitation. [ ---] ZKBMAX = 4000., ! Maximum AGL height in which updrafts can originate. [ m] ZCUTDOWN = 3000., ! Maximum AGL height in which downdrafts can originate. [ m] Z_DETR = 1250., ! Top of the downdraft detrainment layer [ m] MAX_HEAT = 600., ! Maximum heating rate allowed for feedback [ K/day] !---------------------------------------------------------------------------------------! !=======================================================================================! ! The following variables will define the ensemble members. Currently 4 types of ! ! ensemble are available. Odd numbers are preferred over even numbers since everything ! ! is centred around the reference value (except for CLOSURE_TYPE, which is done ! ! differently). These are all (cloud) dependent. ! ! PS1: Use with moderation, these members are multiplicative and can get slow down your ! ! run. ! ! PS2: The perturbation increments are defined inside the code, feel free to take it ! ! out... ! !=======================================================================================! CLOSURE_TYPE = 'KF', ! Closure type to find the reference updraft mass flux: ! PS1. Case insensitive - see table below: !----------------------------------------------------------------------------------! ! Type | Description | Deep | Shal | Reference Paper ! !------+--------------------------------+------+------+----------------------------! ! EN | Ensemble (all closures) | yes | no | Grell and Dévényi, 2002 ! ! NC | Ensemble, excluding LO and MC | yes | yes | Grell and Dévényi, 2002 ! ! GR | Modified quasi-equilibrium ABE | yes | yes | Grell, 1993 ! ! LO | Low level omega | yes | no | Frank and Cohen, 1987 ! ! MC | moisture convergence | yes | no | Krishnamurti et al., 1983 ! ! KF | Instability removal | yes | yes | Kain and Fritsch, 1990 ! ! AS | Quasi equilibrium of ABE | yes | yes | Arakawa and Schubert, 1974 ! !----------------------------------------------------------------------------------! MAXENS_LSF = 1, ! Number of different arbitrary mass fluxes for GR/AS/KF dynamic ! control (it will also create some different numbers for MC) MAXENS_EFF = 1, ! Number of precipitation efficiency members MAXENS_CAP = 1, ! Number of static control members, currently changes in cap_maxs. $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist will allow you to set up and choose several of the BRAMS sub-models. ! !------------------------------------------------------------------------------------------! $MODEL_OPTIONS NADDSC = 0, ! Number of additional scalar species ICORFLG = 1, ! Coriolis flag/2D v-component - 0=off, 1=on IADVEC = 2, ! Advection tendency ! 1- Original ! 2- Freitas et al 2013 IEXEV = 2, ! Exner tendency ! 1- simplest form ! 2- complete, mass conservative, solution ! (works only with sigma z coordinate) IMASSFLX = 1, ! Output advective and convective mass flux for ! Lagrangian Particle Dispersion Modelling? ! 0 - no, 1 - yes IBND = 1, ! Lateral boundary condition flags JBND = 1, ! 1-Klemp/Wilhelmson, 2-Klemp/Lilly, 3-Orlanski ! 4-cyclic CPHAS = 20., ! Phase speed if IBND or JBND = 1 LSFLG = 0, ! Large-scale gradient flag for variables other than ! normal velocity: ! 0 = zero gradient inflow and outflow ! 1 = zero gradient inflow, radiative b.c. outflow ! 2 = constant inflow, radiative b.c. outflow ! 3 = constant inflow and outflow NFPT = 0, ! Rayleigh friction - number of points from the top DISTIM = 200., ! - dissipation time scale ! Radiation parameters ISWRTYP = 4, ! Shortwave radiation type ILWRTYP = 4, ! Longwave radiation type ! 0-none, 1-Chen/Cotton, 2-Mahrer/Pielke ! 3-Harrington, 4-Carma ! ICUMFDBK = 1, ! Add effects of parameterized cumulus in radiation? ! 0 - No , 1 - Yes ! (Used for i?wrtyp = 3 and 4 only) ! ! RADFRQ = 600., ! Freq. of radiation tendency update (s) LONRAD = 1, ! Longitudinal variation of shortwave ! (0-no, 1-yes) ! Surface layer and soil parameterization NPATCH = 3, ! Number of patches per grid cell (min=2) NVEGPAT = 2, ! Number of patches per grid cell to be ! filled from ! vegetation files ! (min of 1, max of NPATCH-1) MIN_PATCH_AREA = 0.01, ! Minimum area for land patches to be accounted ISFCL = 4, ! Surface layer/soil/veg model ! 0 - specified surface layer gradients ! 1 - LEAF-3 submodel ! 2 - LEAF-3/TOPMODEL submodel ! 3 - SIB submodel (not available) ! 4 - "LEAF-4" model. Same as option 0, but vegetation ! radiation is solved by two-stream, and transpiration ! derived from Farquhar model. ! 5 - ED-2 (disabled) DTLEAF = 10., ! 6 1/4 Nominal LEAF-3 time step. This should be positive, and ! less than 30 seconds. If DTLONG is less than DTLEAF, ! then LEAF-3 will use DTLONG instead. NDTVEG = 5, ! Time nesting for vegetation properties. ISTAR = 4, ! Surface layer model (used to compute ustar, tstar,...) ! 1. BRAMS default, based on Louis (1979). ! 2. Oncley and Dudhia (1995) model, based on MM5. ! 3. Beljaars and Holtslag (1991) model. ! 4. CLM-4 (2004). Special function for very unstable and ! very stable cases. IGRNDVAP = 0, ! Ground->canopy conductance. ! 0. Modified Lee Pielke (1992), adding field capacity, ! but using beta factor without the square, like in ! Noilhan and Planton (1989). ! 1. Test # 1 of Mahfouf and Noilhan (1991) ! 2. Test # 2 of Mahfouf and Noilhan (1991) ! 3. Test # 3 of Mahfouf and Noilhan (1991) ! 4. Test # 4 of Mahfouf and Noilhan (1991) ! 5. Combination of test #1 (alpha) and test #2 (soil resistance). UBMIN = 0.25, ! minimum wind speed at the top of the canopy air space UGBMIN = 0.05, ! minimum wind speed at the leaf level USTMIN = 0.01, ! 0.07 -- minimum friction velocity, u*, in m/s. GAMM = 15., ! gamma coefficient for momentum, unstable case (dimensionless) ! Ignored when ISTAR = 1 GAMH = 9., ! gamma coefficient for heat, unstable case (dimensionless) ! Ignored when ISTAR = 1 TPRANDTL = 0.74, ! Turbulent Prandtl number ! Ignored when ISTAR = 1 RIBMAX = 6.00, ! 0.3333333333 -- maximum bulk Richardson number. LEAF_MAXWHC = 0.10, ! Max. water that can be intercepted by leaves, in kg/m2leaf. ICO2 = 1, ! CO2 prognostic level ! 0 - Do not solve it (CO2CON(1) will be used as a time ! invariant CO2 mixing ratio whenever needed); ! 1 - Solve it initialising with a constant value [CO2CON(1)] ! 2 - Solve and initialise with a horizontally homogeneous ! profile, (filling CO2CON with all values from bottom ! to top). ! 3 - Solve and initialise with a RALPH2/SOUND_IN file ! (not available yet). CO2CON = 395., ! Initial constant CO2 value or profile if for all domain. ! If ICO2 = 0 or 1, only the first value will be considered, ! and if ICO2=2, then provide one value per level. ! Unit here is ppm. NVGCON = 7, ! 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 (TEB only) !------------------------------------------------------------------------------- PCTLCON = 1., ! Constant land % if for all domain NSLCON = 16, ! Constant soil type if for all domain ! 1 -- sand | 7 -- silty clay loam | 13 -- bedrock ! ! 2 -- loamy sand | 8 -- clayey loam | 14 -- silt ! ! 3 -- sandy loam | 9 -- sandy clay | 15 -- heavy clay ! ! 4 -- silty loam | 10 -- silty clay | 16 -- clayey sand ! ! 5 -- loam | 11 -- clay | 17 -- clayey silt ! ! 6 -- sandy clay loam | 12 -- peat ! !------------------------------------------------------------------------------------! ! ISOILCOL -- LEAF-3 and ED-2 soil colour classes that the model will use when ! ! ISOILFLG is set to 2. Soil classes are from 1 to 20 (1 = lightest; ! ! 20 = darkest). The values are the same as CLM-4.0. The table is the ! ! albedo for visible and near infra-red. ! !------------------------------------------------------------------------------------! ! ! ! |-----------------------------------------------------------------------| ! ! | | Dry soil | Saturated | | Dry soil | Saturated | ! ! | Class |-------------+-------------| Class +-------------+-------------| ! ! | | VIS | NIR | VIS | NIR | | VIS | NIR | VIS | NIR | ! ! |-------+------+------+------+------+-------+------+------+------+------| ! ! | 1 | 0.36 | 0.61 | 0.25 | 0.50 | 11 | 0.24 | 0.37 | 0.13 | 0.26 | ! ! | 2 | 0.34 | 0.57 | 0.23 | 0.46 | 12 | 0.23 | 0.35 | 0.12 | 0.24 | ! ! | 3 | 0.32 | 0.53 | 0.21 | 0.42 | 13 | 0.22 | 0.33 | 0.11 | 0.22 | ! ! | 4 | 0.31 | 0.51 | 0.20 | 0.40 | 14 | 0.20 | 0.31 | 0.10 | 0.20 | ! ! | 5 | 0.30 | 0.49 | 0.19 | 0.38 | 15 | 0.18 | 0.29 | 0.09 | 0.18 | ! ! | 6 | 0.29 | 0.48 | 0.18 | 0.36 | 16 | 0.16 | 0.27 | 0.08 | 0.16 | ! ! | 7 | 0.28 | 0.45 | 0.17 | 0.34 | 17 | 0.14 | 0.25 | 0.07 | 0.14 | ! ! | 8 | 0.27 | 0.43 | 0.16 | 0.32 | 18 | 0.12 | 0.23 | 0.06 | 0.12 | ! ! | 9 | 0.26 | 0.41 | 0.15 | 0.30 | 19 | 0.10 | 0.21 | 0.05 | 0.10 | ! ! | 10 | 0.25 | 0.39 | 0.14 | 0.28 | 20 | 0.08 | 0.16 | 0.04 | 0.08 | ! ! |-----------------------------------------------------------------------| ! ! ! ! Soil type 21 is a special case in which we use the albedo method that used to be ! ! the default in BRAMS-4.0.6 and ED-2.1. ! !------------------------------------------------------------------------------------! ISOILCOL = 10, ZROUGH = .05, ! Constant roughness if for all domain ALBEDO = .2, ! Constant albedo if not running soil model SEATMP = 298., ! Constant water surface temperature DTHCON = 0., ! Constant sfc layer temp grad for no soil DRTCON = 0., ! Constant sfc layer moist grad for no soil SOIL_MOIST = 'i', ! n => Homogeneous (standard) ! i,h,a => Heterogenous Soil Moist.Init. ! with file. ! i => INITIAL ! h => HISTORY ! a => INITIAL or HISTORY SOIL_MOIST_FAIL = 'l', ! In case that Heterogenous Soil Moist. ! file is not found, what should RAMS do? ! s => STOP the program ! h => Initialize with Homogenous Soil Moist. ! l => Looking for 5 days old files, ! and if not found again then stop. USDATA_IN = '/p1-sonny2/ramsop/dados_brams/soil_moisture_CPTEC/GL_SM.GPNR.', ! Soil Moisture File prefix USMODEL_IN = '/p1-sonny3/marcos/opera_406/model_qe00/surfa/lesmo_qe00', !------------------------------------------------------------------------------------! ! Initial soil moisture is now given as the soil moisture index. Values can be ! ! fraction, in which case they will be linearly interpolated between the special ! ! points (e.g. 0.5 will put soil moisture half way between the wilting point and ! ! field capacity). ! ! -1 = dry air soil moisture ! ! 0 = wilting point ! ! 1 = field capacity ! ! 2 = porosity (saturation) ! !------------------------------------------------------------------------------------! SLZ = -6.000,-4.914,-3.948,-3.100,-2.366,-1.741,-1.222,-0.804,-0.482,-0.249, -0.098,-0.020, ! Soil depth in m SLMSTR = 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, ! Soil moisture index STGOFF = 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, 0.000, ! Soil temperature offset. !------------------------------------------------------------------------------------! !------------------------------------------------------------------------------------! ! ISOILBC -- This controls the soil moisture boundary condition at the bottom. ! ! For regional runs, it is wise to use options 1 or 2, especially for ! ! long runs. ! ! 0. Flat bedrock. Flux from the bottom of the bottommost layer is 0. ! ! 1. Gravitational flow (free drainage). The flux from the bottom of ! ! the bottommost layer is due to gradient of height only. ! ! 2. Lateral drainage. Similar to free drainage, but the gradient is ! ! reduced by the slope not being completely vertical. The reduction ! ! is controlled by variable SLDRAIN. In the future options 0, 1, and ! ! 2 may be combined into a single option. ! ! 3. Aquifer. Soil moisture of the ficticious layer beneath the bottom ! ! is always at saturation. ! !------------------------------------------------------------------------------------! ISOILBC = 1, !------------------------------------------------------------------------------------! !------------------------------------------------------------------------------------! ! SLDRAIN -- This is used only when ISOILBC is set to 2. In this case SLDRAIN is ! ! the equivalent slope that will slow down drainage. If this is set to ! ! zero, then lateral drainage reduces to flat bedrock, and if this is set ! ! to 90, then lateral drainage becomes free drainage. SLDRAIN must be ! ! between 0 and 90. ! !------------------------------------------------------------------------------------! SLDRAIN = 10., !------------------------------------------------------------------------------------! !------------------------------------------------------------------------------------! ! IPERCOL -- This controls percolation and infiltration, and it is used by both ! ! LEAF-3 and ED-2. ! ! 0. Default method. Assumes soil conductivity constant and for the ! ! temporary surface water, it sheds liquid in excess of a 1:9 liquid- ! ! -to-ice ratio through percolation. Temporary surface water exists ! ! only if the top soil layer is at saturation. ! ! 1. Constant soil conductivity, and it uses the percolation model as in ! ! Anderson (1976) NOAA technical report NWS 19. Temporary surface ! ! water may exist after a heavy rain event, even if the soil doesn't ! ! saturate. ! ! 2. Soil conductivity decreases with depth even for constant soil ! ! moisture, otherwise it is the same as 1. ! !------------------------------------------------------------------------------------! IPERCOL = 0, !------------------------------------------------------------------------------------! !------------------------------------------------------------------------------------! ! RUNOFF_TIME -- In case a temporary surface water (TSW) is created, this is the ! ! "e-folding lifetime" of the TSW in seconds due to runoff. If you ! ! don't want runoff to happen, set this to 0. ! !------------------------------------------------------------------------------------! RUNOFF_TIME = 3600.0, !------------------------------------------------------------------------------------! !----- Urban canopy parameterization -------------------------------- IF_URBAN_CANOPY =0, ! Use the urban canopy parameterization (0-no/1-yes)? !------------------------------------------------------------------------- !=======================================================================================! !=======================================================================================! !=======================================================================================! !=======================================================================================! ! This block controls the turbulence and atmospheric boundary layer model. ! !---------------------------------------------------------------------------------------! !--------- Eddy diffusion coefficient parameters ----------------------- IDIFFK = 7,7,7, ! 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,6 - S. Trini Castelli closures ! 7 - Horiz deform/Vert Nakanishi-Niino IBRUVAIS = 1, ! Brunt-Väisälä frequency: ! 1 - Traditional method (g/theta_v d(theta_v)/dz)^˝; ! 2 - Wet enhancement (violates 2nd law of thermo, though). IBOTFLX = 1, ! Horizontal advection flux: ! (this variable will be removed soon) ! 0 - Original code ! 1 - Meester's bug fix IHORGRAD = 1, ! 1 - horiz grad frm decomposed sigma grad ! 2 - true horizontal gradient. ! Non-conserving, but allows small DZ CSX = .32,.32, ! Deformation horiz. K's coefficient CSZ = .35,.35, ! Deformation vert. K's coefficient XKHKM = 3., 3., ! Ratio of horiz K_h to K_m for deformation ZKHKM = 3., 3., ! Ratio of vert K_h to K_m for deformation !---------------------------------------------------------------------------------------! ! The following variables correspond to the "A", "B", and "C" constants for the ! ! Nakanishi and Niino turbulence closure. ! !---------------------------------------------------------------------------------------! NNA = 1.18, 0.665, ! A1 and A2. NNB = 24.0, 15.0 , ! B1 and B2. NNC = 0.137, 0.700, 0.323, 0.00, 0.20, ! C1, C2, C3, C4, and C5. !---------------------------------------------------------------------------! ! The following values can be used to scale the horizontal eddy ! ! viscosity coefficient to the typical value from deformation K (AKSCAL). ! !---------------------------------------------------------------------------! AKMIN = 1.0, 1.0, ! Minimum value for this scaling. AKMAX = 0.0, 0.0, ! Maximum value, for higher elevations. If you ! want AKSCAL to be constant and equal to ! AKMIN, then simply set AKMAX to 0. HGTMIN = 1000.,1000., ! Height above which AKSCAL should start increas- ! ing from AKMIN to AKMAX HGTMAX = 3500.,3500., ! Height above which AKSCAL should remain ! constant and equal to AKMAX. !------------------------------------------------------------------------------ !=======================================================================================! !=======================================================================================! !=======================================================================================! !=======================================================================================! ! The following block controls how BRAMS should solve the water substance, and in ! ! particular the bulk microphysics model. ! !---------------------------------------------------------------------------------------! LEVEL = 3, ! Moisture complexity level ICLOUD = 4, ! Microphysics flags IRAIN = 2, !------------------- IPRIS = 5, ! 1 - diagnostic concen. ISNOW = 2, ! 2 - specified mean diameter IAGGR = 2, ! 3 - specified y-intercept IGRAUP = 2, ! 4 - specified concentration IHAIL = 2, ! 5 - prognostic concentration CPARM = 3.e8, ! cloud droplets, or less frequently, CCN RPARM = 1e-3, ! rain drops PPARM = 0., ! pristine ice (ignored as IPRIS must be 5) SPARM = 1e-3, ! snow APARM = 1e-3, ! aggregates GPARM = 1e-3, ! graupel HPARM = 3e-3, ! hail GNU = 5.,1.,5.,1.,1.,1.,1., ! Gamma shape parms for ! cld rain pris snow aggr graup hail !=======================================================================================! !=======================================================================================! $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls the sounding specification in case INITIAL is set to 0. This ! ! namelist is ignored otherwise. ! !------------------------------------------------------------------------------------------! $MODEL_SOUND ! 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 = 3, ! Specifies what is in RTS array ! 0-dew pnt.(C) 1-dew pnt.(K) ! 2-mix rat(g/kg) ! 3-relative humidity in %, ! 4-dew pnt depression(K) IUSFLG = 1, ! Specifies what is in US and VS arrays ! 0-u,v component(m/s) ! 1-umoms-direction, vmoms-speed HS = 182., PS = 0., 1010., 1000., 925., 850., 700., 600., 500., 400., 300., 200., 100., 70., TS = 25., 18.5, 12.0, 4.5, -24., -11., -37., -56.5, -56.5, -56.5, -56.5, RTS = 70., 70., 70., 70., 20., 20., 20., 20., 10., 10., 10., US = 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, VS = 3., 3., 3., 3., 3., 3., 3., 3., 3., 3., 3., CO2S = 380., 380., 380., 380., 380., 380., 380., 380., 380., 380., 380., $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist is used to define how to plot some of the fields in the standard ! ! output (i.e., on screen). ! !------------------------------------------------------------------------------------------! $MODEL_PRINT NPLT = 0, ! Number of fields printed at each time ! for various cross-sections (limit of 50) IPLFLD = 'UP','THP','THETA','RT','TOTPRE', ! Field names - see table below ! PLFMT(1) = '0PF7.3', ! Format spec. if default is unacceptable IXSCTN = 3,3,3,3,3,3, ! Cross-section type (1=XZ, 2=YZ, 3=XY) ISBVAL = 2,2,2,2,2,2,2, ! Grid-point slab value for third direction ! The following variables can also be set in the namelist: IAA, ! IAB, JOA, JOB, NAAVG, NOAVG, PLTIT, PLCONLO, PLCONHI, and PLCONIN. ! 'UP' - UP(M/S) 'RC' - RC(G/KG) 'PCPT' - TOTPRE ! 'VP' - VP(M/S) 'RR' - RR(G/KG) 'TKE' - TKE ! 'WP' - WP(CM/S) 'RP' - RP(G/KG) 'HSCL' - HL(M) ! 'PP' - PRS(MB) 'RA' - RA(G/KG) 'VSCL' - VL(M) ! 'THP' - THP(K) ! 'THETA'- THETA(K) 'RL' - RL(G/KG) 'TG' - TG (K) ! 'THVP' - THV'(K) 'RI' - RI(G/KG) 'SLM' - SLM (PCT) ! 'TV' - TV(K) 'RCOND'- RD(G/KG) 'CONPR'- CON RATE ! 'RT' - RT(G/KG) 'CP' - NPRIS 'CONP' - CON PCP ! 'RV' - RV(G/KG) 'RTP' - RT'(G/KG) 'CONH' - CON HEAT ! 'CONM' - CON MOIS ! 'THIL' - Theta-il (K) 'TEMP' - temperature (K) ! 'TVP' - Tv' (K) 'THV' - Theta-v (K) ! 'RELHUM'-relative humidity (%) 'SPEED'- wind speed (m/s) ! 'FTHRD'- radiative flux convergence (??) ! 'MICRO'- GASPRC ! 'Z0' - Z0 (M) 'ZI' - ZI (M) 'ZMAT' - ZMAT (M) ! 'USTARL'-USTARL(M/S) 'USTARW'-USTARW(M/S) 'TSTARL'-TSTARL (K) ! 'TSTARW'-TSTARW(K) 'RSTARL'-RSTARL(G/G) 'RSTARW'-RSTARW(G/G) ! 'UW' - UW (M*M/S*S) 'VW' - VW (M*M/S*S) ! 'WFZ' - WFZ (M*M/S*S) 'TFZ' - TFZ (K*M/S) ! 'QFZ' - QFZ (G*M/G*S) 'RLONG'- RLONG ! 'RSHORT'-RSHORT $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls the input of the isentropic analysis. ! !------------------------------------------------------------------------------------------! $ISAN_CONTROL ISZSTAGE = 1, ! Main switches for isentropic-sigz IVRSTAGE = 1, ! "varfile" processing ISAN_INC = 0600, ! ISAN processing increment (hhmm) ! range controlled by TIMMAX, ! IYEAR1,...,ITIME1 GUESS1ST = 'PRESS', ! Type of first guess input- 'PRESS', 'RAMS' I1ST_FLG = 1, ! What to do if first guess file should be ! used but does not exist. ! 1=I know it may not be there, ! skip this data time ! 2=I screwed up, stop the run ! 3=interpolate first guess file from ! nearest surrounding times, stop if unable ! (not yet available) IUPA_FLG = 3, ! UPA-upper air, SFC-surface ISFC_FLG = 3, ! What to do if other data files should be ! uesed, but does not exist. ! 1 = I know it may not be there, ! skip this data time ! 2 = I screwed up, stop the run ! 3 = Try to continue processing anyway ! Input data file prefixes IAPR = '/p1-sonny3/marcos/opera_406/dprep/16071200/dp', IARAWI = '', IASRFCE = '', ! File names and dispose flags VARPFX = '/p1-sonny3/marcos/opera_406/model_qe00/isean/qe00', IOFLGISZ = 0, ! Isen-sigz file flag: 0 = no write, 1 = write IOFLGVAR = 1, ! Var file flag: 0 = no write, 1 = write $END !==========================================================================================! !==========================================================================================! !==========================================================================================! !==========================================================================================! ! This namelist controls the options regarding the execution of the isentropic/ ! ! sigma-z analysis. ! !------------------------------------------------------------------------------------------! $ISAN_ISENTROPIC !----------------------------------- ! Specify isentropic levels !----------------------------------- NISN = 63, ! Number of isentropic levels LEVTH = 248, 250, 252, 254, 256, 258, 260, 262, 264, 266, 268, 270, 272, 274, 276, 278, 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, 850, 900, 950, 1000, NIGRIDS = 1, ! Number of RAMS grids to analyse TOPSIGZ = 25000., ! 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 = 9999, ! 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 IOBSWIN = 1800, STASEP = .1, ! Minimum sfc station separation in degrees. ! Any surface obs within this distance ! of another obs will be thrown out ! unless it has less missing data, ! in which case the other obs will be ! thrown out. IGRIDFL = 3, ! Grid flag=0 if no grid point, only obs ! 1 if all grid point data and obs ! 2 if partial grid point and obs ! 3 if only grid data ! 4 all data... fast GRIDWT = .01,.01, ! Relative weight for the gridded press data ! compared to the observational data in ! the objective analysis GOBSEP = 5., ! Grid-observation separation (degrees) GOBRAD = 5., ! Grid-obs proximity radius (degrees) WVLNTH = 1200., 900., ! Used in S. Barnes objective analysis. ! Wavelength in km to be retained to the ! RESPON % from the data to the upper air ! grids. SWVLNTH = 750., 300., ! Wavelength for surface objective analysis RESPON = .90,.9, ! Percentage of amplitude to be retained. $END !==========================================================================================! !==========================================================================================!