Changeset 13295 for NEMO/trunk/src/ABL
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ABL/ablmod.F90
r13226 r13295 127 127 !! needed for surface boundary condition of TKE 128 128 !! pwndm contains | U10m - U_oce | (see blk_oce_1 in sbcblk) 129 DO_2D _11_11129 DO_2D( 1, 1, 1, 1 ) 130 130 zzoce = pCd_du (ji,jj) * pwndm (ji,jj) 131 131 #if defined key_si3 … … 245 245 ! 246 246 ! Advance u_abl & v_abl to time n+1 247 DO_2D _11_11247 DO_2D( 1, 1, 1, 1 ) 248 248 zcff = ( fft_abl(ji,jj) * rDt_abl )*( fft_abl(ji,jj) * rDt_abl ) ! (f dt)**2 249 249 … … 483 483 DO jk = 2, jpka ! outer loop 484 484 !------------- 485 DO_2D _01_01485 DO_2D( 0, 1, 0, 1 ) 486 486 zcff1 = pblh( ji, jj ) 487 487 zsig = ght_abl(jk) / MAX( jp_pblh_min, MIN( jp_pblh_max, zcff1 ) ) … … 506 506 DO jk = 2, jpka ! outer loop 507 507 !------------- 508 DO_2D _11_11508 DO_2D( 1, 1, 1, 1 ) 509 509 zcff1 = pblh( ji, jj ) 510 510 zsig = ght_abl(jk) / MAX( jp_pblh_min, MIN( jp_pblh_max, zcff1 ) ) … … 580 580 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 581 581 ! 582 DO_2D _11_11582 DO_2D( 1, 1, 1, 1 ) 583 583 ztemp = tq_abl( ji, jj, 2, nt_a, jp_ta ) 584 584 zhumi = tq_abl( ji, jj, 2, nt_a, jp_qa ) … … 589 589 END_2D 590 590 591 DO_2D _01_01591 DO_2D( 0, 1, 0, 1 ) 592 592 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * ( pssu(ji ,jj) + pssu(ji-1,jj) ) 593 593 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * ( pssv(ji,jj ) + pssv(ji,jj-1) ) … … 597 597 ! 598 598 ! ... scalar wind ( = | U10m - U_oce | ) at T-point (masked) 599 DO_2D _11_11599 DO_2D( 1, 1, 1, 1 ) 600 600 zcff = SQRT( zwnd_i(ji,jj) * zwnd_i(ji,jj) & 601 601 & + zwnd_j(ji,jj) * zwnd_j(ji,jj) ) ! * msk_abl(ji,jj) … … 610 610 ! Note the use of 0.5*(2-umask) in order to unmask the stress along coastlines 611 611 ! Note the use of MAX(tmask(i,j),tmask(i+1,j) is to mask tau over ice shelves 612 DO_2D _00_00612 DO_2D( 0, 0, 0, 0 ) 613 613 zcff = 0.5_wp * ( 2._wp - msk_abl(ji,jj)*msk_abl(ji+1,jj) ) 614 614 zztmp = MAX(msk_abl(ji,jj),msk_abl(ji+1,jj)) … … 633 633 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! 634 634 ! ------------------------------------------------------------ ! 635 DO_2D _00_00635 DO_2D( 0, 0, 0, 0 ) 636 636 ptaui_ice(ji,jj) = 0.5_wp * ( rhoa(ji+1,jj) * pCd_du_ice(ji+1,jj) + rhoa(ji,jj) * pCd_du_ice(ji,jj) ) & 637 637 & * ( 0.5_wp * ( u_abl(ji+1,jj,2,nt_a) + u_abl(ji,jj,2,nt_a) ) - pssu_ice(ji,jj) ) … … 646 646 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! 647 647 ! ------------------------------------------------------------ ! 648 DO_2D _00_00648 DO_2D( 0, 0, 0, 0 ) 649 649 650 650 zztmp1 = 0.5_wp * ( u_abl(ji+1,jj ,2,nt_a) + u_abl(ji,jj,2,nt_a) ) … … 1259 1259 smth_b = 1._wp / 4._wp 1260 1260 ! 1261 DO_2D _11_101261 DO_2D( 1, 1, 1, 0 ) 1262 1262 zumsk = msk(ji,jj) * msk(ji+1,jj) 1263 1263 zdX ( ji, jj ) = ( pvar2d( ji+1,jj ) - pvar2d( ji ,jj ) ) * zumsk 1264 1264 END_2D 1265 1265 1266 DO_2D _10_111266 DO_2D( 1, 0, 1, 1 ) 1267 1267 zvmsk = msk(ji,jj) * msk(ji,jj+1) 1268 1268 zdY ( ji, jj ) = ( pvar2d( ji, jj+1 ) - pvar2d( ji ,jj ) ) * zvmsk 1269 1269 END_2D 1270 1270 1271 DO_2D _10_001271 DO_2D( 1, 0, 0, 0 ) 1272 1272 zFY ( ji, jj ) = zdY ( ji, jj ) & 1273 1273 & + smth_a* ( (zdX ( ji, jj+1 ) - zdX( ji-1, jj+1 )) & … … 1275 1275 END_2D 1276 1276 1277 DO_2D _00_101277 DO_2D( 0, 0, 1, 0 ) 1278 1278 zFX( ji, jj ) = zdX( ji, jj ) & 1279 1279 & + smth_a*( (zdY( ji+1, jj ) - zdY( ji+1, jj-1)) & … … 1281 1281 END_2D 1282 1282 1283 DO_2D _00_001283 DO_2D( 0, 0, 0, 0 ) 1284 1284 pvar2d( ji ,jj ) = pvar2d( ji ,jj ) & 1285 1285 & + msk(ji,jj) * smth_b * ( &
Note: See TracChangeset
for help on using the changeset viewer.