Difference between revisions of "Vertical wind interpolation"
| Line 5: | Line 5: | ||
In the code, we thus have two choices: | In the code, we thus have two choices: | ||
| − | # Interpolate | + | # [https://github.com/jbeezley/wrf-fire/blob/master/wrfv2_fire/phys/module_fr_sfire_atm.F#L184 Interpolate] to 6.096m and [https://github.com/jbeezley/wrf-fire/blob/master/wrfv2_fire/phys/module_fr_sfire_atm.F#L622 apply] the wind reduction factors. This is much simpler and faster. |
| − | # Interpolate the wind at each point of the fire mesh to the midflame height at that point separately. | + | # Interpolate the wind at each point of the fire mesh to the midflame height at that point [https://github.com/jbeezley/wrf-fire/blob/master/wrfv2_fire/phys/module_fr_sfire_atm.F#L874 separately]. This is more complicated. |
See Mandel et al. 2011 for more details. | See Mandel et al. 2011 for more details. | ||
| − | The vertical interpolation is | + | The vertical interpolation is somewhat complicated because of the way how WRF represents [[How_to_interpret_WRF-Fire_variables#Wind|the wind speed]] and [[How_to_interpret_WRF-Fire_variables#Location|the vertical coordinate]]. In particular, because the altitude of nodes in WRF changes during the run, the number of vertical levels needed for the interpolation is not known in advance. |
| − | |||
==References== | ==References== | ||
Revision as of 07:35, 1 August 2011
Vertical log interpolation and wind reduction factors
The wind speed that enters the spread rate formula is found by vertical interpolation to a specified height fwh, sometimes called midflame height, using the ideal logarithmic wind profile: the wind speed at height z is assumed to be proportional to log z/z0, where z0 is the roughness height (clearly, this is the height where the wind speed is zero). The special logarithmic interpolation we use preserves such functions just like linear interpolation preserves linear functions. The numbers fwh and z/z0 need to be known at every point of the fire mesh.
In BEHAVE, the wind is measured at the 20ft (6.096m) and log interpolated to the midflame height. The interpolation is the same as multiplying the wind speed by a given wind reduction factor windrf (Baughman and Albini, 1980).
In the code, we thus have two choices:
- Interpolate to 6.096m and apply the wind reduction factors. This is much simpler and faster.
- Interpolate the wind at each point of the fire mesh to the midflame height at that point separately. This is more complicated.
See Mandel et al. 2011 for more details.
The vertical interpolation is somewhat complicated because of the way how WRF represents the wind speed and the vertical coordinate. In particular, because the altitude of nodes in WRF changes during the run, the number of vertical levels needed for the interpolation is not known in advance.
References
- Robert G. Baughman and Frank A. Albini, Estimating Midflame Windspeeds, Sixth Conference on Fire and Forest Meteorology, Seattle, WA April 22-24, 1980, pp. 88-92 pdf
- Jan Mandel, Jonathan D. Beezley, and Adam K. Kochanski, Coupled atmosphere-wildland fire modeling with WRF 3.3 and SFIRE 2011, Geoscientific Model Development (GMD) 4, 591-610, 2011. doi:10.5194/gmd-4-591-2011