Changeset 11631
- Timestamp:
- 2019-10-01T20:22:39+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC
- Files:
-
- 4 edited
-
sbcblk_algo_coare3p0.F90 (modified) (3 diffs)
-
sbcblk_algo_coare3p6.F90 (modified) (1 diff)
-
sbcblk_algo_ecmwf.F90 (modified) (1 diff)
-
sbcblk_algo_ncar.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p0.F90
r11626 r11631 38 38 USE sbc_oce ! Surface boundary condition: ocean fields 39 39 USE sbcblk_phy ! all thermodynamics functions, rho_air, q_sat, etc... !LB 40 USE sbcblk_skin_coare 40 USE sbcblk_skin_coare ! cool-skin/warm layer scheme (CSWL_ECMWF) !LB 41 41 42 42 IMPLICIT NONE … … 323 323 ztmp2 = u_star/vkarmn*LOG(10./z0) ! Neutral wind speed at 10m 324 324 z0 = alfa_charn_3p0(ztmp2)*ztmp1/grav + 0.11_wp*znu_a/u_star ! Roughness length (eq.6) 325 ztmp1 = ( znu_a / (z0*u_star) )**0.6_wp 325 ztmp1 = ( znu_a / (z0*u_star) )**0.6_wp ! (1./Re_r)^0.72 (Re_r: roughness Reynolds number) COARE3.6-specific! 326 326 z0t = MIN( 1.1E-4_wp , 5.5E-5_wp*ztmp1 ) ! Scalar roughness for both theta and q (eq.28) #LOLO: some use 1.15 not 1.1 !!! 327 327 … … 439 439 ! 440 440 ! Charnock's constant, increases with the wind : 441 zgt10 = 0.5 + SIGN(0.5_wp,(zw - 10)) ! If zw<10. --> 0, else --> 1441 zgt10 = 0.5 + SIGN(0.5_wp,(zw - 10)) ! If zw<10. --> 0, else --> 1 442 442 zgt18 = 0.5 + SIGN(0.5_wp,(zw - 18.)) ! If zw<18. --> 0, else --> 1 443 443 ! -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_coare3p6.F90
r11626 r11631 38 38 USE sbc_oce ! Surface boundary condition: ocean fields 39 39 USE sbcblk_phy ! all thermodynamics functions, rho_air, q_sat, etc... !LB 40 USE sbcblk_skin_coare 40 USE sbcblk_skin_coare ! cool-skin/warm layer scheme (CSWL_ECMWF) !LB 41 41 42 42 IMPLICIT NONE -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_ecmwf.F90
r11619 r11631 44 44 PRIVATE 45 45 46 PUBLIC :: 46 PUBLIC :: TURB_ECMWF ! called by sbcblk.F90 47 47 48 48 ! !! ECMWF own values for given constants, taken form IFS documentation... -
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_algo_ncar.F90
r11221 r11631 173 173 END IF 174 174 175 ! Update neutral wind speed at 10m and neutral Cd at 10m (L&Y 2004 eq. 9a)... 176 ! In very rare low-wind conditions, the old way of estimating the 177 ! neutral wind speed at 10m leads to a negative value that causes the code 178 ! to crash. To prevent this a threshold of 0.25m/s is imposed. 175 179 ztmp2 = psi_m(zeta_u) 176 180 IF( ln_cdgw ) THEN ! surface wave case
Note: See TracChangeset
for help on using the changeset viewer.
