Merge pull request #20 from gregory-halverson-jpl/main

sensitivity presentation figures

Author: gregory-halverson-jpl

PMJPL/ECOv002-cal-val-PM-JPL-inputs.csv

@@ -228,6 +228,8 @@ def PMJPL(
     if Ta_C is None:
         raise ValueError("air temperature (Ta_C) not given")
 
+    check_distribution(Ta_C, "Ta_C")
+
     if Tmin_C is None and geometry is not None and time_UTC is not None:
         Tmin_K = GEOS5FP_connection.Tmin_K(
             time_UTC=time_UTC,
@@ -240,6 +242,8 @@ def PMJPL(
     if Tmin_C is None:
         raise ValueError("minimum temperature (Tmin_C) not given")
 
+    check_distribution(Tmin_C, "Tmin_C")
+
     if RH is None and geometry is not None and time_UTC is not None:
         RH = GEOS5FP_connection.RH(
             time_UTC=time_UTC,
@@ -250,9 +254,13 @@ def PMJPL(
     if RH is None:
         raise ValueError("relative humidity (RH) not given")
 
+    check_distribution(RH, "RH")
+
     if elevation_km is None and geometry is not None:
         elevation_km = NASADEM.elevation_km(geometry=geometry)
 
+    check_distribution(elevation_km, "elevation_km")
+
     elevation_m = elevation_km * 1000.0
 
     if IGBP is None and geometry is not None:
@@ -264,6 +272,8 @@ def PMJPL(
             raise ValueError(f"invalid geometry type for IGBP retrieval: {type(geometry)}")
 
         IGBP = load_MCD12C1_IGBP(geometry=IGBP_geometry)
+    
+    check_distribution(np.float32(IGBP), "IGBP")
 
     if regenerate_net_radiation or (Rn_Wm2 is None and albedo is not None and ST_C is not None and emissivity is not None):
         if SWin_Wm2 is None and geometry is not None and time_UTC is not None:
@@ -312,6 +322,8 @@ def PMJPL(
         else:
             raise ValueError("net radiation (Rn_Wm2) not given and cannot be calculated")
 
+    check_distribution(Rn_Wm2, "Rn_Wm2")
+
     if G_Wm2 is None and Rn_Wm2 is not None and ST_C is not None and NDVI is not None and albedo is not None:
         G_Wm2 = calculate_SEBAL_soil_heat_flux(
             Rn=Rn_Wm2,

PMJPL/ECOv002-cal-val-PM-JPL-outputs.csv

@@ -127,6 +127,11 @@ def process_PMJPL_table(
     else:
         elevation_km = None
 
+    if "IGBP" in input_df:
+        IGBP = np.array(input_df.IGBP).astype(np.int8)
+    else:
+        IGBP = None
+
     # --- Handle geometry and time columns ---
     import pandas as pd
     from rasters import MultiPoint, WGS84
@@ -187,6 +192,7 @@ def parse_geom(s):
         SWin_Wm2=SWin_Wm2,
         albedo=albedo,
         Tmin_C=Tmin_C,
+        IGBP=IGBP,
         elevation_km=elevation_km,
         time_UTC=time_UTC,
         upscale_to_daylight=upscale_to_daylight,