Changeset 15113

Timestamp:

2021-07-09T10:04:38+02:00 (3 years ago)

Author:

jchanut

Message:

Correct AGRIF divergence conserving interpolation with cyclic zooms (i.e. nbghostcells_x/nbghostcells_y_s /= nbghostcells), #2638

File:

• NEMO/trunk/src/NST/agrif_oce_interp.F90 (modified) (2 diffs)

NEMO/trunk/src/NST/agrif_oce_interp.F90

@@ -1418 +1418 @@
             DO jj=j1,j2
                IF (utint_stage(ji,jj)==0) THEN 
-                  zx = 2._wp*MOD(ABS(mig0(ji)-nbghostcells), INT(Agrif_Rhox()))/zrhox - 1._wp  
+                  zx = 2._wp*MOD(ABS(mig0(ji)-nbghostcells_x), INT(Agrif_Rhox()))/zrhox - 1._wp  
                   ubdy(ji,jj) = ubdy(ji,jj) + 0.25_wp*(1._wp-zx*zx) * ptab(ji,jj) & 
                               &         / zrhoy *r1_e2u(ji,jj) * umask(ji,jj,1) 
@@ -1522 +1522 @@
             DO jj=j1,j2
                IF (vtint_stage(ji,jj)==0) THEN 
-                  zy = 2._wp*MOD(ABS(mjg0(jj)-nbghostcells), INT(Agrif_Rhoy()))/zrhoy - 1._wp  
+                  zy = 2._wp*MOD(ABS(mjg0(jj)-nbghostcells_y_s), INT(Agrif_Rhoy()))/zrhoy - 1._wp  
                   vbdy(ji,jj) = vbdy(ji,jj) + 0.25_wp*(1._wp-zy*zy) * ptab(ji,jj) & 
                               &         / zrhox * r1_e1v(ji,jj) * vmask(ji,jj,1)