fire_index_haines
Computes the Haines fire index (aka: Lower Atmosphere Severity Index) from a sounding.
Available in version 6.4.0 and later.
Prototype
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/heat_stress.ncl" ; This library is automatically loaded ; from NCL V6.5.0 onward. ; No need for user to explicitly load. function fire_index_haines ( p [*] : numeric, t [*] : numeric, td [*] : numeric, opt [1] : logical ) return_val [3] : float or double
Arguments
pA one-dimensional array containing pressures (hPa).
ttd
One-dimensional arrays containing temperatures (degC) and dew point temperatures (degC), respectively.
optCurrently not used. Set to False.
Return value
The return value is a one-dimensional array of length 3. Elements 0, 1, 2 pertain to low, mid and high elevation.
Description
See: Haines paper and Haines Index.
The Haines_Index is technically the "Lower Atmosphere Severity Index" ==> LASI. It consists of two parts:
- (i) Stability Term: (Tpl - Tp2)
- (ii) Moisture Term : (Tp1 - TDp1)
The Haines Index can range between 2 and 6. The drier and more unstable the lower atmosphere the higher the index. Let ROS mean 'Risk Of Spread'. Then a Haines Severity Index of
- 2-3= a very low ROS
- 4= low ROS
- 5= moderate ROS, and
- 6= high ROS
The Climate Analysis Center uses the low-level Haines for the Coastal and Midlands zones, and the mid-level Haines for the Upstate zones.
The Haines Index is intended to be used all over the United States it is adaptable for three elevation regimes: Low Elevation, Middle Elevation and High Elevation. [See http://www.k3jae.com/wxHaines.php]
- Low Elevation is for fires occurring at or very near sea level.
- Middle Elevation is for fires burning in the 1000-3000 foot elevation range.
- High Elevation is intended for fires burning above 3000 feet elevation.
Examples
; PRESSURE (MB; hPa) p =(/ 1008.,1000.,950.,900.,850.,800.,750.,700.,650.,600., \ 550.,500.,450.,400.,350.,300.,250.,200., \ 175.,150.,125.,100., 80., 70., 60., 50., \ 40., 30., 25., 20. /) ; TEMPERATURE (C) t =(/ 29.3,28.1,23.5,20.9,18.4,15.9,13.1,10.1, 6.7, 3.1, \ -0.5,-4.5,-9.0,-14.8,-21.5,-29.7,-40.0,-52.4, \ -59.2,-66.5,-74.1,-78.5,-76.0,-71.6,-66.7,-61.3, \ -56.3,-51.7,-50.7,-47.5 /) ; MIXING RATIO (g/kg) q =(/ 20.38,19.03,16.14,13.71,11.56,9.80,8.33,6.75,6.06,5.07, \ 3.88, 3.29, 2.39, 1.70,1.00,0.60,0.20,0.00,0.00, \ 0.00, 0.00, 0.00, 0.00,0.00,0.00,0.00,0.00,0.00, \ 0.00, 0.00 /) ; Haines Index (HI) requires dew point temperature so it must be derived ; Change units to those requires by 'relhum' function q = q/1000. ; (kg/kg) t = t+273.15 ; K rh = relhum(t, q, p*100) ; p*100 => mb => Pa td = dewtemp_trh(t,rh) ; dew pt temperature [K] td = td-273.15 ; [C] ; return to units required by 'fire_haines_index' t = t -273.15 ; [C] ;print(p+" "+t+" "+q+" "+rh+" "+td) ; assign metadata ; units [ @ ] p@units = "hPa" t@units = "degC" ; name dimensions [ ! ] p!0 = "p" t!0 = "p" td!0 = "p" ; assign coordinate values [ & ] to named dimensions p&p = p t&p = p td&p = p ; look at data ; printVarSummary(p) ; printVarSummary(t) ; printVarSummary(td) printMinMax(t,0) HI = fire_index_haines(p, t, td, False) print(HI)A sample output
Variable: HI Type: float Total Size: 12 bytes 3 values Number of Dimensions: 1 Dimensions and sizes: [elevation | 3] Coordinates: Number Of Attributes: 5 long_name : Haines Index reference : http://www.nwas.org/digest/papers/1988/Vol13-Issue2-May1988/Pg23-Haines.pdf info : Lower Atmosphere Severity Index (LASI) elevation : ( Low, Medium, High ) details : http://www.erh.noaa.gov/cae/haines.htm (0) 3 (1) 3 (2) 2