Changeset 4681


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

#1345 : correct a bug in CORE bulk formulae in dev_v3_4_STABLE_2012

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r3294 r4681  
    721721           sqrt_Cd=kappa/((kappa/sqrt_Cd_n10) - zpsi_m) ; Cd=sqrt_Cd*sqrt_Cd; 
    722722         ELSE 
    723            !! Shifting the wind speed to 10m and neutral stability : 
    724            U_n10 = dU10*1./(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) !  L & Y eq. (9a) 
     723           !! Shifting the wind speed to 10m and neutral stability : L & Y eq. (9a) 
     724           !   In very rare low-wind conditions, the old way of estimating the 
     725           !   neutral wind speed at 10m leads to a negative value that causes the code 
     726           !   to crash. To prevent this a threshold of 0.25m/s is now imposed. 
     727           U_n10 = MAX( 0.25 , dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) ) 
    725728 
    726729           !! Updating the neutral 10m transfer coefficients : 
     
    857860         zpsi_m  = psi_m(zeta_u) 
    858861         !! 
    859          !! Shifting the wind speed to 10m and neutral stability : (L & Y eq.(9a)) 
    860 !        U_n10 = dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - psi_m(zeta_u))) 
    861          U_n10 = dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) 
     862         !! Shifting the wind speed to 10m and neutral stability : L & Y eq.(9a) 
     863         !   In very rare low-wind conditions, the old way of estimating the 
     864         !   neutral wind speed at 10m leads to a negative value that causes the code 
     865         !   to crash. To prevent this a threshold of 0.25m/s is now imposed. 
     866         U_n10 = MAX( 0.25 , dU10/(1. + sqrt_Cd_n10/kappa*(log(zu/10.) - zpsi_m)) ) 
    862867         !! 
    863868         !! Shifting temperature and humidity at zu :          (L & Y eq. (9b-9c)) 
Note: See TracChangeset for help on using the changeset viewer.