How to diagnose fuel properties in WRF-SFIRE - Revision history

Jmandel: /* Step-by-step instructions */ update for standalone

2022-08-08T20:09:44Z

Step-by-step instructions: update for standalone

← Older revision		Revision as of 20:09, 8 August 2022
Line 17:		Line 17:

	==Step-by-step instructions==		==Step-by-step instructions==
	* Run WRF-SFIRE, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.		===Acquire and visualize ROS===
			* Run WRF-SFIRE, you can kill it right after the the first time step, or follow the instructions for '''fire_ros.exe''' in '''standalone'''
			* There will be file '''fuels.m''' created in the current directory.
	* '''git clone https://github.com/openwfm/wrf-fire-matlab''' and copy the file '''fuels.m''' to the directory created.		* '''git clone https://github.com/openwfm/wrf-fire-matlab''' and copy the file '''fuels.m''' to the directory created.
	* '''cd''' in MATLAB to the directory and type '''startup'''		* '''cd''' in MATLAB to the directory and type '''startup'''
	* ~~Type~~ '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.		* '''cd''' in MATLAB to the directory with the file '''fuels.m''' and type '''fuels'''.
			* This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.
	* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.		* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.
	* ~~Type~~ '''edit fire_ros''' to see the fire spread ~~rate~~ calculation recreated in Matlab.
	* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4) '''Note - at the moment, the checks in plot_fuel are fine, but check_ros is not. '''Opened issue https://github.com/openwfm/wrf-fire-matlab/issues/5 and upstream in WRF-SFIRE.'''		===Compare with Rothermel method===
			* '''edit fire_ros''' to see the Rothermel fire rate of spread calculation recreated in Matlab.
			* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4)
			* '''Note - at the moment, the checks in plot_fuel are fine, but check_ros is not. '''Opened issue https://github.com/openwfm/wrf-fire-matlab/issues/5 and upstream in WRF-SFIRE.'''

	==Diagnostics provided==		==Diagnostics provided==

Jmandel: /* Step-by-step instructions */

2021-10-04T03:19:42Z

Step-by-step instructions

← Older revision		Revision as of 03:19, 4 October 2021
Line 18:		Line 18:
	==Step-by-step instructions==		==Step-by-step instructions==
	* Run WRF-SFIRE, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.		* Run WRF-SFIRE, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.
	* ~~Clone~~ '''https://github.com/openwfm/wrf-fire-matlab''' and copy the file '''fuels.m''' to the directory created.		* '''git clone https://github.com/openwfm/wrf-fire-matlab''' and copy the file '''fuels.m''' to the directory created.
	* '''cd''' in MATLAB to the directory and type '''startup'''		* '''cd''' in MATLAB to the directory and type '''startup'''
	* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.		* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.

Jmandel at 03:18, 4 October 2021

2021-10-04T03:18:42Z

