Changeset 2548


Ignore:
Timestamp:
2011-01-11T17:22:41+01:00 (10 years ago)
Author:
rblod
Message:

Allow reproducibility of GYRE configuration

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90

    r2528 r2548  
    1919   USE lib_mpp         ! distribued memory computing library 
    2020   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     21   USE lib_fortran 
    2122 
    2223   IMPLICIT NONE 
     
    209210      ! Compute the emp flux such as its integration on the whole domain at each time is zero 
    210211      IF( nbench /= 1 ) THEN 
    211          zsumemp = 0.e0   ;   zsurf = 0.e0 
    212          DO jj = 1, jpj 
    213             DO ji = 1, jpi 
    214                zsumemp = zsumemp + emp(ji,jj) * tmask(ji,jj,1) * tmask_i(ji,jj) 
    215                zsurf   = zsurf   +              tmask(ji,jj,1) * tmask_i(ji,jj) 
    216             END DO 
    217          END DO 
    218  
    219          IF( lk_mpp )   CALL mpp_sum( zsumemp  )       ! sum over the global domain 
    220          IF( lk_mpp )   CALL mpp_sum( zsurf    )       ! sum over the global domain 
    221  
     212         zsumemp = GLOB_SUM( emp(:,:) )  
     213         zsurf   = GLOB_SUM( tmask(:,:,1) )  
    222214         ! Default GYRE configuration 
    223215         zsumemp = zsumemp / zsurf 
Note: See TracChangeset for help on using the changeset viewer.