Changeset 11241


Ignore:
Timestamp:
2019-07-10T11:46:48+02:00 (13 months ago)
Author:
girrmann
Message:

dev_r10984_HPC-13 : regroup communications in dyn_spg_ts, does not change the result, see #2285

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DYN/dynspg_ts.F90

    r11240 r11241  
    557557         ! Sea Surface Height at u-,v-points (vvl case only) 
    558558         IF( .NOT.ln_linssh ) THEN                                 
    559             DO jj = 2, jpjm1 
     559            DO jj = 2, jpjm1   ! INNER domain, will be extended to whole domain later 
    560560               DO ji = 2, jpim1      ! NO Vector Opt. 
    561561                  zsshu_a(ji,jj) = r1_2 * ssumask(ji,jj) * r1_e1e2u(ji,jj)    & 
     
    567567               END DO 
    568568            END DO 
    569             CALL lbc_lnk_multi( 'dynspg_ts', zsshu_a, 'U', 1._wp, zsshv_a, 'V', 1._wp ) 
    570569         ENDIF    
    571570         !          
     
    681680         ENDIF 
    682681 
    683           
    684          IF( .NOT.ln_linssh ) THEN                     !* Update ocean depth (variable volume case only) 
    685             hu_e (:,:) = hu_0(:,:) + zsshu_a(:,:) 
     682        
     683         IF( .NOT. ln_linssh ) THEN 
     684            !                                             !* domain lateral boundary 
     685            CALL lbc_lnk_multi( 'dynspg_ts', ua_e   , 'U', -1._wp, va_e   , 'V', -1._wp & 
     686                 &                         , zsshu_a, 'U',  1._wp, zsshv_a, 'V',  1._wp ) 
     687            ! 
     688            hu_e (:,:) = hu_0(:,:) + zsshu_a(:,:)         !* Update ocean depth (variable volume case only) 
    686689            hv_e (:,:) = hv_0(:,:) + zsshv_a(:,:) 
    687690            hur_e(:,:) = ssumask(:,:) / ( hu_e(:,:) + 1._wp - ssumask(:,:) ) 
    688691            hvr_e(:,:) = ssvmask(:,:) / ( hv_e(:,:) + 1._wp - ssvmask(:,:) ) 
    689692            ! 
    690          ENDIF 
    691          !                                             !* domain lateral boundary 
    692          CALL lbc_lnk_multi( 'dynspg_ts', ua_e, 'U', -1._wp, va_e , 'V', -1._wp ) 
     693         ELSE 
     694            !                                             !* domain lateral boundary 
     695            CALL lbc_lnk_multi( 'dynspg_ts', ua_e   , 'U', -1._wp, va_e   , 'V', -1._wp ) 
     696            ! 
     697         ENDIF 
    693698         ! 
    694699         !                                                 ! open boundaries 
Note: See TracChangeset for help on using the changeset viewer.