try to turn on lightbulb for SASW using radice ICs

inputs/standing_accretion_shock.pin

@@ -32,24 +32,27 @@ variables = p.density,  &
             temperature,&
             p.ye,       &
             g.c.coord,  &
-            g.n.coord
+            g.n.coord,  &
+	    flux_divergence, &
+	    src_terms
+
 
 file_type   = hdf5      # Tabular data dump
-dt          = 10        # time increment between outputs
+dt          = 100       # time increment between outputs
 
 <parthenon/time>
 nlim        = -1        # cycle limit
-tlim        = 2000      # time limit
+tlim        = 10000      # time limit
 integrator  = rk2       # time integration algorithm
 ncycle_out  = 100       # interval for stdout summary info
-dt_init_fact = 1.e-6
+dt_init_fact = 1.e-8
 
 <parthenon/mesh>
 nghost = 4
 #refinement = adaptive
 #numlevel = 3
 
-nx1         = 1024           # Number of zones in X1-direction
+nx1         = 512           # Number of zones in X1-direction
 x1min       = 30            # minimum value of X1 | rPNS | x1min = ln R (code units), R = 60 (km) => R = 31.2492 (code units) => x1min = 3.44199
 x1max       = 230           # maximum value of X1 | rmax | x1max => 472.5 (km)
 ix1_bc      = user          # Inner-X1 boundary condition flag
@@ -79,21 +82,21 @@ field = c.c.bulk.rho
 method = derivative_order_1
 max_level = 1
 
-<eos>
-type = IdealGas
-Gamma = 1.33
-Cv = 3.0  
-
 #<eos>
-#type = StellarCollapse
-#filename = /usr/projects/w22_stars/EoS/Hempel_SFHoEOS_rho222_temp180_ye60_version_1.1_20120817.h5
-#use_sp5 = false
+#type = IdealGas
+#Gamma = 1.33
+#Cv = 3.0  
+
+<eos>
+type = StellarCollapse
+filename = /usr/projects/w22_stars/EoS/Hempel_SFHoEOS_rho222_temp180_ye60_version_1.1_20120817.h5
+use_sp5 = false
 
 <fixup>
-enable_floors = true
+enable_floors = false
 enable_mhd_floors = false
-enable_ceilings = true
-enable_flux_fixup = true
+enable_ceilings = false
+enable_flux_fixup = false
 enable_c2p_fixup = true
 floor_type = ConstantRhoSie
 rho0_floor = 1.e-12
@@ -106,7 +109,7 @@ sie_ceiling = 1e-2
 hydro = true
 he = false
 3t = false
-rad = false
+rad = true
 
 <fluid>
 xorder = 2
@@ -131,4 +134,12 @@ fluid_density_cgs = 3.599e17   # 1.e9     # gives mass = 7.07e24 g
 Mdot = 0.2             # Msun / sec
 rShock = 200           # km
 target_mach = -100      # target Mach number for upstream flow in preshock region 
-model_filename = 1d_initial_conditions_radice.txt
+model_filename = 1d_initial_conditions_radice_10k.txt
+
+<radiation>
+method=cooling_function
+do_liebendorfer=false
+do_lightbulb=true
+lum=2.5
+do_gain_reducer=false
+always_gain=false