- Timestamp:
- 2016-11-30T17:56:53+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90
r6140 r7403 57 57 INTEGER :: nadv ! chosen advection scheme 58 58 ! 59 REAL(wp) :: r2dttrc ! vertical profile time-step, = 2 rdt60 ! ! except at nitrrc000 (=rdt) if neuler=061 62 59 !! * Substitutions 63 60 # include "vectopt_loop_substitute.h90" … … 87 84 ! 88 85 CALL wrk_alloc( jpi,jpj,jpk, zun, zvn, zwn ) 89 !90 IF( ( neuler == 0 .AND. kt == nittrc000 ) .OR. ln_top_euler ) THEN ! at nittrc00091 r2dttrc = rdttrc ! = rdttrc (use or restarting with Euler time stepping)92 ELSEIF( kt <= nittrc000 + nn_dttrc ) THEN ! at nittrc000 or nittrc000+193 r2dttrc = 2. * rdttrc ! = 2 rdttrc (leapfrog)94 ENDIF95 86 ! !== effective transport ==! 96 87 DO jk = 1, jpkm1
Note: See TracChangeset
for help on using the changeset viewer.