&problemsize_namelist
 n_r = 64
 n_theta = 192
 nprow = 4
 npcol = 4
 aspect_ratio = 0.5d0
 shell_depth = 1.0d0
/
&numerical_controls_namelist
/
&physical_controls_namelist
 rotation  = .True.
 magnetism = .false.
 viscous_heating = .false.
 ohmic_heating = .false.
 centrifugal_force = .true.
/
&temporal_controls_namelist
 max_time_step = 1.0d-4
 !max_iterations = 40000
 max_simulated_time = 2.5d0
 checkpoint_interval = 10000
 quicksave_interval = 2000
 num_quicksaves=3
 cflmin = 0.4d0
 cflmax = 0.6d0
/
&io_controls_namelist
/
&output_namelist

! New equatorial slices output
equatorial_values    = 1,2,3, 501  ! velocity components
equatorial_frequency = 5000
equatorial_nrec      = 2

!New meridional slices output.  You must also specify the number of phi-indices
!Valid values range from 1 through nphi = 2*n_theta
!The indices below range from phi = 0 up to phi = pi 
meridional_values    = 1,2,3 ! radial and phi components of velocity; temperature
meridional_frequency = 5000
meridional_nrec      = 2
meridional_indices_nrm = 0.7

shellslice_levels_nrm    = 0.1, 0.5, 0.9
shellslice_values    = 1,2,3, 501 ! velocity components
shellslice_frequency = 5000
shellslice_nrec      = 2

shellspectra_levels_nrm    = 0.1, 0.5 , 0.9
shellspectra_values    = 1,2,3  ! velocity and temperature
shellspectra_frequency = 5000
shellspectra_nrec      = 2

azavg_values = 1,2,3,201,202,501,1440,1470, 1201,1202,1203, 1219,1220,1221, 1228,1229,1230, 1237,1238,1239, 1263,1264 ! same as above + r- and theta- mass flux
azavg_frequency = 200
azavg_nrec = 50

point_probe_values = 1,2,3, 1201,1202,1203,1216, 1219,1220, 1221, 1228,1229,1230
point_probe_r_nrm = 0.2, 0.6, 0.9
point_probe_theta_nrm= 0.35
point_probe_phi_nrm = 0.21, 0.72
point_probe_frequency=10000
point_probe_nrec=500
point_probe_cache_size=100


! velocity, temperature, energy fluxes, and Kinetic Energy
shellavg_values = 1,2,3,501, 1440, 1470, 401, 402, 403, 404, 405, 406, 407, 408, 409, 410, 411, 412
shellavg_frequency = 200
shellavg_nrec = 50

! Kinetic energy, Mean KE, Diff-Rot KE, and Convective KE
globalavg_values = 401, 402, 403, 404, 405, 406, 407, 408, 409, 410, 411, 412
globalavg_frequency = 100
globalavg_nrec = 100

sph_mode_ell = 2,4,8
sph_mode_levels_nrm = 0.4, 0.7, 0.9
sph_mode_values = 1,2,3
sph_mode_frequency = 1000
sph_mode_nrec = 50

full3d_values = 501 ! temperature
full3d_frequency = 1000000
/

&Boundary_Conditions_Namelist
no_slip_boundaries = .true.
strict_L_Conservation = .false.
dtdr_bottom = 0.0d0
T_Top    = 0.0d0
T_Bottom = 1.0d0
fix_tvar_top = .true.
fix_tvar_bottom = .true.
/

&Initial_Conditions_Namelist
init_type = 7
temp_amp = 0.01d0
temp_w = 0.01d4
restart_iter = 0
conductive_profile = .true.
/
&Test_Namelist
/
&Reference_Namelist
Ekman_Number = 1.0d-3
Rayleigh_Number = 1.0d6
Prandtl_Number = 1.0d0
Magnetic_Prandtl_Number = 5.0d0
reference_type = 1
heating_type = 0      ! No heating
gravity_power = 0.0d0  ! g ~ radius
froude_number=0.1d0
/
&Transport_Namelist
/