NOAA-OWP / evapotranspiration

Other
3 stars 9 forks source link

Modifying momentum transfer roughness length approximation equation #15

Closed WanruWu-NOAA closed 2 years ago

WanruWu-NOAA commented 2 years ago

To fix the unrealistically large values from PET module options 2, 3 and 5:

  1. Commenting out the zero-plane displacement height (d) function of the vegetation (canopy) height (hc) so to input the d value directly from the config file;
  2. Modifying the momentum-transfer roughness length (zom) from the function of hc to a function of d and keeping the heat-transfer roughness length (zoh) function of zom to accommodate the magnitude logic among the three parameters d, zom and zoh;
  3. Recompile PET module and rebuilt ngen.

Additions

In the config file, after the d value is set, the zom and zoh can be set as "0", and the functions in the code will pick up the d value to compute zom and zoh.

Removals

  1. Commenting out the zero-plane displacement height function of the vegetation (canopy) height //model->pet_params.zero_plane_displacement_height_m=2.0/3.0*model->pet_params.vegetation_height_m;

Changes

  1. Modifying the momentum-transfer roughness length from the function of the vegetation height to a function of the zero-plane displacement height; //model->pet_params.momentum_transfer_roughness_length_m=0.123*model->pet_params.vegetation_height_m; model->pet_params.momentum_transfer_roughness_length_m=0.1845*model->pet_params.zero_plane_displacement_height_m;

Testing

  1. Start with d = 0.001

Screenshots

Notes

-

Todos

-

Checklist

Testing checklist

Target Environment support

Accessibility

Other