Changeset 9939 for NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OFF/nemogcm.F90

Timestamp:

2018-07-13T09:28:50+02:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-04): RK3 branche phased with MLF@9937 branche

Location:

NEMO/branches/2018/dev_r9838_ENHANCE04_RK3

Files:

• . (copied) (copied from NEMO/trunk)
• src/OFF/nemogcm.F90 (modified) (2 diffs)

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

@@ -100 +100 @@
       !                            !==   time stepping   ==!
       !                            !-----------------------!
+      !                                               !== set the model time-step  ==!
+      !
+      IF( l_1st_euler ) THEN   ;   rDt =         rn_Dt   ;   l_1st_euler = .TRUE.    ! start or restart with Euler 1st time-step
+      ELSE                     ;   rDt = 2._wp * rn_Dt   ;   l_1st_euler = .FALSE.   ! restart with leapfrog
+      ENDIF
+      r1_Dt = 1._wp / rDt
+      ! NB: if l_1st_euler=T, rDt will be set to 2*rn_Dt at the end of the 1st time-step (see the DO WHILE below)
+      !
+      !
       istp = nit000
       !
@@ -115 +124 @@
                                 CALL trc_stp    ( istp )         ! time-stepping
                                 CALL stp_ctl    ( istp, indic )  ! Time loop: control and print
+         IF( l_1st_euler ) THEN
+            rDt   = 2._wp * rDt                                  ! recover Leap-frog time-step
+            r1_Dt = 1._wp / rDt
+            l_1st_euler = .FALSE.
+         ENDIF
+         !
          istp = istp + 1
       END DO