search

noaa-oar-arl / canopy-app

Stand-alone/column canopy codes and parameterizations
MIT License
6 stars 6 forks source link

Added vegtype dependent FRP to FlameH calculations. #88

Closed drnimbusrain closed 9 months ago

drnimbusrain commented 10 months ago

This PR address Closes #86 , which stems from Wei-Ting's canopy-wind/WAF paper, and the work of Alexander and Cruz (2012) (Table 1 and S1) that shows different vegetation type dependent relationships for deriving flame height based on fire intensity.

Previously we had used the single relationship based on Byram (1959) (also included in Table 1 of Alexander and Cruz (2012)) to derive flame height from FRP, and here we extend it to a methodology that is approximately related to model input vegetation type (currently supporting MODIS/VIIRS).

For the limited SE example using the FRP to flame height only option in canopy-app:

Byram (1959) Flame Heights image

Alexander and Cruz (2012) - Byram (1959) Flame Heights image

@angehung5 When you return from leave, can you more robustly help test this new feature?

drnimbusrain commented 10 months ago

I convert this to a draft, as when reading Alexander and Cruz, it may be more appropriate to use a vegtype dependent calculations of crown scorch height (hs), which can be related directly to flame height (hf) (see Table 2 and equation 14 in their paper). I will add this option as well.

drnimbusrain commented 10 months ago

I have added the additional option (flameh_frp) to use the vegtype dependent scorch height to flame height relationships from Alexander and Cruz.

Overall, the differences compared to Byram (1959) only are similar as previous for the vegtype dependence, but the different relationship for grasslands/savanna/shublands/crops using scorch heights cause larger flame heights in this case for one grid cell: image

drnimbusrain commented 9 months ago

@angehung5 What do you think about this PR, will you give it a review? Not sure it will have enough impact yet to warrant changing your paper results for midflame wind and WAF. However, you mentioned you had comments on the new NL variable flameh_frp. Also note that I have merged in the latest develop which contains multiple timesteps. Thank you!