Ignore:
Timestamp:
2019-09-19T12:08:31+02:00 (12 months ago)
Author:
jchanut
Message:

#2222, import changes from dev_r10973_AGRIF-01_jchanut_small_jpi_jpj (i.e. #2199)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/OCE/DYN/dynspg_ts.F90

    r11573 r11574  
    483483         !                             ! values of zhup2_e and zhvp2_e on the halo are not needed in bdy_vol2d 
    484484         IF( ln_bdy .AND. ln_vol ) CALL bdy_vol2d( kt, jn, ua_e, va_e, zhup2_e, zhvp2_e ) 
    485          ! 
     485         !       
    486486         !                             ! resulting flux at mid-step (not over the full domain) 
    487487         zhU(1:jpim1,1:jpj  ) = e2u(1:jpim1,1:jpj  ) * ua_e(1:jpim1,1:jpj  ) * zhup2_e(1:jpim1,1:jpj  )   ! not jpi-column 
     
    490490#if defined key_agrif 
    491491         ! Set fluxes during predictor step to ensure volume conservation 
    492          IF( .NOT.Agrif_Root() .AND. ln_bt_fw ) THEN 
    493             IF((nbondi == -1).OR.(nbondi == 2)) THEN 
    494                DO jj = 1, jpj 
    495                   zhU(2:nbghostcells+1,jj) = ubdy_w(1:nbghostcells,jj) * e2u(2:nbghostcells+1,jj) 
    496                   zhV(2:nbghostcells+1,jj) = vbdy_w(1:nbghostcells,jj) * e1v(2:nbghostcells+1,jj) 
    497                END DO 
    498             ENDIF 
    499             IF((nbondi ==  1).OR.(nbondi == 2)) THEN 
    500                DO jj=1,jpj 
    501                   zhU(nlci-nbghostcells-1:nlci-2,jj) = ubdy_e(1:nbghostcells,jj) * e2u(nlci-nbghostcells-1:nlci-2,jj) 
    502                   zhV(nlci-nbghostcells  :nlci-1,jj) = vbdy_e(1:nbghostcells,jj) * e1v(nlci-nbghostcells  :nlci-1,jj) 
    503                END DO 
    504             ENDIF 
    505             IF((nbondj == -1).OR.(nbondj == 2)) THEN 
    506                DO ji=1,jpi 
    507                   zhV(ji,2:nbghostcells+1) = vbdy_s(ji,1:nbghostcells) * e1v(ji,2:nbghostcells+1) 
    508                   zhU(ji,2:nbghostcells+1) = ubdy_s(ji,1:nbghostcells) * e2u(ji,2:nbghostcells+1) 
    509                END DO 
    510             ENDIF 
    511             IF((nbondj ==  1).OR.(nbondj == 2)) THEN 
    512                DO ji=1,jpi 
    513                   zhV(ji,nlcj-nbghostcells-1:nlcj-2) = vbdy_n(ji,1:nbghostcells) * e1v(ji,nlcj-nbghostcells-1:nlcj-2) 
    514                   zhU(ji,nlcj-nbghostcells  :nlcj-1) = ubdy_n(ji,1:nbghostcells) * e2u(ji,nlcj-nbghostcells  :nlcj-1) 
    515                END DO 
    516             ENDIF 
    517          ENDIF 
     492         IF( .NOT.Agrif_Root() .AND. ln_bt_fw ) CALL agrif_dyn_ts_flux( jn, zhU, zhV ) 
    518493#endif 
    519494         IF( ln_wd_il )   CALL wad_lmt_bt(zhU, zhV, sshn_e, zssh_frc, rdtbt)    !!gm wad_lmt_bt use of lbc_lnk on zhU, zhV 
Note: See TracChangeset for help on using the changeset viewer.