Changeset 10297 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/domvvl.F90
- Timestamp:
- 2018-11-12T16:20:57+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/domvvl.F90
r10170 r10297 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 ( dom_vvl_alloc )81 IF( lk_mpp ) 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 ( dom_vvl_alloc )88 IF( lk_mpp ) 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 … … 429 429 END DO 430 430 z_tmax = MAXVAL( ze3t(:,:,:) ) 431 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain431 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 432 432 z_tmin = MINVAL( ze3t(:,:,:) ) 433 IF( lk_mpp ) CALL mpp_min( z_tmin ) ! min over the global domain433 IF( lk_mpp ) CALL mpp_min( 'domvvl', z_tmin ) ! min over the global domain 434 434 ! - ML - test: for the moment, stop simulation for too large e3_t variations 435 435 IF( ( z_tmax > rn_zdef_max ) .OR. ( z_tmin < - rn_zdef_max ) ) THEN … … 493 493 IF ( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 494 494 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( zht(:,:) ) ) 495 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain495 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 496 496 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(SUM(tilde_e3t_a))) =', z_tmax 497 497 END IF … … 502 502 END DO 503 503 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( ht_0(:,:) + sshn(:,:) - zht(:,:) ) ) 504 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain504 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 505 505 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ht_0+sshn-SUM(e3t_n))) =', z_tmax 506 506 ! … … 510 510 END DO 511 511 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( ht_0(:,:) + ssha(:,:) - zht(:,:) ) ) 512 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain512 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 513 513 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ht_0+ssha-SUM(e3t_a))) =', z_tmax 514 514 ! … … 518 518 END DO 519 519 z_tmax = MAXVAL( tmask(:,:,1) * tmask_i(:,:) * ABS( ht_0(:,:) + sshb(:,:) - zht(:,:) ) ) 520 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain520 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 521 521 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ht_0+sshb-SUM(e3t_b))) =', z_tmax 522 522 ! 523 523 z_tmax = MAXVAL( tmask(:,:,1) * ABS( sshb(:,:) ) ) 524 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain524 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 525 525 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(sshb))) =', z_tmax 526 526 ! 527 527 z_tmax = MAXVAL( tmask(:,:,1) * ABS( sshn(:,:) ) ) 528 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain528 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 529 529 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(sshn))) =', z_tmax 530 530 ! 531 531 z_tmax = MAXVAL( tmask(:,:,1) * ABS( ssha(:,:) ) ) 532 IF( lk_mpp ) CALL mpp_max( z_tmax ) ! max over the global domain532 IF( lk_mpp ) CALL mpp_max( 'domvvl', z_tmax ) ! max over the global domain 533 533 IF( lwp ) WRITE(numout, *) kt,' MAXVAL(abs(ssha))) =', z_tmax 534 534 END IF
Note: See TracChangeset
for help on using the changeset viewer.