Changeset 4682
- Timestamp:
- 2014-06-21T11:47:21+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r4624 r4682 820 820 sqrt_Cd=kappa/((kappa/sqrt_Cd_n10) - zpsi_m) ; Cd=sqrt_Cd*sqrt_Cd; 821 821 ELSE 822 !! Shifting the wind speed to 10m and neutral stability : 823 U_n10 = dU10*1./(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) ! L & Y eq. (9a) 822 !! Shifting the wind speed to 10m and neutral stability : L & Y eq. (9a) 823 ! In very rare low-wind conditions, the old way of estimating the 824 ! neutral wind speed at 10m leads to a negative value that causes the code 825 ! to crash. To prevent this a threshold of 0.25m/s is now imposed. 826 U_n10 = MAX( 0.25 , dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) ) 824 827 825 828 !! Updating the neutral 10m transfer coefficients : … … 956 959 zpsi_m = psi_m(zeta_u) 957 960 !! 958 !! Shifting the wind speed to 10m and neutral stability : (L & Y eq.(9a)) 959 ! U_n10 = dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - psi_m(zeta_u))) 960 U_n10 = dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) 961 !! Shifting the wind speed to 10m and neutral stability : L & Y eq.(9a) 962 ! In very rare low-wind conditions, the old way of estimating the 963 ! neutral wind speed at 10m leads to a negative value that causes the code 964 ! to crash. To prevent this a threshold of 0.25m/s is now imposed. 965 U_n10 = MAX( 0.25 , dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) ) 961 966 !! 962 967 !! Shifting temperature and humidity at zu : (L & Y eq. (9b-9c))
Note: See TracChangeset
for help on using the changeset viewer.