aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_Prim2Con.F90
diff options
context:
space:
mode:
authorcott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2011-02-07 02:42:02 +0000
committercott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2011-02-07 02:42:02 +0000
commit54442343375944433c04f86a32370327a4c061d8 (patch)
tree32e0ce2a69e34985ff75536fb071a6952664c4b5 /src/GRHydro_Prim2Con.F90
parent8fbbe803f1a5eb3cb41a76a857eb1d9cb03f6314 (diff)
* add !OMP CRITICAL around CCTK_WARN statements
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@215 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Prim2Con.F90')
-rw-r--r--src/GRHydro_Prim2Con.F908
1 files changed, 7 insertions, 1 deletions
diff --git a/src/GRHydro_Prim2Con.F90 b/src/GRHydro_Prim2Con.F90
index 957a862..0533e27 100644
--- a/src/GRHydro_Prim2Con.F90
+++ b/src/GRHydro_Prim2Con.F90
@@ -46,7 +46,7 @@ subroutine primitive2conservative(CCTK_ARGUMENTS)
integer :: i, j, k
CCTK_REAL :: gxxl,gxyl,gxzl,gyyl,gyzl,gzzl,avg_detl,&
gxxr,gxyr,gxzr,gyyr,gyzr,gzzr,avg_detr
- CCTK_REAL :: xtemp(1)
+ CCTK_REAL :: xtemp
if(evolve_temper.ne.1) then
!$OMP PARALLEL DO PRIVATE(i, j, avg_detl, avg_detr,&
@@ -253,6 +253,7 @@ subroutine prim2con_hot(handle, GRHydro_reflevel, ii, jj, kk, &
drho,deps,temp,ye,dpress,keyerr,anyerr)
keytemp=0
if(anyerr.ne.0) then
+ !OMP CRITICAL
call CCTK_WARN(1,"EOS error in prim2con_hot: lev 2")
write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z
call CCTK_WARN(1,warnline)
@@ -262,11 +263,13 @@ subroutine prim2con_hot(handle, GRHydro_reflevel, ii, jj, kk, &
call CCTK_WARN(1,warnline)
write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel
call CCTK_WARN(1,warnline)
+ !OMP END CRITICAL
endif
else
! This is a way of recovering even on finer refinement levels:
! Use the average temperature at the interface instead of the
! reconstructed specific internal energy.
+ !OMP CRITICAL
call CCTK_WARN(1,"EOS error in prim2con_hot: NOW using averaged temp!")
write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z
call CCTK_WARN(1,warnline)
@@ -276,12 +279,14 @@ subroutine prim2con_hot(handle, GRHydro_reflevel, ii, jj, kk, &
call CCTK_WARN(1,warnline)
write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel
call CCTK_WARN(1,warnline)
+ !OMP END CRITICAL
keytemp=1
temp = temp0
call EOS_Omni_press(handle,keytemp,GRHydro_eos_rf_prec,n,&
drho,deps,temp,ye,dpress,keyerr,anyerr)
keytemp=0
if(anyerr.ne.0) then
+ !OMP CRITICAL
call CCTK_WARN(1,"EOS error in prim2con_hot")
write(warnline,"(3i5,1P10E15.6)") ii,jj,kk,x,y,z
call CCTK_WARN(1,warnline)
@@ -292,6 +297,7 @@ subroutine prim2con_hot(handle, GRHydro_reflevel, ii, jj, kk, &
write(warnline,"(A10,i5)") "reflevel: ", GRHydro_reflevel
call CCTK_WARN(1,warnline)
call CCTK_WARN(0,"Aborting!!!")
+ !OMP END CRITICAL
endif
endif
endif