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 2219 for branches/DEV_r2191_3partymerge2010/NEMO/OPA_SRC/DYN – NEMO

Ignore:
Timestamp:
2010-10-12T15:06:30+02:00 (14 years ago)
Author:
rblod
Message:

Merge branch DEV_1879_mpp_rep into DEV_r2191_3partymerge2010

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r2191_3partymerge2010/NEMO/OPA_SRC/DYN/dynspg_flt.F90

    r2208 r2219  
    4646   USE iom 
    4747   USE restart         ! only for lrst_oce 
     48   USE lib_fortran 
    4849#if defined key_agrif 
    4950   USE agrif_opa_interp 
     
    188189 
    189190#if defined key_obc 
    190       IF( lk_obc )   CALL obc_dyn( kt )   ! Update velocities on each open boundary with the radiation algorithm 
    191       IF( lk_obc )   CALL obc_vol( kt )   ! Correction of the barotropic componant velocity to control the volume of the system 
     191      CALL obc_dyn( kt )      ! Update velocities on each open boundary with the radiation algorithm 
     192      CALL obc_vol( kt )      ! Correction of the barotropic componant velocity to control the volume of the system 
    192193#endif 
    193194#if defined key_bdy 
     
    283284      ! ------------------ 
    284285      rnorme =0.e0 
    285       rnorme = SUM( gcb(1:jpi,1:jpj) * gcdmat(1:jpi,1:jpj) * gcb(1:jpi,1:jpj) * bmask(:,:) ) 
    286       IF( lk_mpp )   CALL mpp_sum( rnorme )   ! sum over the global domain 
     286      rnorme = GLOB_SUM( gcb(1:jpi,1:jpj) * gcdmat(1:jpi,1:jpj) * gcb(1:jpi,1:jpj) * bmask(:,:) ) 
    287287 
    288288      epsr = eps * eps * rnorme 
     
    317317#if defined key_obc 
    318318            ! caution : grad D = 0 along open boundaries 
    319             IF( Agrif_Root() ) THEN 
    320                spgu(ji,jj) = z2dt * ztdgu * obcumask(ji,jj) 
    321                spgv(ji,jj) = z2dt * ztdgv * obcvmask(ji,jj) 
    322             ELSE 
    323                spgu(ji,jj) = z2dt * ztdgu 
    324                spgv(ji,jj) = z2dt * ztdgv 
    325             ENDIF 
     319            spgu(ji,jj) = z2dt * ztdgu * obcumask(ji,jj) 
     320            spgv(ji,jj) = z2dt * ztdgv * obcvmask(ji,jj) 
    326321#elif defined key_bdy 
    327322            ! caution : grad D = 0 along open boundaries 
Note: See TracChangeset for help on using the changeset viewer.