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 13157 – NEMO

Changeset 13157


Ignore:
Timestamp:
2020-06-25T17:31:26+02:00 (4 years ago)
Author:
jchanut
Message:

#2129, swap lbc_lnk and agrif bdy update in split-explicit fs. Overall strategy: agrif first, lateral bcs after

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12973_AGRIF_CMEMS/src/OCE/DYN/dynspg_ts.F90

    r12489 r13157  
    513513         END_2D 
    514514         ! 
     515         ! Duplicate sea level across open boundaries (this is only cosmetic if linssh=T) 
     516         IF( ln_bdy )   CALL bdy_ssh( ssha_e ) 
     517#if defined key_agrif 
     518         IF( .NOT.Agrif_Root() )   CALL agrif_ssh_ts( jn ) 
     519#endif 
     520 
    515521         CALL lbc_lnk_multi( 'dynspg_ts', ssha_e, 'T', 1._wp,  zhU, 'U', -1._wp,  zhV, 'V', -1._wp ) 
    516522         ! 
     
    525531         END IF 
    526532         ! 
    527          ! Duplicate sea level across open boundaries (this is only cosmetic if linssh=T) 
    528          IF( ln_bdy )   CALL bdy_ssh( ssha_e ) 
    529 #if defined key_agrif 
    530          IF( .NOT.Agrif_Root() )   CALL agrif_ssh_ts( jn ) 
    531 #endif 
    532533         !   
    533534         ! Sea Surface Height at u-,v-points (vvl case only) 
     
    643644         ENDIF 
    644645        
    645          IF( .NOT.ln_linssh ) THEN   !* Update ocean depth (variable volume case only) 
     646         IF( .NOT.ln_linssh ) THEN !* Update ocean depth (variable volume case only) 
    646647            hu_e (2:jpim1,2:jpjm1) = hu_0(2:jpim1,2:jpjm1) + zsshu_a(2:jpim1,2:jpjm1) 
    647648            hur_e(2:jpim1,2:jpjm1) = ssumask(2:jpim1,2:jpjm1) / ( hu_e(2:jpim1,2:jpjm1) + 1._wp - ssumask(2:jpim1,2:jpjm1) ) 
    648649            hv_e (2:jpim1,2:jpjm1) = hv_0(2:jpim1,2:jpjm1) + zsshv_a(2:jpim1,2:jpjm1) 
    649650            hvr_e(2:jpim1,2:jpjm1) = ssvmask(2:jpim1,2:jpjm1) / ( hv_e(2:jpim1,2:jpjm1) + 1._wp - ssvmask(2:jpim1,2:jpjm1) ) 
     651         ENDIF 
     652         !                                                 ! open boundaries 
     653         IF( ln_bdy )   CALL bdy_dyn2d( jn, ua_e, va_e, un_e, vn_e, hur_e, hvr_e, ssha_e ) 
     654#if defined key_agrif                                                            
     655         IF( .NOT.Agrif_Root() )  CALL agrif_dyn_ts( jn )  ! Agrif 
     656#endif 
     657         ! 
     658         IF( .NOT.ln_linssh ) THEN   !* Update ocean depth (variable volume case only) 
    650659            CALL lbc_lnk_multi( 'dynspg_ts', ua_e , 'U', -1._wp, va_e , 'V', -1._wp  & 
    651660                 &                         , hu_e , 'U',  1._wp, hv_e , 'V',  1._wp  & 
     
    655664         ENDIF 
    656665         ! 
    657          ! 
    658          !                                                 ! open boundaries 
    659          IF( ln_bdy )   CALL bdy_dyn2d( jn, ua_e, va_e, un_e, vn_e, hur_e, hvr_e, ssha_e ) 
    660 #if defined key_agrif                                                            
    661          IF( .NOT.Agrif_Root() )  CALL agrif_dyn_ts( jn )  ! Agrif 
    662 #endif 
    663666         !                                             !* Swap 
    664667         !                                             !  ---- 
Note: See TracChangeset for help on using the changeset viewer.