New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10372 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/trcstp.F90 – NEMO

Ignore:
Timestamp:
2018-12-04T16:19:46+01:00 (5 years ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 8: remove communications in TOP, see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/trcstp.F90

    r10314 r10372  
    2020   USE trdtrc_oce 
    2121   USE trdmxl_trc 
     22   USE sms_pisces,  ONLY : ln_check_mass 
    2223   ! 
    2324   USE prtctl_trc     ! Print control for debbuging 
     
    7374            cvol(:,:,jk) = e1e2t(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk) 
    7475         END DO 
    75          areatot         = glob_sum( 'trcstp', cvol(:,:,:) ) 
     76         IF ( ln_ctl .OR. kt == nitrst .OR. ( ln_check_mass .AND. kt == nitend )              & 
     77            & .OR. iom_use( "pno3tot" ) .OR. iom_use( "ppo4tot" ) .OR. iom_use( "psiltot" )   & 
     78            & .OR. iom_use( "palktot" ) .OR. iom_use( "pfertot" ) )                           & 
     79            &     areatot = glob_sum( 'trcstp', cvol(:,:,:) ) 
    7680      ENDIF 
    7781      ! 
     
    105109      ENDIF 
    106110      ! 
    107       ztrai = 0._wp                                                   !  content of all tracers 
    108       DO jn = 1, jptra 
    109          ztrai = ztrai + glob_sum( 'trcstp', trn(:,:,:,jn) * cvol(:,:,:)   ) 
    110       END DO 
    111       IF( lwp ) WRITE(numstr,9300) kt,  ztrai / areatot 
     111      IF (ln_ctl ) THEN 
     112         ztrai = 0._wp                                                   !  content of all tracers 
     113         DO jn = 1, jptra 
     114            ztrai = ztrai + glob_sum( 'trcstp', trn(:,:,:,jn) * cvol(:,:,:)   ) 
     115         END DO 
     116         IF( lwm ) WRITE(numstr,9300) kt,  ztrai / areatot 
     117      ENDIF 
    1121189300  FORMAT(i10,D23.16) 
    113119      ! 
Note: See TracChangeset for help on using the changeset viewer.