&gridnl
nelx = 6,
nely = 6,
nelz = 5,
nop = 4,
xdims = 0,400000.0
ydims = 0,400000.0
ztop = 300.0
!Dimension of the domain: (x,y) in (min,max); z in (0,ztop)
x_boundary = 4, 4 !Front and Back (x=-1 and x=+1)
y_boundary = 4, 4 !Left and Right (y=-1 and y=+1)
z_boundary = 4, 4 !Bottom and Top (z=-1 and z=+1)
!These denote constant x,y, or z.
!E.g. z_boundary(1) is the bottom and z_boundary(2) is the top
!BCs are: 0=do nothing; 10=Dirichlet; 3=periodic; 4=NFBC; 6=NRBC
!(see MOD_GLOBAL_GRID for more info)
geometry_type='sphere_hex',
!cube=3D Box;
!sphere_hex=spherical cubed-sphere domain;
!sphere_ico=spherical icosahedral domain
decomp_type='nseam',
!metis3d=3D metis partitioning;
!metis2d=2D metis partitioning;
!nseam=geometric for cubed-sphere;
!geom1=geometric along X for Cube/Cartesian grid
!geom2=geometric along XY for Cube/Cartesian grid
!geom3=geometric along XYZ for Cube/Cartesian grid
!geom4=geometric along XYZ (nproc_z controls Z direction) for Cube/Cartesian grid
!geom5=geometric along Y for Cube/Cartesian grid
!geom6=geometric along YZ for Cube/Cartesian grid
!geom7=geometric along XZ for Cube/Cartesian grid
nproc_z = 2,
!Only used with GEOM4: it controls the decomposition along z. Along XY, it
!will use the square root of nproc/nproc_z
/
&static_database_nl
use_landuse = .true.
use_landmask = .true.
use_terrain = .true.
use_georeference = .true.
/
&physics_nl
mixing_method=0,
!=0; no sub-gridscale mixing
!=1; 1st order closure sub-gridscale mixing (Lilly-Smagorinsky)
surface_fluxes = .false.
/
&convect_nl
run_shallow_q = .false.
/
&radiation_nl
run_shortwave = .false. ! call the shortwave radiation
run_longwave = .false. ! call the longwave radiation
time_sw = 1.0 ! time interval for calling shortwave (s) (not used currently)
time_lw = 1.0 ! time interval for calling longwave (s) (not used currently)
timeofday = 12. ! local time (hours, 0-24, 12=noon)
juldati = 18.0 ! julian date (0-360/361)
reflat = 10.0 ! ref latitude (degrees)
/
&input
dt= 0.01,
time_initial = 0,
time_final = 0.1,
time_restart = 100,
icase=9,
!case0=advection;
!case1 = igw;
!case2=rtb;
!case20=large rtb;
!case200=rtb with tracer;
!case21=3D bubble;
!case3=2d dc;
!case31=3d dc;
!case4=2d lhm (ridge, x-direction);
!case40=2d Nonlinear large Mountain (ridge, x-direction);
!case41=3D Linear Mountain;
!case42=2d lhm (ridge, y-direction);
!case43=3D NonLinear Mountain;
!case5=3D bubble on Sphere;
!case6=Baroclinic Instability on Sphere;
!case7=Baroclinic Instability in Channel: Balanced State
!case71=Baroclinic Instability in Channel: Perturbed State
!case8=Acoustic Wave on Sphere
!case9=Gravity Wave on Sphere
!case1000=Shallow Cumulus Case
ti_method = 'ark2',
!rk; (RK method of order KSTAGES)
!rk35; (Optimized RK35)
!am2; (am2 with ab3)
!bdf22; (bdf2 with 2nd order explicit scheme)
!bdf23; (bdf2 with 3rd order explicit scheme)
!bdf3; (bdf3 with 3rd order explicit scheme)
!lf2;
!ark2a; (2nd order 2-stage method, 1-stage implicit)
!ark2s; (2nd order 2-stage method, 1-stage implicit)
!ark2o; (Original 2nd order 3-stage method, 1-stage implicit)
!ark2n; (New 2nd order 3-stage method, 1-stage implicit)
!ark3;
!ark4;
!ark5
si_method = 'no_schur',
!schur;
!no_schur
si_dimension = '3d',
!1d IMEX in the vertical/radial direction;
!3d IMEX
Iter_Type = 'GMRES', ! = 'GMRES' means run MGS-based GMRES
! = 'PROJ' means run MGS-based GMRES + Projection
! = 'CHEBY' means run Chebyshev (dot-product free)
! = 'BICGS' means run stabilized BiCG method (TBD)
! = 'RICHA' means run Richardson Extrapolation
Proj_Size=30, !=Number of Maximum Projection vectors: Only works for
!Iter_Type=PROJ
Focal_Dist_Adj = 0.981, ! Used only by Chebyshev method and
! adjusts elliptical hull focal distance
! to obtain optimal convergence rate
precon_mode = 'SKIP',
! = 'SKIP' means do not build and do not apply (means precon_order=0)
! = 'BUILD' means:
! 1) compute max/min eigVals and precon coeffs,
! 2) save above data to file, and abort run.
! = 'APPLY' means read above data and run model.
! = 'BOTH' means BUILD and APPLY w/o stopping model
precon_order = 0, ! = 0 means no preconditioner is made/read
! > 0 mean precon K = p(A) where p is a polynomial
! of the order shown and of the type specified below:
precon_type = 'S', ! S = spectrally-optimized polynomial coefficients
! N = Neumann polynomial coefficients
! C = Chebyshev polynomial coefficients
PBNO_pNorm = 4, ! sets p-norm exponent (used in with precon_type = S)
Imag_Wt = 1.0e0, ! Imag Component weight factor (used only with p-norm)
Hmat_Size = 150, ! Size of Hessenberg matrix used to estimate system spectrum
lf2_filter=0., !LF2 Robert-Asselin Filter
lf2_alpha=0.5, !LF2 Off-Centering (0.5 is fully centered)
kstages=5,
solver_type = 'direct',
!=direct uses a direct solver for SI_DIMENSION=1D and SI_METHOD=SCHUR form
!=iterative uses an iterative solver for SI_DIMENSION=1D and SI_METHOD=SCHUR form
solver_tol=1e-5,
!=1e-2 for IMEX-3D
!<=1e-9 for IMEX-1D (e.g., case 21 and case 7)
gamma1=+1,
!-1=relative residual (works best for No Schur);
!+1=absolute residual (works best for Schur)
delta=1,
!-1=No Semi-implicit;
!0=No Semi-implicit but go through algorithm anyway (for debugging);
!1=Semi-Implicit
xmu=0e-2,
ifilter=0,
filter_weight_type='erf',
filter_basis_type='modal',
fname_root='new_case9_visc=0.1',
out_type='nc',
!gks2d=ascii 2d slice;
!gks3d=ascii 3d output;
!nc=netcdf 3d output;
!vtk=Visual ToolKit (Paraview) output
!none=No Output
visc=0.1,
visc_dd=0e4, !divergence damping
nlaplacian=1,
iadapt_timestep=0,
!=1: adapt DT automatically based on stability conditions of each TI
!=0: use given DT
lprint_diagnostics=T,
!T=true=print diagnostics
!F=false=do NOT print diagnostics
lphysics=F,
!T=true=Turn ON PHYSICS
!F=false=Turn OFF PHYSICS
/