Changeset 7427


Ignore:
Timestamp:
2016-12-01T18:51:59+01:00 (5 years ago)
Author:
timgraham
Message:

#1811 Another lk_mpp_rep

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_closea.F90

    r7426 r7427  
    240240      !                                                   !  update emp        ! 
    241241      zfwf = 0.e0_wp                                      !--------------------! 
    242       IF( lk_mpp_rep ) THEN                         ! MPP reproductible calculation 
    243          DO jc = 1, jpncs 
    244             ctmp = CMPLX( 0.e0, 0.e0, wp ) 
    245             DO jj = ncsj1(jc), ncsj2(jc) 
    246                DO ji = ncsi1(jc), ncsi2(jc) 
    247                   ztmp = e1e2t(ji,jj) * ( emp(ji,jj)-rnf(ji,jj) ) * tmask_i(ji,jj) 
    248                   CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
    249                END DO   
    250             END DO  
    251             IF( lk_mpp )   CALL mpp_sum( ctmp ) 
    252             zfwf(jc) = REAL(ctmp,wp) 
    253          END DO 
    254       ELSE                                          ! Standard calculation            
    255          DO jc = 1, jpncs 
    256             DO jj = ncsj1(jc), ncsj2(jc) 
    257                DO ji = ncsi1(jc), ncsi2(jc) 
    258                   zfwf(jc) = zfwf(jc) + e1e2t(ji,jj) * ( emp(ji,jj)-rnf(ji,jj) ) * tmask_i(ji,jj)  
    259                END DO   
    260             END DO  
    261          END DO 
    262          IF( lk_mpp )   CALL mpp_sum ( zfwf(:) , jpncs )       ! mpp: sum over all the global domain 
    263       ENDIF 
     242      DO jc = 1, jpncs 
     243         ctmp = CMPLX( 0.e0, 0.e0, wp ) 
     244         DO jj = ncsj1(jc), ncsj2(jc) 
     245            DO ji = ncsi1(jc), ncsi2(jc) 
     246               ztmp = e1e2t(ji,jj) * ( emp(ji,jj)-rnf(ji,jj) ) * tmask_i(ji,jj) 
     247               CALL DDPDD( CMPLX( ztmp, 0.e0, wp ), ctmp ) 
     248            END DO   
     249         END DO  
     250         IF( lk_mpp )   CALL mpp_sum( ctmp ) 
     251         zfwf(jc) = REAL(ctmp,wp) 
     252      END DO 
    264253 
    265254      IF( cd_cfg == "orca" .AND. kcfg == 2 ) THEN      ! Black Sea case for ORCA_R2 configuration 
Note: See TracChangeset for help on using the changeset viewer.