search

lautenberger / elmfire

Eulerian Level set Model of FIRE spread
https://elmfire.io
Eclipse Public License 2.0
23 stars 11 forks source link

Fireline intensity & wind speed information missing in the initial LIST_BURNED #44

Open yqin123 opened 8 months ago

yqin123 commented 8 months ago

Hi Chris,

I found this problem when initializing a simulation with a customized phi field. I muted the surface fire spreading (by setting Ux=Uy=0) and kept the spotting as the only mechanism. However, my current implementation of the spotting algorithm requires the wind and fireline intensity information from the linked list LIST_BURNT. The initial values LIST_BURNT%WS20_NOW = 0 and LIST_BURNT%FLIN_SRUFACE = 0. These values can not initiate the spotting, thus no fire spreading.

I am wondering whether it is appropriate to initiate the initially burnt cells with the wind speed and assign fire intensity values to some of the cells (cells located at the fire perimeter, for example).

Thanks, Yiren

lautenberger commented 7 months ago

Thank you @yqin123 , this is an edge case that I had not contemplated previously. I do think your suggestion is probably the way to go. If you have an idea for how to implement this, take a crack at it and submit a PR. If not, I have some ideas for how to do it. However, it may take a couple weeks as I'm on travel until after Thanksgiving.