Implemented conversion of Burned Area from m2 to hectares before the …

…fuel load is calculated in fuelload.py. This addresses issue #10.
ecmwf-projects · Apr 23, 2021 · f9306f2 · f9306f2
1 parent 9acf34f
commit f9306f2
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/utils/fuelload.py b/src/utils/fuelload.py
@@ -29,6 +29,11 @@ def fuelload(path_agb: str, path_ba: str, path_fl: str = None) -> xarray.Dataset
     agb_data.values[agb_data.values == AGB_THRESH] = np.nan
     ba_data.values[ba_data.values < BA_THRESH] = np.nan
 
+    # AGB units are Mg/h
+    # BA units are m2, therefore we convert BA to hectares.
+    ba_data = ba_data * 0.0001
+
+    # Now that units are consistent we calculate LOAD = AGB * BA
     fuel_load_dataset = agb_data * ba_data
 
     if path_fl is not None: