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

Timestamp:

Nov 21, 2017, 4:03:44 PM (7 years ago)

Author:

emillour

Message:

Common dynamics:

• enable possiblity to store multiple time steps in the restart.nc file (flag "ecritstart" gives the frequency, in dynamical steps).
• fixed dynredem_mod.F90 to correctly write multiple time steps.
• fixed computation of JH_cur in the mars case where "hour_ini" contains the initial time of day read from the start.nc file
• minor fix in dynetat0.F90

RY

File:

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

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

@@ -19 +19 @@
      &                       dissip_period,offline,ip_ebil_dyn,
      &                       ok_dynzon,periodav,ok_dyn_ave,iecri,
-     &                       ok_dyn_ins,output_grads_dyn
+     &                       ok_dyn_ins,output_grads_dyn,ecritstart
       use exner_hyb_m, only: exner_hyb
       use exner_milieu_m, only: exner_milieu
@@ -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
+     .                  start_time,dt,hour_ini
 
       IMPLICIT NONE
@@ -131 +131 @@
 c   variables pour le fichier histoire
 !      REAL dtav      ! intervalle de temps elementaire
+      LOGICAL lrestart
 
       REAL tppn(iim),tpps(iim),tpn,tps
@@ -289 +290 @@
 c   ----------------------------------
 
+c     RMBY: check that hour_ini and start_time are not both non-zero
+      if ((hour_ini.ne.0.0).and.(start_time.ne.0.0)) then
+        write(*,*) "ERROR: hour_ini = ", hour_ini, 
+     &             "start_time = ", start_time
+        abort_message = 'hour_ini and start_time both nonzero'
+        call abort_gcm(modname,abort_message,1)
+      endif
+
    1  CONTINUE ! Matsuno Forward step begins here
 
@@ -296 +305 @@
       jD_cur = jD_ref + day_ini - day_ref +                             &
      &          (itau+1)/day_step 
-      jH_cur = jH_ref + start_time +                                    &
-     &          mod(itau+1,day_step)/float(day_step) 
+      IF (planet_type .eq. "mars") THEN
+        jH_cur = jH_ref + hour_ini +                                    &
+     &           mod(itau+1,day_step)/float(day_step) 
+      ELSE
+        jH_cur = jH_ref + start_time +                                  &
+     &           mod(itau+1,day_step)/float(day_step)
+      ENDIF
       jD_cur = jD_cur + int(jH_cur)
       jH_cur = jH_cur - int(jH_cur)
@@ -351 +365 @@
         jD_cur = jD_ref + day_ini - day_ref +
      &            (itau+1)/day_step
-        jH_cur = jH_ref + start_time +
-     &            mod(itau+1,day_step)/float(day_step) 
+        IF (planet_type .eq. "mars") THEN
+          jH_cur = jH_ref + hour_ini + 
+     &             mod(itau+1,day_step)/float(day_step) 
+        ELSE
+          jH_cur = jH_ref + start_time +
+     &             mod(itau+1,day_step)/float(day_step)
+        ENDIF
         jD_cur = jD_cur + int(jH_cur)
         jH_cur = jH_cur - int(jH_cur)
@@ -529 +548 @@
            ENDIF
 
-           jH_cur = jH_ref + start_time +                               &
-     &          mod(itau+1,day_step)/float(day_step) 
-           IF ((planet_type .eq."generic").or.
-     &         (planet_type .eq."mars")) THEN
+           IF (planet_type .eq. "mars") THEN
+             jH_cur = jH_ref + hour_ini +                                 &
+     &                mod(itau,day_step)/float(day_step)
+           ELSE IF (planet_type .eq. "generic") THEN
              jH_cur = jH_ref + start_time +                               &
-     &          mod(itau,day_step)/float(day_step)
+     &                mod(itau,day_step)/float(day_step)
+           ELSE
+             jH_cur = jH_ref + start_time +                               &
+     &                mod(itau+1,day_step)/float(day_step)
            ENDIF
            jD_cur = jD_cur + int(jH_cur)
@@ -861 +883 @@
             ENDIF ! of IF(MOD(itau,iecri).EQ.0)
 
-            IF(itau.EQ.itaufin) THEN
-
+c           Determine whether to write to the restart.nc file
+c           Decision can't be made in one IF statement as if
+c           ecritstart==0 there will be a divide-by-zero error
+            lrestart = .false.
+            IF (itau.EQ.itaufin) THEN
+              lrestart = .true.
+            ELSE IF (ecritstart.GT.0) THEN
+              IF (MOD(itau,ecritstart).EQ.0) lrestart  = .true.
+            ENDIF
+
+c           Write to the restart.nc file
+            IF (lrestart) THEN
               if (planet_type=="mars") then
                 CALL dynredem1("restart.nc",REAL(itau)/REAL(day_step),
@@ -872 +904 @@
               CLOSE(99)
               !!! Ehouarn: Why not stop here and now?
-            ENDIF ! of IF (itau.EQ.itaufin)
+            ENDIF ! of IF (lrestart)
 
 c-----------------------------------------------------------------------
@@ -992 +1024 @@
               ENDIF ! of IF(MOD(itau,iecri         ).EQ.0) 
 
-              IF(itau.EQ.itaufin) THEN
+c             Determine whether to write to the restart.nc file
+c             Decision can't be made in one IF statement as if
+c             ecritstart==0 there will be a divide-by-zero error
+              lrestart = .false.
+              IF (itau.EQ.itaufin) THEN
+                lrestart = .true.
+              ELSE IF (ecritstart.GT.0) THEN
+                IF (MOD(itau,ecritstart).EQ.0) lrestart = .true.
+              ENDIF
+
+c             Write to the restart.nc file
+              IF (lrestart) THEN
                 if (planet_type=="mars") then
                   CALL dynredem1("restart.nc",REAL(itau)/REAL(day_step),
@@ -1000 +1043 @@
      &                         vcov,ucov,teta,q,masse,ps)
                 endif
-              ENDIF ! of IF(itau.EQ.itaufin)
+              ENDIF ! of IF (lrestart)
 
               forward = .TRUE.