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 10288 for NEMO/branches/2018/dev_r9866_HPC_03_globcom/src/OCE/SBC/sbcfwb.F90 – NEMO

Ignore:
Timestamp:
2018-11-07T18:25:49+01:00 (5 years ago)
Author:
francesca
Message:

reduce global communications, see #2010

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9866_HPC_03_globcom/src/OCE/SBC/sbcfwb.F90

    r9598 r10288  
    4444   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    4545   !! $Id$ 
    46    !! Software governed by the CeCILL licence     (./LICENSE) 
     46   !! Software governed by the CeCILL license (see ./LICENSE) 
    4747   !!---------------------------------------------------------------------- 
    4848CONTAINS 
     
    155155            ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 
    156156            ! 
    157             zsurf_neg = glob_sum( e1e2t(:,:)*ztmsk_neg(:,:) )  ! Area filled by <0 and >0 erp  
    158             zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) ) 
    159157            !                                                  ! fwf global mean (excluding ocean to ice/snow exchanges)  
    160158            z_fwf     = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 
    161159            !             
    162160            IF( z_fwf < 0._wp ) THEN         ! spread out over >0 erp area to increase evaporation 
     161                zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) ) 
    163162                zsurf_tospread      = zsurf_pos 
    164163                ztmsk_tospread(:,:) = ztmsk_pos(:,:) 
    165164            ELSE                             ! spread out over <0 erp area to increase precipitation 
     165                zsurf_neg = glob_sum( e1e2t(:,:)*ztmsk_neg(:,:) )  ! Area filled by <0 and >0 erp  
    166166                zsurf_tospread      = zsurf_neg 
    167167                ztmsk_tospread(:,:) = ztmsk_neg(:,:) 
Note: See TracChangeset for help on using the changeset viewer.