diff options
Diffstat (limited to 'src/GRHydro_Source.F90')
| -rw-r--r-- | src/GRHydro_Source.F90 | 18 |
1 files changed, 8 insertions, 10 deletions
diff --git a/src/GRHydro_Source.F90 b/src/GRHydro_Source.F90 index 553c916..1923c2d 100644 --- a/src/GRHydro_Source.F90 +++ b/src/GRHydro_Source.F90 @@ -71,7 +71,7 @@ subroutine SourceTerms(CCTK_ARGUMENTS) CCTK_INT :: i, j, k, na, nb, nc CCTK_REAL :: one, two, half CCTK_REAL :: t00, t0a, t0b, t0c, taa, tab, tac, tbb, tbc, tcc - CCTK_REAL :: sqrtdet, det, uaa, uab, uac, ubb, ubc, ucc, rhoenthalpyW2 + CCTK_REAL :: uaa, uab, uac, ubb, ubc, ucc, rhoenthalpyW2 CCTK_REAL :: shifta, shiftb, shiftc, velashift, velbshift, velcshift CCTK_REAL :: vlowa, vlowb, vlowc CCTK_REAL :: da_betaa, da_betab, da_betac, db_betaa, db_betab,& @@ -191,7 +191,7 @@ subroutine SourceTerms(CCTK_ARGUMENTS) !$OMP PARALLEL DO PRIVATE(i, j, k, local_spatial_order,& !$OMP localgaa,localgab,localgac,localgbb,localgbc,localgcc,& - !$OMP det,sqrtdet,rhoenthalpyW2,shifta,shiftb,shiftc,& + !$OMP rhoenthalpyW2,shifta,shiftb,shiftc,& !$OMP da_betaa,da_betab,da_betac,db_betaa,db_betab,db_betac,& !$OMP dc_betaa,dc_betab,dc_betac,velashift,velbshift,velcshift,& !$OMP vlowa,vlowb,vlowc,t00,t0a,t0b,t0c,taa,tab,tac,tbb,tbc,tcc,& @@ -221,10 +221,8 @@ subroutine SourceTerms(CCTK_ARGUMENTS) localgbc = g23(i,j,k) localgcc = g33(i,j,k) - det = SPATIAL_DETERMINANT(localgaa, localgab, localgac,\ - localgbb, localgbc, localgcc) - sqrtdet = sqrt(det) - call UpperMetric(uaa, uab, uac, ubb, ubc, ucc, det, localgaa,& + call UpperMetric(uaa, uab, uac, ubb, ubc, ucc, & + sdetg(i,j,k)*sdetg(i,j,k), localgaa,& localgab, localgac, localgbb, localgbc, localgcc) !!$ All the matter variables (except velocity) always appear @@ -442,10 +440,10 @@ subroutine SourceTerms(CCTK_ARGUMENTS) invalp densrhs(i,j,k) = 0.d0 - srhs(i,j,k,1) = alp(i,j,k)*sqrtdet*sa_source - srhs(i,j,k,2) = alp(i,j,k)*sqrtdet*sb_source - srhs(i,j,k,3) = alp(i,j,k)*sqrtdet*sc_source - taurhs(i,j,k) = alp(i,j,k)*sqrtdet*tau_source + srhs(i,j,k,1) = alp(i,j,k)*sdetg(i,j,k)*sa_source + srhs(i,j,k,2) = alp(i,j,k)*sdetg(i,j,k)*sb_source + srhs(i,j,k,3) = alp(i,j,k)*sdetg(i,j,k)*sc_source + taurhs(i,j,k) = alp(i,j,k)*sdetg(i,j,k)*tau_source enddo enddo |