adjust calculation of the above-canopy Richardson number

CH-Earth / summa

Structure for Unifying Multiple Modeling Alternatives:

http://www.ral.ucar.edu/projects/summa

GNU General Public License v3.0

80 stars 105 forks source link

Open martynpclark opened 9 years ago

martynpclark commented 9 years ago

The Richardson number is computed based on a the gradients in wind speed and height. Above the canopy this should be given as

(g(h_ref - d)/(u_air - u_zpd)²) * (T_air - T_cas)/(0.5_(T_air + T_cas)) instead of (g_h_ref/u_air²) * (T_air - T_cas)/(0.5*(T_air + T_cas))

This may also require adjustments to the derivative calculations because of u_zpd (need to check)

bartnijssen commented 7 years ago

Bumped milestone since 2.0 has been released

martynpclark commented 7 years ago

@arbennett -- was this fixed as part of your changes to canopy snow?