ZETA_SETUP
&meso dt = 120. ; time step for this segment [sec] mstart = 1 ; starting time index (set to 1 for 1st segment) mend = 1000 ; ending time index mcn = 1000 ; frequency of restart writes [timesteps] mhn = 200 ; frequency of history writes [timesteps] cfil = 0.05 ; Robert (Asselin) filter coefficient ht = 18.e+3 ; height of rigid lid [m] epsi = 0.1 ; geometric scale factor in vertical coordinate transformation rlonmin = -30.0 ; longitude of western boundary mass point [deg] -180. < rlonmin < rlonmax < 180. rlonmax = 30.0 ; longitude of eastern boundary mass point [deg] -180. < rlonmin < rlonmax < 180. rlatmin = 42.5 ; latitude of south wall [deg] -90. < rlatmin < rlatmax < 90. rlatmax = 47.5 ; latitude of north wall [deg] -90. < rlatmin < rlatmax < 90. shear = 0.e-3 ; vertical shear in basic state zonal wind [/sec] bvfsq = 1.0e-4 ; Brunt-Vaisala frequency squared of tropospheric basic state [/sec**2] tho = 300. ; potential temperature at the ground [K] thy = 0.0 ; unused rkh = 5.0e+4 ; horizontal diffusion coefficient [m**2/sec] rkv = 20. ; vertical diffusion coefficient [m**2/sec] fnorm = 4.5 ; unused rkh4 = 0.00 ; fourth-order horizontal diffusion coefficient [m**4/sec] &end &mesomnt hmnt = 1000. ; maximum mountain height [m] efwx = 50.e+3 ; mountain e-folding width in zonal direction [m] efwy = 100.e+3 ; mountain e-folding width in merid direction [m] ywid = 5000.e+3 ; width of constant-height plateau in merid direction [m] imnt = 0. ; longitude of mountain peak [deg] jmnt = 45. ; latitude of mountain peak [deg] angle = 0. ; clockwise rotation angle of mountain [degrees] &end &mesoinit U0 = 20.0 ; surface wind [m/sec] eadyamp = 0.0 ; amplitude of initial Eady wave [mb] rsshramp = 0.0 ; amplitude of Rossby-Haurwitz wave [m/sec] lread = F ; read in initial condition? lmntint = F ; interpolate fields back to heights without topography? masscont = F ; apply mass continuity over topography initially? jetset = F ; set a jet structure for the basic flow? jetbeg = 20. ; latitude of jet beginning [deg] jetend = 70. ; latitude of jet end [deg] KSTRATOS = 99 ; z-index of tropopause KINTERF = 99 ; z-index of shear change &end &tropopause trpmid = 14.e+3 ; height of tropopause at middle of jet [m] trpamp = 2.0e+3 ; amplitude of tropopause variation at middle of jet [m] trpyscldeg = 7. ; y-scale of the tropopause variation [degrees] uvscl = 3.0e+3 ; vertical scale of u-variation [m] (d(Cor)/dzee ~ tanh((z-ztrop)/uvscl)) &end &mesobc ispwid = 11 ; width of sponge in zonal direction [index units] (0 => no sponge) jspwid = 0 ; width of polar sponge [index units] (0 => no sponge) bndwrit = F ; write out zonal boundaries for subsequent input? bndread = F ; read in zonal boundary conditions? &end &mesovmix rnumin = 0.1 ; minimum value of viscosity coefficient [m**2/sec] rnumax = 10.0 ; maximum value of viscosity coefficient [m**2/sec] richcrit = 2.0 ; critical Richardson number where rnuv=0.5(rnuvmax+rnuvmin) delrich = 2.0 ; half-width of richardson number transition zone &end &mesovort vortamp = 0.0e-4 ; amplitude of initial barotropic vortex [/s] vefwx = 2000.e+3 ; e-folding width in x of barotropic vortex [m] vefwy = 2000.e+3 ; e-folding width in x of barotropic vortex [m] ivort = 15 ; x-index of barotropic vortex center jvort = 34 ; y-index of barotropic vortex center &end &upforc upvamp = 1.0 ; amplitude of upstream meridional velocity (Eady wave) forcing [m/s] upvper = 100000. ; period of upstream meridional velocity forcing [s] &end &densprof rho0 = 1.0 ; density at the ground [kg/m**3] h_rho_inv = 1. ; anelastic switch (0 => Bousinesq, 1 => anelastic) &end &tracepnt ITRACE = 370 ; x-index of trace point JTRACE = 97 ; y-index of trace point KTRACE = 3 ; z-index of trace point &end &mesotraj mtn = 0 ; timesteps between subdomain writes ixtb = 49 ; x-index of western subdomain boundary ixte = 97 ; x-index of eastern subdomain boundary iytb = 56 ; y-index of southern subdomain boundary iyte = 168 ; y-index of northern subdomain boundary iztb = 1 ; z-index of bottom of subdomain boundary izte = 15 ; z-index of top of subdomain boundary &end &sclzonp mszp = 0 ; time between renormalizations [sec] (0 => no renormalization) pamp = 0.01 ; max zonal pert. amp of v after renormalization [m/sec] &end