Difference between revisions of "SFIRE variables"

From openwfm
Jump to navigation Jump to search
(Created page with '400px The fire mesh is 2D. Every atmosphere cell is divided into '''sr_x''' by '''sr_y''' fire cells. The submesh ratios '''sr_x''' and '''sr_y''' are spec…')
 
 
(21 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
{{users guide}}
 +
{{Hangon}}
 +
 +
==Introduction==
 
[[File:firemesh.png|400px]]
 
[[File:firemesh.png|400px]]
  
Line 10: Line 14:
 
  from netCDF4 import Dataset
 
  from netCDF4 import Dataset
 
  f=Dataset('wrfout','r')
 
  f=Dataset('wrfout','r')
−
  sr_x=len(f.dimensions['west_east_subgrid'])/(len(f.dimensions['west_east'])+1)
+
  sr_x=len(f.dimensions['west_east_subgrid'])//(len(f.dimensions['west_east'])+1)
−
  sr_y=len(f.dimensions['south_north_subgrid'])/(len(f.dimensions['south_north'])+1)
+
  sr_y=len(f.dimensions['south_north_subgrid'])//(len(f.dimensions['south_north'])+1)
−
  fgrnhfx=f.variables['FGRNHFX']
+
  fgrnhfx=f.variables['FGRNHFX'][:,:-sr_y,:-sr_x]
−
fgrnhfx=fgrnhfx[...,:-sr_y,:-sr_x]
+
 
 +
Variables on the fire mesh are located at the centers of the fire mesh cells of a 2D fire mesh.
 +
 
 +
==List of variables==
 +
The authoritative information can be always found in {{WRF-Fire-file|wrfv2_fire/Registry/registry.fire|the Registry}}.
 +
 
 +
{| border="1"
 +
|+ Fire variables on fire grid
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | NFUEL_CAT
 +
|fuel data
 +
|
 +
|-
 +
| scope="row" | ZSF
 +
|height of surface above sea level
 +
|m
 +
|-
 +
| scope="row" | DZDXF
 +
|surface gradient x
 +
|1
 +
|-
 +
| scope="row" | DZDYF
 +
|surface gradient y
 +
|1
 +
|}
 +
 
 +
 
 +
{| border="1"
 +
|+Fire variables on atm grid
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | RTHFRTEN
 +
|temperature tendency
 +
|K/s
 +
|-
 +
| scope="row" | RQVFRTEN
 +
|humidity tendency
 +
|
 +
|}
 +
 
 +
 
 +
{| border="1"
 +
|+Diagnostics
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | AVG_FUEL_FRAC
 +
|fuel remaining averaged to atmospheric grid
 +
|1
 +
|-
 +
| scope="row" | GRNHFX
 +
|heat flux from ground fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | GRNQFX
 +
|moisture flux from ground fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | CANHFX
 +
|heat flux from crown fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | CANQFX
 +
|moisture flux from crown fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | UAH
 +
|wind at fire_wind_heigh
 +
|m/s
 +
|-
 +
| scope="row" | VAH
 +
|wind at fire_wind_heigh
 +
|m/s
 +
|}
 +
 
 +
 
 +
{| border="1"
 +
|+Sfire variables on fire grid
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | TIGN_G
 +
|ignition time on ground
 +
|s
 +
|-
 +
| scope="row" | LFN
 +
|level function
 +
|1
 +
|-
 +
| scope="row" | FUEL_FRAC
 +
|fuel remaining
 +
|1
 +
|-
 +
| scope="row" | FMC_G
 +
|fuel moisture contents
 +
|1
 +
|-
 +
| scope="row" | FIRE_AREA
 +
|fraction of cell area on fire
 +
|1
 +
|-
 +
| scope="row" | UF
 +
|fire wind
 +
|m/s
 +
|-
 +
| scope="row" | VF
 +
|fire wind
 +
|m/s
 +
|-
 +
| scope="row" | FGRNHFX
 +
|heat flux from ground fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | FGRNQFX
 +
|moisture flux from ground fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | FCANHFX
 +
|heat flux from crown fire
 +
|W/m<sup>2
 +
|-
 +
| scope="row" | FCANQFX
 +
|moisture flux from crown fire
 +
|W/m<sup>2
 +
|}
 +
 
 +
 
 +
{| border="1"
 +
|+Diagnostics for the actual modeled fire
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | ROS
 +
|rate of spread in the normal direction to the fireline
 +
|m/s
 +
|-
 +
| scope="row" | FLINEINT
 +
|fireline intensity
 +
|W/m
 +
|-
 +
| scope="row" | FLINEINT2
 +
|alternative fireline intensity
 +
|J/m/s<sup>2
 +
|}
 +
 
 +
 
 +
{| border="1"
 +
|+Diagnostics for fire risk rating - independent on any actual fire going on
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | F_ROS0
 +
|base rate of spread in all directions
 +
|m/s
 +
|-
 +
| scope="row" | F_ROSX
 +
|X component of the spread vector driven by wind and slope
 +
|m/s
 +
|-
 +
| scope="row" | F_ROSY
 +
|Y component of the spread vector driven by wind and slope
 +
|m/s
 +
|-
 +
| scope="row" | F_ROS
 +
|max spread rate in any direction
 +
|m/s
 +
|-
 +
| scope="row" | F_INT
 +
|fire reaction intensity for risk rating, without fire
 +
|J/m<sup>2</sup>/s
 +
|-
 +
| scope="row" | F_LINEINT
 +
|Byram fireline intensity for risk rating, without fire
 +
|J/m/s
 +
|-
 +
| scope="row" | F_LINEINT2
 +
|alternative fireline intensity for risk rating, without fire"
 +
|J/m/s<sup>2
 +
|}
 +
 
 +
 
 +
{| border="1"
 +
|+Constant data arrays
 +
| scope="col" | Variable name
 +
| scope="col" | Description
 +
| scope="col" | Unit
 +
|-
 +
| scope="row" | FXLONG
 +
|longitude of midpoints of fire cells
 +
|degrees
 +
|-
 +
| scope="row" |FXLAT
 +
|latitude of midpoints of fire cells
 +
|degrees
 +
|-
 +
| scope="row" |FUEL_TIME
 +
|fuel
 +
|
 +
|-
 +
| scope="row" |BBB
 +
|fuel
 +
|
 +
|-
 +
| scope="row" |PHISC
 +
|fuel
 +
|
 +
|-
 +
| scope="row" |PHIWC
 +
|fuel
 +
|
 +
|-
 +
| scope="row" |R_0
 +
|fuel
 +
|
 +
|-
 +
| scope="row" |FGIP
 +
|fuel
 +
|
 +
|-
 +
| scope="row" |FZ0
 +
|fuel roughness height
 +
|
 +
|-
 +
| scope="row" |FWH
 +
|fuel fire wind height
 +
|
 +
|-
 +
| scope="row" |ISCHAP
 +
|fuel
 +
|
 +
|}
 +
 
 +
== Fire intensities ==
 +
We introduce three fire intensities:
 +
 
 +
1. Reaction : The heat release rate per unit area of the front.
 +
 
 +
2. Byram :  The heat produced per unit length of the fireline in unit time (J/m/s) in the so-called flaming zone behind the fireline.
 +
 
 +
3. Fireline, new : The amount of heat generated by the advancing fireline from the newly burning fuel only, in the unit of time.
 +
 
 +
Definitions and implementations of fire intensities include variables :
 +
 
 +
* H<sub>F</sub> : Heat contents of the fuel (J/kg)
 +
 
 +
* R<sub>S</sub> : Fire spread rate (m/s)
 +
 
 +
* F : Fuel remaining at time t (kg/m<sup>2</sup>)
 +
 
 +
* F<sub>0</sub> : Initial fuel load (kg/m<sup>2</sup>)
 +
 
 +
* F<sub>n</sub> : Fuel load when reaction ends (kg/m<sup>2</sup>)
 +
 
 +
* F<sub>r</sub> : Fuel fraction burnt in the flaming (or reaction) zone (1)
 +
 
 +
* T<sub>f</sub> : fuel burn time (1-1/e  63% of the fuel burned) (s)
  
−
Variables on the fire mesh are located at the centers of the fire mesh cells of a 2D fire mesh. The position of the centers is stored in the arrays '''FXLONG''' and '''FXLAT'''.
+
* <var>&tau;</var><sub>R</sub> : Reaction time (s)
  
−
Fire variables important for visualization:
+
{| border="1"
 +
|+
 +
| scope="col" | Intensity
 +
| scope="col" | Definition
 +
| scope="col" | Implementation
 +
| scope="col" | Unit
 +
| scope="col" | Depends on rates spread R<sub>S</sub>
 +
| scope="col" | Depends on rates burn 1/T<sub>f</sub>
 +
|-
 +
| scope="row" | Reaction
 +
|I<sub>R</sub> = H<sub>F</sub>(F<sub>0</sub>-F<sub>n</sub>)/<var>&tau;</var><sub>R</sub>
 +
|I<sub>R</sub> = H<sub>F</sub>F<sub>0</sub>F<sub>r</sub>/(-T<sub>f</sub>ln(1-F<sub>r</sub>)
 +
| J/m<sup>2</sup>/s
 +
| no
 +
| yes
 +
|-
 +
| scope="row" | Byram
 +
| I = H<sub>F</sub>R<sub>S</sub>F
 +
| I = H<sub>F</sub>R<sub>S</sub>F<sub>0</sub>F<sub>r</sub>
 +
| J/m/s
 +
| yes
 +
| no
 +
|-
 +
| scope="row" | Fireline, new
 +
| J(<var>&Delta;</var>t<sup>2</sup>) &asymp; H<sub>F</sub>R<sub>S</sub>F<sub>0</sub> &int;<sub>a</sub><sup><var>&Delta;</var>t</sup><var>&tau;</var><sub>R</sub>/T<sub>f</sub> d<var>&tau;</var>
 +
| J = H<sub>F</sub>R<sub>S</sub>F<sub>0</sub>/2T<sub>f</sub>
 +
| J/m/s<sup>2</sup>
 +
| yes
 +
| yes
 +
|}
  
−
* '''FGRNHFX''' - heat flux from the ground fire. Higher value means a more intense fire, so it should be brighter.
+
==See also==
−
* '''FUEL_FRAC''' - the fuel fraction remaning. 1 means no fire yet, so one can display the original ground colors. 0 or small means burned, so the ground can have a different color, e.g. it could be blackened.
 
−
* '''FIRE_AREA''' - the proportion of the cell area that is on fire. 0 means no fire yet, 1 means all cell is on fire. It can be also used for ground colors.
 
  
−
See the [[Media:Users guide chap-wrf-fire.pdf|WRF-Fire user's guide]] for the description of all fire model variables.
+
* [[How to interpret WRF variables]]

Latest revision as of 08:22, 27 January 2024

Back to the WRF-SFIRE user guide.
This page is under construction. More contents should be coming soon.

Introduction

Firemesh.png

The fire mesh is 2D. Every atmosphere cell is divided into sr_x by sr_y fire cells. The submesh ratios sr_x and sr_y are specified in the file namelist.input. These values are not stored in the NetCDF files, but can be computed from the dimension sizes. 

sr_x=west_east_subgrid/west_east_stag = west_east_subgrid/(west_east + 1)
sr_y=south_north_subgrid/south_north_stag = south_north_subgrid/(south_north + 1)

For historical reasons, the fire grid dimensions in the output files are larger than what is actually used by the code internally. The extra space is at the end of the variables of size sr_x in x and sr_y in y. The following is an example of python code for correctly reading FGRNHFX from the file wrfout. 

from netCDF4 import Dataset
f=Dataset('wrfout','r')
sr_x=len(f.dimensions['west_east_subgrid'])//(len(f.dimensions['west_east'])+1)
sr_y=len(f.dimensions['south_north_subgrid'])//(len(f.dimensions['south_north'])+1)
fgrnhfx=f.variables['FGRNHFX'][:,:-sr_y,:-sr_x]

Variables on the fire mesh are located at the centers of the fire mesh cells of a 2D fire mesh.

List of variables

The authoritative information can be always found in the Registry.

Fire variables on fire grid
Variable name Description Unit
NFUEL_CAT fuel data
ZSF height of surface above sea level m
DZDXF surface gradient x 1
DZDYF surface gradient y 1


Fire variables on atm grid
Variable name Description Unit
RTHFRTEN temperature tendency K/s
RQVFRTEN humidity tendency


Diagnostics
Variable name Description Unit
AVG_FUEL_FRAC fuel remaining averaged to atmospheric grid 1
GRNHFX heat flux from ground fire W/m2
GRNQFX moisture flux from ground fire W/m2
CANHFX heat flux from crown fire W/m2
CANQFX moisture flux from crown fire W/m2
UAH wind at fire_wind_heigh m/s
VAH wind at fire_wind_heigh m/s


Sfire variables on fire grid
Variable name Description Unit
TIGN_G ignition time on ground s
LFN level function 1
FUEL_FRAC fuel remaining 1
FMC_G fuel moisture contents 1
FIRE_AREA fraction of cell area on fire 1
UF fire wind m/s
VF fire wind m/s
FGRNHFX heat flux from ground fire W/m2
FGRNQFX moisture flux from ground fire W/m2
FCANHFX heat flux from crown fire W/m2
FCANQFX moisture flux from crown fire W/m2


Diagnostics for the actual modeled fire
Variable name Description Unit
ROS rate of spread in the normal direction to the fireline m/s
FLINEINT fireline intensity W/m
FLINEINT2 alternative fireline intensity J/m/s2


Diagnostics for fire risk rating - independent on any actual fire going on
Variable name Description Unit
F_ROS0 base rate of spread in all directions m/s
F_ROSX X component of the spread vector driven by wind and slope m/s
F_ROSY Y component of the spread vector driven by wind and slope m/s
F_ROS max spread rate in any direction m/s
F_INT fire reaction intensity for risk rating, without fire J/m2/s
F_LINEINT Byram fireline intensity for risk rating, without fire J/m/s
F_LINEINT2 alternative fireline intensity for risk rating, without fire" J/m/s2


Constant data arrays
Variable name Description Unit
FXLONG longitude of midpoints of fire cells degrees
FXLAT latitude of midpoints of fire cells degrees
FUEL_TIME fuel
BBB fuel
PHISC fuel
PHIWC fuel
R_0 fuel
FGIP fuel
FZ0 fuel roughness height
FWH fuel fire wind height
ISCHAP fuel

Fire intensities

We introduce three fire intensities:

1. Reaction : The heat release rate per unit area of the front.

2. Byram : The heat produced per unit length of the fireline in unit time (J/m/s) in the so-called flaming zone behind the fireline.

3. Fireline, new : The amount of heat generated by the advancing fireline from the newly burning fuel only, in the unit of time.

Definitions and implementations of fire intensities include variables :

  • HF : Heat contents of the fuel (J/kg)
  • RS : Fire spread rate (m/s)
  • F : Fuel remaining at time t (kg/m2)
  • F0 : Initial fuel load (kg/m2)
  • Fn : Fuel load when reaction ends (kg/m2)
  • Fr : Fuel fraction burnt in the flaming (or reaction) zone (1)
  • Tf : fuel burn time (1-1/e 63% of the fuel burned) (s)
  • τR : Reaction time (s)
Intensity Definition Implementation Unit Depends on rates spread RS Depends on rates burn 1/Tf
Reaction IR = HF(F0-Fn)/τR IR = HFF0Fr/(-Tfln(1-Fr) J/m2/s no yes
Byram I = HFRSF I = HFRSF0Fr J/m/s yes no
Fireline, new J(Δt2) ≈ HFRSF0 ∫aΔtτR/Tf dτ J = HFRSF0/2Tf J/m/s2 yes yes

See also

Retrieved from "http://wiki.openwfm.org/mediawiki/index.php?title=SFIRE_variables&oldid=4536"