- Timestamp:
- 2021-12-03T20:32:50+01:00 (3 years ago)
- Location:
- NEMO/branches/2021/dev_r14318_RK3_stage1
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14318_RK3_stage1
- Property svn:externals
-
old new 9 9 10 10 # SETTE 11 ^/utils/CI/sette@14244 sette 11 ^/utils/CI/sette@HEAD sette 12
-
- Property svn:externals
-
NEMO/branches/2021/dev_r14318_RK3_stage1/src/ABL/ablmod.F90
r14239 r15574 128 128 !! needed for surface boundary condition of TKE 129 129 !! pwndm contains | U10m - U_oce | (see blk_oce_1 in sbcblk) 130 DO_2D( 1, 1, 1, 1)130 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 131 131 zzoce = pCd_du (ji,jj) * pwndm (ji,jj) 132 132 #if defined key_si3 … … 250 250 ! 251 251 ! Advance u_abl & v_abl to time n+1 252 DO_2D( 1, 1, 1, 1)252 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 253 253 zcff = ( fft_abl(ji,jj) * rDt_abl )*( fft_abl(ji,jj) * rDt_abl ) ! (f dt)**2 254 254 … … 488 488 DO jk = 2, jpka ! outer loop 489 489 !------------- 490 DO_2D( 0, 1, 0, 1)490 DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls ) 491 491 zcff1 = pblh( ji, jj ) 492 492 zsig = ght_abl(jk) / MAX( jp_pblh_min, MIN( jp_pblh_max, zcff1 ) ) … … 511 511 DO jk = 2, jpka ! outer loop 512 512 !------------- 513 DO_2D( 1, 1, 1, 1)513 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 514 514 zcff1 = pblh( ji, jj ) 515 515 zsig = ght_abl(jk) / MAX( jp_pblh_min, MIN( jp_pblh_max, zcff1 ) ) … … 534 534 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 535 535 ! 536 CALL lbc_lnk _multi( 'ablmod', u_abl(:,:,:,nt_a ), 'T', -1._wp, v_abl(:,:,:,nt_a) , 'T', -1._wp )537 CALL lbc_lnk _multi( 'ablmod', tq_abl(:,:,:,nt_a,jp_ta), 'T', 1._wp , tq_abl(:,:,:,nt_a,jp_qa), 'T', 1._wp , kfillmode = jpfillnothing ) ! ++++ this should not be needed...536 CALL lbc_lnk( 'ablmod', u_abl(:,:,:,nt_a ), 'T', -1._wp, v_abl(:,:,:,nt_a) , 'T', -1._wp ) 537 CALL lbc_lnk( 'ablmod', tq_abl(:,:,:,nt_a,jp_ta), 'T', 1._wp , tq_abl(:,:,:,nt_a,jp_qa), 'T', 1._wp , kfillmode = jpfillnothing ) ! ++++ this should not be needed... 538 538 ! 539 539 #if defined key_xios … … 585 585 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 586 586 ! 587 DO_2D( 1, 1, 1, 1)587 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 588 588 ztemp = tq_abl( ji, jj, 2, nt_a, jp_ta ) 589 589 zhumi = tq_abl( ji, jj, 2, nt_a, jp_qa ) … … 595 595 END_2D 596 596 597 DO_2D( 0, 1, 0, 1)597 DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls ) 598 598 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * ( pssu(ji ,jj) + pssu(ji-1,jj) ) 599 599 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * ( pssv(ji,jj ) + pssv(ji,jj-1) ) 600 600 END_2D 601 601 ! 602 CALL lbc_lnk _multi( 'ablmod', zwnd_i(:,:) , 'T', -1.0_wp, zwnd_j(:,:) , 'T', -1.0_wp )602 CALL lbc_lnk( 'ablmod', zwnd_i(:,:) , 'T', -1.0_wp, zwnd_j(:,:) , 'T', -1.0_wp ) 603 603 ! 604 604 ! ... scalar wind ( = | U10m - U_oce | ) at T-point (masked) 605 DO_2D( 1, 1, 1, 1)605 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 606 606 zcff = SQRT( zwnd_i(ji,jj) * zwnd_i(ji,jj) & 607 607 & + zwnd_j(ji,jj) * zwnd_j(ji,jj) ) ! * msk_abl(ji,jj) … … 625 625 END_2D 626 626 ! 627 CALL lbc_lnk _multi( 'ablmod', ptaui(:,:), 'U', -1.0_wp, ptauj(:,:), 'V', -1.0_wp )627 CALL lbc_lnk( 'ablmod', ptaui(:,:), 'U', -1.0_wp, ptauj(:,:), 'V', -1.0_wp ) 628 628 629 629 CALL iom_put( "taum_oce", ptaum ) … … 645 645 & * ( 0.5_wp * ( v_abl(ji,jj+1,2,nt_a) + v_abl(ji,jj,2,nt_a) ) - pssv_ice(ji,jj) ) 646 646 END_2D 647 CALL lbc_lnk _multi( 'ablmod', ptaui_ice, 'U', -1.0_wp, ptauj_ice, 'V', -1.0_wp )647 CALL lbc_lnk( 'ablmod', ptaui_ice, 'U', -1.0_wp, ptauj_ice, 'V', -1.0_wp ) 648 648 ! 649 649 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab2d_1=ptaui_ice , clinfo1=' abl_stp: putaui : ' & … … 664 664 & * ( zztmp2 - pssv_ice(ji,jj) ) 665 665 END_2D 666 CALL lbc_lnk _multi( 'ablmod', ptaui_ice, 'U', -1.0_wp, ptauj_ice,'V', -1.0_wp )666 CALL lbc_lnk( 'ablmod', ptaui_ice, 'U', -1.0_wp, ptauj_ice,'V', -1.0_wp ) 667 667 ! 668 668 IF(sn_cfctl%l_prtctl) THEN … … 1265 1265 smth_b = 1._wp / 4._wp 1266 1266 ! 1267 DO_2D( 1, 0, 1, 1)1267 DO_2D( nn_hls, nn_hls-1, nn_hls, nn_hls ) 1268 1268 zumsk = msk(ji,jj) * msk(ji+1,jj) 1269 1269 zdX ( ji, jj ) = ( pvar2d( ji+1,jj ) - pvar2d( ji ,jj ) ) * zumsk 1270 1270 END_2D 1271 1271 1272 DO_2D( 1, 1, 1, 0)1272 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls-1 ) 1273 1273 zvmsk = msk(ji,jj) * msk(ji,jj+1) 1274 1274 zdY ( ji, jj ) = ( pvar2d( ji, jj+1 ) - pvar2d( ji ,jj ) ) * zvmsk 1275 1275 END_2D 1276 1276 1277 DO_2D( 0, 0, 1, 0)1277 DO_2D( nn_hls-1, nn_hls-1, nn_hls, nn_hls-1 ) 1278 1278 zFY ( ji, jj ) = zdY ( ji, jj ) & 1279 1279 & + smth_a* ( (zdX ( ji, jj+1 ) - zdX( ji-1, jj+1 )) & … … 1281 1281 END_2D 1282 1282 1283 DO_2D( 1, 0, 0, 0)1283 DO_2D( nn_hls, nn_hls-1, nn_hls-1, nn_hls-1 ) 1284 1284 zFX( ji, jj ) = zdX( ji, jj ) & 1285 1285 & + smth_a*( (zdY( ji+1, jj ) - zdY( ji+1, jj-1)) & … … 1287 1287 END_2D 1288 1288 1289 DO_2D( 0, 0, 0, 0)1289 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 1290 1290 pvar2d( ji ,jj ) = pvar2d( ji ,jj ) & 1291 1291 & + msk(ji,jj) * smth_b * ( &
Note: See TracChangeset
for help on using the changeset viewer.