Changeset 11205 for NEMO/branches/2019/dev_r10973_AGRIF-01_jchanut_small_jpi_jpj/src/OCE/DYN/dynspg_ts.F90
- Timestamp:
- 2019-07-02T14:25:46+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10973_AGRIF-01_jchanut_small_jpi_jpj/src/OCE/DYN/dynspg_ts.F90
r10742 r11205 797 797 IF( ln_bdy .AND. ln_vol ) CALL bdy_vol2d( kt, jn, ua_e, va_e, zhup2_e, zhvp2_e ) 798 798 ! 799 #if defined key_agrif 800 ! Set fluxes during predictor step to ensure volume conservation 801 IF( .NOT.Agrif_Root() .AND. ln_bt_fw ) CALL agrif_dyn_ts( jn ) 802 #endif 799 803 zwx(:,:) = e2u(:,:) * ua_e(:,:) * zhup2_e(:,:) ! fluxes at jn+0.5 800 804 zwy(:,:) = e1v(:,:) * va_e(:,:) * zhvp2_e(:,:) 801 805 ! 802 #if defined key_agrif803 ! Set fluxes during predictor step to ensure volume conservation804 IF( .NOT.Agrif_Root() .AND. ln_bt_fw ) THEN805 IF((nbondi == -1).OR.(nbondi == 2)) THEN806 DO jj = 1, jpj807 zwx(2:nbghostcells+1,jj) = ubdy_w(1:nbghostcells,jj) * e2u(2:nbghostcells+1,jj)808 zwy(2:nbghostcells+1,jj) = vbdy_w(1:nbghostcells,jj) * e1v(2:nbghostcells+1,jj)809 END DO810 ENDIF811 IF((nbondi == 1).OR.(nbondi == 2)) THEN812 DO jj=1,jpj813 zwx(nlci-nbghostcells-1:nlci-2,jj) = ubdy_e(1:nbghostcells,jj) * e2u(nlci-nbghostcells-1:nlci-2,jj)814 zwy(nlci-nbghostcells :nlci-1,jj) = vbdy_e(1:nbghostcells,jj) * e1v(nlci-nbghostcells :nlci-1,jj)815 END DO816 ENDIF817 IF((nbondj == -1).OR.(nbondj == 2)) THEN818 DO ji=1,jpi819 zwy(ji,2:nbghostcells+1) = vbdy_s(ji,1:nbghostcells) * e1v(ji,2:nbghostcells+1)820 zwx(ji,2:nbghostcells+1) = ubdy_s(ji,1:nbghostcells) * e2u(ji,2:nbghostcells+1)821 END DO822 ENDIF823 IF((nbondj == 1).OR.(nbondj == 2)) THEN824 DO ji=1,jpi825 zwy(ji,nlcj-nbghostcells-1:nlcj-2) = vbdy_n(ji,1:nbghostcells) * e1v(ji,nlcj-nbghostcells-1:nlcj-2)826 zwx(ji,nlcj-nbghostcells :nlcj-1) = ubdy_n(ji,1:nbghostcells) * e2u(ji,nlcj-nbghostcells :nlcj-1)827 END DO828 ENDIF829 ENDIF830 #endif831 806 IF( ln_wd_il ) CALL wad_lmt_bt(zwx, zwy, sshn_e, zssh_frc, rdtbt) 832 807
Note: See TracChangeset
for help on using the changeset viewer.