Ignore:
Timestamp:
2020-02-28T16:55:11+01:00 (14 months ago)
Author:
davestorkey
Message:

Preparation for new timestepping scheme #2390.
Main changes:

  1. Initial euler timestep now handled in stp and not in TRA/DYN routines.
  2. Renaming of all timestep parameters. In summary, the namelist parameter is now rn_Dt and the current timestep is rDt (and rDt_ice, rDt_trc etc).
  3. Renaming of a few miscellaneous parameters, eg. atfp → rn_atfp (namelist parameter used everywhere) and rau0 → rho0.

This version gives bit-comparable results to the previous version of the trunk.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/step.F90

    r12377 r12489  
    102102      ! 
    103103      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     104      ! model timestep 
     105      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
     106      ! 
     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 
     112      ! 
     113      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    104114      ! update I/O and calendar  
    105115      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
     
    335345#endif 
    336346      ! 
     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 
     352      ! 
    337353      IF( ln_timing )   CALL timing_stop('stp') 
    338354      ! 
Note: See TracChangeset for help on using the changeset viewer.