Stability upgrades

AsterAnalysis/interface.ccl

@@ -8,14 +8,10 @@ USES INCLUDE HEADER: aster_fd.hxx aster_interp.hxx aster_utils.hxx
 
 CCTK_REAL mhd_analysis TYPE=gf CENTERING={ccc} tags='checkpoint="no"'
 {
-divA, A_norm, divB, B_norm, b2small, mean_B_numerator, mean_B_denominator, mhd_press_ratio, mhd_energy_ratio, magnetic_energy, kinetic_energy, poynting_flux_x, poynting_flux_y, poynting_flux_z, poynting_flux_norm, poynting_vector_x, poynting_vector_y, poynting_vector_z, poynting_scalar, lambda_MRI_x, lambda_MRI_y, lambda_MRI_z
+divA, A_norm, divB, mean_B_numerator, mean_B_denominator, mhd_press_ratio, mhd_energy_ratio, magnetic_energy, kinetic_energy, poynting_flux_x, poynting_flux_y, poynting_flux_z, poynting_flux_norm, poynting_vector_x, poynting_vector_y, poynting_vector_z, poynting_scalar, lambda_MRI_x, lambda_MRI_y, lambda_MRI_z
 } "Variables used for analysis purposes when magnetic fields are active"
 
 CCTK_REAL mhd_spherical TYPE=gf CENTERING={ccc} tags='checkpoint="no"'
 {
 Brho, Btheta, Bphi, B_Poloidal_norm, B_Toroidal_norm, mean_B_Tor_numerator, mean_B_Tor_denominator, mean_B_Pol_numerator
 } "For decomposition of B field into toroidal (phi) and poloidal parts"
-
-public:
-
-CCTK_REAL w_lorentz TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Lorentz factor"

AsterAnalysis/schedule.ccl

@@ -13,15 +13,16 @@ if (compute_MHD) {
     READS: HydroBaseX::vel(everywhere)
     READS: HydroBaseX::Bvec(everywhere)
     READS: AsterX::zvec(everywhere)
+    READS: AsterX::B_norm(everywhere)
+    READS: AsterX::b2small(everywhere)
+    READS: AsterX::w_lorentz(everywhere)
     READS: AsterX::Avec_x(everywhere)
     READS: AsterX::Avec_y(everywhere)
     READS: AsterX::Avec_z(everywhere)
     READS: AsterX::dBx_stag(everywhere)
     READS: AsterX::dBy_stag(everywhere)
     READS: AsterX::dBz_stag(everywhere)
-    WRITES: w_lorentz(interior)
     WRITES: mhd_analysis(interior)
-    SYNC: w_lorentz
   }
   "Compute MHD analysis quantities"
 }
