- Timestamp:
- 2021-06-18T12:26:11+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DYN/dynspg_ts.F90
r15004 r15018 499 499 ! 500 500 ! ! resulting flux at mid-step (not over the full domain) 501 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 502 zhV(1:jpi ,1:jpjm1) = e1v(1:jpi ,1:jpjm1) * va_e(1:jpi ,1:jpjm1) * zhvp2_e(1:jpi ,1:jpjm1) ! not jpj-row 501 DO_2D( 1, 0, 1, 1 ) ! not jpi-column 502 zhU(ji,jj) = e2u(ji,jj) * ua_e(ji,jj) * zhup2_e(ji,jj) 503 END_2D 504 DO_2D( 1, 1, 1, 0 ) ! not jpj-row 505 zhV(ji,jj) = e1v(ji,jj) * va_e(ji,jj) * zhvp2_e(ji,jj) 506 END_2D 503 507 ! 504 508 #if defined key_agrif … … 538 542 ! sum over sub-time-steps to decide which baroclinic velocities to set to zero (zuwdav2 is only used when ln_wd_dl_bc=True) 539 543 IF ( ln_wd_dl_bc ) THEN 540 zuwdav2(1:jpim1,1:jpj ) = zuwdav2(1:jpim1,1:jpj ) + za2 * zuwdmask(1:jpim1,1:jpj ) ! not jpi-column 541 zvwdav2(1:jpi ,1:jpjm1) = zvwdav2(1:jpi ,1:jpjm1) + za2 * zvwdmask(1:jpi ,1:jpjm1) ! not jpj-row 544 DO_2D( 1, 0, 1, 1 ) ! not jpi-column 545 zuwdav2(ji,jj) = zuwdav2(ji,jj) + za2 * zuwdmask(ji,jj) 546 END_2D 547 DO_2D( 1, 1, 1, 0 ) ! not jpj-row 548 zvwdav2(ji,jj) = zvwdav2(ji,jj) + za2 * zvwdmask(ji,jj) 549 END_2D 542 550 END IF 543 551 ! … … 579 587 IF( ln_wd_il ) THEN ! W/D : gravity filters applied on pressure gradient 580 588 CALL wad_spg( zsshp2_e, zcpx, zcpy ) ! Calculating W/D gravity filters 581 zu_spg(2:jpim1,2:jpjm1) = zu_spg(2:jpim1,2:jpjm1) * zcpx(2:jpim1,2:jpjm1) 582 zv_spg(2:jpim1,2:jpjm1) = zv_spg(2:jpim1,2:jpjm1) * zcpy(2:jpim1,2:jpjm1) 589 DO_2D( 0, 0, 0, 0 ) 590 zu_spg(ji,jj) = zu_spg(ji,jj) * zcpx(ji,jj) 591 zv_spg(ji,jj) = zv_spg(ji,jj) * zcpy(ji,jj) 592 END_2D 583 593 ENDIF 584 594 ! … … 670 680 671 681 IF( .NOT.ln_linssh ) THEN !* Update ocean depth (variable volume case only) 672 hu_e (2:jpim1,2:jpjm1) = hu_0(2:jpim1,2:jpjm1) + zsshu_a(2:jpim1,2:jpjm1) 673 hur_e(2:jpim1,2:jpjm1) = ssumask(2:jpim1,2:jpjm1) / ( hu_e(2:jpim1,2:jpjm1) + 1._wp - ssumask(2:jpim1,2:jpjm1) ) 674 hv_e (2:jpim1,2:jpjm1) = hv_0(2:jpim1,2:jpjm1) + zsshv_a(2:jpim1,2:jpjm1) 675 hvr_e(2:jpim1,2:jpjm1) = ssvmask(2:jpim1,2:jpjm1) / ( hv_e(2:jpim1,2:jpjm1) + 1._wp - ssvmask(2:jpim1,2:jpjm1) ) 682 DO_2D( 0, 0, 0, 0 ) 683 hu_e (ji,jj) = hu_0(ji,jj) + zsshu_a(ji,jj) 684 hur_e(ji,jj) = ssumask(ji,jj) / ( hu_e(ji,jj) + 1._wp - ssumask(ji,jj) ) 685 hv_e (ji,jj) = hv_0(ji,jj) + zsshv_a(ji,jj) 686 hvr_e(ji,jj) = ssvmask(ji,jj) / ( hv_e(ji,jj) + 1._wp - ssvmask(ji,jj) ) 687 END_2D 676 688 ENDIF 677 689 !
Note: See TracChangeset
for help on using the changeset viewer.