Changeset 420 for palm/trunk/SOURCE/disturb_field.f90

Timestamp:

Jan 13, 2010 3:10:53 PM (14 years ago)

Author:

franke

Message:

collision of cloud droplets has changed in advec_particles
bugfixes in advec_particles and collision_efficiency
change in disturb_field to make runs reproducible on HLRN

File:

• palm/trunk/SOURCE/disturb_field.f90 (modified) (5 diffs)

palm/trunk/SOURCE/disturb_field.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! A loop has been splitted to make runs reproducible on HLRN systems
 !
 ! Former revisions:
@@ -98 +98 @@
 !
 !-- Exchange of ghost points for the random perturbation
+
     CALL exchange_horiz( dist1 )
 
@@ -104 +105 @@
 !-- Neighboured grid points in all three directions are used for the
 !-- filter operation.
-    DO  i = nxl, nxr
-       DO  j = nys, nyn
+!-- Loop has been splitted to make runs reproducible on HLRN systems using
+!-- compiler option -O3
+     DO  i = nxl, nxr
+        DO  j = nys, nyn
           DO  k = disturbance_level_ind_b-1, disturbance_level_ind_t+1
-             dist2(k,j,i) = ( dist1(k,j,i-1) + dist1(k,j,i+1) + dist1(k,j-1,i) &
-                            + dist1(k,j+1,i) + dist1(k+1,j,i) + dist1(k-1,j,i) &
+             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
+          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
           ENDDO
-       ENDDO
-    ENDDO
+        ENDDO
+     ENDDO
 
 !
@@ -119 +126 @@
 !-- Afterwards, filter operation and exchange of ghost points are repeated.
     CALL exchange_horiz( dist2 )
+
     DO  i = nxl, nxr
        DO  j = nys, nyn
@@ -129 +137 @@
        ENDDO
     ENDDO
+
     CALL exchange_horiz( dist1 )