Changeset 1322 for palm/trunk/SOURCE/fft_xy.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/fft_xy.f90 (modified) (18 diffs)

palm/trunk/SOURCE/fft_xy.f90

@@ -16 +16 @@
 !
 ! Copyright 1997-2014 Leibniz Universitaet Hannover
-!--------------------------------------------------------------------------------!
+!------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! -----------------
-! 
+! REAL functions provided with KIND-attribute
 !
 ! Former revisions:
@@ -413 +413 @@
 
                    DO  i = 0, (nx+1)/2
-                      ar(i,j,k) = REAL( cwork(i) )
+                      ar(i,j,k) = REAL( cwork(i), KIND=wp )
                    ENDDO
                    DO  i = 1, (nx+1)/2 - 1
@@ -441 +441 @@
 
                    DO  i = 0, nx
-                      ar(i,j,k) = REAL( cwork(i) )
+                      ar(i,j,k) = REAL( cwork(i), KIND=wp )
                    ENDDO
 
@@ -519 +519 @@
 
                       DO  i = 0, (nx+1)/2
-                         ar_2d(i,j) = REAL( x_out(i) ) / ( nx+1 )
+                         ar_2d(i,j) = REAL( x_out(i), KIND=wp ) / ( nx+1 )
                       ENDDO
                       DO  i = 1, (nx+1)/2 - 1
@@ -528 +528 @@
 
                       DO  i = 0, (nx+1)/2
-                         ar(i,j,k) = REAL( x_out(i) ) / ( nx+1 )
+                         ar(i,j,k) = REAL( x_out(i), KIND=wp ) / ( nx+1 )
                       ENDDO
                       DO  i = 1, (nx+1)/2 - 1
@@ -689 +689 @@
 
                    DO  i = 0, (nx+1)/2
-                      ar(i,j,k)      = REAL( ar_tmp(i,j,k) )  * dnx
+                      ar(i,j,k)      = REAL( ar_tmp(i,j,k), KIND=wp )  * dnx
                    ENDDO
 
@@ -794 +794 @@
              CALL FFTN( cwork, ishape )
              DO  i = 0, (nx+1)/2
-                ar(i) = REAL( cwork(i) )
+                ar(i) = REAL( cwork(i), KIND=wp )
              ENDDO
              DO  i = 1, (nx+1)/2 - 1
@@ -813 +813 @@
 
              DO  i = 0, nx
-                ar(i) = REAL( cwork(i) )
+                ar(i) = REAL( cwork(i), KIND=wp )
              ENDDO
 
@@ -862 +862 @@
 
              DO  i = 0, (nx+1)/2
-                ar(i) = REAL( x_out(i) ) / ( nx+1 )
+                ar(i) = REAL( x_out(i), KIND=wp ) / ( nx+1 )
              ENDDO
              DO  i = 1, (nx+1)/2 - 1
@@ -1056 +1056 @@
 
                    DO  j = 0, (ny+1)/2
-                      ar_tr(j,i,k) = REAL( cwork(j) )
+                      ar_tr(j,i,k) = REAL( cwork(j), KIND=wp )
                    ENDDO
                    DO  j = 1, (ny+1)/2 - 1
@@ -1084 +1084 @@
 
                    DO  j = 0, ny
-                      ar(j,i,k) = REAL( cwork(j) )
+                      ar(j,i,k) = REAL( cwork(j), KIND=wp )
                    ENDDO
 
@@ -1160 +1160 @@
 
                    DO  j = 0, (ny+1)/2
-                      ar_tr(j,i,k) = REAL( y_out(j) ) / (ny+1)
+                      ar_tr(j,i,k) = REAL( y_out(j), KIND=wp ) / (ny+1)
                    ENDDO
                    DO  j = 1, (ny+1)/2 - 1
@@ -1305 +1305 @@
 
                    DO  j = 0, (ny+1)/2
-                      ar(j,i,k)      = REAL( ar_tmp(j,i,k) )  * dny
+                      ar(j,i,k)      = REAL( ar_tmp(j,i,k), KIND=wp )  * dny
                    ENDDO
 
@@ -1412 +1412 @@
 
              DO  j = 0, (ny+1)/2
-                ar(j) = REAL( cwork(j) )
+                ar(j) = REAL( cwork(j), KIND=wp )
              ENDDO
              DO  j = 1, (ny+1)/2 - 1
@@ -1431 +1431 @@
 
              DO  j = 0, ny
-                ar(j) = REAL( cwork(j) )
+                ar(j) = REAL( cwork(j), KIND=wp )
              ENDDO
 
@@ -1480 +1480 @@
 
              DO  j = 0, (ny+1)/2
-                ar(j) = REAL( y_out(j) ) / (ny+1)
+                ar(j) = REAL( y_out(j), KIND=wp ) / (ny+1)
              ENDDO
              DO  j = 1, (ny+1)/2 - 1
@@ -1674 +1674 @@
              DO  k = 1, nz
                 DO  i = 0, (nx+1)/2
-                   ar(i,k) = REAL( work(i+1,k) )
+                   ar(i,k) = REAL( work(i+1,k), KIND=wp )
                 ENDDO
                 DO  i = 1, (nx+1)/2 - 1
@@ -1817 +1817 @@
              DO  k = 1, nz
                 DO  j = 0, (ny+1)/2
-                   ar(j,k) = REAL( work(j+1,k) )
+                   ar(j,k) = REAL( work(j+1,k), KIND=wp )
                 ENDDO
                 DO  j = 1, (ny+1)/2 - 1