Changeset 11574 for NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/OCE/DYN/dynspg_ts.F90
- Timestamp:
- 2019-09-19T12:08:31+02:00 (5 years ago)
- 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 483 483 ! ! values of zhup2_e and zhvp2_e on the halo are not needed in bdy_vol2d 484 484 IF( ln_bdy .AND. ln_vol ) CALL bdy_vol2d( kt, jn, ua_e, va_e, zhup2_e, zhvp2_e ) 485 ! 485 ! 486 486 ! ! resulting flux at mid-step (not over the full domain) 487 487 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 … … 490 490 #if defined key_agrif 491 491 ! 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 ) 518 493 #endif 519 494 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.