c===================================================================
c
c The experiment namelist sets the casename for naming the 
c  GrADS output.  Max 72 characters and avoid using underscores.
c  Enter filename within single quotes (e.g., 'test.number.1')
c
c===================================================================

 &experiment
	casename = 'rolls',
 $

c===================================================================
c
c The grid_run namelist specifies model dimensions, integration
c  length, and plotting output interval
c
c  nx - number of horizontal grid points - max NXM in storage.txt
c     default is 101
c  nz - number of vertical grid points - max NZM in storage.txt
c     default is 84
c  dx - horizontal grid spacing (m; default = 1000.)
c  dz - vertical grid spacing (m; default = 250.)
c  dt - time step (s; default = 1.0)
c  timend - integration stop time (s; default = 7200.)
c  plot - plotting interval (s; default = 300.)
c     * if interval < 60 sec, GrADS will report time incorrectly
c===================================================================

 &grid_run
	nx = 301,
	nz = 52,
	dx = 500.,
	dz = 200.,
	dt = 0.5,
	timend = 19200.,
	plot = 300.,
 $

c===================================================================
c 
c The framework namelist specifies whether the model is compressible
c  or anelastic, the sound speed and whether pressure decompositions
c  are performed
c
c  ipressure - output pressure decompositions if 1 (default is 0)
c     * can make model integrations much longer in some cases
c  ianelastic (1 = anelastic model, 0 - compressible [default])
c  csnd - sound speed (m$s; default = 50.)
c     * 50 m$s is very unrealistic but permits fast integration
c
c===================================================================

 &framework
	ipressure = 1,
	ianelastic = 0,
	csnd = 50.,
 $

c===================================================================
c 
c The numerics namelist sets up numerical aspects of the
c  integration, including diffusion and boundary handling
c
c  byteswap (1 = byteswapping needed; default is 0)
c  cstar - lateral boundary wave speed (m$s; default = 50.)
c  dkx - horizontal diffusion coefficient (m^2$s; default = 500.)
c   * rethink
c  dkz - vertical diffusion coefficient (m^2$s; default = 100.)
c   * rethink
c  eps - time smoothing coefficient (ndim; default = 0.005)
c
c===================================================================

 &numerics
	byteswap = 0,
	cstar = 100.,
	dkx = 1250.,
	dkz = 20.,
	eps = 0.005,
 $

c===================================================================
c 
c The environ namelist sets up the model initial state
c
c   bvpbl - PBL tropospheric BV freq (1$s; default = 0.01)
c   pbld - PBL depth (m; default = 2000.)
c   bvtropo - free tropospheric BV freq (1$s; default = 0.01)
c   tropo - tropopause height (m; default = 12000.)
c   bvstrat - stratospheric BV freq (1$s; default = 0.02)
c   psurf - surface pressure (mb; default = 965.)
c   usurf - surface wind speed (m$s; default = 0.)
c   shear1 - vertical shear for first layer (1$s; default = 0.)
c   depth1 - thickness of first layer (m; default = 3000.)
c   shear2 - vertical shear for second layer (1$s; default = 0.)
c   depth2 - thickness of second layer (m; default = 1500.)
c   shear3 - vertical shear above second layer (1$s; default = 0.)
c    * this layer extends to model top
c
c===================================================================

 &environ
	bvpbl = 0.01,
	pbld = 2000.,
	bvtropo = 0.01,
	tropo = 12000.,
	bvstrat = 0.02,
	psurf = 965.,
	usurf = 0.,
	shear1 = 0.,
	sdepth1 = 1800.,
	shear2 = 0.,
	sdepth2 = 2500.,
	shear3 = 0.,
 $

c===================================================================
c
c The thermal namelist specifies characteristics of an impulsive
c  initial thermal bubble
c
c  ithermal (1 = turn thermal on; default is 0)
c  delt - amplitude of thermal (K)
c  radx - horizontal radius of thermal (m; default = 8000.)
c  radz - vertical radius of thermal (m; default = 4000.)
c  zcnt - center height of thermal (m; default = 3000.)
c
c===================================================================

 &thermal
	ithermal = 0,
	delt = 3.,
	radx = 8000.,
	radz = 4000.,
	zcnt = 3000.,
 $

