Changeset 15564
- Timestamp:
- 2021-12-01T17:35:59+01:00 (3 years ago)
- Location:
- NEMO/releases/r4.0/r4.0-HEAD/src/NST
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/NST/agrif_oce_update.F90
r10068 r15564 404 404 !! 405 405 INTEGER :: ji,jj,jk,jn 406 REAL(wp) :: z tb, ztnu, ztno406 REAL(wp) :: ze3b, ztb, ztnu, ztno 407 407 !!--------------------------------------------- 408 408 ! … … 434 434 DO ji = i1, i2 435 435 IF( tabres(ji,jj,jk,jn) /= 0._wp ) THEN 436 ztb = tsb(ji,jj,jk,jn) * e3t_b(ji,jj,jk) ! fse3t_b prior update should be used 436 ze3b = e3t_b(ji,jj,jk) & ! Recover e3t_b before update 437 & - rn_atfp * ( e3t_n(ji,jj,jk) - e3t_a(ji,jj,jk) ) 438 ztb = tsb(ji,jj,jk,jn) * ze3b 437 439 ztnu = tabres(ji,jj,jk,jn) 438 440 ztno = tsn(ji,jj,jk,jn) * e3t_a(ji,jj,jk) … … 573 575 ! 574 576 INTEGER :: ji, jj, jk 575 REAL(wp) :: zrhoy, z ub, zunu, zuno577 REAL(wp) :: zrhoy, ze3b, zub, zunu, zuno 576 578 !!--------------------------------------------- 577 579 ! … … 588 590 ! 589 591 IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN ! Add asselin part 590 zub = ub(ji,jj,jk) * e3u_b(ji,jj,jk) ! fse3t_b prior update should be used 592 ze3b = e3u_b(ji,jj,jk) & ! Recover e3u_b before update 593 & - rn_atfp * ( e3u_n(ji,jj,jk) - e3u_a(ji,jj,jk) ) 594 zub = ub(ji,jj,jk) * ze3b 591 595 zuno = un(ji,jj,jk) * e3u_a(ji,jj,jk) 592 596 zunu = tabres(ji,jj,jk,1) … … 759 763 ! 760 764 INTEGER :: ji, jj, jk 761 REAL(wp) :: zrhox, z vb, zvnu, zvno765 REAL(wp) :: zrhox, ze3b, zvb, zvnu, zvno 762 766 !!--------------------------------------------- 763 767 ! … … 778 782 ! 779 783 IF (.NOT.(lk_agrif_fstep.AND.(neuler==0))) THEN ! Add asselin part 780 zvb = vb(ji,jj,jk) * e3v_b(ji,jj,jk) ! fse3t_b prior update should be used 784 ze3b = e3v_b(ji,jj,jk) & ! Recover e3v_b before update 785 & - rn_atfp * ( e3v_n(ji,jj,jk) - e3v_a(ji,jj,jk) ) 786 zvb = vb(ji,jj,jk) * ze3b 781 787 zvno = vn(ji,jj,jk) * e3v_a(ji,jj,jk) 782 788 zvnu = tabres(ji,jj,jk,1) -
NEMO/releases/r4.0/r4.0-HEAD/src/NST/agrif_top_update.F90
r11078 r15564 173 173 !! 174 174 INTEGER :: ji,jj,jk,jn 175 REAL(wp) :: ztb, z tnu, ztno175 REAL(wp) :: ztb, ze3b, ztnu, ztno 176 176 !!---------------------------------------------------------------------- 177 177 ! … … 204 204 DO ji=i1,i2 205 205 IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN 206 ztb = trb(ji,jj,jk,jn) * e3t_b(ji,jj,jk) ! fse3t_b prior update should be used 206 207 ze3b = e3t_b(ji,jj,jk) & ! Recover e3t_b before update 208 & - rn_atfp * ( e3t_n(ji,jj,jk) - e3t_a(ji,jj,jk) ) 209 ztb = trb(ji,jj,jk,jn) * ze3b 207 210 ztnu = tabres(ji,jj,jk,jn) 208 211 ztno = trn(ji,jj,jk,jn) * e3t_a(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.