[Physics] Adding Bonan Wind Profiles with Unified RSL Approach

[drnimbusrain]

@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