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

Timestamp:

May 9, 2014 2:03:54 PM (10 years ago)

Author:

knoop

Message:

Parallel random number generator added (preliminary version).

File:

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

palm/trunk/SOURCE/disturb_field.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Parallel random number generator added
 ! 
 ! Former revisions:
@@ -72 +72 @@
     USE random_function_mod,                                                   &
         ONLY: random_function
+        
+    USE random_generator_parallel,                                             &
+        ONLY:  random_number_parallel, random_seed_parallel, random_dummy,     &
+               id_random_array, seq_random_array
 
     IMPLICIT NONE
@@ -112 +116 @@
                 ENDIF
              ENDDO
+          ENDDO
+       ENDDO
+    ELSEIF ( random_generator == 'random-parallel' )  THEN
+       DO  i = dist_nxl(dist_range), dist_nxr(dist_range)
+          DO  j = dist_nys(dist_range), dist_nyn(dist_range)
+             CALL random_seed_parallel( put=seq_random_array(:, j, i) )
+             DO  k = disturbance_level_ind_b, disturbance_level_ind_t
+                CALL random_number_parallel( random_dummy )
+                randomnumber = 3.0_wp * disturbance_amplitude *                &
+                               ( random_dummy - 0.5_wp )
+                IF ( nxl <= i  .AND.  nxr >= i  .AND.  nys <= j  .AND.         &
+                     nyn >= j )                                                &
+                THEN
+                   dist1(k,j,i) = randomnumber
+                ENDIF
+             ENDDO
+!              IF (myid == 0) PRINT*, nxl, i, nxr, i, nys, j, nyn, j
+             CALL random_seed_parallel( get=seq_random_array(:, j, i) )
           ENDDO
        ENDDO