User: rhaas Date: 2013/03/27 08:47 PM
Modified: /trunk/src/ GRHydro_PPM.F90, GRHydro_PPMM.F90
Log: GRHydro: explicitly cast from CCTK_REAL to int to avoid warning
From: Roland Haas rhaas@tapir.caltech.edu
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_PPM.F90 Delta lines: +7 -7 =================================================================== --- trunk/src/GRHydro_PPM.F90 2013-03-28 01:46:46 UTC (rev 497) +++ trunk/src/GRHydro_PPM.F90 2013-03-28 01:47:00 UTC (rev 498) @@ -494,7 +494,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) if (abs(1.d0 - flatten) > 0.d0) then trivial_rp(i-1) = .false. @@ -829,7 +829,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) if (abs(1.d0 - flatten) > 0.d0) then trivial_rp(i-1) = .false. @@ -1196,7 +1196,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) tempplus(i) = flatten * tempplus(i) + (1.d0 - flatten) * temperature(i) tempminus(i) = flatten * tempminus(i) + (1.d0 - flatten) * temperature(i) @@ -1459,7 +1459,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) tempplus(i) = flatten * tempplus(i) + (1.d0 - flatten) * temperature(i) tempminus(i) = flatten * tempminus(i) + (1.d0 - flatten) * temperature(i) @@ -1679,7 +1679,7 @@ else !!$ Really implement C&W, page 197; which requires stencil 4. do itracer=1,number_of_tracers do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) tracerplus(i,itracer) = flatten * tracerplus(i,itracer) + & (1.d0 - flatten) * tracer(i,itracer) @@ -1910,7 +1910,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) Y_e_plus(i) = flatten * Y_e_plus(i) + & (1.d0 - flatten) * Y_e(i) @@ -2040,7 +2040,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(one, -dpress(i)) + s=int(sign(one, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) Y_e_plus(i) = flatten * Y_e_plus(i) + & (1.d0 - flatten) * Y_e(i)
File [modified]: GRHydro_PPMM.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PPMM.F90 2013-03-28 01:46:46 UTC (rev 497) +++ trunk/src/GRHydro_PPMM.F90 2013-03-28 01:47:00 UTC (rev 498) @@ -411,7 +411,7 @@ end do else !!$ Really implement C&W, page 197; which requires stencil 4. do i = 4, nx - 3 - s=sign(1.d0, -dpress(i)) + s=int(sign(1.d0, -dpress(i))) flatten = max(tilde_flatten(i), tilde_flatten(i+s)) if (abs(1.d0 - flatten) > 0.d0) then trivial_rp(i-1) = .false.
commits@lists.einsteintoolkit.org