Difference between revisions of "WRF-Fire development notes"

From openwfm
Jump to navigation Jump to search
(WPS)
(→‎In progress: add instructions, rm names - already identified by branches)
Line 15: Line 15:
  
 
==In progress==
 
==In progress==
 +
 +
:''If you start on an item from the wish list, please move the item here and add links to the work branches.''
  
 
* Better quadrature for fuel left
 
* Better quadrature for fuel left
** The quadrature is done separately on 4 smaller cells per fire cell. The interpolation of time from ignition to the smaller cells should take into consideration the case when some of the nodes interpolated from are ignited and some are not. Volodymyr Kondratenko, [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/jm2/fuel branch jm2/fuel]
+
** The quadrature is done separately on 4 smaller cells per fire cell. The interpolation of time from ignition to the smaller cells should take into consideration the case when some of the nodes interpolated from are ignited and some are not. [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/jm2/fuel branch jm2/fuel]
** The current scheme is second order accurate when all 4 corners of a fire mesh cell are on fire, and exact when all 4 corners are not on fire. In the case when only some of the 4 corners are on fire, the only requirements are that the transition when nodes ignite is continuous and monotonous. In that case, the scheme may not be very accurate. A better scheme would be accurate in more cases and have natural invariance properties. Minjeong Kim, [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/jm2/mkim branch jm2/mkim] [http://repo.or.cz/w/wrffire.git/commit/4cf71f8d9e4309281321ad7d38e4f028a005fa8a commit]
+
** The current scheme is second order accurate when all 4 corners of a fire mesh cell are on fire, and exact when all 4 corners are not on fire. In the case when only some of the 4 corners are on fire, the only requirements are that the transition when nodes ignite is continuous and monotonous. In that case, the scheme may not be very accurate. A better scheme would be accurate in more cases and have natural invariance properties. [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/jm2/mkim branch jm2/mkim] [http://repo.or.cz/w/wrffire.git/commit/4cf71f8d9e4309281321ad7d38e4f028a005fa8a commit]
 
:This affects the amount of heat output from the cell in the timesteps when the fireline crosses the cell.
 
:This affects the amount of heat output from the cell in the timesteps when the fireline crosses the cell.
 
* Input to the fire model the gradient rather than terrain height.  This will allow the gradient to be smoothly interpolated in preprocessing in the WPS or in setting up ideal cases. No artificial terrain smoothing is needed. The current scheme interpolates terrain height linearly, and if the mesh the terrain is interpolated from is much coarser than the fire mesh, this results in gradient jumps along the terrain mesh lines. Then, the  fire spread rate jumps along the terrain mesh lines. ''See [[Talk:WRF-Fire wish list#Proposed computation of smooth terrain gradient|discussion]].'' [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/jm2/grad Branch jm2/grad].  
 
* Input to the fire model the gradient rather than terrain height.  This will allow the gradient to be smoothly interpolated in preprocessing in the WPS or in setting up ideal cases. No artificial terrain smoothing is needed. The current scheme interpolates terrain height linearly, and if the mesh the terrain is interpolated from is much coarser than the fire mesh, this results in gradient jumps along the terrain mesh lines. Then, the  fire spread rate jumps along the terrain mesh lines. ''See [[Talk:WRF-Fire wish list#Proposed computation of smooth terrain gradient|discussion]].'' [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/jm2/grad Branch jm2/grad].  
* Read the topography from a file in an ideal run.  Fix/extend surface initialization in ideal cases to allow for defining surface properties like roughness, ground temperature, albedo etc., required by soil models and radiation codes. See [[Talk:WRF-Fire wish list#Proposed input of additional data in ideal run|discussion]]. Adam Kochanski, [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/adamko/surface branch adamko/surface]
+
* Read the topography from a file in an ideal run.  Fix/extend surface initialization in ideal cases to allow for defining surface properties like roughness, ground temperature, albedo etc., required by soil models and radiation codes. See [[Talk:WRF-Fire wish list#Proposed input of additional data in ideal run|discussion]]. [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/adamko/surface branch adamko/surface]
* Support hyperbolic vertical mesh grading in ideal run not only exponential. Adam Kochanski, [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/adamko/surface branch adamko/surface]
+
* Support hyperbolic vertical mesh grading in ideal run not only exponential. [http://repo.or.cz/w/wrffire.git/shortlog/refs/heads/adamko/surface branch adamko/surface]
  
 
==Done==
 
==Done==

Revision as of 15:20, 30 May 2010

This page tracks activity in the WRF-Fire software development.

Wish list

  • Get rid of the ugly includes in the argument lists to get the coefficients down the call chain to the spread rate calculations. Pass the coefficients by pointers defined in a module instead. Declare the pointers in the module, or pass them around in a single derived type argument. This is a necessary step before including the weather and fire components into VisTrails separately.
  • A standalone driver for the fire code that is independent of WRF. See Talk:WRF-Fire development notes#Standalone driver
  • Merge with WRF 3.2
  • Modify WPS to compute terrain gradient at data resolution and interpolate it to the fire mesh.
  • Gradual ignition from a point or a line, without limitation on the initial radius.
    • The current scheme requires that the ignition area is several fire mesh steps large. When such area ignites all at once, all the heat is output to the atmosphere at once, instead of being output gradually and allowing ground layer circulation to develop over time.
  • Canopy fire
    • Canopy fire collocated with ground fire is easy. A separate canopy fire would be harder.
  • Output to WRF-Chem
    • Output of the combustion products into WRF is also needed to visualize smoke.

In progress

If you start on an item from the wish list, please move the item here and add links to the work branches.
  • Better quadrature for fuel left
    • The quadrature is done separately on 4 smaller cells per fire cell. The interpolation of time from ignition to the smaller cells should take into consideration the case when some of the nodes interpolated from are ignited and some are not. branch jm2/fuel
    • The current scheme is second order accurate when all 4 corners of a fire mesh cell are on fire, and exact when all 4 corners are not on fire. In the case when only some of the 4 corners are on fire, the only requirements are that the transition when nodes ignite is continuous and monotonous. In that case, the scheme may not be very accurate. A better scheme would be accurate in more cases and have natural invariance properties. branch jm2/mkim commit
This affects the amount of heat output from the cell in the timesteps when the fireline crosses the cell.
  • Input to the fire model the gradient rather than terrain height. This will allow the gradient to be smoothly interpolated in preprocessing in the WPS or in setting up ideal cases. No artificial terrain smoothing is needed. The current scheme interpolates terrain height linearly, and if the mesh the terrain is interpolated from is much coarser than the fire mesh, this results in gradient jumps along the terrain mesh lines. Then, the fire spread rate jumps along the terrain mesh lines. See discussion. Branch jm2/grad.
  • Read the topography from a file in an ideal run. Fix/extend surface initialization in ideal cases to allow for defining surface properties like roughness, ground temperature, albedo etc., required by soil models and radiation codes. See discussion. branch adamko/surface
  • Support hyperbolic vertical mesh grading in ideal run not only exponential. branch adamko/surface

Done

  • Restart commit
  • Walking line ignition (as the fireman walks igniting the fire) commit
  • Remove (by default) the warm perturbation (bubble) from atmospheric initializations in idealized case. Modify the code and add namelist variables allowing for ingesting the cold/warm perturbation of a specified size (for people interested in fire propagation in downburst/convective environments) (Adam Kochanski) commit

External links