Changeset 10359 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/ICE/icedyn_adv_umx.F90

Timestamp:

2018-11-25T22:33:50+01:00 (5 years ago)

Author:

smasson

Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 4c: introduce a better version of the non-blocking mpp_max, see #2133

File:

• NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/ICE/icedyn_adv_umx.F90 (modified) (2 diffs)

NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/ICE/icedyn_adv_umx.F90

@@ -73 +73 @@
       INTEGER  ::   ipl                     ! third dimention of tracer array
 
-      REAL(wp) ::   zcfl , zusnit, zdt      !   -      -
+      REAL(wp) ::   zcfl(2), zusnit, zdt      !   -      -
       REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zudy, zvdx, zcu_box, zcv_box
       REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   zpato
@@ -86 +86 @@
       !     When needed, the advection split is applied at the next time-step in order to avoid blocking global comm.
       !     ...this should not affect too much the stability... Was ok on the tests we did...
-      zcfl =            MAXVAL( ABS( pu_ice(:,:) ) * rdt_ice * r1_e1u(:,:) )
-      zcfl = MAX( zcfl, MAXVAL( ABS( pv_ice(:,:) ) * rdt_ice * r1_e2v(:,:) ) )
-      IF( zcfl > 0.5 ) THEN   ;   l_split_advumx(1) = .TRUE.    ! split advection time-step if CFL violated
-      ELSE                    ;   l_split_advumx(1) = .FALSE.
-      ENDIF
+      zcfl(1) =               MAXVAL( ABS( pu_ice(:,:) ) * rdt_ice * r1_e1u(:,:) )
+      zcfl(1) = MAX( zcfl(1), MAXVAL( ABS( pv_ice(:,:) ) * rdt_ice * r1_e2v(:,:) ) )
       
-      ! non-blocking global communication send l_split_advumx(1) and receive l_split_advumx(2)
-      IF( lk_mpp )   CALL mpp_ilor( l_split_advumx, ldlast = kt == nitend - nn_fsbc + 1 )
-
-      IF( l_split_advumx(2) ) THEN   ;   initad = 2   ;   zusnit = 0.5_wp    ! split defined at the previous time-step
-      ELSE                           ;   initad = 1   ;   zusnit = 1.0_wp
+      ! non-blocking global communication send zcfl(1) and receive zcfl(2)
+      IF( lk_mpp )   CALL mpp_delay_max( 'icedyn_adv_umx', 'advumx_delay', zcfl, kt == nitend - nn_fsbc + 1 )
+
+      IF( zcfl(2) > .5 ) THEN   ;   initad = 2   ;   zusnit = 0.5_wp   ! zcfl(2) corresponding to zcfl(1) of the previous time-step
+      ELSE                      ;   initad = 1   ;   zusnit = 1.0_wp
       ENDIF