Changeset 13109 for utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_dom_update.F90

Timestamp:

2020-06-16T12:07:28+02:00 (4 years ago)

Author:

rblod

Message:

ticket #2129 : major corrections in domcfg

File:

• utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_dom_update.F90 (modified) (6 diffs)

utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_dom_update.F90

@@ -22 +22 @@
       IF( Agrif_Root() ) return
 
-      CALL agrif_update_variable(bottom_level_id,locupdate=(/npt_copy,0/),procname = update_bottom_level)
+      CALL agrif_update_variable(bottom_level_id,procname = update_bottom_level)
       !
       Agrif_UseSpecialValueInUpdate = .TRUE.
@@ -41 +41 @@
       !
       !!---------------------------------------------------------------------- 
+      REAL(WP),DIMENSION(jpi,jpj) :: zk
       !
       IF( before) THEN
@@ -49 +50 @@
          WHERE ( mbkt(i1:i2,j1:j2) .EQ. 0 )
             ssmask(i1:i2,j1:j2) = 0.
+            mbkt(i1:i2,j1:j2)   = 1
          ELSEWHERE
             ssmask(i1:i2,j1:j2) = 1.
-         END WHERE           
+         END WHERE 
+         zk(:,:) = REAL(mbkt(:,:),wp); CALL lbc_lnk('update_bottom',zk,'T',1.); mbkt(:,:) = MAX(NINT(zk(:,:)),1)
+         CALL lbc_lnk('update_bottom',ssmask,'T',1.)          
       ENDIF
       !
@@ -71 +75 @@
             DO jj=j1,j2
                DO ji=i1,i2
-                   IF( mbkt(ji,jj) .LE. jk ) THEN
+                   IF( mbkt(ji,jj) .GE. jk ) THEN
                       tabres(ji,jj,jk) = e3t_0(ji,jj,jk)
                    ELSE
@@ -83 +87 @@
             DO jj=j1,j2
                DO ji=i1,i2
-                   IF( mbkt(ji,jj) .LE.jk ) THEN
+                   IF( mbkt(ji,jj) .GE. jk ) THEN
                       e3t_0(ji,jj,jk) = MAX(tabres(ji,jj,jk),MIN(e3zps_min,e3t_1d(jk)*e3zps_rat))
                    ELSE
@@ -91 +95 @@
             END DO
          END DO
+
+         CALL lbc_lnk('update_e3t',e3t_0,'T',1.)
          !
       ENDIF