c===================================================================
c
c The surface_flux namelist sets up at least part of a domain
c  to represent a heated surface, to heat atmosphere from below
c
c  ishflux (1 = turn surface heat flux on; default is 0)
c  tdelt - initial ground-air T difference (K; default = 12)
c  icoast - gridpt location of coastline (0 for all land; default = 30)
c  cdh - effective heat flux coefficient (ndim; default = 7.2e-3)
c  irand - (1 = impose randomness on surface heat flux; default is 0)
c
c===================================================================

 &surface_flux
	ishflux = 1,
	tdelt = 12.,
	icoast = 0,
	cdh = 7.2e-3,
	irand = 1,
 $

c===================================================================
c
c The streamfunction namelist sets up a momentum source used
c  to excite gravity waves.  The source may be isolated or repeated,
c  oscillatory or steady
c
c  istrfcn (1 = turn streamfunction forcing on; default is 0)
c  s_repeat(1 for repeated source, 0 for single source; default 0)
c  s_ampl - amplitude of momentum source (kg$m/s/s; default = 40.)
c  s_naught - height of source center (m; default = 6000.)
c  s_hwavel - horizontal wavelength of source (m; default = 40000.)
c  s_vwavel - vertical wavelength of source (m; default = 18000.)
c  s_period - oscillation period (s; default = 1200.)
c  * s_period = 0 makes source steady
c
c===================================================================

 &streamfunction
	istrfcn = 0,
	s_repeat = 1,
	s_ampl = .40,
	s_znaught = 1000.,
	s_hwavel = 10000.,
	s_vwavel = 6000.,
	s_period = 0.,
 $

c===================================================================
c
c The atmos_heat_source namelist sets up a free tropospheric heat source
c  used to excite gravity waves.  The source may be oscillatory or steady
c
c  ihsrc (1 = turn tropospheric heat source on; default is 0)
c  h_ampl - amplitude of heat source (K$s; default = 0.075)
c  h_radius_x - horizontal radius of heat source (m; default = 3000.)
c  h_radius_z - vertical radius of heat source (m; default = 3000.)
c  h_center_z - height of heat source center (m; default = 3000.)
c  h_freq - frequency for heat source oscillation (1$s; default = 0.005)
c  h_modes - number of vertical modes (ndim, max 2; default = 2)
c
c===================================================================

 &atmos_heat_source
	ihsrc = 0,
	h_ampl = 0.025,
	h_radius_x = 3000.,
	h_radius_z = 3000.,
	h_center_z = 3000.,
	h_freq = 0.000,
	h_modes = 2,
 $

c===================================================================
c
c The rotunno_seabreeze namelist implements a lower tropospheric
c  heat source following Rotunno (1983), useful for long-term
c  integrations of the sea-land-breeze circulation
c 
c  iseabreeze (1 = turn Rotunno heat source on; default is 0)
c  sb_ampl - amplitude of heat source (K$s; default = 0.000175)
c  sb_x0 - controls heat source shape at coastline (m; default = 1000.)
c  sb_z0 - controls heat source shape at coastline (m; default = 1000.)
c  sb_period - period of heating, in days (default = 1.0)
c  sb_latitude - latitude for experiment (degrees; default = 60.)
c  sb_linear (1 = linearize model; default = 1)
c
c===================================================================

 &rotunno_seabreeze
	iseabreeze = 0,
	sb_ampl = 0.000175,
	sb_x0 = 1000.,
	sb_z0 = 1000.,
	sb_period = 1.0,
	sb_latitude = 60.,
	sb_linear = 1,
 $

c===================================================================
c
c The cooling_zone namelist implements a lower tropospheric
c  cooling source following Fovell and Tan (2000) or an impulsive
c  block of cold air of specified dimensions
c 
c  icoolzone (1 = Fovell-Tan storm-adaptive cooling zone,
c             2 = impulsive cold block; default is 0)
c  cz_ampl - cooling zone amplitude (K; default = 4.5)
c  * cz_ampl is a positive number for cooling
c  cz_rightedge - zone right (east) edge (%; default = 0.60)
c  * expressed as a domain fraction; e.g., 0.60 = 60% of way across domain
c  cz_depth - cooling zone depth (m; default = 2500.)
c  cz_width - width of cooling zone (m; default = 46000.)
c  cz_coolrate - cooling rate (K$s; default = 600.)
c
c===================================================================

 &cooling_zone
	icoolzone = 0,
	cz_ampl = 9.0,
	cz_rightedge = 0.30,
	cz_depth = 2000.,
	cz_width = 40000.,
	cz_coolrate = 600.,
 $