Changeset 8535 for branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

Timestamp:

2017-09-18T18:36:54+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- trdtrc -- found ATF trend bug relative to TRC case

File:

• branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 (modified) (3 diffs)

branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

@@ -49 +49 @@
    USE agrif_opa_interp
 #endif
+
 
    IMPLICIT NONE
@@ -314 +315 @@
       IF( kt == kit000 )  THEN
          IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) 'tra_nxt_vvl : time stepping', cdtype
+         IF(lwp) WRITE(numout,*) 'tra_nxt_vvl : time stepping ', cdtype
          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
       ENDIF
@@ -340 +341 @@
       DO jn = 1, kjpt      
          DO jk = 1, jpkm1
-            zfact = 1._wp / r2dtra(jk)
+            !! zfact = 1._wp / r2dtra(jk)
+            zfact = 0.5_wp *p2dt(jk)
             zfact1 = atfp * p2dt(jk)
             zfact2 = zfact1 / rau0