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

palm/trunk/SOURCE/init_3d_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Parallel random number generator added
 ! 
 ! Former revisions:
@@ -198 +198 @@
     
     USE random_function_mod 
+    
+    USE random_generator_parallel,                                             &
+        ONLY:  random_number_parallel, random_seed_parallel, random_dummy,     &
+               id_random_array, seq_random_array
     
     USE statistics,                                                            &
@@ -481 +485 @@
 
 !
+!-- ID-array and state-space-array for the parallel random number generator
+    IF ( random_generator == 'random-parallel' )  THEN
+       ALLOCATE ( seq_random_array(5,nysg:nyng,nxlg:nxrg) )
+       ALLOCATE ( id_random_array(0:ny,0:nx) )
+       seq_random_array = 0
+       id_random_array  = 0
+    ENDIF
+    
+!
 !-- 3D-arrays for the leaf area density and the canopy drag coefficient
     IF ( plant_canopy ) THEN
@@ -913 +926 @@
 
 !
+!--    Initialize the random number generators (from numerical recipes)
+       CALL random_function_ini
+       IF ( random_generator == 'random-parallel' )  THEN
+          CALL random_function_ini
+          !- asigning an ID to every vertical gridpoint column
+          DO j=0,ny
+             DO i=0,nx
+                id_random_array(j,i) = 1000000*i + j
+             END DO
+          ENDDO
+          !- executing the initial ran_parallel
+          DO j = nysg, nyng
+             DO i = nxlg, nxrg
+                CALL random_seed_parallel (random_sequence=id_random_array(j, i))
+                CALL random_number_parallel (random_dummy)
+                CALL random_seed_parallel (get=seq_random_array(:, j, i))
+             END DO
+          ENDDO
+       ENDIF
+
+!
 !--    Initialize fluxes at bottom surface
        IF ( use_surface_fluxes )  THEN
@@ -1098 +1132 @@
           q(nzb,:,:) = q(nzb,:,:) + q_surface_initial_change
        ENDIF
-!
-!--    Initialize the random number generator (from numerical recipes)
-       CALL random_function_ini
 
 !