c arps4.5.2.1_C115 c c ################################################################## c ################################################################## c ###### ###### c ###### INPUT FILE FOR ARPS IN NAMELIST FORMAT ###### c ###### ###### c ###### Developed by ###### c ###### Center for Analysis and Prediction of Storms ###### c ###### University of Oklahoma ###### c ###### ###### c ################################################################## c ################################################################## c c c####################################################################### c c This file contains the input parameters in the NAMELIST format. c for ARPS version 4.0 or later and a few utility programs. c c For ARPS official release, this file is configured for a c supercell storm simulation, described in the validation chapter c of ARPS User's Guide. c c####################################################################### c c NAMELIST is not a standard FORTRAN 77 feature but is supported by c almost all computer system vendors. c c The format for NAMELIST statement inside the program is: c c NAMELIST /name/ variable_list c c "name" is a NAMELIST name that is enclosed in slashes. It must not c be the same as a variable or array name, a named constant, or c a procedure. c c variable_list is list of variables or arrays separated by commas. c c The format of NAMELIST data: c c 1. NAMELIST input data c Input data must be in a special form to be read using a NAMELIST c list. The first character in each record to be read must be blank. c The second character in the first record of a group of data records c must be an ampersand (&) immediately followed by NAMELIST name. c This name is followed by (with a blank in between) data items c separated by commas (a comma after the last item is optional). c The end of a data group is signaled by &END. c c The form of the data items in an input record is: c c var_name = constant (ends with a comma) c c - The var_name can be an array element name or a variable name. c - The constant can be integer, real, complex, logical or character. c - Subscripts must be integer constants. c c array_name = Set of Constant (separated by commas) c c - The set of constants consists of the type integer, real, c complex, logical ot character. c - The number of constants must be less than or equal to the number c of element in the array. c - Successive occurrences of the same constant can be represented c in the form c*constant, where c is a nonzero integer constant c specifying the number of times the constant is to occur. c c 2. NAMELIST output data c When output data is written using a NAMELIST list, it is written in a c form that can be read using a NAMELIST list. All variable and array c names specified in the NAMELIST list and their values are written out, c each according to its type. Character data is included between c apostrophes. The fields for the data are made large enough to contain c all significant digits. The values of a complete array are written c out in column-major order. c c c Note that only lines between &NAMELIST_NAME and &END are read as the c input data, and there must be a blank space in front of the '&' sign. c c Comments can be written between these data blocks. We are using 'c' c in the first column of comment line only to distinguish them from the c data statement. c c####################################################################### c c Author: c c Ming Xue (10/1/1990) c c Modification history: c c 10/1/1993 (Ming Xue & Adwait Sathye) c Converted to namelist format. c c 04/07/94 (Yuhe Liu) c c Added surface model flag, data input flags for soil and vegetation c data and initial values, length of time step for surface model c integration, and user specified surface data and variables to the c namelist &soil_ebm. c c 04/26/94 (Ming Xue) c c Added comments and a number of additional parameters. c c 01/13/95 (Alan Shapiro, Steven Lazarus, Yvette Richardson) c c Documentation clean-up. c c 01/28/95 (Gene Bassett) c c Added namelist input for arpsr2h (gridinit). c c 02/07/1995 (Yuhe Liu) c c Added a new input parameter, veg0, to the namelist, &soil_ebm. c c 05/25/1995 (Alan Shapiro) c c Documentation clean-up. c c 01/31/1996 (V. Wong and X. Song) c c Added a new input parameter, qpfgfrq, to the namelist, µphysics. c c 02/05/1996 (Donghai Wang and Yuhe Liu) c c Added two parameters to control the calculation related to map c projection factor. c c 03/26/1996 (Yuhe Liu) c c Added a namelist &radiation c c 04/02/1996 (Donghai Wang, X. Song and M. Xue) c c Added two parameters to control the implicit treatment for c the vertical mixing. c c 05/07/1996 (Donghai Wang and M. Xue) c c Added a new option for Rayleigh damping. c c 03/27/1997 (Yuhe Liu) c c A new namelist, arpsagri, was added for ARPS Adaptive Grid c Refinement Interface (AGRI). c c 07/29/97 (Dan Weber) c c Added fftopt for use with the tbc=4 upper radiation condition. c c 10/21/97 (Donghai Wang) c Added two options,rhofctopt for using total density (rho) in the c calculation of the pressure gradient force terms, and buoy2nd c for the second order terms in the linerized buoyancy terms. c c 04/15/98 (Donghai Wang) c Added a new option for Kain-Fritsch cumulus scheme to feed back c the convectively generated rainwater into grid-resolved rainwater c (or snow) fields. c c####################################################################### c####################################################################### c c ARPSAGRI Namelist for parameters used in AGRI c c To use AGRI (Adaptive Grid Refinement Interface), the regular ARPS c driver program in arps##.f is replaced by ARPSAGRI driver in c src/agri. c c To use AGRI, compile using c c makearps arpsagri c c instead of `makearps arps`. The files in src/agri are then used. c c For regular ARPS run with executable arps built by c 'makearps arps', this namelist block is not use. It will only c be used when arpsagri is used which is built by 'makearps arpsagri'. c c The executable will be arpsagri instead of arps. To run AGRI under c the ARPS directory, enter c c bin/arpsagri < input/arps.input >! arps.output c c The base grid in the following refers to the grid at the lowest c level. All other grids sit on top of the base grid. There can only c be one base grid. c c ATTENTION: Before compile arpsagri, you need to set parameter lstore c in file include/alloc.h. c c lstore is the size of the main storage array 'a' used by ARPSAGRI c (i.e., the entire program). Array 'a' has to be large enough to c store all permanant arrays of all grids plus work arrays shared by c all grids. c c It has to be predetermined before compilation. Statistics on its use c are output on restart and exit of the interface. c Variable ihighwater is the maximum amount of a(*) used. c c####################################################################### c c nxc Dimension size in X-direction for base coarse grid c nyc Dimension size in Y-direction for base coarse grid c nzc Dimension size in Z-direction for base coarse grid c nx,ny,nz in include/dims.inc will not be used. c c levfix Level above which to regrid at start (or restart) time. c intrat Spatial refinement ratio between parent (coarse_ and c child (fine) grids, i.e., dx_course/dx_fine. c c intratt Temporal refinement ratio, i.e., dtbig_parent/dtbig_child. c Currently, the refinement ratio has to be the same between c all levels. c c kcheck Number of steps between automatic regridding. Currently c not implemented. c c intrpodr Order of spatial interpolation between coarse and fine c grids. c c = 1, for linear interpolation c = 2, for quadratic interpolation c c rstart Is this run a restart from a previous AGRI restart data? c = .true., yes c = .false., no c c runold The runname of the old run the current job restarts from, c when rstart=.true., i.e., when this run is a restart run. c c rstime Time at which to restart, if rstart=.true. c runold and rstime are used for determine the name of the c restart file. c c rstdump Dump data for later restart at the end of the run? c = .true., yes c = .fasle., no c c grdsrt regrid (i.e., place new fine grids and replace old ones if c exist at the same level) at start? c = .true., yes c = .false., no c c verbose1 Diagnostic output flag - grid location, size, etc... c = .true., output; c = .false., no output c c verbose2 Diagnostic output flag - grid variable information c = .true., output; c = .false., no output c c c verbose3 Diagnostic output flag - io information c = .true., output; c = .false., no output c c verbose4 Diagnostic output flag - storage information c = .true., output; c = .false., no output c c verbose5 Diagnostic output flag - solver diagnostics c = .true., output; c = .false., no output c c verbose6 Diagnostic output flag - random information c = .true., output; c = .false., no output c c nfinelv Number of new grid levels to be placed c c ngrdnew(nfinelv) c Number of new grids to be placed for each level c c ixc(ngrdnew,nfinelv) A real array. c The coordinate of each new grid's center in X-direction c in the unit of grid point index i of the base-grid. c c jyc(ngrdnew,nfinelv) - A real array. c The coordinate of each new grid's center in Y-direction c in the unit of grid point index j of the base-grid. c c ixln(ngrdnew,nfinelv) - A real array. c The length of each new grid in X-direction in terms of c the number of grid intervals of the base grid. c E.g., if a new grid's length is n*dx_base, then ixln=n. c c jyln(ngrdnew,nfinelv) - A real array. c The width of each new grid in Y-direction in terms of the c number of grid intervals of the base grid. c E.g., if a new grid's width is n*dy_base, then iyln=n. c c gangle(ngrdnew,nfinelv) - A real array. c The angle between the x-axis of the new grid and the c x-axis of the base grid. c c####################################################################### &arpsagri nxc = 67, nyc = 67, nzc = 35, levfix = 1, intrat = 3, intratt = 3, intrpodr = 1, kcheck = 50000, verbose1 = .false., verbose2 = .false., verbose3 = .false., verbose4 = .false., verbose5 = .false., verbose6 = .false., rstart = .false., runold = 'may20, old run name if restart', rstime = 3600.0, rstdump = .false., grdsrt = .true., nfinelv = 0, ngrdnew(1) = 1, ixc(1,1) = 34.,jyc(1,1)=34.,ixln(1,1)=22.,jyln(1,1)=22.,gangle(1,1)=0., ngrdnew(2) = 1, ixc(1,2) = 34.,jyc(1,2)=34.,ixln(1,2)=7.,jyln(1,2)=7.,gangle(1,2)=0., &END c####################################################################### c c COMMENT_LINES Comments c c nocmnt Number of comment lines c cmnt Comments c c####################################################################### &comment_lines nocmnt = 2, cmnt(1) = 'ARPS 4.5.2', cmnt(2) = 'roll test', &END c####################################################################### c c runname, a string of up to 80 characters long, is used to identify c this job. c c The first 6 characters, or the characters before either a blank space or c comma, will be used to construct output file names. Not more than 6 c characters are used to define runname. This character string will be c printed on plots produced by ARPSPLT. c c####################################################################### &jobname runname = 'C1', &END c####################################################################### c c The model can be run in 3D, 2D x-z plane, 2D y-z plane or 1D vertical c column mode. Please set c c runmod = 1 for 3-D run; c 2 for 2-D xz plane run; c 3 for 2-D yz plane run; c 4 for vertical 1-D run. c c####################################################################### &model_configuration runmod = 1, &END c####################################################################### c c initime = 'yyyy-mm-dd.hr:mn:se', UTC(GMT) date/time. yyyy is a c 4-digit integer for year, and mm, dd, c hr, mn, and se are 2-digit integers c for month, day, hour, minute, and c second, respectively. For example, c 20:30Z, Dec. 19, 1994 would be c represented as: 1994-12-19.20:30:00 c c initopt Model initialization option. c = 1, initialize using analytic functions; c = 2, initialize from a restart file; c = 3, initialize from an external data set. See inifmt below c for available history data formats. c c timeopt Options to check the consistency of user specified time c (initime and tstart) with the time of history data which is c used to start up ARPS (initopt=3) c = 1, warning on inconsistence and continue using initime and tstart c = 2, warning on inconsistence and continue using data time c = else, warning on inconsistence and stop, default c c inibasopt Initialization option for base state fields. c = 1 external sounding; c = 2, isentropic atmosphere; c = 3, isothermal atmosphere; c = 4, constant static stability atmosphere; c = 5, analytic thermodynamic sounding c (Weisman and Klemp 1982, MWR). c = 6, constant density, pot. tem. and hydrostatic base state c c For options 2, 3, 4, 5 and 6 the wind profile is specified c using option viniopt. c Note: base state fields initialized this way will be c overwritten for initopt = 2 or 3, the base state c variables are contained in the restart or external data file. c Note: if inibasopt = 4, then user must specify the static c stability in subroutine INIBASE. c c viniopt Initialization option for base state wind fields. c = 1, user specified constant ubar0 and vbar0; c = 2, user specified wind profile. c Option viniopt will be used if option inibasopt = 2, 3, 4, or 5. c Note: if option 2 is chosen, the user must specify the c desired wind profile in subroutine INIBASE. c c ubar0 Constant base state wind in x-direction. Effective when c inibasopt .ne. 1 and viniopt = 1. c c vbar0 Constant base state wind in y-direction. Effective when c inibasopt .ne. 1 and viniopt = 1. c c soilinitopt Iterative soil initialization option c = 0, no iteration c = 1, integrate soil model using the initial atmospheric c forcing. c c soiltintv Time interval for the initial integration of soil c model c c pt0opt Initial potential temperature perturbation option for initopt=1. c = 0, no initial perturbation; c = 1, bubble-shaped initial perturbation; c = 2, random initial perturbation; c = 3, symmetric random initial perturbation; c = 4, Skamarock and Klemp (1994) initial perturbation. c = 5, Soup can shaped perturbation (when used to test advection c one should also set buoyopt=0 to turn off buoyancy). c = 6, a bubble specified in terms of temperature perturbation c instead of PT. The amplitude of T pert. is prpert0. c c ptpert0 Magnitude of initial potential temperature perturbation (K). c pt0radx Bubble radius in x-direction. c pt0rady Bubble radius in y-direction. c pt0radz Bubble radius in z-direction. c pt0ctrx x coordinate of bubble center. c pt0ctry y coordinate of bubble center. c pt0ctrz z coordinate of bubble center. c c sndfile Name of the sounding file. c rstinf Name of restart file for initopt = 2. c c inifmt Data format of external data files (inifile, inigbf) for initopt = 3. c = 1, unformatted binary data; c = 2, formatted ascii data; c = 3, NCSA HDF format data; c = 4, Packed binary data; c = 7, NetCDF format (no longer support starting from 4.1.5); c = 8, Packed NetCDF format (no longer support since from 4.1.5). c = 10,GRIB format c c inifile Name of external data history file for initopt = 3. c inigbf Name of base-state/grid file for initopt = 3. c c####################################################################### &initialization initime = '2003-06-20.19:31:00', initopt = 1, timeopt = 0, pt0opt = 0, ptpert0 = 1.0, pt0radx = 4000.0, pt0rady = -2000.0, pt0radz = 1500.0, pt0ctrx = 10000.0, pt0ctry = 4750.0, pt0ctrz = 1500.0, rstinf = 'C1.rst010800.04', inifmt = 1, inifile = 'may20.grb003600', inigbf = 'may20.grbgrdbas', inibasopt = 1, viniopt = 1, ubar0 = 0.0, vbar0 = 0.0, sndfile = 'C1_quarter_circle_RUN06.snd', soilinitopt = 0, soiltintv = 1800.0, &END c####################################################################### c Options and parameters related to nudging data assimilation c c nudgopt Analysis increment nudging option. c = 0, no nudging c = 1, nudging with uniform time weight over nudging window c c ndstart Time (sec) of beginning of nudging window. c Due to complications with the first leap-frog time c step, it is recommended that ndstart >= dtbig c c ndstop Time (sec) of end of nudging window c c ndintvl Time interval (sec) to apply nudging c Due to the nature of leap-frog integration it is recommended c that ndintvl be an odd multiple of dtbig. c c ndgain Multiplier to apply to nudging at each step c (typically 1.0 =< ndgain =< 1.2) c c incrfnam File containing analysis increments (from ADAS) c c nudgu Option to apply nudging to u wind component (0:no, 1:yes) c c nudgv Option to apply nudging to v wind component (0:no, 1:yes) c c nudgw Option to apply nudging to w wind component (0:no, 1:yes) c c nudgp Option to apply nudging to pressure (0:no, 1:yes) c c nudgpt Option to apply nudging to potential temperature (0:no, 1:yes) c c nudgqv Option to apply nudging to specific humidity (0:no, 1:yes) c c nudgqc Option to apply nudging to cloud water (0:no, 1:yes) c c nudgqr Option to apply nudging to rain (0:no, 1:yes) c c nudgqi Option to apply nudging to ice (0:no, 1:yes) c c nudgqs Option to apply nudging to snow (0:no, 1:yes) c c nudgqh Option to apply nudging to hail (0:no, 1:yes) C c####################################################################### c &nudging nudgopt = 0, ndstart = 6.0, ndstop = 300.0, ndintvl = 12.0, ndgain = 1.0, incrfnam = 'dummy', nudgu = 1, nudgv = 1, nudgw = 0, nudgp = 1, nudgpt = 1, nudgqv = 1, nudgqc = 1, nudgqr = 1, nudgqi = 0, nudgqs = 0, nudgqh = 0, &END c####################################################################### c c Options and parameters related to terrain specification. c c ternopt Model terrain option. c = 0, no terrain, flat ground; c = 1, analytic mountain profile; c = 2, terrain data read in from file terndta (defined later) c mntopt Option for choosing idealized mountain type. c = 1, Bell-shaped mountain, default; c = 2, user specified (in code). c mntopt used only for ternopt = 1. c Note: For mntopt = 2, the user must specify the c desired terrain in subroutine INIGRD. c c The following options are used if ternopt = 1: c hmount Mountain height (m). c mntwidx Half-width of bell-shaped mountain in x-dir. c mntwidy Half-width of bell-shaped mountain in y-dir. c mntctrx x coordinate of the bell-shaped mountain center. c mntctry y coordinate of the bell-shaped mountain center. c c The following option is used if ternopt = 2: c terndta Name of the terrain data file for ternopt=2. c c####################################################################### &terrain ternopt = 0, mntopt = 1, hmount = 0.000, mntwidx = 10000.000, mntwidy = 10000.000, mntctrx = 10000.000, mntctry = 10000.000, terndta = 'arpstern.dat', &END c####################################################################### c c dx Grid spacing in x-direction in computational c and physical space (m). c dy Grid spacing in y-direction in computational c and physical space (m). c dz Averaged vertical grid spacing in transformed c computational space (m). c c strhopt Grid stretching option. c = 0, no vertical stretching; c = 1, vertical stretching with f=z**3 function for dz; c = 2, vertical stretching with hyperbolic tangent (see User's Guide). c dzmin Minimum vertical grid spacing in physical space (m). Used c if strhopt = 1 or 2. c zrefsfc Reference height of the surface (ground level) (m). c c dlayer1 Height (m) of the layer beneath which stretching is not applied. c 0.0 =< dlayer1 < (nz-3)*dz c c dlayer2 Depth of the mid-layer with stretched vertical spacing (m) c 0.0 =< dlayer2 < (nz-3)*dz and 0.0 =< dlayer1+dlayer2 < (nz-3)*dz c For consistency, dlayer2 is reset to: min(dlayer2,ztop-dlayer1). c c strhtune Tuning parameter used when strhopt = 2. c A value between 0.2 and 5.0 is recommended. Stretching c becomes more linear as strhtune increases. Default value is 1.0. c c zflat Height at which the grid becomes flat in the c terrain-following coordinate transformation (m). c c ctrlat Latitude of the model physical domain center (deg. N). c ctrlon Longitude of the model physical domain center (deg. E). c c Parameters ctrlat and ctrlon do not have to be set when c initopt=3, since the values in the header of the input data c file will be used in this case. c c####################################################################### &grid dx = 500.000, dy = 500.000, dz = 275.000, strhopt = 1, dzmin = 50.000, zrefsfc = 0.0, dlayer1 = 0.0, dlayer2 = 1.0e5, strhtune = 1.0, zflat = 1.0e5, ctrlat = 36.62, ctrlon = -97.50, &END c####################################################################### c c projection parameters: c c mapproj Map projection option. c = 0, no map projection; c = 1, polar projection; c = 2, Lambert projection; c = 3, Mercator projection. c trulat1 1st true latitude of map projection. c trulat2 2nd true latitude of map projection (used only by mapproj = 2). c trulon True longitude of map projection. c sclfct Map scale factor (default is 1.0). c c The above five parameters do not have to be set when c initopt=3, since the values in the header of the input data c file will be used in this case. c c mpfctopt Option parameter for map factor c = 0, map factor set to 1 c = 1, map factor calculated according to mapproj c c mptrmopt Option parameter for map factor terms in momentum advection c = 0, ignore the terms c = 1, include the terms c c####################################################################### &projection mapproj = 0, trulat1 = 30.0, trulat2 = 60.0, trulon = -80.75, sclfct = 1.0, mpfctopt = 1, mptrmopt = 1, maptest = 0, &END c####################################################################### c c dtbig Large time step (s) for model integration. c tstart Model start time. In the restart case (initopt=2), c this value is reset to the time in the restart data. c tstop Stop time for the model integration. c c####################################################################### ×tep dtbig = 6.0, tstart= 0.0, tstop = 3600.0, &END c####################################################################### c c vimplct Vertically implicit option for the w and p equations. c = 0, explicit solution; c = 1, implicit solution. c ptsmlstp Option for integrating potential temperature equation. c = 0, solve potential temperature eq. outside small time steps; c = 1, solve potential temperature eq. within small time steps. c csopt Sound wave speed option used in the pressure equation. c = 1, csound = cp/cv * rd * T; c = 2, csound = cp/cv * rd * T * csfactr; c = 3, csound = specfied constant. c Option 1 should be used whenever possible. Reduced sound c wave speed may result in inaccurate solution. c c csfactr Multiplication factor for the sound speed if csopt=2. c csound User specified constant sound speed if csopt=3. c tacoef Weighting coefficient for time average in the vertically c implicit solver. (see User's Guide) c dtsml Small time step (s) for integrating acoustic wave modes. c c####################################################################### &acoustic_wave vimplct = 1, tacoef = 0.6, csopt = 1, csfactr = 0.5, csound = 150.0, ptsmlstp = 0, dtsml = 1.0, &END c####################################################################### c c buoyopt Flag for turning buoyancy terms on or off. c = 1, buoyancy terms included; c = 0, buoyancy terms turned off. c buoy2nd Option for the second order terms in the linerized buoyancy terms. c = 1, including the 2nd-order terms; c = 0, only the 1st-order terms. c c rhofctopt Option for removing the density approximation in the c pressure gradient force(PGF) terms. c = 1, removing the approximation,using total density in PGF terms; c = 0, using the base state density(rhobar) in PGF terms. c c bsnesq Bousinessq approximation (used by inibaseopt=3 or 4 only)? c = 1, yes c = 0, no c c peqopt Option for an alternative formulation for pressure equation c = 1, Original formulation as described in ARPS 4.0 User's Guide. c = 2, An alternative formulation for special applications. c Option 1 recommended. c c####################################################################### &equation_formulation buoyopt = 1, buoy2nd = 1, rhofctopt = 1, bsnesq = 0, peqopt = 1, &END c####################################################################### c c madvopt Momentum advection option. c = 1, second order advection; c = 2, fourth order horizontal and second order vertical advection. c = 3, fourth order advection in both the horizontal and vertical. c c sadvopt Scalar advection option. c = 1, second order advection; c = 2, fourth order horizontal and second order vertical advection. c = 3, fourth order advection in both the horizontal and vertical. c = 4, Zalesak's multi-dimensional version of FCT based on c second-order centered and first-order upstream schemes. c FCT is applied to potential temperature, water variables c and TKE, while either 2nd or 4th order advection is used c for pressure depending option fctorderopt. c = 5, simple positive definite advection (MPDCD) scheme c based on flux correction/limiting on leapfrog-centered c advective fluxes. c With this option, positive definite water variables and TKE c are advected using this scheme while potential temperature c and pressure are advected by either 2nd or 4th-order centered c scheme (i.e., sadvopt=1 or 3) depending option fctorderopt. c c fctorderopt Option of the spactial order of accuracy of the FCT advection c (sadvopt=4) and MPDCD advection schemes (sadvopt=5) c = 1 2nd order c = 2 4th order c c fctadvptprt Option for FCT advection for potential temperature. c Used only when sadvopt=4. c = 0, FCT scheme is applied to ptbar+ptprt. Not recommended. c = 1, FCT scheme is applied to ptprt only. This option is RECOMMENDED! c = 2, FCT scheme is applied to ptbar+ptprt-min(ptbar+ptprt). c c The most accurate (also most expensive) choices are: c madvopt=3, sadvopt=4, fctorderopt=2 with fctadvptprt=1. c The most econimical choices are: c madvopt=1, sadvopt=1. c c####################################################################### &numerics madvopt = 1, sadvopt = 1, fctorderopt = 2, fctadvptprt = 1, &END c####################################################################### c c lbcopt Lateral boundary condition option. c = 1, All boundary condition options except externally-forced c (option 5 for ebc, wbc, nbc, or sbc is not allowed); c = 2, Externally-forced lateral boundary conditions. In this c case, ebc, wbc, nbc, and sbc will be set to 5. c c wbc West boundary condition option. c ebc East boundary condition option. c sbc South boundary condition option. c nbc North boundary condition option. c = 1 Rigid wall; c = 2 Periodic; c = 3 Zero gradient; c = 4 Radiation (open) lateral boundary; c = 5 Externally-forced lateral boundary; c = 6 Nested grid lateral boundary. c tbc Top boundary condition option. c = 1 Rigid wall; c = 2 Periodic; c = 3 Zero gradient: c = 4 Linear hydrostatic radiation top boundary: c References: Klemp and Durran MWR, 1983 and Chen MWR, 1991 c c ********THIS OPTION CANNOT BE USED ON MPP MACHINES********* c c This condition requires a statically stable base state c at scalar nz-2. It will run with a neutral environment c but the accuracy (and application of the condition) is c questionable. In addition, zflat must be set c to a level at or below the scalar point nz-3. c c fftopt Fast Fourier Transform method for use with the upper c boundary tbc=4. c c = 1, periodic transform used, edges are assumed to be c equal in value. c c Requires special dimensions for (nx,ny) given by c c nx-1 = 2**P * 3**Q * 5**R c ny-1 = 2**P * 3**Q * 5**R where c c P .GE. 1 , Q .GE. 0 , and R .GE. 0 . (nx,ny must be odd!) c For a xz run nx should be odd and ny = 4. c For a yz run ny should be odd and nx = 4. c For a xyz run nx and ny should be odd. c c = 2, even Cosine transform used, edges are NOT assumed to be c equal in value. To determine the nx and ny required for c this fft choice, ADD 1 to the value obtained from the c above equation: c c nx-1 = 2**P * 3**Q * 5**R + 1 c ny-1 = 2**P * 3**Q * 5**R + 1 c c P .GE. 1 , Q .GE. 0 , and R .GE. 0 . (nx,ny SHOULD be even!) c For a xz run nx should be even and ny = 4. c For a yz run ny should be even and nx = 4. c For a xyz run nx and ny should be even. c NOTE: The simulation will NOT stop if an incorrect even value c is selected. In this case the transform will be the c slower simple fourier transform, NOT a FAST fourier transform. c c bbc Bottom boundary condition option. c = 1 Rigid wall; c = 2 Periodic; c = 3 Zero gradient; c c rbcopt Radiation lateral boundary condition option (used if c radiation condition is chosen for wbc,ebc,sbc,nbc). c (Note: These condition are applied to horizontal c velocities u and v ONLY.) c = 1, Klemp & Wilhelmson type with constant phase speed, c; c computed AND applied on the SMALL time step. c = 2, Klemp & Wilhelmson type with constant phase speed, c; c computed on the BIG time step and applied on the c SMALL time step. c = 3, Orlanski (1976) condition computed on the BIG time step c and applied on the SMALL time step. c = 4, Klemp-Lilly /Durran (1983) condition for u and v. c Computes the Orlanski phase speed on the BIG time step c and vertically averages the phase speed. The phase c speed is then applied on the small time step. c c c_phase Constant phase speed for rbcopt=1 (only). c rlxlbc Relaxation coeff. used by RBC option c 0.0=< rlxlbc =< 0.5 c pdetrnd Option switch for detrending the pressure field. c With the option on, the domain averaged perturbation c Exner function is reset to zero every time step to remove c domain-wide pressure drift/trend sometimes seen when c open boundary condition is used. c The detrending SHOULD NOT be used when the model is c initialized with 3D fields (real data). c = 0, no detrend; c = 1, with detrend. c c####################################################################### &boundary_condition_options lbcopt = 1, wbc = 2, ebc = 2, sbc = 2, nbc = 2, c_phase = 50.0, rlxlbc = 0.0, rbcopt = 1, tbc = 4, fftopt = 1, bbc = 1, pdetrnd = 0, &END c####################################################################### c c EXBCPARA Parameters by the external boundary conditions (lbcopt=2). c c exbcname The prefix of the input external boundary file names. c c tinitebd = The time in 'yyyy-mm-dd.hh:mm:ss' format for the c first external boundary data file. The file must be c be named in format exbcname//'.yyyymmdd.hhmmss'. c c tintvebd = Time interval (s) at which external boundary data files c will be searched. c c ngbrz = Number of grid zones in the boundary relaxation zone. c brlxhw = Half-width of the boundary relaxation function in term c of the number of grid zones (a real number). c cbcdmp = The coefficient of relaxation in the relaxation zone (1/s). c cbcmix = The coefficient (1/s) of additional second-order horizontal c computational mixing in the boundary zone. c c####################################################################### &exbcpara exbcname = 'arpsexbc', tinitebd = '1977-05-20.21:00:00', tintvebd = 3600, ngbrz = 5, brlxhw = 2.3, cbcdmp = 0.0033333333, cbcmix = 1.0e-3, &END c####################################################################### c c coriopt Coriolis term option. c = 0, no coriolis terms. Default; c = 1, coriolis terms involving horizontal wind; c = 2, coriolis terms involving both horizontal and vertical wind; c = 3, as 1, but the Coriolis parameters are latitude dependent; c = 4, as 2, but the Coriolis parameters are latitude dependent; c c coriotrm An option for imposing an approximate geostrophic initial c balance between the base state winds and the pressure gradient c force. If coriotrm=2 the Coriolis terms in the momentum c equations are modified from their standard formulation so that c f(u-ubar) and f(v-vbar) are used in place of fu and fv. Here c f*ubar and f*vbar represent the geostrophic pressure gradient c forces associated with ubar and vbar. It may be desirable c to impose this balance if the model is initialized from a c single sounding (the base-state pressure gradient being zero) c and the user wishes to redefine the pressure gradient to be c approximately consistent with a geostrophic balance. This c option is not used if coriopt=0 (no Coriolis force). c c = 1, No balancing step. Total u and v are used in the Coriolis terms; c = 2, An approximate geostrophic balance imposed initially. c u-ubar and v-vbar are used in the place of u and v in the c Coriolis formulation. c c####################################################################### &coriolis_force coriopt = 0, coriotrm = 1, &END c####################################################################### c c Subgrid-scale turbulent mixing parameters. c c tmixopt Control parameter for turbulent mixing options. c = 0, zero turbulent mixing; c = 1, constant mixing coefficient, km = tmixcst; c = 2, Smagorinsky mixing coefficient; c = 3, Smagorinsky mixing coefficient c plus a constant coeffcient, tmixcst; c = 4, 1.5 TKE turbulent mixing. c c trbisotp Option for isotropic subgrid scale turbulence. c = 0, the turbulence is assumed to be anisotropic, c Use when dx>>dz. c = 1, the turbulence is assumed to be isotropic (default). c Use when dx ~ dz. c c tkeopt Option for 1.5 order TKE formulation used by tmixopt=4 c = 1, after Moeng and Wyngaard (default); c = 2, after Deardroff; c = 3, after Sun and Chang (1986) J.Climate Appl. Meteor. c c trbvimp Option for implicit treatment of vertical mixing c = 0, vertical explicit (default); c = 1, vertical implicit (always use this option c when tmixopt=4 AND tkeopt=3) c c tmixvert Option for computing only the vertical mixing terms. c c Use for synoptic scale forecasts (where dx>>dz). c Time intensive horizontal components in the turbulence c mixing are neglected. c c = 0, full turbulence formulation. c = 1, only vertical components are retained. c c alfcoef Time average weighting coefficient (for past time level) used c in the vertically implicit mixing. When trbvimp=0, c it is reset to 1. Otherwise, 0.25 is recommended. c c prantl Constant turbulent prandtl number used by Smagorinsky option c c tmixcst Value of the constant mixing coefficient (m**2/s) when c tmixopt=1 or 3 c c kmlimit Nondimensional upper limit on mixing coefficient. Upper c limit for stability is 1. c c####################################################################### &turbulence tmixopt = 4, trbisotp = 1, tkeopt = 1, tmixcst = 0.0, tmixvert = 0, prantl = 1.0, trbvimp = 1, alfcoef = 0.25, kmlimit = 1.0, &END c####################################################################### c c Computational mixing parameters. c c c cmix2nd 2nd order computational mixing option. c = 0, mixing off; c = 1, mixing on. c cfcm2h 2nd order horizontal computational mixing coefficient scaled c by horizontal grid spacing (1/s). c cfcm2v 2nd order vertical computational mixing coefficient scaled c by vertical grid spacing (1/s). c c cmix4th 4th order computational mixing option. c = 0, mixing off; c = 1, mixing on. c cfcm4h 4th order horizontal computational mixing coefficient scaled c by horizontal grid spacing (1/s). c cfcm4v 4th order vertical computational mixing coefficient scaled c by vertical grid spacing (1/s). c scmixfctr Reduction factor of the computational mixing coefficients c for scalars relative to those of velocities, the c-mixing c coefficients are multiplied by a factor of scmixfctr for scalars. c Default is 1. c c####################################################################### &computational_mixing cmix2nd = 1, cfcm2h = 0.0, cfcm2v = 1.0e-3, cmix4th = 1, cfcm4h = 1.0e-3, cfcm4v = 0.0e-4, scmixfctr = 1.0, &END c####################################################################### c c Acoustic wave divergence damping parameters. c c divdmp Acoustic wave divergence damping option. c = 0, divergence damping off; c = 1, isotropic divergence damping on. c = 2, anisotripic divergence damping on. c divdmpndh Non-dimensional divergence damping coefficient in c horizontal direction c divdmpndv Non-dimensional divergence damping coefficient in c vertical direction c c####################################################################### &divergence_damping divdmp = 1, divdmpndh = 0.05, divdmpndv = 0.05, &END c####################################################################### c c Upper level Rayleigh damping parameters. c c raydmp Rayleigh damping option. c = 0, Damping off; c = 1, Damping difference between total and base state fields; c = 2, Damping difference between total and external fields c defined in the EXBC file. In this case, lbcopt must c be set to 2. c c zbrdmp Height of the bottom of Rayleigh damping (sponge) layer (m). c cfrdmp Rayleigh damping coefficient (1/second). c c####################################################################### &rayleigh_damping raydmp = 1, cfrdmp = 0.00333, zbrdmp = 12000.0, &END c####################################################################### c c Robert-Asselin time filter coefficient for leapfrog time c c flteps Robert-Asselin time filter coefficient (non-dimensional). c c####################################################################### &asselin_time_filter flteps = 0.1, &END c####################################################################### c c Moisture/microphysics parameters: c c moist Moist processes option. c = 0, Dry run, all processes related to moisture are turned off; c = 1, Moist processes are activated. If mphyopt.eq.0, c water (qc and qr) or ice equations will not be solved. c c mphyopt Microphysics option. c = 0, No microphysics process. Warm (liquid) saturation c adjustment is performed; c = 1, Kessler warm rain microphysics; c = 2, Ice microphysics; c = 3, Schultz NEM ice microphysics. c c cnvctopt Option for convective cumulus parameterizations. c = 0, no convective parameterization and grid-scale c condensation; c = 1, Kuo scheme with its own grid-scale condensation; c = 2, Kuo scheme and Kessler warm rain microphysics; c = 3, Kain and Fritsch cumulus parameterization c c kffbfct Factor for Kain-Fritsch scheme,to feed convectively c generated rainwater into grid-resolved rainwater c (or snow) field. kffbfct is the fraction of available c precipitation to be fed back (0.0 - 1.0). c =0.0, no feed back; c =1.0, all convective rainwater feed back, c so no cumulus rainfall in this case. c 0.0 < kffbfct <= 1.0 recommended when horizontal grid c spacing is less than 25km. c c The following four parameters are used by Kuo scheme only. c c wcldbs Vertical motion needed at cloud base for convection. c confrq Frequency of convection parameters' updates in seconds c qpfgfrq Frequency of grid parameters' updates in seconds c idownd Downdraft flag. c = 0, no downdrafts; c = 1, simple downdraft model. c c####################################################################### c µphysics moist = 1, mphyopt = 1, cnvctopt = 0, kffbfct = 0.0, wcldbs = 0.005, confrq = 600.0, qpfgfrq = 120.0, idownd = 1, &END c c####################################################################### c c Radiation physics input parameters: c c radopt Option to switch on/off radiation physics c = 0, No radiation physics; c = 1, Simplified surface radiation physics; c = 2, Atmospheric radiation transfer parameterization. c c Notes: 1) When radopt=2, the dimension parameters, nx_rad, ny_rad and c nz_rad in file dims.inc must be set equal to nx, ny and c nz in the same file. Otherwise they can and should be c set to 1 to save significant memory usage. c c 2) When sfcphy is chosen to 3 or 4, radopt=0 will be reset c to 1 in order to compute the surface energy balance for c soil model. c c radstgr Option for radiation computing at staggering points (used by c radopt = 2 only). c c = 0, No staggering; Radiation calculation on x-y plane is at c all points; c = 1, staggering; Radiation calculation on x-y plane is at c (even,even) and (odd,odd) points. The values at c (even,odd) (odd,even) points are averaged from the c surrounding four points. For example for nx=ny=9, the c directly calculation are performed at the "x" points, c then calculate radiation variables at "o" by averaging c from their surrounding "x" points. (The "." points are c not updated since they are unused for scalar variables). c This scheme can reduce ALMOST HALF of radiation calculation. c c c j c c 9 | . . . . . . . . . c 8 | o x o x o x o x . c 7 | x o x o x o x o . c 6 | o x o x o x o x . c 5 | x o x o x o x o . c 4 | o x o x o x o x . c 3 | x o x o x o x o . c 2 | o x o x o x o x . c 1 | x o x o x o x o . c +------------------- i c 1 2 3 4 5 6 7 8 9 c c On boundary, the zero-gradient is assumed. c c NOTE: if ny is not odd the message passing version will not produce c results which are exactly the same as a single processor run for c radstgr=1. c c rlwopt Option to choose the longwave schemes. c = 0, transmission functions are c computed using the k-distribution method with linear c pressure scaling. cooling rates are not calculated c accurately for pressures less than 20 mb. The c computation is faster with this option. c = 1, transmission functions in the c co2, o3 in the co2, o3, and the three water vapor bands c with strong absorption are computed using table look-up. c cooling rates are computed accurately from the surface c up to 0.01 mb. c c dtrad Time interval (seconds) to update the radiation forcing c (used by radopt = 2 only). c c raddiag Option to dump radiation variables to a file in GrADS c format for diagnostic review. The frequency is controled by c dtrad (used by radopt = 2 only). c c = 0, no such dump c = 1, dump to a file with a name like 'runname.radout' c and its control file has a name like 'runname.radctl' c c####################################################################### c &radiation radopt = 2, radstgr = 1, rlwopt = 1, dtrad = 600.0, raddiag = 0, &END c c####################################################################### c c sfcphy Surface physics options. c = 0, No surface physics; c = 1, Surface fluxes are calculated from constant surface c drag coefficients, and user-specified values of surface c potential temperature and relative humidity; c = 2, Surface fluxes are calculated from the c stability-dependent surface drag coefficients, and user- c specified values of surface potential temperature and c relative humidity; c = 3, Surface fluxes are calculated from constant surface c drag coefficients, and predicted surface temperature c and surface volumetric water content; c = 4, Surface fluxes are calculated from the stability-dependent c surface drag coefficients, and predicted surface c temperature and surface volumetric water content. c c landwtr Flag indicating whether or not a distinction is made between c land and water surfaces in the surface physics calculations. c = 0 No distinction between land and water c = 1 Land and water are treated differently. c cdhwtropt Option to use cdhwtr instead of calculated values for cdh c (and cdq) over water even for sfcphy=2 or 4. c = 0 Use calculated values c = 1 used specified value cdhwtr. c cdmlnd Land surface momentum drag coefficient. c cdmwtr Water surface momentum drag coefficient. c cdhlnd Land surface heat exchange coefficient. c cdhwtr Water surface heat exchange coefficient. c cdqlnd Land surface moisture exchange coefficient. c cdqwtr Water surface moisture exchange coefficient. c c pbldopt The option for PBL depth determination. c = 1, PBL depth is user-specified (as pbldpth0); c = 2, PBL depth is diagnosed; c c pbldpth0 Specified PBL depth for option 1 and 2. c lsclpbl0 PBL length scale used for tkeopt=3 c (0.25 recommended by Sun and Chang 1986). c c tqflxdis Option for distributing heat and moisture fluxes qudratically c in a specified depth dtqflxdis. Use only when near surface c vertical resolution is high (<50m). c = 0, no distribution; c = 1, with distribution c = 2, with distribution over a depth according to similarity c dtqflxdis Depth of flux distribution for tqflxdis=1, 200 m recommended. c c smthflx Option to smooth surface fluxes c = 0, no smoothing c = 1, smoothing c numsmth Number of smooth passes (>=1 if smthflx=1) c c sfcdiag Flag controlling output of surface diagnostic calculations. c c####################################################################### &surface_physics sfcphy = 3, landwtr = 0, cdhwtropt= 0, cdmlnd = 3.0e-3, cdmwtr = 1.0e-3, cdhlnd = 2.0e-3, cdhwtr = 1.0e-3, cdqlnd = 1.2e-3, cdqwtr = 0.7e-3, pbldopt = 2, pbldpth0 = 1400.0, lsclpbl0 = 0.25, tqflxdis = 2, dtqflxdis= 200.0, smthflx = 0, numsmth = 1, sfcdiag = 0, &END c####################################################################### c c The following surface parameters are valid for sfcphy = 3 and 4: c c sfcdat Option for defining the surface characteristics when initopt.ne.2. c c = 1, Surface characteristics are defined using input parameters. c = 2, Surface characteristics are defined using file sfcdtfl; c = 3, Same as sfcdat =2, except when initopt=3 and c the variables are present in the grid/base state c file inibgf, the values in inibgf will be used instead. c c This option is not used when initopt=2, i.e., for restart runs. c In this case, data in the restart file will be used. c c styp Soil type (an integer). Used if sfcdat=1. c The soil type is based on USDA definitions along with c categories for ice and water. c c 01 Sand c 02 Loamy sand 11 17 23 c 03 Sandy loam 14 20 26 27 c 04 Silt loam c 05 Loam 12 18 24 c 06 Sandy clay loam 15 21 28 c 07 Silty clay loam c 08 Clay loam 13 c 09 Sandy clay 19 25 c 10 Silty clay 16 22 c 11 Clay 29 30 31 c 12 Ice 34 c 13 Water 00 c c Note: The numbers on the right hand side above represent c Mylne and Henderson-Sellers soil classes. c c Default: 3 for Norman, Oklahoma c c vtyp Vegetation type (an integer). Used if sfcdat=1. c c 01 Desert c 02 Tundra c 03 Grassland c 04 Grassland with shrub cover c 05 Grassland with tree cover c 06 Deciduous forest c 07 Evergreen forest c 08 Rain forest c 09 Ice c 10 Cultivation c 11 Bog or marsh c 12 Dwarf shrub c 13 Semidesert c 14 Water c c Default: 10 for Norman, Oklahoma c c lai Leaf Area Index. Used if sfcdat=1. Default: 0.31 c c roufns0 Surface roughness. Used if sfcdat=1. Default: 0.01 c c veg0 Vegetation fraction. Used if sfcdat=1. Default: 0.3 c c sfcdtfl Data file containing the surface characteristics c (soil and vegetation type, leaf area index and surface roughness). c c c soilinit Soil model variable initialization option used when initopt.ne.2. c c = 1, Soil model variables are initialized using input parameters; c c = 2, Soil model variables are initialized using values found c in file soilinfl. c For variables missing in soilinfl, the values in initial c file inifile will be used when initopt=3. In another word, c the values in soilinfl take precedence over those in inifile. c c = 3, As soilinit=2, except that the values found in inifile take c precedence over those found in soilinfl. c c = 4, Soil temperature variables are initialized by adding c offsets to the surface air temperatue, while soil moisture c variables are initialized from given saturation rates. c The canopy water amount is initialized from its default c value, wetcanp0 though. c c This option is not used when initopt=2. c When initopt=2, data in the restart file will be used. c c ptslnd0 Initial land surface potential temperature (K). c Used by option soilinit=1. c c ptswtr0 Initial water surface potential temperature over water (K). c Used by option soilinit=1. c c tsiol0 Initial deep soil temperature in K. Used by option soilinit=1. c c wetsfc0 Initial ground surface soil moisture. Used by option soilinit=1. c c wetdp0 Initial deep soil moisture. Used by option soilinit=1. c c wetcanp0 Initial canopy moisture. Used by option soilinit=1. c c snowdpth0 Initial snow depth (m). Used by option soilinit=1. c c tsprt Offset of tsfc from sfc air temperature c c t2prt Offset of tsoil from sfc air temperature c c wgrat Saturation rate of sfc soil moisture c c w2rat Saturation rate of deep soil moisture c c soilinfl Data file containing the initial values of soil model variables c (ground surface temperature, deep soil temperature, c ground surface soil moisture, deep soil moisture, and c canopy moisture) c c dtsfc Time step for surface (soil) model integration. dtsfc =< dtbig. c c Note: The above options are effective only when sfcphy = 3 or 4. c c####################################################################### &soil_ebm sfcdat = 1, styp = 3, vtyp = 1, lai0 = 0.31, roufns0 = 0.01, veg0 = 0.3, sfcdtfl = 'test.sfcdata', soilinit = 4, ptslnd0 = 293.0, ptswtr0 = 288.0, tsoil0 = 297.0, wetsfc0 = 0.25, wetdp0 = 0.25, wetcanp0 = 0.00, snowdpth0 = 0.0, tsprt = 0.0, t2prt = 0.0, wgrat = 0.7, w2rat = 0.7, soilinfl = 'may20.soilinit', dtsfc = 60.0, &END c####################################################################### c c Options for grid translation and cell-tracking. c c cltkopt Option for performing cell-tracking (tracking positions of cells). c = 0, no cell-tracking, c = 1, cell-tracking on. If grdtrns=2, cltkopt will be forced to 1. c c c c NOTE: Cell-tracking is NOT supported at this time. c c c c grdtrns = 0, No grid translation. c = 1, Grid translation based on user-specified domain translation c speed (umove,vmove). The speed does not change during the c simulation. c = 2, Grid translation based on results from the cell-tracking c algorithm. Grid motion is such that center of mass of c cells is kept near center of grid. Grid motion changes c during the simulation. c = 3, Grid translation based on results from the optimal scalar c pattern-translation algorithm. Grid moves at "optimal" c patterm-translation speed. Grid motion changes during c the simulation. c c umove User-specified initial domain translation speed in x-dir. c vmove User-specified initial domain translation speed in y-dir. c They remain unchanged for grdtrns=1, and are adjusted c during the run for grdtrns=2 or 3. c They are not used when grdtrns=0. c c chkdpth The domain depth over which scalar pattern translation is c computed when grdtrns = 3. c twindow The time window within which the average domain translation c speed is calculated. Used by option grdtrns=3. c c tceltrk time interval between calls to the cell tracking routine. c tcrestr time required for the cell center to be restored to the c domain center. Used when grdtrns=2. c c####################################################################### &grdtrans cltkopt = 0, tceltrk = 120.0, tcrestr = 1800.0, grdtrns = 0, chkdpth = 2500.0, twindow = 300.0, umove = 0.0, vmove = 0.0, &END c####################################################################### c c history_dump parameters. c c hdmpopt = History data dump option. c = 1, linear dump, start from tstrdmp; c = 2, dump at model time specified by user. c c hdmpfmt History data dump format option. c = 0, no data dump is produced; c = 1, unformatted binary data dump; c = 2, formatted ascii data dump; c = 3, NCSA HDF4 format data dump; c = 4, Packed binary data dump; c = 5, dump for Savi3D visualization package; c = 6, binary allowing data point skip; c = 7, NetCDF format (no longer support starting from 4.1.5); c = 8, Packed NetCDF format (no longer support since from 4.1.5); c = 9, GrADS data dump; c = 10, GRIB data dump; c = 11, Vis5D data dump. c c grbpkbit Number of bits in packing GRIB data c = 16 (default) c c hdfcompr HDF4 compression option (for hdmpfmt=3). c = 0 (default), no compression; c = 1, RLE compression; c = 2, high gzip compression; c = 3, fast gzip compression; c = 4, adaptive or skipping Huffman compression. c c thisdmp time interval (s) between history data dumps when hdmpopt=1. c Choose 0.0 if no history data dump is desired. c tstrtdmp time at which history dumps start. c c numhdmp number of history dumps specified by user for hdmpopt=2. c Choose 0 if no history data dump is desired. c hdmptim array of maximum size 100 where the user specified history c dumping times are stored. c c####################################################################### &history_dump hdmpopt = 1, hdmpfmt = 9, grbpkbit = 16, hdfcompr = 0, thisdmp = 180.0, tstrtdmp = 0.0, numhdmp = 3, hdmptim(1) = 0., hdmptim(2) = 3600., hdmptim(3) = 7200., &END c####################################################################### c c output control parameters. c c dirname Name of directory into which output files are written. c c exbcdmp Flag to dump ARPS external boundary data files. c = 0, no EXBC dumping; c = 1, EXBC dumping. c c extdadmp Flag to dump the fields that contains external data array c interpolated to the current model time. When lbcopt.ne.2, c reset extdadmp to 0. c c filcmprs Option to compress the history dumping files. c = 0, history files not compressed. c = 1, history files compressed; c readyfl Option to create a marker file (same name is the history dump c but with "_ready" appended to the end) to indicate that c writing of the history dump has completed. c = 0, do not create a ready file. c = 1, create a ready file; c grdout Grid output option. c = 0, no grid array output in time-dependent history files; c = 1, grid arrays written in time-dependent history files. c basout Base state field output option. c = 0, no base state arrays time-dependent history files; c = 1, base state arrays written in time-dependent history files. c varout Perturbation fields output option in history dump. c = 0, no perturbation wind, pressure or pot. temperature output; c = 1, perturbation wind, pressure and pot. temperature output. c mstout Moist variable output option in history dump. c = 0, no moisture variables output; c = 1, qv, qc, qr, qi, qs and qh output. c iceout Ice variable outout option in history dump. c = 0, no ice variables output; c = 1, qi, qs and qh output. c tkeout TKE output option in history dump. c = 0, no c = 1, yes c trbout Turbulence field (km) output option in history dump. c = 0, no km output; c = 1, km output. c rainout Option for surface accumulated rainfall array output. c = 0, no, c = 1, yes. c sfcout Soil model variable output option in history dump. c = 0, no, c = 1, yes. c landout Option for soil and vegetation property array output in c history dump. c = 0, no, c = 1, yes. c prcout Precipitation rates output option in history dump. c = 0, no, c = 1, yes. c radout Radiation arrays output option in history dump. c = 0, no, c = 1, yes. c flxout Surface fluxes output option in history dump. c = 0, no, c = 1, yes. c c qcexout Option for qc array output in EXBC file dump when exbcdmp=1 c = 0, no, c = 1, yes. c c qrexout Option for qr array output in EXBC file dump when exbcdmp=1 c = 0, no, c = 1, yes. c c qiexout Option for qi array output in EXBC file dump when exbcdmp=1 c = 0, no, c = 1, yes. c c qsexout Option for qs array output in EXBC file dump when exbcdmp=1 c = 0, no, c = 1, yes. c c qhexout Option for qh array output in EXBC file dump when exbcdmp=1 c = 0, no, c = 1, yes. c c tfmtprt time interval (s) between formatted print out. Choose c 0.0 if no print out is desired. c trstout time interval between restart data dumps. c tmaxmin time interval between max/min statistics calc. c tenergy time interval between energy statistics calc. c c imgopt HDF image dumping option, 0 or 1. c timgdmp time interval between HDF image dumps. c c pltopt Graphic plotting option, 0 or 1. c tplots time interval between graphic plotting calls. c c####################################################################### &output dirname = './', filcmprs = 0, readyfl = 0, basout = 0, grdout = 0, varout = 1, mstout = 1, iceout = 1, tkeout = 1, trbout = 1, rainout = 1, sfcout = 1, landout = 0, prcout = 0, radout = 0, flxout = 0, exbcdmp = 0, extdadmp = 0, qcexout = 0, qrexout = 0, qiexout = 0, qsexout = 0, qhexout = 0, tfmtprt = 0.0, trstout = 999999.0, tmaxmin = 6.0, tenergy = 0.0, imgopt = 0, timgdmp = 60.0, pltopt = 0, tplots = 600.0, &END c####################################################################### c c Debug information printing level. c c lvldbg Level of debug information printing. c =0, no printing; c =1, Print variables in big t-step; c =2, Print forcings in big t-step; c =3, Print variables in small t-step; c =4, Print forcings in small t-step; c =5, Print individual forcing terms and misc. info. c c####################################################################### &debug lvldbg = 0, &END c####################################################################### c c EXTDFILE Namelist used by EXT2ARPS program (not by c ARPS) for the external boundary data pre-processing. c c iorder Order of polynomial used for interpolation c = 1 Linear c = 2 Quadratic c = 3 Cubic (default) c Option "3" is recommended c c nsmooth Number of smoothing passes after interpolation c 1 (default) or 2 recommended. c c obropt O'Brien adjustment option. Determines c distribution of mean divergence error used to c enforce upper boundary condition of w=0 c = 1 Linearly in computational z. (default) c = 2 Linearly in physical z. c = 3 Linearly in potential temperature. c c Adding 10 to the previously listed options makes c w=0 at the bottom of the Rayleigh damping layer c (zbrdmp) rather that at the physcial top of c the domain. c c obropt=11 is recommended for most c ext2arps applications to large-scale models c c obrzero Level above which w is set to zero for the O'Brien c adjustment. c c hydradj Option for adjusting pressure to balance gravity c and buoyancy to force dw/dt=0 for output field. c = 0 no correction c = 1 correction applied beginning at surface c = 2 correction applied beginning at top c = 3 hydrostatic equation integrated from c surface for total-p, pbar subtracted. c hydradj=0 suggested c c wndadj Option for adjusting interpolated winds. c = 0 no adjustment from read-in values c if w is not read-in, it is set to zero c = 1 w set so that wcont = 0, no adjustment to u or v c = 2 w set using integrated divergence with OBrien c correction to satisfy top and bottom bc c = 3 as in 2, but u and v adjusted for c to remove 3-D divergence. c c dir_extd = The directory that contains the external data files c that are to be converted/processed. c c extdname = Prefix string of external file name c extdopt Options for external data sources c = 0, ARPS, default c = 1, NMC RUC (Hybrid-B) data in GRIB (grid #87) c = 2, NMC ETA data in GRIB (grid #212) c = 3, OLAPS data c = 4, GEMPAK RUC data c = 5, GEMPAK ETA data (see also extdopt=10) c = 6, COAMPS data c = 7, NMC RUC AWIPS data in GRIB, grid #211 c = 8, NMC T62 Gaussian grid reanalysis data in GRIB c = 9, GEMPAK RUC-2 data c = 10, GEMPAK ETA data, grid #104 c = 11, NCEP native coordinate RUC2 data in GRIB (grid #236) c = 12, NCEP isobaric RUC2 data in GRIB (grid #236) c >= 12, User defined c NCEP AVN GRIB global data, grid #3 (not implemented yet). c c extdfmt = Flag for ARPS (external) input data format used when c extdopt=0. For definition, see description for hdmpfmt. c c nextdfil = The number of external data files for EXT2ARPS to c process. c c extdtime = The time corresponding to the external data files. c The format is a concatenation of the inidate and c initime formats, plus forecast time information, namely: c extdtime(1) = 'yyyy-mm-dd.hh:mm:ss+hhh:mm:ss', c For example: c extdtime(1) = '1977-05-20.21:00:00+003:00:00' c represents 3 hr forecast fields starting from c 21z, May 20 1977 c c####################################################################### &extdfile extdopt = 0, extdfmt = 10, dir_extd = './', extdname = 'may20', nextdfil = 1, extdtime(1) = '1977-05-20.21:00:00+000:00:00', extdtime(2) = '1977-05-20.21:00:00+000:00:00', extdtime(3) = '1977-05-20.21:00:00+000:00:00', iorder = 2, nsmooth= 1, hydradj= 0, wndadj = 2, obropt = 11, obrzero = 12000., &END c####################################################################### c c SOIL_VEG_DATA Additional namelist used by ARPSSFC program (not by c ARPS) for the surface data pre-processing. c c schmopt Option of schemes to generate the distribution of soil c and vegetation data set in ARPS domain. c = 0, Constant for the entire domain; c = 1, One constant value in a user specified rectangular c region (foreground) and another in the rest area c (background); c = 2, Constant for the foreground region and c real data for the background; c = 3, Real data for the entire domain files. c = 4, Use real data for stype and vtype (retrieved from files c stypfl and vtypfl). Lai, veg & roufns are derived c from vegetation type using a look-up table (in c file lkupfl). c c sdatopt Option to select soil texture data sets c = 1, select STATSGO soil texture data set. The c compressed data set can be obtained from CAPS ftp c directory: c c ftp://ftp.caps.ou.edu/ftp/pub/ARPS/arps40.data/arpssfc.data/soil_1km.data.gz c c = 2, select GED soil type data set. The compressed data c sets can be obtained from CAPS ftp directory: c c ftp://ftp.caps.ou.edu/ftp/pub/ARPS/arps40.data/arpssfc.data/whsoil.data.Z c owe14d.data.Z c ndvi90##.data.Z c where ## represents two-digit of month number. c c vdatopt Option to select vegetation type data sets c = 1, select North American 1km data set. c = 2, select GED vegetation type data set. c c ndatopt Option to select NDVI data sets c = 1, select North American 1km data set. c = 2, select GED NDVI data set. c c vfrcopt Option to select the vegetation fraction data set c c = 0, DeFault built-in table conversion. c c = 1, select NESDIS green vegetation fraction data sets. c 0.144 x 0.144 degree resolution. c The compressed data set can be obtained from the CAPS c ftp directory: c c ftp://ftp.caps.ou.edu/pub/ARPS/arps40.data/arpssfc.data/gvegf.data.tar.gz c c nstyp Number of soil types for each grid point (only used c for sdatopt=1. See sdatopt for detail). c 0 < nstyp <= nstyps, which is the maximum number c defined in the program. c c nsmthsl Number of smoothing passes applied to veg, ln(roufns), c and lai. c c fgbgni = Beginning index (i) in x-dir of the foreground region. c fgendi = Ending index (i) in x-dir of the foreground region. c fgbgnj = Beginning index (j) in y-dir of the foreground region. c fgendj = Ending index (j) in y-dir of the foreground region. c c fgstyp = Soil type for the foreground. c fgvtyp = Vegetation type for the foreground. c fglai = Leaf Area Index for the foreground. c fgrfns = Surface roughness for the foreground. c fgveg = Vegetation fraction for the foreground. c c bgstyp = Soil type for the background. c bgvtyp = Vegetation type for the background. c bglai = Leaf Area Index for the background. c bgrfns = Surface roughness for the background. c bgveg = Vegetation fraction for the background. c c stypout = Flag for output of soil type. c vtypout = Flag for output of vegetation type. c vfrcdr = directory name for the vegetataion fraction data file. c laiout = Flag for output of Leaf Area Index. c rfnsout = Flag for output of surface roughness. c vegout = Flag for output of vegetation fraction. c ndviout = Flag for output of NDVI. c c drawval = Option for NCARG to draw integer value at each grid point c c stypfl = File name of soil classes. c vtypfl = File name of vegetation classes. c ndvifl = File name of NDVI. c lkupfl = File name of look-up table (for schmopt=4). c c Note: The soil and vegetation data are stored in CAPS ftp site, c ftp://ftp.caps.ou.edu/pub/ARPS/arps40.data/arpssfc.data/ c c####################################################################### &soil_veg_data schmopt = 1, sdatopt = 1, fstypfl = '/model/arpssfc.data/soil_1km.data', bstypfl = '/model/arpssfc.data/whsoil.data', vdatopt = 1, fvtypfl = '/model/arpssfc.data/naoge1_01l.img', bvtypfl = '/model/arpssfc.data/owe14d.data', ndatopt = 1, fndvifl = '/model/arpssfc.data/namay92ndl_1km.img', bndvifl = '/model/arpssfc.data/ndvi9005.data', vfrcopt = 1, vfrcdr = './', lkupfl = 'sfc_winter.tbl', nstyp = 3, nsmthsl = 3, fgbgni = 1, fgendi = 45, fgbgnj = 1, fgendj = 1, fgstyp = 13, bgstyp = 3, fgvtyp = 14, bgvtyp = 10, fglai = 0.00, bglai = 0.31, fgrfns = 0.001, bgrfns = 0.02, fgveg = 0.000, bgveg = 0.3, stypout = 1, vtypout = 1, laiout = 1, rfnsout = 1, vegout = 1, ndviout = 1, drawval = 0, &end c####################################################################### c c End of input file c c#######################################################################