- Timestamp:
- 2020-10-14T18:31:28+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/stpMLF.F90
r13427 r13608 1 MODULE st epMLF1 MODULE stpMLF 2 2 !!====================================================================== 3 !! *** MODULE st ep***3 !! *** MODULE stpMLF *** 4 4 !! Time-stepping : manager of the ocean, tracer and ice time stepping 5 !! using Modified Leap Frog for OCE 5 6 !!====================================================================== 6 7 !! History : OPA ! 1991-03 (G. Madec) Original code … … 32 33 !! 4.0 ! 2017-05 (G. Madec) introduction of the vertical physics manager (zdfphy) 33 34 !! 4.1 ! 2019-08 (A. Coward, D. Storkey) rewrite in preparation for new timestepping scheme 34 !! 4.x ! 2020-08 (S. Techene, G. Madec) quasi eulerian coordinate time stepping 35 !! 4.x ! 2020-08 (S. Techene, G. Madec) quasi eulerian coordinate time stepping 35 36 !!---------------------------------------------------------------------- 36 37 … … 45 46 USE step_oce ! time stepping definition modules 46 47 ! 47 USE iom ! xIOs server 48 USE domqco 48 USE domqco ! quasi-eulerian coordinate 49 49 USE traatfqco ! time filtering (tra_atf_qco routine) 50 50 USE dynatfqco ! time filtering (dyn_atf_qco routine) … … 57 57 PUBLIC stp_MLF ! called by nemogcm.F90 58 58 59 ! !----------------------------------------------------------------------60 !! time level indices61 !!---------------------------------------------------------------------- 62 INTEGER, PUBLIC :: Nbb, Nnn, Naa, Nrhs !! used by nemo_init59 ! !** time level indices **! 60 INTEGER, PUBLIC :: Nbb, Nnn, Naa, Nrhs !: used by nemo_init 61 62 !! * Substitutions 63 63 # include "domzgr_substitute.h90" 64 64 !!---------------------------------------------------------------------- … … 94 94 INTEGER :: ji, jj, jk ! dummy loop indice 95 95 INTEGER :: indic ! error indicator if < 0 96 !!gm kcall can be removed, I guess97 INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt)98 !!st patch99 96 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgdept 100 97 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zssh_f … … 119 116 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 120 117 ! 121 IF( l_1st_euler ) THEN 122 ! start or restart with Euler 1st time-step 123 rDt = rn_Dt 118 IF( l_1st_euler ) THEN ! start or restart with Euler 1st time-step 119 rDt = rn_Dt 124 120 r1_Dt = 1._wp / rDt 125 121 ENDIF 126 122 ! 127 123 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 128 129 124 ! update I/O and calendar 130 125 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< … … 374 369 ! 375 370 IF( l_1st_euler ) THEN ! recover Leap-frog timestep 376 rDt = 2._wp * rn_Dt371 rDt = 2._wp * rn_Dt 377 372 r1_Dt = 1._wp / rDt 378 373 l_1st_euler = .FALSE. … … 477 472 478 473 !!====================================================================== 479 END MODULE st epMLF474 END MODULE stpMLF
Note: See TracChangeset
for help on using the changeset viewer.