← Older revision		Revision as of 03:18, 4 October 2021
Line 1:		Line 1:
	{{users guide}}		{{users guide}}
	The fuel properties in [[WRF-~~Fire~~]] are given in file [https://github.com/openwfm/WRF-SFIRE/blob/master/test/em_fire/namelist.fire namelist.fire]. These serve to create coefficients in the Rothermel's formula. To aid diagnostics it is useful to graph the resulting fire spread rate as a function of wind and slope. The graphs are available in metric units, as well as in English units for direct comparison with the published spread rate graphs for the Scott-Burgan fuel categories.		The fuel properties in [[WRF-SFIRE]] are given in file [https://github.com/openwfm/WRF-SFIRE/blob/master/test/em_fire/namelist.fire namelist.fire]. These serve to create coefficients in the Rothermel's formula. To aid diagnostics it is useful to graph the resulting fire spread rate as a function of wind and slope. The graphs are available in metric units, as well as in English units for direct comparison with the published spread rate graphs for the Scott-Burgan fuel categories.

	==Prerequisites==		==Prerequisites==
Line 33:		Line 33:

	==Works with==		==Works with==
			''To be updated''
	* Version {{WRF-Fire-commit\|334bf9417b5397ed3692b8f469a4914ed433e051\|16 Sep 2010}}		* Version {{WRF-Fire-commit\|334bf9417b5397ed3692b8f469a4914ed433e051\|16 Sep 2010}}
	* Version {{WRF-Fire-commit\|bce8561325b4ced49303e941931d78f63b3c843b\|29 Aug 2012}} and hopefully later		* Version {{WRF-Fire-commit\|bce8561325b4ced49303e941931d78f63b3c843b\|29 Aug 2012}} and hopefully later

Jmandel: Jmandel moved page How to diagnose fuel properties in WRF-Fire to How to diagnose fuel properties in WRF-SFIRE

2021-10-04T03:15:42Z

Jmandel moved page How to diagnose fuel properties in WRF-Fire to How to diagnose fuel properties in WRF-SFIRE

← Older revision	Revision as of 03:15, 4 October 2021
(No difference)

Jmandel: /* Step-by-step instructions */

2021-10-03T22:40:56Z

Step-by-step instructions

← Older revision		Revision as of 22:40, 3 October 2021
Line 23:		Line 23:
	* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.		* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.
	* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.		* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.
	* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4) '''Note - the checks in plot_fuel are fine, but check_ros is not. '''Opened issue https://github.com/openwfm/wrf-fire-matlab/issues/5 and upstream in WRF-SFIRE.'''		* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4) '''Note - at the moment, the checks in plot_fuel are fine, but check_ros is not. '''Opened issue https://github.com/openwfm/wrf-fire-matlab/issues/5 and upstream in WRF-SFIRE.'''

	==Diagnostics provided==		==Diagnostics provided==

Jmandel: /* Step-by-step instructions */

2021-10-03T22:26:17Z

Step-by-step instructions

← Older revision		Revision as of 22:26, 3 October 2021
Line 18:		Line 18:
	==Step-by-step instructions==		==Step-by-step instructions==
	* Run WRF-SFIRE, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.		* Run WRF-SFIRE, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.
	* Clone '''[https://github.com/openwfm/wrf-fire-matlab]''' and copy the file '''fuels.m''' to the directory created.		* Clone '''https://github.com/openwfm/wrf-fire-matlab''' and copy the file '''fuels.m''' to the directory created.
	* '''cd''' in MATLAB to the directory and type '''startup'''		* '''cd''' in MATLAB to the directory and type '''startup'''
	* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.		* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.

Jmandel: /* Step-by-step instructions */ opened issue

2021-10-03T03:47:48Z

Step-by-step instructions: opened issue

← Older revision		Revision as of 03:47, 3 October 2021
Line 23:		Line 23:
	* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.		* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.
	* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.		* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.
	* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4) '''Note - the checks in plot_fuel are fine, but check_ros is not.		* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4) '''Note - the checks in plot_fuel are fine, but check_ros is not. '''Opened issue https://github.com/openwfm/wrf-fire-matlab/issues/5 and upstream in WRF-SFIRE.'''

	==Diagnostics provided==		==Diagnostics provided==

Jmandel: /* Step-by-step instructions */

2021-10-03T03:35:49Z

Step-by-step instructions

