[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 498)

rhaas at tapir.caltech.edu rhaas at tapir.caltech.edu
Wed Mar 27 20:47:06 CDT 2013


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 at 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.



More information about the Commits mailing list