Changeset 541 for codes/icosagcm/trunk/src

Timestamp:

06/21/17 17:48:15 (7 years ago)

Author:

millour

Message:

Fix dissipation initialization to be identical when changing number of processes.

File:

• codes/icosagcm/trunk/src/dissip_gcm.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/dissip_gcm.f90

@@ -87 +87 @@
   
             
-  INTEGER :: i,j,ij,ind,it,iter
+  INTEGER :: i,j,ij,ind,it,iter,M
 
    rayleigh_friction_key='none'
@@ -159 +159 @@
       u=f_u(ind)
 
+! set random seed to get reproductibility when using a different number of process
+      CALL RANDOM_SEED(size=M)
+      CALL RANDOM_SEED(put=(/(i,i=1,M)/))
+
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end
@@ -247 +251 @@
       u=f_u(ind)
 
-      DO j=jj_begin,jj_end
+! set random seed to get reproductibility when using a different number of process
+      CALL RANDOM_SEED(size=M)
+      CALL RANDOM_SEED(put=(/(i,i=1,M)/))
+ 
+       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end
           ij=(j-1)*iim+i   
@@ -335 +343 @@
       CALL swap_geometry(ind)
       theta=f_theta(ind)
- 
+
+! set random seed to get reproductibility when using a different number of process
+      CALL RANDOM_SEED(size=M)
+      CALL RANDOM_SEED(put=(/(i,i=1,M)/))
+
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end