# Heathland

The Anderson (2015) rate of spread model was developed for predicting fire spread rates in Heathland.

Heathland

### Fuel inputs

• temp
• rel_hum
• WF
• fuel_height
• R0
• delta
• Temperature (degrees C)
• Relative humidity (%)
• Wind reduction factor (0 – 1)
• Average vegetation height (m)
• Rate of fire spread for zero wind (m/min)
• Solar radiation factor (0 or 1)

### Code

```// Heathland model - Anderson et al. (2015)

// -------------------------------------------
// Model parameters
// These must be defined below, or included as a user-defined layer
//
// 1. Temperature, 'temp' (input)
// 2. Relative humidity, 'rel_hum' (input)
// 3. Wind reduction factor, 'WF' (input), 0.67 for heath-shrublands and 0.35 for woodlands.
// 4. Average fuel height, 'fuel_height' (input)
// 5. Rate of spread for zero wind, 'R0' (input), recommended as 5 m/min.
// 6. Solar radiation factor, 'delta' (input), {0,1}, delta = 1 for sunny days from 12:00 - 17:00 from October to March (high solar radiation) and 0 otherwise.
// -------------------------------------------

// Backing and flanking coefficients compared to head fire ROS
REAL wind_speed = length(wind_vector);

// Calculate length-to-breadth ratio (LBR) for grassland
REAL LBR = 1.0;
if (wind_speed < 5){
LBR = 1.0;
} else {
LBR = 1.1*pow(wind_speed, 0.464);
}

// Determine coefficient for flank rank of spread
REAL cc = sqrt(1.0-pow(LBR, -2.0));
REAL cb = (1.0-cc)/(1.0+cc);
REAL a_LBR = 0.5*(cb+1.0);
REAL cf = a_LBR/LBR;

// Determine shape parameters
REAL f = 0.5*(1.0+cb);
REAL g = 0.5*(1.0-cb);
REAL h = cf;

// Now calculate a speed coefficient using normal flow formula
REAL speed_fraction = (g*wdot+sqrt(h*h+(f*f-h*h)*wdot*wdot));

// Initialising the solar radiation variable
REAL delta;
if (hour > 11 && hour < 17){
delta = 1;}
else {
delta = 0;}

// Calculating moisture content
REAL MC = 4.37 + 0.161*rel_hum - 0.1*(temp - 25) - delta * 0.027 * rel_hum;