diff options
-rw-r--r-- | interface.ccl | 8 | ||||
-rw-r--r-- | schedule.ccl | 6 |
2 files changed, 10 insertions, 4 deletions
diff --git a/interface.ccl b/interface.ccl index ef1cc27..2b58ca5 100644 --- a/interface.ccl +++ b/interface.ccl @@ -345,7 +345,7 @@ real tau type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolo real scon[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="D" tensorweight=+1.0 jacobian="inverse_jacobian" interpolator="matter"' "generalized momenta" -real Bcons[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="D" tensorweight=+1.0 jacobian="inverse_jacobian" interpolator="matter"' "B-field conservative variable" +real Bcons[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="D" tensorparity=-1 tensorweight=+1.0 jacobian="inverse_jacobian" interpolator="matter"' "B-field conservative variable" real Y_e_con type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" tensorweight=+1.0 jacobian="inverse_jacobian" interpolator="matter"' "Conserved electron fraction" @@ -367,9 +367,9 @@ real psidcrhs type = GF Timelevels = 1 tags='Prolongation="None" checkpoint="no" real divB type = GF Timelevels = 1 tags='Prolongation="None" checkpoint="no"' "Magnetic field constraint" -real bcom[3] type = GF Timelevels = 1 tags='Prolongation="none" tensortypealias="U"' "b^i: comoving contravariant magnetic field 4-vector spatial components" -real bcom0 type = GF Timelevels = 1 tags='Prolongation="none" tensortypealias="Scalar"' "b^0 component of the comoving contravariant magnetic field 4-vector" -real bcom_sq type = GF Timelevels = 1 tags='Prolongation="none" tensortypealias="Scalar"' "half of magnectic pressure: contraction of b_a b^a " +real bcom[3] type = GF Timelevels = 1 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" tensorparity=-1 interpolator="matter"' "b^i: comoving contravariant magnetic field 4-vector spatial components" +real bcom0 type = GF Timelevels = 1 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "b^0 component of the comoving contravariant magnetic field 4-vector" +real bcom_sq type = GF Timelevels = 1 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "half of magnectic pressure: contraction of b_a b^a " ################################################## ### variables in the local tensor basis diff --git a/schedule.ccl b/schedule.ccl index 8005d7e..05d031b 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -82,6 +82,12 @@ if(CCTK_Equals(Bvec_evolution_method,"GRHydro")) { STORAGE:divB } + if (calculate_bcom) + { + STORAGE:GRHydro::bcom + STORAGE:GRHydro::bcom0 + STORAGE:GRHydro::bcom_sq + } } STORAGE:GRHydro_eos_scalars STORAGE:GRHydro_minima |