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

palm/trunk/SOURCE/disturb_field.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! REAL constants defined as wp-kind
 !
 ! Former revisions:
@@ -98 +98 @@
           DO  j = dist_nys(dist_range), dist_nyn(dist_range)
              DO  k = disturbance_level_ind_b, disturbance_level_ind_t
-                randomnumber = 3.0 * disturbance_amplitude *                   &
-                               ( random_function( iran ) - 0.5 )
+                randomnumber = 3.0_wp * disturbance_amplitude *                &
+                               ( random_function( iran ) - 0.5_wp )
                 IF ( nxl <= i  .AND.  nxr >= i  .AND.  nys <= j  .AND.         &
                      nyn >= j )                                                &
@@ -113 +113 @@
              DO  k = disturbance_level_ind_b, disturbance_level_ind_t
 #if defined( __nec )
-                randomnumber = 3.0 * disturbance_amplitude *                   &
-                               ( RANDOM( 0 ) - 0.5 )
+                randomnumber = 3.0_wp * disturbance_amplitude *                &
+                               ( RANDOM( 0 ) - 0.5_wp )
 #else
                 CALL RANDOM_NUMBER( randomnumber )
-                randomnumber = 3.0 * disturbance_amplitude *                   &
-                                ( randomnumber - 0.5 )
+                randomnumber = 3.0_wp * disturbance_amplitude *                &
+                                ( randomnumber - 0.5_wp )
 #endif
                 IF ( nxl <= i .AND. nxr >= i .AND. nys <= j .AND. nyn >= j )   &
@@ -145 +145 @@
              dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1)                  &
                             + dist1(k,j+1,i) + dist1(k+1,j,i)                  &
-                            ) / 12.0
+                            ) / 12.0_wp
           ENDDO
           DO  k = disturbance_level_ind_b-1, disturbance_level_ind_t+1
               dist2(k,j,i) = dist2(k,j,i) + ( dist1(k,j-1,i) + dist1(k-1,j,i)  &
-                            + 6.0 * dist1(k,j,i)                               &
-                            ) / 12.0
+                            + 6.0_wp * dist1(k,j,i)                            &
+                            ) / 12.0_wp
           ENDDO
         ENDDO
@@ -165 +165 @@
              dist1(k,j,i) = ( dist2(k,j,i-1) + dist2(k,j,i+1) + dist2(k,j-1,i) &
                             + dist2(k,j+1,i) + dist2(k+1,j,i) + dist2(k-1,j,i) &
-                            + 6.0 * dist2(k,j,i)                               &
-                            ) / 12.0
+                            + 6.0_wp * dist2(k,j,i)                            &
+                            ) / 12.0_wp
           ENDDO
        ENDDO