Changeset 2507 for trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

Timestamp:

Apr 28, 2021, 1:00:13 PM (4 years ago)

Author:

romain.vande

Message:

For LMDZ MARS: Update of day_ini, time and hour_ini in restart and retstartfi.
Hour_ini is obsolete. If we write one restart file: day_ini is the last day
of the simulation and the remaining time is in Time (often=0), if we write
multiple restart nothing changes

File:

• trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F (modified) (7 diffs)

trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

@@ -31 +31 @@
      .                  statcl,conser,apdiss,purmats,tidal,ok_strato
       USE temps_mod, ONLY: jD_ref,jH_ref,itaufin,day_ini,day_ref,
-     .                  start_time,dt,hour_ini
+     .                  start_time,dt,hour_ini,day_end
 
       IMPLICIT NONE
@@ -278 +278 @@
             abort_message = 'PK non strictement decroissante'
             call abort_gcm(modname,abort_message,1)
-c           write(*,*) "ATTENTION, Test PK deconnecté..."
+c           write(*,*) "ATTENTION, Test PK deconnecte..."
           endif
         enddo
@@ -500 +500 @@
        IF ((planet_type.eq."titan").and.(tidal)) then
 c-----------------------------------------------------------------------
-c   Marées gravitationnelles causées par Saturne
+c   Marees gravitationnelles causees par Saturne
 c   B. Charnay (28/10/2010)
 c   ----------------------------------------------------------
@@ -579 +579 @@
 c+jld
 
-c  Diagnostique de conservation de l'énergie : initialisation
+c  Diagnostique de conservation de l'energie : initialisation
          IF (ip_ebil_dyn.ge.1 ) THEN 
           ztit='bil dyn'
@@ -635 +635 @@
          ENDIF
 
-c  Diagnostique de conservation de l'énergie : difference
+c  Diagnostique de conservation de l'energie : difference
          IF (ip_ebil_dyn.ge.1 ) THEN 
           ztit='bil phys'
@@ -907 +907 @@
             IF (lrestart) THEN
               if (planet_type=="mars") then
-                CALL dynredem1("restart.nc",REAL(itau)/REAL(day_step),
-     &                         vcov,ucov,teta,q,masse,ps)
+                if(ecritstart.GT.0) then
+                 CALL dynredem1("restart.nc",
+     &                        REAL(itau)/REAL(day_step),
+     &                        vcov,ucov,teta,q,masse,ps)
+                else
+                 CALL dynredem1("restart.nc",
+     &             REAL(itau)/REAL(day_step)-(day_end-day_ini),
+     &                        vcov,ucov,teta,q,masse,ps)
+                endif
               else
                 CALL dynredem1("restart.nc",start_time,
@@ -1053 +1060 @@
               IF (lrestart) THEN
                 if (planet_type=="mars") then
-                  CALL dynredem1("restart.nc",REAL(itau)/REAL(day_step),
+                  if(ecritstart.GT.0) then
+                        CALL dynredem1("restart.nc",
+     &                         REAL(itau)/REAL(day_step),
      &                         vcov,ucov,teta,q,masse,ps)
+                                  else
+                        CALL dynredem1("restart.nc",
+     &          REAL(itau)/REAL(day_step)-(day_end-day_ini),
+     &                         vcov,ucov,teta,q,masse,ps)
+                                  endif
                 else
                   CALL dynredem1("restart.nc",start_time,