diff options
Diffstat (limited to 'schedule.ccl')
-rw-r--r-- | schedule.ccl | 42 |
1 files changed, 34 insertions, 8 deletions
diff --git a/schedule.ccl b/schedule.ccl index 2e32506..2b72ab4 100644 --- a/schedule.ccl +++ b/schedule.ccl @@ -43,20 +43,46 @@ if (CCTK_Equals(initial_hydro,"shocktube")) { if(CCTK_Equals(shocktube_type,"diagshock")) { - schedule GRHydro_Diagshock_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection - { - LANG: Fortran - } "Diagonal shock boundary conditions" + if(clean_divergence){ + schedule GRHydro_Diagshock_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection + { + LANG: Fortran + SYNC: GRHydro::dens, GRHydro::tau, GRHydro::scon, GRHydro::Bcons, GRHydro::psidc + } "Diagonal shock boundary conditions" + } else { + schedule GRHydro_Diagshock_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection + { + LANG: Fortran + SYNC: GRHydro::dens, GRHydro::tau, GRHydro::scon, GRHydro::Bcons + } "Diagonal shock boundary conditions" + } } if(CCTK_Equals(shocktube_type,"diagshock2d")) { - schedule GRHydro_Diagshock2D_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection - { - LANG: Fortran - } "2-D Diagonal shock boundary conditions" + schedule GRHydro_Diagshock2D_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection + { + LANG: Fortran + } "2-D Diagonal shock boundary conditions" } +# if(CCTK_Equals(shocktube_type,"diagshock2d")) +# { +# if(clean_divergence){ +# schedule GRHydro_Diagshock2D_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection +# { +# LANG: Fortran +# SYNC: GRHydro::dens, GRHydro::tau, GRHydro::scon, GRHydro::Bcons, GRHydro::psidc +# } "2-D Diagonal shock boundary conditions" +# } else { +# schedule GRHydro_Diagshock2D_BoundaryM IN ApplyBCs AFTER BoundaryConditions AFTER Boundary::Boundary_ClearSelection +# { +# LANG: Fortran +# SYNC: GRHydro::dens, GRHydro::tau, GRHydro::scon, GRHydro::Bcons +# } "2-D Diagonal shock boundary conditions" +# } +# } + } else { schedule GRHydro_shocktube in HydroBase_Initial { |