@@ -30,7 +31,7 @@ if (compute_MHD_spherical) {
   SCHEDULE AsterAnalysis_MHD_Spherical AT analysis AFTER AsterAnalysis_MHD
   {
     LANG: C
-    READS: w_lorentz(interior)
+    READS: AsterX::w_lorentz(interior)
     READS: ADMBaseX::metric(interior)
     READS: HydroBaseX::Bvec(interior)
     READS: HydroBaseX::rho(interior)

AsterAnalysis/src/analysis_mhd.cxx

@@ -56,13 +56,11 @@ extern "C" void AsterAnalysis_MHD(CCTK_ARGUMENTS) {
           v_up(1) = vely(p.I);
           v_up(2) = velz(p.I);
           v_low = calc_contraction(g_avg, v_up);
-          w_lorentz(p.I) = calc_wlorentz(v_low, v_up);
 
         } else {
 
           const vec<CCTK_REAL, 3> z_up{zvec_x(p.I), zvec_y(p.I), zvec_z(p.I)};
           const vec<CCTK_REAL, 3> z_low = calc_contraction(g_avg, z_up);
-          w_lorentz(p.I) = calc_wlorentz_zvec(z_low, z_up);
           v_up = z_up / w_lorentz(p.I);
           v_low = z_low / w_lorentz(p.I);
         }
@@ -76,7 +74,6 @@ extern "C" void AsterAnalysis_MHD(CCTK_ARGUMENTS) {
         const vec<CCTK_REAL, 3> B_up{Bvecx(p.I), Bvecy(p.I), Bvecz(p.I)};
         const vec<CCTK_REAL, 3> B_low = calc_contraction(g_avg, B_up);
         const CCTK_REAL B2big = calc_contraction(B_up, B_low);
-        B_norm(p.I) = sqrt(B2big);
 
         /* cell-centered A^i and A_i */
         const vec<CCTK_REAL, 3> A_low{calc_avg_e2c<0>(gf_Avec(0), p),
@@ -93,10 +90,6 @@ extern "C" void AsterAnalysis_MHD(CCTK_ARGUMENTS) {
         const vec<CCTK_REAL, 3> bs =
             (B_up / w_lorentz(p.I) + alp_avg * bs0 * velshift);
 
-        /* b^2 */
-        b2small(p.I) = (calc_contraction(B_up, B_low) + pow2(alp_avg * bs0)) /
-                       pow2(w_lorentz(p.I));
-
         /* inverse Beta and magnetization */
         mhd_press_ratio(p.I) = b2small(p.I) / (2.0 * press(p.I));
         mhd_energy_ratio(p.I) = b2small(p.I) / (2.0 * rho(p.I));

AsterX/interface.ccl

@@ -13,16 +13,19 @@ USES INCLUDE HEADER: aster_fd.hxx aster_interp.hxx aster_utils.hxx
 
 PUBLIC:
 
-CCTK_REAL dens TYPE=gf CENTERING={ccc} TAGS='rhs="densrhs" dependents="HydroBaseX::rho TmunuBaseX::eTtt AsterX::flux_x AsterX::flux_y AsterX::flux_z"' "Conserved rest mass density"
+# ----- evolved state variables
 
-CCTK_REAL mom TYPE=gf CENTERING={ccc} TAGS='rhs="momrhs" dependents="HydroBaseX::vel TmunuBaseX::eTti"'
+# Note that the prolongation method and order are hard-coded and will
+# supersede parfile choices
+CCTK_REAL cons_vector TYPE=gf CENTERING={ccc} TAGS='rhs="rhs_vector" dependents="HydroBaseX::rho HydroBaseX::eps HydroBaseX::press HydroBaseX::temperature HydroBaseX::Ye HydroBaseX::Bvec HydroBaseX::vel HydroBaseX::entropy TmunuBaseX::eTtt TmunuBaseX::eTti TmunuBaseX::eTij AsterX::flux_x AsterX::flux_y AsterX::flux_z" prolongation_type="minmod" prolongation_order=1' "Conserved quantities during AsterX evolution"
 {
-  momx, momy, momz
-} "Conserved momentum density"
+  dens, momx, momy, momz, tau, DYe, DEnt
+} "Conserved density, momenta, internal energy, electron fraction and entropy"
 
-CCTK_REAL tau TYPE=gf CENTERING={ccc} TAGS='rhs="taurhs" dependents="HydroBaseX::eps HydroBaseX::press HydroBaseX::temperature HydroBaseX::Ye HydroBaseX::Bvec TmunuBaseX::eTij"' "Conserved internal energy density"
-
-CCTK_REAL DYe TYPE=gf CENTERING={ccc} TAGS='rhs="DYe_rhs" dependents="HydroBaseX::Ye AsterX::flux_x AsterX::flux_y AsterX::flux_z"' "Conserved electron fraction"
+CCTK_REAL rhs_vector TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "RHSs of conserved variables"
+{
+  densrhs, momxrhs, momyrhs, momzrhs, taurhs, DYe_rhs, DEntrhs
+} "Conserved density, momenta, internal energy, electron fraction and entropy RHS"
 
 # DEnt is the conserved "evolved" entropy.
 # The corresponding primitive is stored in HydroBaseX::entropy.
@@ -34,7 +37,8 @@ CCTK_REAL DYe TYPE=gf CENTERING={ccc} TAGS='rhs="DYe_rhs" dependents="HydroBaseX
 # refer to "physical" and "evolved" entropy, respectively. In this 
 # evolution thorn we always refer to DEnt and entropy to describe 
 # the "evolved" quantities. 
-CCTK_REAL DEnt TYPE=gf CENTERING={ccc} TAGS='rhs="DEntrhs" dependents="HydroBaseX::entropy AsterX::flux_x AsterX::flux_y AsterX::flux_z"' "Advected entropy density"
+
+# ----- fluxes
 
 CCTK_REAL flux_x TYPE=gf CENTERING={vcc} TAGS='checkpoint="no"'
 {
@@ -66,42 +70,47 @@ CCTK_REAL flux_z TYPE=gf CENTERING={ccv} TAGS='checkpoint="no"'
   fzBx fzBy fzBz
 } "Fluxes in z direction"
 
+# ----- flux limiter
+
 CCTK_REAL theta_tot TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "min(thetas) taken over all cell faces"
 
 CCTK_REAL theta_x TYPE=gf CENTERING={vcc} TAGS='checkpoint="no"' "flux limiter in x"
 CCTK_REAL theta_y TYPE=gf CENTERING={cvc} TAGS='checkpoint="no"' "flux limiter in y"
 CCTK_REAL theta_z TYPE=gf CENTERING={ccv} TAGS='checkpoint="no"' "flux limiter in z"
 
-CCTK_REAL densrhs TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Conserved rest mass density RHS"
+# ----- shock detector
 
-CCTK_REAL momrhs TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"'
-{
-  momxrhs, momyrhs, momzrhs
-} "Conserved momentum density RHS"
-
-CCTK_REAL taurhs TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Conserved internal energy density RHS" 
+CCTK_REAL LOflag TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Shock detector switch"
+CCTK_REAL etac TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Shock detector measure"
 
-CCTK_REAL DYe_rhs TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Conserved electron fraction RHS"
-CCTK_REAL DEntrhs TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Advected entropy density RHS"
+# ----- staggered vector potential (edges)
 
-CCTK_REAL Avec_x TYPE=gf CENTERING={cvv} TAGS='rhs="Avec_x_rhs"' "x-component of vector potential"
-CCTK_REAL Avec_y TYPE=gf CENTERING={vcv} TAGS='rhs="Avec_y_rhs"' "y-component of vector potential"
-CCTK_REAL Avec_z TYPE=gf CENTERING={vvc} TAGS='rhs="Avec_z_rhs"' "z-component of vector potential"
+# Note that the prolongation method and order are hard-coded and will
+# supersede parfile choices
+CCTK_REAL Avec_x TYPE=gf CENTERING={cvv} TAGS='rhs="Avec_x_rhs" prolongation_type="hermite" prolongation_order=5' "x-component of vector potential"
+CCTK_REAL Avec_y TYPE=gf CENTERING={vcv} TAGS='rhs="Avec_y_rhs" prolongation_type="hermite" prolongation_order=5' "y-component of vector potential"
+CCTK_REAL Avec_z TYPE=gf CENTERING={vvc} TAGS='rhs="Avec_z_rhs" prolongation_type="hermite" prolongation_order=5' "z-component of vector potential"
 
 CCTK_REAL Avec_x_rhs TYPE=gf CENTERING={cvv} TAGS='checkpoint="no"' "x-component of vector potential RHS"
 CCTK_REAL Avec_y_rhs TYPE=gf CENTERING={vcv} TAGS='checkpoint="no"' "y-component of vector potential RHS"
 CCTK_REAL Avec_z_rhs TYPE=gf CENTERING={vvc} TAGS='checkpoint="no"' "z-component of vector potential RHS"
 
+# ----- cell-centered magnetic field
+
 CCTK_REAL dB TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"'
 {
   dBx, dBy, dBz
 } "Densitized magnetic field components: sqrt(det(gamma_ij))*B^k"
 
+# ----- staggered magnetic field (faces)
+
 CCTK_REAL dBx_stag TYPE=gf CENTERING={vcc} TAGS='checkpoint="no"' "x-component of staggered densitized magnetic field"
 CCTK_REAL dBy_stag TYPE=gf CENTERING={cvc} TAGS='checkpoint="no"' "y-component of staggered densitized magnetic field"
 CCTK_REAL dBz_stag TYPE=gf CENTERING={ccv} TAGS='checkpoint="no"' "z-component of staggered densitized magnetic field"
 
-CCTK_REAL Psi TYPE=gf CENTERING={vvv} TAGS='rhs="Psi_rhs"' "t-component of vector potential: sqrt(det(gamma_ij))*Phi"
+# ----- staggered gauge functions
+
+CCTK_REAL Psi TYPE=gf CENTERING={vvv} TAGS='rhs="Psi_rhs" prolongation_type="hermite" prolongation_order=5' "t-component of vector potential: sqrt(det(gamma_ij))*Phi"
 
 CCTK_REAL Psi_rhs TYPE=gf CENTERING={vvv} TAGS='checkpoint="no"' "t-component of vector potential RHS"
 
@@ -110,21 +119,29 @@ CCTK_REAL Fbeta TYPE=gf CENTERING={vvv} TAGS='checkpoint="no"'
   Fbetax, Fbetay, Fbetaz
 } "Auxiliary variables which appear in the rhs of Psi"
 
-CCTK_REAL LOflag TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Shock detector switch"
-CCTK_REAL etac TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Shock detector measure"
-
-CCTK_REAL Ex TYPE=gf CENTERING={cvv} TAGS='checkpoint="no"' "x-component of E"
-CCTK_REAL Ey TYPE=gf CENTERING={vcv} TAGS='checkpoint="no"' "y-component of E"
-CCTK_REAL Ez TYPE=gf CENTERING={vvc} TAGS='checkpoint="no"' "z-component of E"
-
 CCTK_REAL G TYPE=gf CENTERING={vvv} TAGS='checkpoint="no"' "Auxiliary variable which appears in the rhs of Avec"
 
 CCTK_REAL Fx_stag TYPE=gf CENTERING={cvv} TAGS='checkpoint="no"' "x-component of F"
 CCTK_REAL Fy_stag TYPE=gf CENTERING={vcv} TAGS='checkpoint="no"' "y-component of F"
 CCTK_REAL Fz_stag TYPE=gf CENTERING={vvc} TAGS='checkpoint="no"' "z-component of F"
 
+# ----- staggered electric fields (edges)
+
+CCTK_REAL Ex TYPE=gf CENTERING={cvv} TAGS='checkpoint="no"' "x-component of E"
+CCTK_REAL Ey TYPE=gf CENTERING={vcv} TAGS='checkpoint="no"' "y-component of E"
+CCTK_REAL Ez TYPE=gf CENTERING={vvc} TAGS='checkpoint="no"' "z-component of E"
+
+# ----- auxiliary variables
+
 CCTK_REAL con2prim_flag TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "Flag indicating failure of con2prim"
 
+CCTK_REAL diagnostics TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"'
+{
+  w_lorentz
+  B_norm
+  b2small
+} "Diagnostic grid functions"
+
 # saved_prims are computed only once at the beginning of the simulation
 # at the outer physical boundary because they are not needed afterwards 
 CCTK_REAL saved_prims TYPE=gf CENTERING={ccc}
@@ -145,7 +162,7 @@ CCTK_REAL svec TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"'
   svec_x, svec_y, svec_z
 } ""
 
-#grid functions required for upwindCT computation of electric field
+# ----- grid functions required for upwindCT computation of electric field
 
 CCTK_REAL vbar_xface TYPE=gf CENTERING={vcc} TAGS='checkpoint="no"'
 {

AsterX/param.ccl

@@ -21,7 +21,7 @@ BOOLEAN use_entropy_fix "Shall we use inversion based on entropy as a backup?" S
 
 BOOLEAN prolongate_dbstag "Shall we prolongate dBstag?" STEERABLE=always
 {
-} "yes"
+} "no"
 
 KEYWORD flux_type "Flux solver" STEERABLE=always
 {
@@ -184,7 +184,8 @@ USES CCTK_REAL alp_thresh
 USES CCTK_REAL rho_BH
 USES CCTK_REAL eps_BH
 USES CCTK_REAL vwlim_BH
-USES CCTK_REAL cons_error_limit
+USES CCTK_REAL sigma_max
+USES CCTK_REAL inv_beta_max
 USES BOOLEAN use_z
 USES BOOLEAN use_temperature
 

AsterX/schedule.ccl

@@ -1,12 +1,13 @@
 # Schedule definitions for thorn AsterX
 
-STORAGE: dens mom tau DEnt DYe dB Psi HydroBaseX::Bvec dBx_stag dBy_stag dBz_stag
+STORAGE: cons_vector dB Psi HydroBaseX::Bvec dBx_stag dBy_stag dBz_stag
 STORAGE: flux_x flux_y flux_z
-STORAGE: densrhs momrhs taurhs DEntrhs Avec_x_rhs Avec_y_rhs Avec_z_rhs Psi_rhs
+STORAGE: rhs_vector Avec_x_rhs Avec_y_rhs Avec_z_rhs Psi_rhs
 STORAGE: ADMBaseX::metric ADMBaseX::lapse ADMBaseX::shift ADMBaseX::curv
 STORAGE: G Fbeta Fx_stag Fy_stag Fz_stag
 STORAGE: TmunuBaseX::eTtt TmunuBaseX::eTti TmunuBaseX::eTij
-
+STORAGE: Ex Ey Ez
+STORAGE: con2prim_flag diagnostics
 
 
 if(unit_test)
@@ -104,8 +105,8 @@ SCHEDULE AsterX_Prim2Con_Initial IN AsterX_InitialGroup AFTER AsterX_CheckPrims
   READS: ADMBaseX::metric(interior)
   READS: HydroBaseX::rho(interior) HydroBaseX::vel(interior) HydroBaseX::eps(interior) HydroBaseX::press(interior) HydroBaseX::Bvec(interior)
   READS: HydroBaseX::entropy(interior) HydroBaseX::Ye(interior)
-  WRITES: dens(interior) tau(interior) DEnt(interior) mom(interior) DYe(interior) dB(interior)
-  SYNC: dens tau DEnt mom dYe dB 
+  WRITES: cons_vector(interior) dB(interior)
+  SYNC: cons_vector dB 
 } "Compute conserved variables from primitive variables at initial"
 
 if (zero_psi) {
@@ -160,15 +161,15 @@ SCHEDULE AsterX_Sync AT postregrid
 {
   LANG: C
   OPTIONS: global
-  SYNC: dens tau mom DEnt DYe Avec_x Avec_y Avec_z Psi
+  SYNC: cons_vector Avec_x Avec_y Avec_z Psi
   SYNC: saved_prims
 } "Synchronize"
 
 SCHEDULE AsterX_Sync IN ODESolvers_PostStep
 {
   LANG: C
   OPTIONS: global
-  SYNC: dens tau mom DEnt DYe Avec_x Avec_y Avec_z Psi
+  SYNC: cons_vector Avec_x Avec_y Avec_z Psi
 } "Synchronize"
 
 SCHEDULE AsterX_Sync AT post_recover_variables BEFORE ODESolvers_PostStep AFTER AsterX_TouchNonEvolVarsGroup
@@ -226,17 +227,19 @@ SCHEDULE AsterX_Con2Prim IN AsterX_Con2PrimGroup AFTER AsterX_ComputedBFromdBsta
   READS: AsterMasks::aster_mask_vc(everywhere)
   READS: ADMBaseX::metric(everywhere)
   READS: ADMBaseX::lapse(everywhere)
-  READS: dens(everywhere) tau(everywhere) DEnt(everywhere) mom(everywhere) DYe(everywhere) dB(everywhere)
+  READS: cons_vector(everywhere) dB(everywhere)
+  READS: ADMBaseX::shift(everywhere)
   READS: saved_prims(everywhere)
   READS: Avec_x(everywhere) Avec_y(everywhere) Avec_z(everywhere)
   WRITES: con2prim_flag(everywhere)
+  WRITES: diagnostics(everywhere)
   WRITES: AsterMasks::aster_mask_cc(everywhere)
   WRITES: HydroBaseX::rho(everywhere) HydroBaseX::vel(everywhere) HydroBaseX::eps(everywhere) HydroBaseX::press(everywhere) HydroBaseX::Bvec(everywhere)
   WRITES: HydroBaseX::temperature(everywhere) HydroBaseX::entropy(everywhere) HydroBaseX::Ye(everywhere)
   WRITES: saved_prims(everywhere)
   WRITES: zvec(everywhere)
   WRITES: svec(everywhere)
-  WRITES: dens(everywhere) tau(everywhere) DEnt(everywhere) mom(everywhere) DYe(everywhere) dB(everywhere)
+  WRITES: cons_vector(everywhere) dB(everywhere)
 } "Calculate primitive variables from conservative variables"
 
 SCHEDULE AsterX_ApplyOuterBCOnPrim IN AsterX_Con2PrimGroup AFTER AsterX_Con2Prim
@@ -281,7 +284,7 @@ SCHEDULE AsterX_Fluxes IN AsterX_RHSGroup
   READS: zvec_x(everywhere) zvec_y(everywhere) zvec_z(everywhere)
   READS: svec_x(everywhere) svec_y(everywhere) svec_z(everywhere)
   READS: Avec_x(everywhere) Avec_y(everywhere) Avec_z(everywhere) Psi(everywhere)
-  READS: dens(everywhere) DYe(everywhere)
+  READS: cons_vector(everywhere)
   READS: LOflag(everywhere)
   WRITES: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)
   WRITES: theta_x(everywhere) theta_y(everywhere) theta_z(everywhere)
@@ -332,23 +335,23 @@ SCHEDULE AsterX_SourceTerms IN AsterX_RHSGroup AFTER AsterX_Fluxes
   READS: zvec_x(everywhere), zvec_y(everywhere), zvec_z(everywhere)
   READS: svec_x(everywhere), svec_y(everywhere), svec_z(everywhere)
   READS: HydroBaseX::Bvec(everywhere)
-  WRITES: densrhs(interior) taurhs(interior) momrhs(interior) DEntrhs(interior) DYe_rhs(interior)
+  WRITES: rhs_vector(interior)
 } "Calculate the source terms and compute the RHS of the hydro equations"
 
 SCHEDULE AsterX_RHS IN AsterX_RHSGroup AFTER AsterX_SourceTerms
 {
   LANG: C
   READS: ADMBaseX::metric(everywhere) ADMBaseX::lapse(everywhere) ADMBaseX::shift(everywhere)
   READS: flux_x(everywhere) flux_y(everywhere) flux_z(everywhere)
-  READS: densrhs(interior) taurhs(interior) momrhs(interior) DEntrhs(interior) DYe_rhs(interior)
+  READS: rhs_vector(interior)
   READS: Psi(everywhere)
   READS: Ex(interior) Ey(interior) Ez(interior)
   READS: G(everywhere)
   READS: Fbeta(everywhere)
   READS: Fx_stag(everywhere) Fy_stag(everywhere) Fz_stag(everywhere)
   READS: theta_x(everywhere) theta_y(everywhere) theta_z(everywhere)
   READS: LOflag(everywhere) 
-  WRITES: densrhs(interior) taurhs(interior) momrhs(interior) DEntrhs(interior) DYe_rhs(interior)
+  WRITES: rhs_vector(interior)
   WRITES: theta_tot(interior)
   WRITES: Avec_x_rhs(interior) Avec_y_rhs(interior) Avec_z_rhs(interior) Psi_rhs(interior)
 } "Update the RHS of the hydro equations with the flux contributions"
@@ -364,15 +367,15 @@ if(interpolate_failed_c2p)
     READS: con2prim_flag(interior)
     READS: HydroBaseX::rho(interior) HydroBaseX::vel(interior) HydroBaseX::eps(interior) HydroBaseX::press(interior) HydroBaseX::Bvec(interior)
     READS: saved_prims(interior)
-    READS: dens(interior) tau(interior) mom(interior)
+    READS: cons_vector(interior)
     WRITES: con2prim_flag(interior)
     WRITES: HydroBaseX::rho(interior) HydroBaseX::vel(interior) HydroBaseX::eps(interior) HydroBaseX::press(interior)
     WRITES: saved_prims(interior)
-    WRITES: dens(interior) tau(interior) mom(interior)
+    WRITES: cons_vector(interior)
     SYNC: con2prim_flag
     SYNC: HydroBaseX::rho HydroBaseX::vel HydroBaseX::eps HydroBaseX::press
     SYNC: saved_prims
-    SYNC: dens tau mom
+    SYNC: cons_vector
   } "Interpolate those grid points where con2prim failed"
 }