diff --git a/Exec/DevTests/TemperatureSource/input_sounding b/Exec/DevTests/TemperatureSource/input_sounding
new file mode 100644
index 000000000..52b74b123
--- /dev/null
+++ b/Exec/DevTests/TemperatureSource/input_sounding
@@ -0,0 +1,3 @@
+1000.0 300.0 0.0
+   0.0 300.0 0.0 10.0 0.0
+1000.0 300.0 0.0 10.0 0.0
diff --git a/Exec/DevTests/TemperatureSource/neutral_column.i b/Exec/DevTests/TemperatureSource/neutral_column.i
new file mode 100644
index 000000000..35a2e131a
--- /dev/null
+++ b/Exec/DevTests/TemperatureSource/neutral_column.i
@@ -0,0 +1,49 @@
+# ------------------  INPUTS TO MAIN PROGRAM  -------------------
+stop_time = 999.9
+max_step = 20
+
+amrex.fpe_trap_invalid = 1
+
+fabarray.mfiter_tile_size = 1024 1024 1024
+
+# PROBLEM SIZE & GEOMETRY
+geometry.prob_extent    =   40.  40.  640.
+amr.n_cell              =   4    4    64
+
+geometry.is_periodic = 1 1 0
+
+zlo.type = "SlipWall"
+zhi.type = "SlipWall"
+
+# INITIALIZATION
+erf.init_type = input_sounding
+erf.init_sounding_ideal = true
+
+# TIME STEP CONTROL
+erf.fixed_dt                    = 0.1
+
+# DIAGNOSTICS & VERBOSITY
+amr.v               = 1     # verbosity in Amr.cpp
+erf.v               = 1     # verbosity in ERF.cpp -- needs to be 1 to write out data_log files
+erf.sum_interval    = 1     # timesteps between computing mass
+
+# REFINEMENT / REGRIDDING
+amr.max_level       = 0       # maximum level number allowed
+
+# CHECKPOINT FILES
+erf.check_file      = chk        # root name of checkpoint file
+erf.check_int       = -1         # number of timesteps between checkpoints
+
+# PLOTFILES
+erf.plot_file_1     = plt       # prefix of plotfile name
+erf.plot_int_1      = 5         # number of timesteps between plotfiles (DEBUG)
+erf.plot_vars_1     = density x_velocity y_velocity z_velocity pressure theta
+
+# SOLVER CHOICES
+erf.use_gravity = true
+erf.use_coriolis = false
+
+erf.molec_diff_type = "None"
+erf.les_type = "None"
+
+erf.add_custom_rhotheta_forcing = true