Balbi model

[janmandel]

Add Balbi fire model as an alternative to Rothermel model.

[janmandel]

Run by @Fergui , comparing master 681dc5a and balbi 9e70428 with ibeh=1, hill test.

Differences in diagnostic quantities, F_LINEINT, F_LINEINT2, F_ROS, F_ROS0, F_ROSX, F_ROSY.
The fire progress itself, LFN, FIRE_AREA and the heat fluxes, GRNHFX FGRNHFX are the same.

[janmandel]

From 9e70428
ibeh in namelist fire selects spread formula 0 = MacArthur, 1 = Rothermel, 2 = Balbi
If if ibeh=2, call fire_ros_balbi otherwise fire_ros_cawfe (Rothermel).

Subroutine fire_ros_balbi receives the number of fuel category, and windspeed and slope, both projected on the propagation direction, if fire_advection=0. fire_advection not 0 was something else but is not supported by the current code.
The fuel category is currently ignored and it uses constants for the fuel properties instead.

[janmandel]

Do do:

In #3 standalone:

• add balbi main program to call fire_ros_balbi, read namelist.fire, print result. Done by call call init_fuel_cats in fire_ros_main.F 6e0e54a
• read and pass on also air temperature, surface air density, fire front width.
• merge here 691d84e

In module_fr_sfire_phys:

• the fuel properties in the grid cell are already available indexed by nfuel
• add any new tables such as live fuel (from Anderson 1982), add adjusted fuel loading if needed - do not change existing fuel table
• pass through surface air temperature, surface air density
• Add fire max (33m/50m) front width as a fuel property.

Above module_fr_sfire_phys:

• bring in surface air temperature and surface air density.
• process the level set function to compute fireline curvature, then assign fire front width as const/curvature

Future:

• treat wind vector, slope vector, and normal vector in some proper way.
• Scott-Burgan categories

Inside fire_ros_balbi:

• translate the inputs into what Balbi formula needs
• upgrade the algorithm to current version

[janmandel]

@jstu42 @adamk0 @Fergui
I finished the merge and advanced master, develop-45, and develop-3 to the same commit
d9b313e. Please let me know if you see any problems.
Please git pull and continue using develop-45 to continue in Balbi #45. Subroutine fire_ros_balbi is not getting the values of tempf, rhof, and ffwidth from WRF yet (previously t2f, psfcf, ffw) but you can keep using the values handed from standalone fire_ros.exe for testing, with the new names of course.
This also brings standalone #3 to v4.4 #39

[janmandel]

@jstu42 @adamk0 @Jeremy-Benik @Fergui

Feeding fire_ros_balbi from WRF arrays 62ce9e0 when ibeh=2 in namelist.fire. It also creates fuels.m but using fuels.m with ROS diagnostics and visualization needs to be tested.

fire_ros_balbi iterations result in negative ROS (and therefore do not converge) when the wind speed in the propagation direction is negative.

The code is in branch develop-45

[Fergui]

Branch develop-45 is in develop-60 at 4993bf1
Closing as duplicate issue continued in #60