Changeset 1013 for palm/trunk/SOURCE/poisfft.f90

Timestamp:

Sep 21, 2012 7:03:55 AM (12 years ago)

Author:

raasch

Message:

FLOAT type conversion replaced by REAL

File:

• palm/trunk/SOURCE/poisfft.f90 (modified) (4 diffs)

palm/trunk/SOURCE/poisfft.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! FLOAT type conversion replaced by REAL
 !
 ! Former revisions:
@@ -406 +406 @@
              IF ( i >= 0 .AND. i <= nnxh )  THEN
                 ll(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
-                                          FLOAT( nx+1 ) ) ) / ( dx * dx ) + &
+                                          REAL( nx+1 ) ) ) / ( dx * dx ) + &
                         2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                          FLOAT( ny+1 ) ) ) / ( dy * dy )
+                                          REAL( ny+1 ) ) ) / ( dy * dy )
              ELSE
                 ll(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
-                                          FLOAT( nx+1 ) ) ) / ( dx * dx ) + &
+                                          REAL( nx+1 ) ) ) / ( dx * dx ) + &
                         2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                          FLOAT( ny+1 ) ) ) / ( dy * dy )
+                                          REAL( ny+1 ) ) ) / ( dy * dy )
              ENDIF
              DO  k = 0,nz-1
@@ -423 +423 @@
 #else
           DO  i = 0, nnxh
-             ll(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / FLOAT( nx+1 ) ) ) / &
+             ll(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / REAL( nx+1 ) ) ) / &
                            ( dx * dx ) + &
-                     2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / FLOAT( ny+1 ) ) ) / &
+                     2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / REAL( ny+1 ) ) ) / &
                            ( dy * dy )
              DO  k = 0, nz-1
@@ -1489 +1489 @@
              IF ( i >= 0 .AND. i <= nnxh ) THEN
                 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
-                                         FLOAT( nx+1 ) ) ) * ddx2 + &
+                                         REAL( nx+1 ) ) ) * ddx2 + &
                        2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                         FLOAT( ny+1 ) ) ) * ddy2
+                                         REAL( ny+1 ) ) ) * ddy2
              ELSE
                 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
-                                         FLOAT( nx+1 ) ) ) * ddx2 + &
+                                         REAL( nx+1 ) ) ) * ddx2 + &
                        2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                         FLOAT( ny+1 ) ) ) * ddy2
+                                         REAL( ny+1 ) ) ) * ddy2
              ENDIF
           ENDDO