Difference between revisions of "Talk:WRF-Fire"

From openwfm
Jump to navigation Jump to search
m (typo)
 
(19 intermediate revisions by 3 users not shown)
Line 1: Line 1:
This is a talk page in blog format. Enter your questions, wishes, or anything else as separate sections at the bottom (click the + tab above). Start your response with : on a new line to offset, and sign your contribution by <nowiki>~~~~</nowiki>. See [[Help:Talk]] for more editing conventions for talk pages.
+
{{talk}}
  
==Current wish list==
+
:''The WRF-Fire wish list has been moved to [[WRF-Fire wish list]] and the discussion to [[Talk:WRF-Fire wish list]]. This page now serves for a discussion about the WRF-Fire page itself.''
* <s>Restart</s> done
 
* Read topography from file in ideal run
 
* Smoothing the terrain gradients for fire spread
 
* Canopy fire
 
* Output to WRF-Chem
 
 
 
Anything more? [[User:Jmandel|Jmandel]] 18:59, 19 January 2010 (UTC)
 
 
 
==Proposed input of additional data in ideal run==
 
 
 
* '''Single topography for both the atmosphere and the fire models''' - from an arbitrary uniform mesh aligned with coordinate axes
 
* '''Fuel data for the fire model''' - at the centers of the fire mesh
 
* '''Surface data for WRF''' -  see e.g. <code>/dyn_em/module_initialize_seabreeze2d_x.F</code> around <code> call ns_set_iswater</code>, at centers of the atmosphere mesh ("mass points")
 
 
 
The input will be done similarly as in <code>dyn_em/module_initialize_b_wave.F</code> per suggestion of the WRF developers. See <code>subroutine read_input_jet</code>.
 
 
 
This is as far as I want to go in ideal.exe, more complicated data input should be done through real.exe.
 
 
 
[[User:Jmandel|Jmandel]] 01:29, 1 February 2010 (UTC)
 
 
 
== Proposed computation of smooth terrain gradient ==
 
 
 
Currently the terrain height zsf is given at the centers of the fire grid cells. The terrain gradient is then computed in Fire by central differences, and this is the gradient that enters in the Rothermel's formula. When the terrain height is interpolated from coarser data (by bilinear interpolation), the gradient is discontinuous at data cell edges, which causes jumps in fire propagation. Instead, I propose:
 
# The gradient components will be computed on a staggered grid directly from the data, on the mesh where the data is given. That is, the x-component of the gradient will be computed at midpoints in the x-direction between the data points by central differencing, and similarly the y-component of  the gradient will be computed at midpoint the y-direction between the data points.
 
# The two components of the gradient will be interpolated to the centers of the fire grid nodes by bilinear interpolation
 
# Terrain height for the atmospheric model will be computed by bilinear interpolation from the data
 
The staggered mesh where the gradient is computed in 1. is reminiscent of the staggered mesh WRF uses for the wind; in fact, the interpolation in 2. is exactly the same that I already use for the interpolation of the wind to the fire mesh. This should result in
 
# Gradient as accurate as the resolution of the topography data allows, without an unnecessary loss of accuracy
 
# Continous gradient, without any jumps
 
# No artificial smoothing needed
 
The terrain height data can/should be given in its native resolution, unrelated to the any model mesh sizes. This will be easy to handle when terrain data is loaded from a file in an ideal run. In a real run, it might be more appropriate to do this computation in WPS.
 
[[User:Jmandel|Jmandel]] 00:48, 8 February 2010 (UTC)
 

Latest revision as of 21:01, 28 May 2010

This is a talk page in blog format. Enter your questions, wishes, or anything else as separate sections at the bottom (click the + tab above). Start your response with : on a new line to offset, and sign your contribution by ~~~~. Please reply to entries here not on the author's talk page so that the discussion stays in one piece; you may want to put a short message on the author's talk page by placing the code {{message|Talk:Name of the talk page you want to bring to their attention}} ~~~~ there. Please see Help:Talk for more on talk pages.
The WRF-Fire wish list has been moved to WRF-Fire wish list and the discussion to Talk:WRF-Fire wish list. This page now serves for a discussion about the WRF-Fire page itself.