BESS-JPL 1.13.0__py3-none-any.whl → 1.15.0__py3-none-any.whl

BESS_JPL/BESS_JPL.py

@@ -6,7 +6,7 @@ import numpy as np
 import rasters as rt
 from rasters import Raster, RasterGeometry
 
-from check_distribution import check_distribution
+from .check_distribution import check_distribution
 
 from sun_angles import calculate_SZA_from_DOY_and_hour
 from solar_apparent_time import solar_day_of_year_for_area, solar_hour_of_day_for_area

BESS_JPL/canopy_shortwave_radiation.py

@@ -1,6 +1,6 @@
 import numpy as np
 
-from check_distribution import check_distribution
+from .check_distribution import check_distribution
 
 
 def canopy_shortwave_radiation(

BESS_JPL/check_distribution.py

@@ -0,0 +1,118 @@
+"""
+Replacement for check_distribution package to avoid pycksum dependency.
+
+This module provides data validation functionality that was previously
+provided by the check-distribution package, which depends on pycksum
+that doesn't compile on Python 3.12.
+"""
+
+import logging
+from typing import Union, Optional
+from datetime import datetime
+
+# Try to import numpy, but make it optional
+try:
+    import numpy as np
+    HAS_NUMPY = True
+except ImportError:
+    np = None
+    HAS_NUMPY = False
+
+
+logger = logging.getLogger(__name__)
+
+
+def check_distribution(
+    data: Union["np.ndarray", float, int, None], 
+    name: str, 
+    time_UTC: Optional[datetime] = None
+) -> None:
+    """
+    Check and validate data distribution, logging warnings for invalid values.
+    
+    This function validates input data and logs information about the distribution,
+    including warnings for NaN values, infinite values, and basic statistics.
+    
+    Args:
+        data: The data to validate (numpy array, scalar, or None)
+        name: Name of the variable for logging purposes
+        time_UTC: Optional timestamp for logging context
+        
+    Returns:
+        None
+        
+    Raises:
+        Nothing - this function only logs warnings and information
+    """
+    
+    # Create log prefix with time if provided
+    time_prefix = f"[{time_UTC}] " if time_UTC else ""
+    
+    if data is None:
+        logger.warning(f"{time_prefix}{name}: data is None")
+        return
+    
+    # Handle numpy arrays or array-like objects
+    if hasattr(data, 'shape') and hasattr(data, 'dtype'):
+        # This is likely a numpy array or array-like object
+        
+        if hasattr(data, 'size') and data.size == 0:
+            logger.warning(f"{time_prefix}{name}: array is empty")
+            return
+        
+        if HAS_NUMPY:
+            # Use numpy functions if available
+            # Check for NaN values
+            if np.any(np.isnan(data)):
+                nan_count = np.sum(np.isnan(data))
+                nan_percentage = (nan_count / data.size) * 100
+                logger.warning(f"{time_prefix}{name}: contains {nan_count} NaN values ({nan_percentage:.1f}%)")
+            
+            # Check for infinite values  
+            if np.any(np.isinf(data)):
+                inf_count = np.sum(np.isinf(data))
+                inf_percentage = (inf_count / data.size) * 100
+                logger.warning(f"{time_prefix}{name}: contains {inf_count} infinite values ({inf_percentage:.1f}%)")
+                
+            # Calculate statistics for finite values only
+            finite_data = data[np.isfinite(data)]
+            if len(finite_data) > 0:
+                min_val = np.min(finite_data)
+                max_val = np.max(finite_data)
+                mean_val = np.mean(finite_data)
+                std_val = np.std(finite_data)
+                
+                logger.debug(f"{time_prefix}{name}: shape={data.shape}, "
+                            f"min={min_val:.6f}, max={max_val:.6f}, "
+                            f"mean={mean_val:.6f}, std={std_val:.6f}")
+            else:
+                logger.warning(f"{time_prefix}{name}: no finite values found")
+        else:
+            # Fallback when numpy is not available
+            logger.debug(f"{time_prefix}{name}: array-like data validated (numpy not available for detailed analysis)")
+            
+    else:
+        # Handle scalar values
+        try:
+            # Convert to float to check for NaN/inf
+            float_val = float(data)
+            
+            if HAS_NUMPY:
+                if np.isnan(float_val):
+                    logger.warning(f"{time_prefix}{name}: scalar value is NaN")
+                elif np.isinf(float_val):
+                    logger.warning(f"{time_prefix}{name}: scalar value is infinite")
+                else:
+                    logger.debug(f"{time_prefix}{name}: scalar value = {float_val:.6f}")
+            else:
+                # Basic validation without numpy
+                import math
+                if math.isnan(float_val):
+                    logger.warning(f"{time_prefix}{name}: scalar value is NaN")
+                elif math.isinf(float_val):
+                    logger.warning(f"{time_prefix}{name}: scalar value is infinite")
+                else:
+                    logger.debug(f"{time_prefix}{name}: scalar value = {float_val:.6f}")
+                
+        except (ValueError, TypeError) as e:
+            logger.warning(f"{time_prefix}{name}: cannot convert to float for validation: {e}")

BESS_JPL/model.py

@@ -6,7 +6,7 @@ import numpy as np
 import rasters as rt
 from rasters import Raster, RasterGeometry
 
-from check_distribution import check_distribution
+from .check_distribution import check_distribution
 
 from sun_angles import calculate_SZA_from_DOY_and_hour
 from solar_apparent_time import solar_day_of_year_for_area, solar_hour_of_day_for_area
@@ -16,7 +16,7 @@ from gedi_canopy_height import load_canopy_height, GEDI_DOWNLOAD_DIRECTORY
 from FLiESANN import FLiESANN
 from GEOS5FP import GEOS5FP
 from MODISCI import MODISCI
-from NASADEM import NASADEM
+from NASADEM import NASADEMConnection
 
 from .constants import *
 from .C3_photosynthesis import *
@@ -87,6 +87,8 @@ def BESS_JPL(
         peakVCmax_C3: np.ndarray = None,  # peak maximum carboxylation rate for C3 plants
         peakVCmax_C4: np.ndarray = None,  # peak maximum carboxylation rate for C4 plants
         CI: Union[Raster, np.ndarray] = None,
+        MODISCI_connection: MODISCI = None,
+        NASADEM_connection: NASADEMConnection = None,
         resampling: str = RESAMPLING,
         GEDI_download_directory: str = GEDI_DOWNLOAD_DIRECTORY):  # clumping index
     if geometry is None and isinstance(ST_C, Raster):
@@ -103,7 +105,11 @@ def BESS_JPL(
         raise ValueError("no time given between time_UTC, day_of_year, and hour_of_day")
 
     if elevation_km is None and geometry is not None:
-        elevation_km = NASADEM.elevation_km(geometry=geometry)
+        if NASADEM_connection is None:
+            from NASADEM import NASADEMConnection
+            NASADEM_connection = NASADEMConnection()
+
+        elevation_km = NASADEM_connection.elevation_km(geometry=geometry)
 
     # load air temperature in Celsius if not provided
     if Ta_C is None:
@@ -269,9 +275,11 @@ def BESS_JPL(
     if SZA is None:
         SZA = calculate_SZA_from_DOY_and_hour(geometry.lat, geometry.lon, day_of_year, hour_of_day)
 
+    if MODISCI_connection is None:
+        MODISCI_connection = MODISCI()
+
     if CI is None and geometry is not None:
-        modisci = MODISCI()
-        CI = modisci.CI(geometry=geometry, resampling=resampling)
+        CI = MODISCI_connection.CI(geometry=geometry, resampling=resampling)
 
     # canopy height defaults to zero
     canopy_height_meters = np.where(np.isnan(canopy_height_meters), 0, canopy_height_meters)

BESS_JPL/version.txt

@@ -1 +1 @@
-1.13.0
+1.15.0

{bess_jpl-1.13.0.dist-info → bess_jpl-1.15.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: BESS-JPL
-Version: 1.13.0
+Version: 1.15.0
 Summary: Breathing Earth System Simulator (BESS) Gross Primary Production (GPP) and Evapotranspiration (ET) Model Python
 Author-email: Gregory Halverson <gregory.h.halverson@jpl.nasa.gov>
 Project-URL: Homepage, https://github.com/JPL-Evapotranspiration-Algorithms/BESS-JPL
@@ -9,13 +9,12 @@ Classifier: Operating System :: OS Independent
 Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: check-distribution
 Requires-Dist: FLiESANN>=1.5.0
 Requires-Dist: gedi-canopy-height>=1.1.0
 Requires-Dist: GEOS5FP>=1.1.1
 Requires-Dist: koppengeiger>=1.0.4
 Requires-Dist: MODISCI>=1.3.0
-Requires-Dist: NASADEM
+Requires-Dist: NASADEM>=1.3.0
 Requires-Dist: numpy
 Requires-Dist: rasters
 Requires-Dist: solar-apparent-time>=1.3.2

{bess_jpl-1.13.0.dist-info → bess_jpl-1.15.0.dist-info}/RECORD

@@ -1,4 +1,4 @@
-BESS_JPL/BESS_JPL.py,sha256=CTby93_HwhONe2NOlt-8DuqDH7jL9mAc3euuzbouANg,1365
+BESS_JPL/BESS_JPL.py,sha256=-OieL_tZ04yhlv0z3_s9XmrEuAPueIvGON1d1Aj-F2M,1366
 BESS_JPL/C3_photosynthesis.py,sha256=-ormDycaWttCOoYJdOoOV_zlM6n0nxkgXxj7fyH_mIs,2760
 BESS_JPL/C4_fraction.jpeg,sha256=ECaEYWA8MnNd5z0Yq3beeUCR93KtTDxM8MrH-YYon3Y,2746
 BESS_JPL/C4_fraction.tif,sha256=DFS-J08JCjnf0hi_T_dFridAegVTiU3FGyMBjS9yc1w,25446
@@ -19,10 +19,11 @@ BESS_JPL/ball_berry_slope_C4.tif,sha256=UedoUOnOjBxQ9pAaf5nXgDa8wvKOHVBD3sVToV12
 BESS_JPL/calculate_VCmax.py,sha256=DDZk7v-ME_UtupiPPvTsGqI55-3JUv6i_jDYQozcqzQ,1959
 BESS_JPL/canopy_energy_balance.py,sha256=k7uEm15iv7cG11grvb3hZSgZ45nLbPkeniWJFk7K9rE,6452
 BESS_JPL/canopy_longwave_radiation.py,sha256=TH5YaxyVgVP58moZ0aDGkMukJfcPQQgiYbuCHO6TxPc,3407
-BESS_JPL/canopy_shortwave_radiation.py,sha256=TzDGFFBVHcZ3s5v-R9pQYeW9PPvYDQ168VdoqtCXh4U,9767
+BESS_JPL/canopy_shortwave_radiation.py,sha256=jkMK3iYq1WsebfQWLrESd5ZwZ31vW-EKRpyn92_IxJU,9768
 BESS_JPL/carbon_uptake_efficiency.jpeg,sha256=wz2cLk_4-DRdF71xR2c-cx1MsE0CH-EDWjs2pDJcRF0,704540
 BESS_JPL/carbon_uptake_efficiency.tif,sha256=Nnc8h2ImU6_sip0WzDh3q0BRoBe-CfCbwzDtYo5BSog,819269
 BESS_JPL/carbon_water_fluxes.py,sha256=qkCq_B-_0F9XQz9IDWQBaF-OpiDKIuZmhS7vyYKtYkI,13746
+BESS_JPL/check_distribution.py,sha256=T6lYN_vBCb1duoJ1BZ0XLRXttEI81S4BwK0DKb8KDeI,4508
 BESS_JPL/constants.py,sha256=tOYubq17dnLVKhtTI0CJUQAEeROE4KeSBpOVBqZprUU,131
 BESS_JPL/interpolate_C3_C4.py,sha256=adIyk03Yt6SJ_jFCUm6ivK9QpZYw10C3wJKpx1O56Pk,394
 BESS_JPL/kn.jpeg,sha256=wsI40jabT9_iPpQmOxd3GZzfAN7oZ2DubiARmULFUS0,854097
@@ -38,16 +39,16 @@ BESS_JPL/load_kn.py,sha256=_IE7w_VfahubfqjvgE6F6HCO0SizcpMu104_7ZXOiaw,345
 BESS_JPL/load_peakVCmax_C3.py,sha256=ISaCyaN81F3CyYbYwUqvUfQ_vjn4_jCpt1_x_b4fY5o,365
 BESS_JPL/load_peakVCmax_C4.py,sha256=8vtjEiJKCqSK04DfTs8Zgn-rWwGucmeeLMzblLT1U-Y,365
 BESS_JPL/meteorology.py,sha256=xIPQnIk9wGJPy8esJyxn35d6c_7UyiNuWTJXDnuSwf8,6995
-BESS_JPL/model.py,sha256=T4vjlFsX1BXA_784DJ5RHklOHDGZJ_IxomwT2S3V9wI,22490
+BESS_JPL/model.py,sha256=ioCuQXPO3tbhF-PXlCm1I7oo1Apno2FRjAANm0ShLG8,22811
 BESS_JPL/peakVCmax_C3.jpeg,sha256=nVvwLx8JyRhtm5lWRW93HLz0mInshEdOCJ1tAdcFqa8,1006133
 BESS_JPL/peakVCmax_C3.tif,sha256=ax6wCOPc_ovgJJl9YRjPWNY13OCGlzs2djXqwob_U7A,1583598
 BESS_JPL/peakVCmax_C4.jpeg,sha256=s7dhpcD573KW9Se4mejDeSzbSHqPtQY2EL6KJKt7ZIo,961497
 BESS_JPL/peakVCmax_C4.tif,sha256=EST4_yy-HHYmazIv--RePL_bhLejMWql6hoV9EE38ok,1556928
 BESS_JPL/process_BESS_table.py,sha256=TYsDEZNELGdUc0xkB1MRQIDkVonNL5d49q7Zl1iaJvo,1472
 BESS_JPL/soil_energy_balance.py,sha256=5TKDVkKYJ8jhjuiwbRva_03fi_0gTM0IAKbSm4WhWnY,944
-BESS_JPL/version.txt,sha256=dIz5gc5PjwOt3_4pgn-AK1ryzAJ9xRlX-XHBsjuLTL8,6
-bess_jpl-1.13.0.dist-info/licenses/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
-bess_jpl-1.13.0.dist-info/METADATA,sha256=nIdcT3g6u2s-MEdvhE2aI1ttR1lKyI3-in_bvPEBNmw,6340
-bess_jpl-1.13.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-bess_jpl-1.13.0.dist-info/top_level.txt,sha256=GaKnzt-BBktYn1o-w4Qzh_jHxse4Y3ACOxZQrB2ufhc,9
-bess_jpl-1.13.0.dist-info/RECORD,,
+BESS_JPL/version.txt,sha256=-Q2F9kDoBjnO1g1dZcNwbI4pu99GLJDvdHkxmPdU4JU,6
+bess_jpl-1.15.0.dist-info/licenses/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
+bess_jpl-1.15.0.dist-info/METADATA,sha256=N6XioN8Coziy-RjHi-pXGw_CRN4KWqhwyCzBpG7HPec,6313
+bess_jpl-1.15.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+bess_jpl-1.15.0.dist-info/top_level.txt,sha256=GaKnzt-BBktYn1o-w4Qzh_jHxse4Y3ACOxZQrB2ufhc,9
+bess_jpl-1.15.0.dist-info/RECORD,,