Changeset 14952 for utils/tools/DOMAINcfg/src/agrif_recompute_scales.F90

Timestamp:

2021-06-04T17:42:50+02:00 (3 years ago)

Author:

jchanut

Message:

#2638, closed domains AGRIF new convention + various add ons

File:

• utils/tools/DOMAINcfg/src/agrif_recompute_scales.F90 (modified) (1 diff)

utils/tools/DOMAINcfg/src/agrif_recompute_scales.F90

@@ -39 +39 @@
          END DO
       END DO
-      IF ( ln_isfcav ) THEN
-      ! (ISF) define e3uw (adapted for 2 cells in the water column)
-      print *,'NOT READY SINCE:'
-      print *,'MBATHY HAS NOT BEEN CORRECTED / UPDATED'
-      print *,'EVEN NOT COMPUTED IN THE CASE ln_read_cfg = .TRUE.'
-      STOP
-         DO jj = 2, jpjm1 
-            DO ji = 2, jpim1   ! vector opt. 
-               ikb = MAX(mbathy (ji,jj),mbathy (ji+1,jj))
-               ikt = MAX(misfdep(ji,jj),misfdep(ji+1,jj))
-               IF (ikb == ikt+1) e3uw_0(ji,jj,ikb) =  MIN( gdept_0(ji,jj,ikb  ), gdept_0(ji+1,jj  ,ikb  ) ) &
-                                       &            - MAX( gdept_0(ji,jj,ikb-1), gdept_0(ji+1,jj  ,ikb-1) )
-               ikb = MAX(mbathy (ji,jj),mbathy (ji,jj+1))
-               ikt = MAX(misfdep(ji,jj),misfdep(ji,jj+1))
-               IF (ikb == ikt+1) e3vw_0(ji,jj,ikb) =  MIN( gdept_0(ji,jj,ikb  ), gdept_0(ji  ,jj+1,ikb  ) ) &
-                                       &            - MAX( gdept_0(ji,jj,ikb-1), gdept_0(ji  ,jj+1,ikb-1) )
-            END DO
-         END DO
-      END IF
+!      IF ( ln_isfcav ) THEN
+!      ! (ISF) define e3uw (adapted for 2 cells in the water column)
+!      print *,'NOT READY SINCE:'
+!      print *,'MBATHY HAS NOT BEEN CORRECTED / UPDATED'
+!      print *,'EVEN NOT COMPUTED IN THE CASE ln_read_cfg = .TRUE.'
+!      STOP
+!         DO jj = 2, jpjm1 
+!            DO ji = 2, jpim1   ! vector opt. 
+!               ikb = MAX(mbathy (ji,jj),mbathy (ji+1,jj))
+!               ikt = MAX(misfdep(ji,jj),misfdep(ji+1,jj))
+!               IF (ikb == ikt+1) e3uw_0(ji,jj,ikb) =  MIN( gdept_0(ji,jj,ikb  ), gdept_0(ji+1,jj  ,ikb  ) ) &
+!                                       &            - MAX( gdept_0(ji,jj,ikb-1), gdept_0(ji+1,jj  ,ikb-1) )
+!               ikb = MAX(mbathy (ji,jj),mbathy (ji,jj+1))
+!               ikt = MAX(misfdep(ji,jj),misfdep(ji,jj+1))
+!               IF (ikb == ikt+1) e3vw_0(ji,jj,ikb) =  MIN( gdept_0(ji,jj,ikb  ), gdept_0(ji  ,jj+1,ikb  ) ) &
+!                                       &            - MAX( gdept_0(ji,jj,ikb-1), gdept_0(ji  ,jj+1,ikb-1) )
+!            END DO
+!         END DO
+!      END IF
 
       CALL lbc_lnk('toto', e3u_0 , 'U', 1._wp )   ;   CALL lbc_lnk('toto', e3uw_0, 'U', 1._wp )   ! lateral boundary conditions