Changeset 262 for codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/timeloop_gcm.f90

Timestamp:

07/29/14 18:30:42 (10 years ago)

Author:

ymipsl

Message:

Save starting iteration itau0 in dynamico restart file.
Timeloop now begin from itau0 which can be read from start file or set to 0 by default.

YM

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/timeloop_gcm.f90 (modified) (2 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/timeloop_gcm.f90

@@ -197 +197 @@
 !$OMP END MASTER   
 
+  CALL check_conserve(f_ps,f_dps,f_u,f_theta_rhodz,f_phis,itau0)  
+
   CALL trace_on
   
-  DO it=0,itaumax
+  DO it=itau0+1,itau0+itaumax
     
     IF (xios_output) CALL xios_update_calendar(it)
-    IF (MOD(it,itau_sync)==0) THEN
+    IF (it==itau0+1 .OR. MOD(it,itau_sync)==0) THEN
       CALL send_message(f_ps,req_ps0)
       CALL wait_message(req_ps0)
@@ -294 +296 @@
   ENDDO
 
-  CALL write_etat0(f_ps, f_phis,f_theta_rhodz,f_u,f_q) 
+  CALL write_etat0(itau0+itaumax,f_ps, f_phis,f_theta_rhodz,f_u,f_q) 
 
   CALL check_conserve(f_ps,f_dps,f_u,f_theta_rhodz,f_phis,it)