diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-29 17:17:01 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-29 17:17:01 +0000 |
commit | 63edce19a45f1ecf92680655d30505ed199df466 (patch) | |
tree | e19eafbf53a5e61af11d9f851902ffa309283c76 /src | |
parent | eb8aadfb83add7cffda819a262022aee5c5e4595 (diff) |
fix error checking in GRHydro_Boundaries
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@305 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_Boundaries.F90 | 76 |
1 files changed, 38 insertions, 38 deletions
diff --git a/src/GRHydro_Boundaries.F90 b/src/GRHydro_Boundaries.F90 index a432d09..5fe3f85 100644 --- a/src/GRHydro_Boundaries.F90 +++ b/src/GRHydro_Boundaries.F90 @@ -170,7 +170,7 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS) DECLARE_CCTK_FUNCTIONS integer, dimension(3) :: sw - integer :: ierr + integer :: ierr = 0 CCTK_REAL :: pi = 3.141569d0 integer :: i,j,k CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates @@ -185,61 +185,61 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS) !!$Flat boundaries if required if (CCTK_EQUALS(bound,"flat")) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::dens", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::tau", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::w_lorentz", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::rho", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::press", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::eps", "Flat") if (general_coordinates .ne. 0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::lvel", "Flat") else - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + 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 = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::lBvec", "Flat") else - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::Bvec", "Flat") endif - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::Bcons", "Flat") if(clean_divergence.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::psidc", "Flat") endif endif if(evolve_tracer.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::GRHydro_tracers", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::GRHydro_cons_tracers", "Flat") endif if(evolve_y_e.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::Y_e", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::Y_e_con", "Flat") endif if(evolve_temper.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::temperature", "Flat") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::entropy", "Flat") endif @@ -247,56 +247,56 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS) endif if (CCTK_EQUALS(bound,"none")) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::dens", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::tau", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::scon", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::w_lorentz", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::rho", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::press", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::eps", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::vel", "None") if(evolve_mhd.ne.0) then if (general_coordinates .ne. 0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::lBvec", "None") else - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::Bvec", "None") endif - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::Bcons", "None") if(clean_divergence.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::psidc", "None") endif endif if(evolve_tracer.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::GRHydro_tracers", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::GRHydro_cons_tracers", "None") endif if(evolve_y_e.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::Y_e", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "GRHydro::Y_e_con", "None") endif if(evolve_temper.ne.0) then - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::temperature", "None") - ierr = Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & + ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, & "HydroBase::entropy", "None") endif |