Changeset 1322 for palm/trunk/SOURCE/tridia_solver.f90

Timestamp:

Mar 20, 2014 4:38:49 PM (10 years ago)

Author:

raasch

Message:

REAL functions and a lot of REAL constants provided with KIND-attribute,
some small bugfixes

File:

• palm/trunk/SOURCE/tridia_solver.f90 (modified) (3 diffs)

palm/trunk/SOURCE/tridia_solver.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! REAL functions provided with KIND-attribute
 ! 
 ! Former revisions:
@@ -170 +170 @@
                    IF ( i >= 0  .AND.  i <= nnxh )  THEN
                       ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
-                                                  REAL( nx+1 ) ) ) / ( dx * dx ) + &
+                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
                                 2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                                  REAL( ny+1 ) ) ) / ( dy * dy )
+                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
                    ELSE
                       ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
-                                                  REAL( nx+1 ) ) ) / ( dx * dx ) + &
+                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
                                 2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                                  REAL( ny+1 ) ) ) / ( dy * dy )
+                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
                    ENDIF
                 ELSE
                    IF ( i >= 0  .AND.  i <= nnxh )  THEN
                       ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
-                                                  REAL( nx+1 ) ) ) / ( dx * dx ) + &
+                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
                                 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( ny+1-j ) ) / &
-                                                  REAL( ny+1 ) ) ) / ( dy * dy )
+                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
                    ELSE
                       ll(i,j) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
-                                                  REAL( nx+1 ) ) ) / ( dx * dx ) + &
+                                                  REAL( nx+1, KIND=wp ) ) ) / ( dx * dx ) + &
                                 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( ny+1-j ) ) / &
-                                                  REAL( ny+1 ) ) ) / ( dy * dy )
+                                                  REAL( ny+1, KIND=wp ) ) ) / ( dy * dy )
                    ENDIF
                 ENDIF
@@ -558 +558 @@
              IF ( i >= 0 .AND. i <= nnxh ) THEN
                 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * i ) / &
-                                         REAL( nx+1 ) ) ) * ddx2 + &
+                                         REAL( nx+1, KIND=wp ) ) ) * ddx2 + &
                        2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                         REAL( ny+1 ) ) ) * ddy2
+                                         REAL( ny+1, KIND=wp ) ) ) * ddy2
              ELSE
                 l(i) = 2.0 * ( 1.0 - COS( ( 2.0 * pi * ( nx+1-i ) ) / &
-                                         REAL( nx+1 ) ) ) * ddx2 + &
+                                         REAL( nx+1, KIND=wp ) ) ) * ddx2 + &
                        2.0 * ( 1.0 - COS( ( 2.0 * pi * j ) / &
-                                         REAL( ny+1 ) ) ) * ddy2
+                                         REAL( ny+1, KIND=wp ) ) ) * ddy2
              ENDIF
           ENDDO