@zmoon @angehung5 Decided to make this a separate issue to add the Bonan wind profiles with Unified RSL approach (eventual RSL_OPT = 1) to canopy-app, so we can be more robust and fully compare in-canopy wind profiles to Massman/MOST approach (RSL_OPT=0) for Wei-Ting's paper as well.
@zmoon @angehung5 Decided to make this a separate issue to add the Bonan wind profiles with Unified RSL approach (eventual RSL_OPT = 1) to canopy-app, so we can be more robust and fully compare in-canopy wind profiles to Massman/MOST approach (RSL_OPT=0) for Wei-Ting's paper as well.
The code for the 2018 paper can be found at https://github.com/gbonan/CLM-ml_v0
This isn't all of the code you need to run it; if interested I have a full setup at https://github.com/zmoon/CLM-ml-mb
Wind profile computation starts with the above-canopy part (Eq. 19 in the paper): https://github.com/gbonan/CLM-ml_v0/blob/881a5c44b86f92c4ce65b3268b190297fbb0bd65/CanopyTurbulenceMod.F90#L432 This uses utility functions here: https://github.com/gbonan/CLM-ml_v0/blob/881a5c44b86f92c4ce65b3268b190297fbb0bd65/CanopyTurbulenceMod.F90#L921 Obukhov length needed but we could make it an input.
Followed by within-canopy (Eq. 21 in the paper): https://github.com/gbonan/CLM-ml_v0/blob/881a5c44b86f92c4ce65b3268b190297fbb0bd65/CanopyTurbulenceMod.F90#L478