Changeset 1023 for codes/icosagcm/trunk/src/time

Timestamp:

03/16/20 23:56:25 (4 years ago)

Author:

ymipsl

Message:

• 1+1=2 : missing sync for rhodz
• add "sync" parameter in run.def to define synchronisation period (in time step) ; default value is define by model (close to every 10 timestep)

YM

File:

• codes/icosagcm/trunk/src/time/timeloop_gcm.F90 (modified) (5 diffs)

codes/icosagcm/trunk/src/time/timeloop_gcm.F90

@@ -12 +12 @@
   
   INTEGER, PARAMETER :: sync_it=10
-  TYPE(t_message),SAVE :: req_ps0, req_mass0, req_theta_rhodz0, req_u0, req_q0, req_W0, req_geopot0
+  TYPE(t_message),SAVE :: req_ps0, req_rhodz0, req_mass0, req_theta_rhodz0, req_u0, req_q0, req_W0, req_geopot0
   LOGICAL, SAVE :: positive_theta
   INTEGER :: itau_prof, id_timeloop, id_dyn, id_phys, id_dissip, id_adv, id_diags
@@ -167 +167 @@
     CALL init_message(f_ps,req_i0,req_ps0)
     CALL init_message(f_mass,req_i0,req_mass0)
+    CALL init_message(f_rhodz,req_i0,req_rhodz0)
     CALL init_message(f_theta_rhodz,req_i0,req_theta_rhodz0)
     CALL init_message(f_u,req_e0_vect,req_u0)
@@ -213 +214 @@
       IF (MOD(86400,INT(i*dt))==0 .AND. ABS((sync_it-itau_sync)*1./sync_it )/sync_it < (sync_it-itau_sync)*1./sync_it) itau_sync=i
     ENDDO
+    CALL getin("itau_sync",itau_sync)
     IF (is_master) PRINT*,"Synchronize frontier every itau_sync =",itau_sync
       
@@ -233 +235 @@
     END DO
     !$OMP BARRIER
+
     fluxt_zero=.TRUE.
 
@@ -287 +290 @@
           CALL send_message(f_mass,req_mass0)
           CALL wait_message(req_mass0)
+          CALL send_message(f_rhodz,req_rhodz0)
+          CALL wait_message(req_rhodz0)
           CALL send_message(f_theta_rhodz,req_theta_rhodz0) 
           CALL wait_message(req_theta_rhodz0)