Changeset 15437

Timestamp:

2021-10-22T14:21:20+02:00 (2 years ago)

Author:

jchanut

Message:

agrif fixes: i) Add masking prior interpolation on "before" tracers: these seem to be strangly unmasked which induces extrapolation issues. ii) Add flux limiter on tracer sponge (recall sponge is applied on anomalies) to prevent from generating spurious extrema.

Location:

NEMO/trunk/src/NST

Files:

• agrif_oce_interp.F90 (modified) (1 diff)
• agrif_oce_sponge.F90 (modified) (6 diffs)
• agrif_top_sponge.F90 (modified) (1 diff)

NEMO/trunk/src/NST/agrif_oce_interp.F90

@@ -1578 +1578 @@
          DO jj = j1, j2
             DO ji = i1, i2
-               zh = 0._wp
-               DO jk = 1, mbkt_parent(ji, jj)-1
-                  zh = zh + e3t0_parent(ji,jj,jk)
-               END DO
-               e3t0_parent(ji,jj,mbkt_parent(ji,jj)) = ht0_parent(ji, jj) - zh
+               IF ( mbkt_parent(ji,jj) > 1 ) THEN
+                  zh = 0._wp
+                  DO jk = 1, mbkt_parent(ji, jj)-1
+                     zh = zh + e3t0_parent(ji,jj,jk)
+                  END DO
+                  e3t0_parent(ji,jj,mbkt_parent(ji,jj)) = ht0_parent(ji, jj) - zh
+               ENDIF  
             END DO
          END DO

NEMO/trunk/src/NST/agrif_oce_sponge.F90

@@ -386 +386 @@
       INTEGER  ::   iku, ikv
       REAL(wp) :: ztsa, zabe1, zabe2, zbtr, zhtot
+      REAl(wp) :: zflag, zdmod, zdtot
       REAL(wp), DIMENSION(i1-1:i2,j1-1:j2,jpk) :: ztu, ztv
       REAL(wp), DIMENSION(i1:i2,j1:j2,jpk,n1:n2) ::tsbdiff
@@ -401 +402 @@
                DO jj=j1,j2
                   DO ji=i1,i2
-                     tabres(ji,jj,jk,jn) = ts(ji,jj,jk,jn,Kbb_a)
+                     ! JC: masking is mandatory here: before tracer field seems 
+                     !     to hold non zero values where tmask=0
+                     tabres(ji,jj,jk,jn) = ts(ji,jj,jk,jn,Kbb_a) * tmask(ji,jj,jk)
                   END DO
                END DO
@@ -536 +539 @@
                   DO ji = i1,i2-1
                      zabe1 = rn_sponge_tra * r1_Dt * umask(ji,jj,jk) * e1e2u(ji,jj) * e3u(ji,jj,jk,Kmm_a)
-                     ztu(ji,jj,jk) = zabe1 * fspu(ji,jj) * ( tsbdiff(ji+1,jj  ,jk,jn) - tsbdiff(ji,jj,jk,jn) ) 
+                     zdtot =  tsbdiff(ji+1,jj,jk,jn) -  tsbdiff(ji,jj,jk,jn) 
+                     zdmod =       ts(ji+1,jj,jk,jn,Kbb_a) - ts(ji,jj,jk,jn,Kbb_a)
+                     zflag = 0.5_wp + SIGN(0.5_wp, zdtot*zdmod)
+                     ztu(ji,jj,jk) = zabe1 * fspu(ji,jj) * ( zflag * zdtot + (1._wp - zflag) * zdmod ) 
                   END DO
                END DO
@@ -544 +550 @@
                      zabe2 = rn_sponge_tra * r1_Dt * vmask(ji,jj,jk) * e1e2v(ji,jj) * e3v(ji,jj,jk,Kmm_a)
                      ztv(ji,jj,jk) = zabe2 * fspv(ji,jj) * ( tsbdiff(ji  ,jj+1,jk,jn) - tsbdiff(ji,jj,jk,jn) )
+                     zdtot =  tsbdiff(ji,jj+1,jk,jn) -  tsbdiff(ji,jj,jk,jn) 
+                     zdmod =       ts(ji,jj+1,jk,jn,Kbb_a) - ts(ji,jj,jk,jn,Kbb_a)
+                     zflag = 0.5_wp + SIGN(0.5_wp, zdtot*zdmod)
+                     ztv(ji,jj,jk) = zabe2 * fspv(ji,jj) * ( zflag * zdtot + (1._wp - zflag) * zdmod ) 
                   END DO
                END DO
@@ -612 +622 @@
             DO jj=j1,j2
                DO ji=i1,i2
-                  tabres(ji,jj,jk,m1) = uu(ji,jj,jk,Kbb_a)
+                  tabres(ji,jj,jk,m1) = uu(ji,jj,jk,Kbb_a) * umask(ji,jj,jk)
                END DO
             END DO
@@ -797 +807 @@
             DO jj=j1,j2
                DO ji=i1,i2
-                  tabres(ji,jj,jk,m1) = vv(ji,jj,jk,Kbb_a)
+                  tabres(ji,jj,jk,m1) = vv(ji,jj,jk,Kbb_a) * vmask(ji,jj,jk)
                END DO
             END DO

NEMO/trunk/src/NST/agrif_top_sponge.F90

@@ -90 +90 @@
                DO jj=j1,j2
                   DO ji=i1,i2
-                     tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kbb_a)
+                     tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kbb_a) * tmask(ji,jj,jk) 
                   END DO
                END DO