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 – NEMO

Changeset 10372


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

Location:
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP
Files:
3 edited

Legend:

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

    r10170 r10372  
    113113 
    114114      ! 0. initialization 
    115       DO jn = 1, jptra 
    116          ! 
    117          IF( l_trdtrc )   ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends 
    118          ! 
    119          IF( nn_ice_tr == -1 ) THEN    ! No tracers in sea ice (null concentration in sea ice) 
     115      IF( nn_ice_tr == -1 ) THEN    ! No tracers in sea ice (null concentration in sea ice) 
     116         ! 
     117         DO jn = 1, jptra 
    120118            DO jj = 2, jpj 
    121119               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    123121               END DO 
    124122            END DO 
    125          ELSE 
     123         END DO 
     124         ! 
     125       ELSE 
     126         ! 
     127         DO jn = 1, jptra 
    126128            DO jj = 2, jpj 
    127129               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    143145               END DO 
    144146            END DO 
    145          ENDIF 
    146          ! 
    147          CALL lbc_lnk( 'trcsbc', sbc_trc(:,:,jn), 'T', 1. ) 
    148          !                                       Concentration dilution effect on tracers due to evaporation & precipitation  
     147         END DO 
     148      ENDIF 
     149      ! 
     150      CALL lbc_lnk( 'trcsbc', sbc_trc(:,:,:), 'T', 1. ) 
     151      !                                       Concentration dilution effect on tracers due to evaporation & precipitation  
     152      DO jn = 1, jptra 
     153         ! 
     154         IF( l_trdtrc )   ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends 
     155         ! 
    149156         DO jj = 2, jpj 
    150157            DO ji = fs_2, fs_jpim1   ! vector opt. 
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/trcini.F90

    r10345 r10372  
    7272      CALL trc_ice_ini   ! Tracers in sea ice 
    7373      ! 
    74       IF(lwp) CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea ) 
     74      IF(lwm) CALL ctl_opn( numstr, 'tracer.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp , narea ) 
    7575      ! 
    7676      CALL trc_ini_state  !  passive tracers initialisation : from a restart or from clim 
  • 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.