Changeset 13498 for NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90

Timestamp:

2020-09-21T15:56:43+02:00 (4 years ago)

Author:

jchanut

Message:

#2222 1) Define agrif sponge at f-points 2) Correct bathymetry check with vertical refinement just over needed points

File:

• NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90 (modified) (4 diffs)

NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce_interp.F90

@@ -1589 +1589 @@
       IF(lk_west) THEN
          ispon  = nn_sponge_len * Agrif_irhox()
-         istart = nn_hls + 2                                 ! halo + land + 1
-         iend   = nn_hls + 1 + nbghostcells + ispon          ! halo + land + nbghostcells + sponge
+         istart = nn_hls + 2                                  ! halo + land + 1
+         iend   = nn_hls + 1 + nbghostcells + ispon           ! halo + land + nbghostcells + sponge
          jstart = nn_hls + 2
          jend   = jpjglo - nn_hls - 1
@@ -1611 +1611 @@
       IF(lk_east) THEN
          ispon  = nn_sponge_len * Agrif_irhox() 
-         istart = jpiglo - ( nn_hls + nbghostcells + ispon ) ! halo + land + nbghostcells + sponge - 1
-         iend   = jpiglo - ( nn_hls + 1 )                    ! halo + land + 1                     - 1
+         istart = jpiglo - ( nn_hls + nbghostcells + ispon )  ! halo + land + nbghostcells + sponge - 1
+         iend   = jpiglo - ( nn_hls + 1 )                     ! halo + land + 1                     - 1
          jstart = nn_hls + 2
-         jend   = jpjglo - nn_hls - 1
+         jend   = jpjglo - nn_hls - 1 
          DO ji = mi0(istart), mi1(iend)
             DO jj = mj0(jstart), mj1(jend)
@@ -1623 +1623 @@
             END DO
          END DO
-         DO ji = mi0(istart), mi1(iend-1)
+         DO ji = mi0(istart+1), mi1(iend-1)
             DO jj = mj0(jstart), mj1(jend)
                IF ( ABS(hu0_parent(ji,jj)-hu_0(ji,jj)) > 1.e-3 ) iindic = iindic + 1
@@ -1667 +1667 @@
             END DO
          END DO
-         DO jj = mj0(jstart), mj1(jend-1)
+         DO jj = mj0(jstart+1), mj1(jend-1)
             DO ji = mi0(istart), mi1(iend)
                IF ( ABS(hv0_parent(ji,jj)-hv_0(ji,jj)) > 1.e-3 ) iindic = iindic + 1