New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12424 for NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/step.F90 – NEMO

Ignore:
Timestamp:
2020-02-20T16:06:41+01:00 (4 years ago)
Author:
davestorkey
Message:
  1. Rename r2dt -> rDt
  2. Rename r1_2dt -> r1_Dt
  3. Reorganise management of initial Euler timestep for leapfrogging.

This version passes all SETTE tests and bit-compares with the trunk @ 12377

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/step.F90

    r12406 r12424  
    105105      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    106106      ! 
    107       IF( l_1st_euler ) THEN   ;   r2dt =         rn_Dt   ! start or restart with Euler 1st time-step 
    108       ELSE                     ;   r2dt = 2._wp * rn_Dt   ! restart with leapfrog 
    109       ENDIF 
    110       r1_2dt = 1._wp / r2dt 
     107      IF( l_1st_euler ) THEN   
     108         ! start or restart with Euler 1st time-step 
     109         rDt =  rn_Dt    
     110         r1_Dt = 1._wp / rDt 
     111      ENDIF 
    111112      ! 
    112113      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    344345#endif 
    345346      ! 
    346       IF( l_1st_euler ) l_1st_euler = .FALSE.      ! recover Leap-frog timestep 
     347      IF( l_1st_euler ) THEN         ! recover Leap-frog timestep 
     348         rDt = 2._wp * rn_Dt    
     349         r1_Dt = 1._wp / rDt 
     350         l_1st_euler = .FALSE.       
     351      ENDIF 
    347352      ! 
    348353      IF( ln_timing )   CALL timing_stop('stp') 
Note: See TracChangeset for help on using the changeset viewer.