- Timestamp:
- 2018-12-17T17:07:54+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/DOM/domvvl.F90
r10345 r10402 79 79 & dtilde_e3t_a(jpi,jpj,jpk) , un_td (jpi,jpj,jpk) , vn_td (jpi,jpj,jpk) , & 80 80 & STAT = dom_vvl_alloc ) 81 IF( lk_mpp )CALL mpp_sum ( 'domvvl', dom_vvl_alloc )81 CALL mpp_sum ( 'domvvl', dom_vvl_alloc ) 82 82 IF( dom_vvl_alloc /= 0 ) CALL ctl_warn('dom_vvl_alloc: failed to allocate arrays') 83 83 un_td = 0._wp … … 86 86 IF( ln_vvl_ztilde ) THEN 87 87 ALLOCATE( frq_rst_e3t(jpi,jpj) , frq_rst_hdv(jpi,jpj) , hdiv_lf(jpi,jpj,jpk) , STAT= dom_vvl_alloc ) 88 IF( lk_mpp )CALL mpp_sum ( 'domvvl', dom_vvl_alloc )88 CALL mpp_sum ( 'domvvl', dom_vvl_alloc ) 89 89 IF( dom_vvl_alloc /= 0 ) CALL ctl_warn('dom_vvl_alloc: failed to allocate arrays') 90 90 ENDIF … … 431 431 END DO 432 432 z_tmax = MAXVAL( ze3t(:,:,:) ) 433 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain433 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 434 434 z_tmin = MINVAL( ze3t(:,:,:) ) 435 IF( lk_mpp )CALL mpp_min( 'domvvl', z_tmin ) ! min over the global domain435 CALL mpp_min( 'domvvl', z_tmin ) ! min over the global domain 436 436 ! - ML - test: for the moment, stop simulation for too large e3_t variations 437 437 IF( ( z_tmax > rn_zdef_max ) .OR. ( z_tmin < - rn_zdef_max ) ) THEN … … 495 495 IF ( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 496 496 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( zht(:,:) ) ) 497 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain497 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 498 498 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(SUM(tilde_e3t_a))) =', z_tmax 499 499 END IF … … 504 504 END DO 505 505 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( ht_0(:,:) + sshn(:,:) - zht(:,:) ) ) 506 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain506 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 507 507 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ht_0+sshn-SUM(e3t_n))) =', z_tmax 508 508 ! … … 512 512 END DO 513 513 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( ht_0(:,:) + ssha(:,:) - zht(:,:) ) ) 514 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain514 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 515 515 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ht_0+ssha-SUM(e3t_a))) =', z_tmax 516 516 ! … … 520 520 END DO 521 521 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( ht_0(:,:) + sshb(:,:) - zht(:,:) ) ) 522 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain522 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 523 523 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ht_0+sshb-SUM(e3t_b))) =', z_tmax 524 524 ! 525 525 z_tmax = MAXVAL( tmask(:,:,1) * ABS( sshb(:,:) ) ) 526 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain526 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 527 527 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(sshb))) =', z_tmax 528 528 ! 529 529 z_tmax = MAXVAL( tmask(:,:,1) * ABS( sshn(:,:) ) ) 530 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain530 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 531 531 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(sshn))) =', z_tmax 532 532 ! 533 533 z_tmax = MAXVAL( tmask(:,:,1) * ABS( ssha(:,:) ) ) 534 IF( lk_mpp )CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain534 CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 535 535 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax 536 536 END IF
Note: See TracChangeset
for help on using the changeset viewer.