[Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 458)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Mon Jan 14 08:23:40 CST 2013
User: rhaas
Date: 2013/01/14 08:23 AM
Modified:
/trunk/
schedule.ccl
/trunk/src/
GRHydro_Boundaries.F90, GRHydro_Con2PrimAM.F90
Log:
GRHydro: Starting to push vector potential changes
* incl boundary calls for Avec
* disable call to BvecfromAvec from within loop, this one is a schedule
function call now.
From: Philipp Moesta <pmoesta at tapir.caltech.edu>
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_Boundaries.F90
Delta lines: +120 -134
===================================================================
--- trunk/src/GRHydro_Boundaries.F90 2013-01-14 14:23:38 UTC (rev 457)
+++ trunk/src/GRHydro_Boundaries.F90 2013-01-14 14:23:40 UTC (rev 458)
@@ -85,11 +85,6 @@
call SetCartSymGN(ierr, cctkGH, sym, "HydroBase::entropy")
!endif
endif
-
- if(evolve_entropy.ne.0) then
- call SetCartSymGN(ierr, cctkGH, sym, "HydroBase::entropy")
- call SetCartSymGN(ierr, cctkGH, sym, "GRHydro::entropycons")
- endif
call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::atmosphere_mask_real")
@@ -99,9 +94,9 @@
!if (sync_conserved_only .eq. 0) then
if (general_coordinates .ne. 0) then
- call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[0]")
+ call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[0]")
else
- call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]")
+ call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]")
endif
!endif
call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[0]")
@@ -118,15 +113,15 @@
!if (sync_conserved_only .eq. 0) then
if (general_coordinates .ne. 0) then
- call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[1]")
+ call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[1]")
else
- call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]")
+ call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]")
endif
!endif
call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[1]")
if(evolve_mhd.ne.0) then
!if (sync_conserved_only .eq. 0) then
- call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]")
+ call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]")
!endif
call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::Bcons[1]")
endif
@@ -137,9 +132,9 @@
!if (sync_conserved_only .eq. 0) then
if (general_coordinates .ne. 0) then
- call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[2]")
+ call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[2]")
else
- call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]")
+ call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]")
endif
!endif
call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[2]")
@@ -206,31 +201,31 @@
"GRHydro::scon", "Flat")
if (sync_conserved_only .eq. 0) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::w_lorentz", "Flat")
+ "HydroBase::w_lorentz", "Flat")
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::rho", "Flat")
+ "HydroBase::rho", "Flat")
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::press", "Flat")
+ "HydroBase::press", "Flat")
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::eps", "Flat")
+ "HydroBase::eps", "Flat")
if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lvel", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lvel", "Flat")
else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::vel", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::vel", "Flat")
endif
endif
if(evolve_mhd.ne.0) then
- if (sync_conserved_only .eq. 0) then
- if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lBvec", "Flat")
- else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Bvec", "Flat")
- endif
- endif
+ if (sync_conserved_only .eq. 0) then
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "Flat")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "Flat")
+ endif
+ endif
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"GRHydro::Bcons", "Flat")
if(clean_divergence.ne.0) then
@@ -239,10 +234,22 @@
endif
endif
+ if (CCTK_Equals(Bvec_evolution_method, "GRHydro_Avec")) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Avec", "Flat")
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "Flat")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "Flat")
+ endif
+ endif
+
if(evolve_tracer.ne.0) then
if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::GRHydro_tracers", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::GRHydro_tracers", "Flat")
endif
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"GRHydro::GRHydro_cons_tracers", "Flat")
@@ -250,8 +257,8 @@
if(evolve_y_e.ne.0) then
if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Y_e", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Y_e", "Flat")
endif
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"GRHydro::Y_e_con", "Flat")
@@ -259,21 +266,13 @@
if(evolve_temper.ne.0) then
if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::temperature", "Flat")
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::temperature", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::entropy", "Flat")
endif
endif
- if(evolve_entropy.ne.0) then
- if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "Flat")
- endif
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::entropycons", "Flat")
- endif
endif
@@ -287,44 +286,56 @@
if (sync_conserved_only .eq. 0) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::w_lorentz", "None")
+ "HydroBase::w_lorentz", "None")
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::rho", "None")
+ "HydroBase::rho", "None")
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::press", "None")
+ "HydroBase::press", "None")
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::eps", "None")
+ "HydroBase::eps", "None")
if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lvel", "None")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lvel", "None")
else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::vel", "None")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::vel", "None")
endif
endif
if(evolve_mhd.ne.0) then
- if (sync_conserved_only .eq. 0) then
- if (general_coordinates .ne. 0) then
+ if (sync_conserved_only .eq. 0) then
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "None")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "None")
+ endif
+ endif
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::Bcons", "None")
+ if(clean_divergence.ne.0) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lBvec", "None")
- else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Bvec", "None")
- endif
- endif
+ "GRHydro::psidc", "None")
+ endif
+ endif
+
+ if (CCTK_Equals(Bvec_evolution_method, "GRHydro_Avec")) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::Bcons", "None")
- if(clean_divergence.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::psidc", "None")
+ "HydroBase::Avec", "None")
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "Flat")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "Flat")
endif
- endif
+ endif
if(evolve_tracer.ne.0) then
if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::GRHydro_tracers", "None")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::GRHydro_tracers", "None")
endif
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"GRHydro::GRHydro_cons_tracers", "None")
@@ -332,8 +343,8 @@
if(evolve_y_e.ne.0) then
if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Y_e", "None")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Y_e", "None")
endif
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"GRHydro::Y_e_con", "None")
@@ -341,22 +352,13 @@
if(evolve_temper.ne.0) then
if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::temperature", "None")
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "None")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::temperature", "None")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::entropy", "None")
endif
endif
- if(evolve_entropy.ne.0) then
- if (sync_conserved_only .eq. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "None")
- endif
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::entropycons", "None")
- endif
-
endif
if (CCTK_EQUALS(bound,"scalar")) then
@@ -415,13 +417,13 @@
"HydroBase::vel", "Flat")
endif
if(evolve_mhd.ne.0) then
- if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lBvec", "Flat")
- else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Bvec", "Flat")
- endif
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "Flat")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "Flat")
+ endif
endif
!if(evolve_tracer.ne.0) then
@@ -441,10 +443,6 @@
! "HydroBase::entropy", "Flat")
endif
- if(evolve_entropy.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "Flat")
- endif
endif
@@ -466,13 +464,14 @@
endif
if(evolve_mhd.ne.0) then
- if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lBvec", "None")
- else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Bvec", "None")
- endif
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "None")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "None")
+ endif
+
endif
!if(evolve_tracer.ne.0) then
@@ -494,11 +493,6 @@
! "HydroBase::entropy", "None")
endif
- if(evolve_entropy.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "None")
- endif
-
end if
if (CCTK_EQUALS(bound,"scalar")) then
@@ -558,25 +552,24 @@
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"HydroBase::vel", "Flat")
endif
-
if(evolve_mhd.ne.0) then
- if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lBvec", "Flat")
- else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Bvec", "Flat")
- endif
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "Flat")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "Flat")
+ endif
endif
if(evolve_tracer.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::GRHydro_tracers", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::GRHydro_tracers", "Flat")
endif
if(evolve_y_e.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Y_e", "Flat")
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Y_e", "Flat")
endif
if(evolve_temper.ne.0) then
@@ -586,10 +579,6 @@
"HydroBase::entropy", "Flat")
endif
- if(evolve_entropy.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "Flat")
- endif
endif
@@ -611,18 +600,20 @@
endif
if(evolve_mhd.ne.0) then
- if (general_coordinates .ne. 0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "GRHydro::lBvec", "None")
- else
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::Bvec", "None")
- endif
+ if (general_coordinates .ne. 0) then
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "GRHydro::lBvec", "None")
+ else
+ ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+ "HydroBase::Bvec", "None")
+ endif
+
endif
if(evolve_tracer.ne.0) then
ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
"GRHydro::GRHydro_tracers", "None")
+
endif
if(evolve_y_e.ne.0) then
@@ -638,11 +629,6 @@
"HydroBase::entropy", "None")
endif
- if(evolve_entropy.ne.0) then
- ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
- "HydroBase::entropy", "None")
- endif
-
end if
if (CCTK_EQUALS(bound,"scalar")) then
File [modified]: GRHydro_Con2PrimAM.F90
Delta lines: +0 -3
===================================================================
--- trunk/src/GRHydro_Con2PrimAM.F90 2013-01-14 14:23:38 UTC (rev 457)
+++ trunk/src/GRHydro_Con2PrimAM.F90 2013-01-14 14:23:40 UTC (rev 458)
@@ -102,9 +102,6 @@
!as conformally flat only for con2prim inversion purposes.
posdef = .true.
- !! Get Bvec from Avec before rest, so we can reuse pointwise MHD recovery
- call GRHydro_BvecfromAvec(CCTK_PASS_FTOF)
-
if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then
g11 => gaa
g12 => gab
Directory: /trunk/
==================
File [modified]: schedule.ccl
Delta lines: +11 -6
===================================================================
--- trunk/schedule.ccl 2013-01-14 14:23:38 UTC (rev 457)
+++ trunk/schedule.ccl 2013-01-14 14:23:40 UTC (rev 458)
@@ -1066,21 +1066,26 @@
{
LANG: Fortran
} "Convert initial data given in primive variables to conserved variables - MHD version"
- }
- else if(CCTK_Equals(Bvec_evolution_method,"GRHydro_Avec"))
+ } else if(CCTK_Equals(Bvec_evolution_method,"GRHydro_Avec"))
{
- schedule Conservative2PrimitiveAM IN HydroBase_Con2Prim AS Con2Prim IF GRHydro::execute_MoL_PostStep
+ schedule GRHydro_BvecfromAvec IN HydroBase_Con2Prim
{
LANG: Fortran
+ } "Populate Bvec from Avec"
+# SCHEDULE group HydroBase_Boundaries AT CCTK_POSTINITIAL AFTER GRHydro_PoloidalMagFieldM AFTER GRHydro_BvecfromAvec
+ SCHEDULE group HydroBase_Boundaries IN HydroBase_Con2Prim AFTER GRHydro_BvecfromAvec
+ {
+ } "Call boundary conditions after magnetic field initial data setup"
+ schedule Conservative2PrimitiveAM IN HydroBase_Con2Prim AS Con2Prim AFTER HydroBase_Boundaries IF GRHydro::execute_MoL_PostStep
+ {
+ LANG: Fortran
} "Convert back to primitive variables (general) - MHD with Avec version"
schedule Primitive2ConservativeCellsAM IN HydroBase_Prim2ConInitial
{
LANG: Fortran
} "Convert initial data given in primive variables to conserved variables - MHD with Avec version"
- }
- else
- {
+ } else {
schedule Conservative2Primitive IN HydroBase_Con2Prim AS Con2Prim IF GRHydro::execute_MoL_PostStep
{
LANG: Fortran
More information about the Commits
mailing list