Changeset 9863 for NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OFF/nemogcm.F90

Timestamp:

2018-06-30T12:51:02+02:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-04): simplified implementation of the Euler stepping at nit000

File:

• NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OFF/nemogcm.F90 (modified) (2 diffs)

NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OFF/nemogcm.F90

@@ -100 +100 @@
       !                            !==   time stepping   ==!
       !                            !-----------------------!
+      !                                               !== set the model time-step  ==!
+      !
+      IF( l_1st_euler ) THEN   ;   r2dt =         rn_rdt   ;   l_1st_euler = .TRUE.    ! start or restart with Euler 1st time-step
+      ELSE                     ;   r2dt = 2._wp * rn_rdt   ;   l_1st_euler = .FALSE.   ! restart with leapfrog
+      ENDIF
+      r1_2dt = 1._wp / r2dt
+      ! NB: if l_1st_euler=T, r2dt will be set to 2*rdt at the end of the 1st time-step (in step.F90)
+      !     Done here (not in domain.F90) as in ASM initialization an Euler 1st time step can be forced
+      !
+      !
       istp = nit000
       !
@@ -115 +125 @@
                                 CALL trc_stp    ( istp )         ! time-stepping
                                 CALL stp_ctl    ( istp, indic )  ! Time loop: control and print
+         IF( l_1st_euler ) THEN
+            r2dt   = 2._wp * rn_rdt                              ! recover Leap-frog time-step
+            r1_2dt = 1._wp / r2dt
+            l_1st_euler = .FALSE.
+         ENDIF
+         !
          istp = istp + 1
       END DO