Changeset 4682


Ignore:
Timestamp:
2014-06-21T11:47:21+02:00 (7 years ago)
Author:
gm
Message:

#1345 : correct a bug in CORE bulk formulae in the v3.6alpha (trunk)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r4624 r4682  
    820820           sqrt_Cd=kappa/((kappa/sqrt_Cd_n10) - zpsi_m) ; Cd=sqrt_Cd*sqrt_Cd; 
    821821         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)) ) 
    824827 
    825828           !! Updating the neutral 10m transfer coefficients : 
     
    956959         zpsi_m  = psi_m(zeta_u) 
    957960         !! 
    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)) ) 
    961966         !! 
    962967         !! Shifting temperature and humidity at zu :          (L & Y eq. (9b-9c)) 
Note: See TracChangeset for help on using the changeset viewer.