make this an option

Author: zingale

Diff for: Exec/science/massive_star/_prob_params

@@ -1,7 +1,10 @@
 model_name      character       ""        y
 
+interpolate_pres   int          0         y
+
 perturb_model   int             0         y
 velpert_amplitude  real         1.e5      y
 velpert_scale   real            1.e7      y
 
 isi28           int            -1
+

Diff for: Exec/science/massive_star/problem_initialize_state_data.H

@@ -33,7 +33,7 @@ void problem_initialize_state_data (int i, int j, int k,
 
     state(i,j,k,URHO) = interpolate(dist, model::idens);
     state(i,j,k,UTEMP) = interpolate(dist, model::itemp);
-    state(i,j,k,UPRES) = interpolate(dist, model::ipres);
+    Real pres = interpolate(dist, model::ipres);
     for (int n = 0; n < NumSpec; n++) {
         state(i,j,k,UFS+n) = amrex::max(interpolate(dist, model::ispec+n), small_x);
     }
@@ -123,17 +123,21 @@ void problem_initialize_state_data (int i, int j, int k,
     eos_t eos_state;
     eos_state.rho = state(i,j,k,URHO);
     eos_state.T = state(i,j,k,UTEMP);
-    eos_state.p = state(i,j,k,UPRES);
+    eos_state.p = pres;
     for (int n = 0; n < NumSpec; n++) {
         eos_state.xn[n] = state(i,j,k,UFS+n);
     }
     for (int n = 0; n < NumAux; n++) {
         eos_state.aux[n] = state(i,j,k,UFX+n);
     }
 
-    eos(eos_input_rp, eos_state);
+    if (problem::interpolate_pres == 1) {
+        eos(eos_input_rp, eos_state);
+        state(i,j,k,UTEMP) = eos_state.T;
+    } else {
+        eos(eos_input_rt, eos_state);
+    }
 
-    state(i,j,k,UTEMP) = eos_state.T;
     state(i,j,k,UEINT) = state(i,j,k,URHO) * eos_state.e;
     state(i,j,k,UEDEN) = state(i,j,k,URHO) * eos_state.e;