Changeset 14072 for NEMO/trunk/src/ABL/ablmod.F90
- Timestamp:
- 2020-12-04T08:48:38+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ABL/ablmod.F90
r13295 r14072 18 18 USE dom_oce, ONLY : tmask 19 19 USE sbc_oce, ONLY : ght_abl, ghw_abl, e3t_abl, e3w_abl, jpka, jpkam1, rhoa 20 USE sbcblk ! use rn_efac , cdn_oce21 USE sbc blk_phy ! use some physical constants for flux computation20 USE sbcblk ! use rn_efac 21 USE sbc_phy ! Catalog of functions for physical/meteorological parameters in the marine boundary layer 22 22 ! 23 23 USE prtctl ! Print control (prt_ctl routine) … … 135 135 ustar2(ji,jj) = zzoce 136 136 #endif 137 zrough(ji,jj) = ght_abl(2) * EXP( - vkarmn / SQRT( MAX( Cdn_oce(ji,jj), 1.e-4 ) ) ) !<-- recover the value of z0 from Cdn_oce 137 !#LB: sorry Cdn_oce is gone: 138 !zrough(ji,jj) = ght_abl(2) * EXP( - vkarmn / SQRT( MAX( Cdn_oce(ji,jj), 1.e-4 ) ) ) !<-- recover the value of z0 from Cdn_oce 138 139 END_2D 140 141 zrough(:,:) = z0_from_Cd( ght_abl(2), pCd_du(:,:) / MAX( pwndm(:,:), 0.5_wp ) ) ! #LB: z0_from_Cd is define in sbc_phy.F90... 142 139 143 ! 140 144 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< … … 252 256 & + rDt_abl * fft_abl(ji, jj) * v_abl( ji, jj, jk, nt_n ) ) & 253 257 & / (1._wp + gamma_Cor*gamma_Cor*zcff) 254 258 255 259 v_abl( ji, jj, jk, nt_a ) = e3t_abl(jk) *( & 256 260 & (1._wp-gamma_Cor*(1._wp-gamma_Cor)*zcff) * v_abl( ji, jj, jk, nt_n ) & … … 590 594 591 595 DO_2D( 0, 1, 0, 1 ) 592 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * ( pssu(ji ,jj) + pssu(ji-1,jj) ) 593 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * ( pssv(ji,jj ) + pssv(ji,jj-1) ) 596 zwnd_i(ji,jj) = u_abl(ji ,jj,2,nt_a) - 0.5_wp * ( pssu(ji ,jj) + pssu(ji-1,jj) ) 597 zwnd_j(ji,jj) = v_abl(ji,jj ,2,nt_a) - 0.5_wp * ( pssv(ji,jj ) + pssv(ji,jj-1) ) 594 598 END_2D 595 599 ! … … 633 637 ! Wind stress relative to the moving ice ( U10m - U_ice ) ! 634 638 ! ------------------------------------------------------------ ! 635 DO_2D( 0, 0, 0, 0 ) 639 DO_2D( 0, 0, 0, 0 ) 636 640 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 641 & * ( 0.5_wp * ( u_abl(ji+1,jj,2,nt_a) + u_abl(ji,jj,2,nt_a) ) - pssu_ice(ji,jj) )
Note: See TracChangeset
for help on using the changeset viewer.