Changeset 13204 for utils/tools/DOMAINcfg/src/agrif_domzgr.F90

Timestamp:

2020-07-02T10:38:35+02:00 (4 years ago)

Author:

smasson

Message:

tools: update with tools_dev_r12970_AGRIF_CMEMS

File:

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

utils/tools/DOMAINcfg/src/agrif_domzgr.F90

@@ +1 @@
+MODULE agrif_domzgr
+
+   USE agrif_profiles
+   USE dom_oce
+
+   IMPLICIT NONE
+   PRIVATE
+  
+   PUBLIC :: agrif_create_bathy_meter
+ 
+
+CONTAINS
+
 #if defined key_agrif
-subroutine agrif_domzgr
-end subroutine agrif_domzgr
 
-subroutine agrif_create_bathy_meter
-use agrif_profiles
-external :: init_bathy
+   SUBROUTINE agrif_create_bathy_meter
 
-call Agrif_Init_variable(bathy_id, procname = init_bathy)
+      CALL Agrif_Init_variable(bathy_id, procname = init_bathy)
 
-end subroutine agrif_create_bathy_meter
+   END SUBROUTINE agrif_create_bathy_meter
 
-    SUBROUTINE init_bathy( ptab, i1, i2, j1, j2, before, nb,ndir)
-    use dom_oce
+   SUBROUTINE init_bathy( ptab, i1, i2, j1, j2, before, nb,ndir)
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE interpsshn  ***
@@ -20 +28 @@
       LOGICAL                         , INTENT(in   ) ::   before
       INTEGER                         , INTENT(in   ) ::   nb , ndir
-      LOGICAL  ::   western_side, eastern_side,northern_side,southern_side
       !
       !!----------------------------------------------------------------------  
-      !
-         western_side  = (nb == 1).AND.(ndir == 1)
-         eastern_side  = (nb == 1).AND.(ndir == 2)
-         southern_side = (nb == 2).AND.(ndir == 1)
-         northern_side = (nb == 2).AND.(ndir == 2)
+      INTEGER :: ji,jj
+
       IF( before) THEN
          ptab(i1:i2,j1:j2) = bathy(i1:i2,j1:j2)
+         DO jj=j1,j2
+            DO ji=i1,i2
+               ptab(ji,jj) = SUM( e3t_0(ji,jj, 1:mbkt(ji,jj) ) ) * ssmask(ji,jj)
+            END DO
+         END DO
       ELSE
          bathy(i1:i2,j1:j2)=ptab
@@ -35 +44 @@
       !
    END SUBROUTINE init_bathy
+
 #else
-subroutine agrif_domzgr_empty
-end subroutine agrif_domzgr_empty
+   SUBROUTINE agrif_create_bathy_meter
+   END SUBROUTINE agrif_create_bathy_meter
 #endif
+END MODULE agrif_domzgr