Changeset 14721

Timestamp:

2021-04-16T16:11:55+02:00 (3 years ago)

Author:

jchanut

Message:

#2638, (partially) correct bathymetry transitioning in corners + add lk_north=.false. case

File:

• utils/tools/DOMAINcfg/src/agrif_connect.F90 (modified) (3 diffs)

utils/tools/DOMAINcfg/src/agrif_connect.F90

@@ -33 +33 @@
       CALL Agrif_Bc_variable(e3t_connect_id, procname = connect_e3t_connect)
       Agrif_UseSpecialValue = .FALSE.
-      !    
+      !
    END SUBROUTINE agrif_boundary_connections
 
@@ -118 +118 @@
 
         ! Connected bathymetry
-               bathy_local(ji,jj)=(1.-ztabramp(ji,jj))*bathy_local(ji,jj)+ztabramp(ji,jj)*bathy_interp(ji,jj)
+               IF( e3t_interp(ji,jj,1) == -10 ) THEN
+                  bathy_local(ji,jj)=(1.-ztabramp(ji,jj))*bathy_local(ji,jj)+ztabramp(ji,jj)*bathy_interp(ji,jj)
+               ENDIF
             END DO
          END DO
@@ -236 +238 @@
 
       ! --- North --- !
-      IF( (nbondj == 1) .OR. (nbondj == 2) ) THEN
+      IF(( (nbondj == 1) .OR. (nbondj == 2) ).AND.(lk_north)) THEN
          ind2 = jpjglo - (nn_hls + nbghostcells - 1) - istart
          ind1 = ind2 -ispongearea