Changeset 6349 for branches/UKMO/product_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

Timestamp:

2016-02-24T14:18:49+01:00 (8 years ago)

Author:

davestorkey

Message:

Temporary fix for bug in NEMO whereby grid cell thicknesses were written
to diagnostics at wrong time level. NB. This has now been fixed at the
latest version of 3.6_stable (ticket #1682). Only putting the fix in this
branch in order to rerun some old experiments.

File:

• branches/UKMO/product_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90 (modified) (2 diffs)

branches/UKMO/product_diagnostics/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

@@ -594 +594 @@
          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   - interpolate scale factors and compute depths for next time step'
       ENDIF
+
+      ! Write outputs
+      ! =============
+      CALL iom_put(     "e3t" , fse3t_n  (:,:,:) )
+      CALL iom_put(     "e3u" , fse3u_n  (:,:,:) )
+      CALL iom_put(     "e3v" , fse3v_n  (:,:,:) )
+      CALL iom_put(     "e3w" , fse3w_n  (:,:,:) )
+      CALL iom_put( "tpt_dep" , fsde3w_n (:,:,:) )
+      IF( iom_use("e3tdef") )   &
+         CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 )
+
       !
       ! Time filter and swap of scale factors
@@ -666 +677 @@
       END DO
 
-      ! Write outputs
-      ! =============
-      CALL iom_put(     "e3t" , fse3t_n  (:,:,:) )
-      CALL iom_put(     "e3u" , fse3u_n  (:,:,:) )
-      CALL iom_put(     "e3v" , fse3v_n  (:,:,:) )
-      CALL iom_put(     "e3w" , fse3w_n  (:,:,:) )
-      CALL iom_put( "tpt_dep" , fsde3w_n (:,:,:) )
-      IF( iom_use("e3tdef") )   &
-         CALL iom_put( "e3tdef"  , ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 )
 
       ! write restart file