← Older revision		Revision as of 03:35, 3 October 2021
Line 17:		Line 17:

	==Step-by-step instructions==		==Step-by-step instructions==
	* Run WRF-~~Fire~~, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.		* Run WRF-SFIRE, you can kill it right after the the first time step. There will be file '''fuels.m''' created in the current directory.
	* ~~Start Matlab in~~ the ~~directory~~ '''~~test/em_fire~~''' to ~~set up~~ the ~~search path properly~~.		* Clone '''[https://github.com/openwfm/wrf-fire-matlab]''' and copy the file '''fuels.m''' to the directory created.
	* ~~In Matlab, navigate~~ to the directory ~~with the~~ '''~~fuels.m~~''' ~~file~~		* '''cd''' in MATLAB to the directory and type '''startup'''
	* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.		* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.
	* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.		* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.
	* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.		* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.
	* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4)		* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4) '''Note - the checks in plot_fuel are fine, but check_ros is not.

	==Diagnostics provided==		==Diagnostics provided==

Jmandel: /* Step-by-step instructions */

2021-10-03T03:24:29Z

Step-by-step instructions

← Older revision		Revision as of 03:24, 3 October 2021
Line 21:		Line 21:
	* In Matlab, navigate to the directory with the '''fuels.m''' file		* In Matlab, navigate to the directory with the '''fuels.m''' file
	* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.		* Type '''fuels'''. This will create variable '''fuel'''. Say you want to check fuel 3. Type '''fuel(3)''' to see what is there.
	* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-~~Fire~~ and stored in the '''fuel''' variable.		* Type '''plot_fuel(fuel(3))''' to create the graph below in metric units, and '''plot_fuel(fuel(3),'sb')''' for the same units as in Scott and Burgan (2005). Type '''help plot_fuel''' for further information. Note: The graph data were created in WRF-SFIRE and stored in the '''fuel''' variable.
	* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.		* Type '''edit fire_ros''' to see the fire spread rate calculation recreated in Matlab.
	* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-~~Fire~~, it should return rounding error only (less than 1e-4)		* Type '''big(check_ros(fuel))''' to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file '''fuels.m'''. It should return rounding error only (less than 1e-4)

	==Diagnostics provided==		==Diagnostics provided==

Jmandel: /* Location of the codes */ update

2021-10-03T03:21:18Z

Location of the codes: update

← Older revision		Revision as of 03:21, 3 October 2021
Line 9:		Line 9:
	==Location of the codes==		==Location of the codes==
	* The fuel properties of every category of the fuel are defined in file [https://github.com/openwfm/WRF-SFIRE/blob/master/test/em_fire/namelist.fire '''namelist.fire'''], which needs to be in the current directory.		* The fuel properties of every category of the fuel are defined in file [https://github.com/openwfm/WRF-SFIRE/blob/master/test/em_fire/namelist.fire '''namelist.fire'''], which needs to be in the current directory.
	* The fire rate of spread computation is done in file '''[~~http~~://github.com/~~jbeezley~~/~~wrf~~-~~fire~~/blob/~~334bf9417b5397ed3692b8f469a4914ed433e051/wrfv2_fire~~/phys/module_fr_sfire_phys.F ~~WRFV3/phys/~~module_fr_sfire_phys.F]''' split between two subroutines. First, the coefficients at every point are precomputed from fuel properties in ~~[http://github.com/jbeezley/wrf-fire/blob/334bf9417b5397ed3692b8f469a4914ed433e051/wrfv2_fire/phys/module_fr_sfire_phys.F#L413~~ subroutine '''set_fire_params'''] at initialization. In every time step, the spread rate at a point is then computed from those coefficients in ~~[http://github.com/jbeezley/wrf-fire/blob/334bf9417b5397ed3692b8f469a4914ed433e051/wrfv2_fire/phys/module_fr_sfire_phys.F#L672~~ subroutine '''fire_ros'''].		* The fire rate of spread computation is done in file '''[https://github.com/openwfm/WRF-SFIRE/blob/master/phys/module_fr_sfire_phys.F module_fr_sfire_phys.F]''', split between two subroutines. First, the coefficients at every point are precomputed from fuel properties in subroutine '''set_fire_params''' at initialization. In every time step, the spread rate at a point is then computed from those coefficients in subroutine '''fire_ros'''.
	* The computation has been duplicated in Matlab for diagnostic purposes in file '''[~~http~~://github.com/~~jbeezley~~/wrf-fire/blob/~~334bf9417b5397ed3692b8f469a4914ed433e051/other/Matlab~~/vis3d/fire_ros.m ~~other/Matlab/vis3d/~~fire_ros.m]'''		* The computation has been duplicated in Matlab for diagnostic purposes in file '''[https://github.com/openwfm/wrf-fire-matlab/blob/master/vis3d/fire_ros.m fire_ros.m]'''
	* The plot routine to create graphs like shown here is '''[~~http~~://github.com/~~jbeezley~~/wrf-fire/blob/~~334bf9417b5397ed3692b8f469a4914ed433e051/other/Matlab~~/vis3d/plot_fuel.m ~~other/Matlab/vis3d/~~plot_fuel.m]'''		* The plot routine to create graphs like shown here is '''[https://github.com/openwfm/wrf-fire-matlab/blob/master/vis3d/plot_fuel.m plot_fuel.m]'''

	==Heat flux diagnostics==		==Heat flux diagnostics==