Changeset 15564 for NEMO/releases/r4.0/r4.0-HEAD/src/NST/agrif_top_update.F90

Timestamp:

2021-12-01T17:35:59+01:00 (3 years ago)

Author:

jchanut

Message:

#2751, fixes small constancy preservation issue with AGRIF

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/NST/agrif_top_update.F90 (modified) (2 diffs)

NEMO/releases/r4.0/r4.0-HEAD/src/NST/agrif_top_update.F90

@@ -173 +173 @@
       !!
       INTEGER :: ji,jj,jk,jn
-      REAL(wp) :: ztb, ztnu, ztno
+      REAL(wp) :: ztb, ze3b, ztnu, ztno
       !!----------------------------------------------------------------------
       !
@@ -204 +204 @@
                      DO ji=i1,i2
                         IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN
-                           ztb  = trb(ji,jj,jk,jn) * e3t_b(ji,jj,jk) ! fse3t_b prior update should be used
+
+                           ze3b = e3t_b(ji,jj,jk) & ! Recover e3t_b before update
+                                & - rn_atfp * ( e3t_n(ji,jj,jk) - e3t_a(ji,jj,jk) )
+                           ztb  = trb(ji,jj,jk,jn) * ze3b 
                            ztnu = tabres(ji,jj,jk,jn)
                            ztno = trn(ji,jj,jk,jn) * e3t_a(ji,jj,jk)