BESS-JPL 1.26.0__py3-none-any.whl

BESS_JPL/retrieve_BESS_inputs.py

@@ -0,0 +1,279 @@
+from typing import Union
+from datetime import datetime
+import logging
+import numpy as np
+
+import rasters as rt
+from rasters import Raster, RasterGeometry
+
+from check_distribution import check_distribution
+
+from sun_angles import calculate_SZA_from_DOY_and_hour
+from solar_apparent_time import calculate_solar_day_of_year, calculate_solar_hour_of_day
+from koppengeiger import load_koppen_geiger
+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 NASADEMConnection
+
+from .constants import *
+from .exceptions import *
+from .colors import *
+from .C3_photosynthesis import *
+from .C4_photosynthesis import *
+from .canopy_energy_balance import *
+from .canopy_longwave_radiation import *
+from .canopy_shortwave_radiation import *
+from .carbon_water_fluxes import *
+from .FVC_from_NDVI import *
+from .interpolate_C3_C4 import *
+from .LAI_from_NDVI import *
+from .load_C4_fraction import *
+from .load_carbon_uptake_efficiency import *
+from .load_kn import *
+from .load_NDVI_minimum import *
+from .load_NDVI_maximum import *
+from .load_peakVCmax_C3 import *
+from .load_peakVCmax_C4 import *
+from .load_ball_berry_intercept_C3 import *
+from .load_ball_berry_slope_C3 import *
+from .load_ball_berry_slope_C4 import *
+from .calculate_VCmax import *
+from .meteorology import *
+from .soil_energy_balance import *
+from .retrieve_BESS_JPL_GEOS5FP_inputs import retrieve_BESS_JPL_GEOS5FP_inputs
+
+logger = logging.getLogger(__name__)
+
+def retrieve_BESS_inputs(ST_C: Union[Raster, np.ndarray],  # surface temperature in Celsius
+        NDVI: Union[Raster, np.ndarray],  # NDVI
+        albedo: Union[Raster, np.ndarray],  # surface albedo
+        geometry: RasterGeometry = None,
+        time_UTC: datetime = None,
+        hour_of_day: np.ndarray = None,
+        day_of_year: np.ndarray = None,
+        GEOS5FP_connection: GEOS5FP = None,
+        elevation_m: Union[Raster, np.ndarray] = None,  # elevation in meters
+        Ta_C: Union[Raster, np.ndarray] = None,  # air temperature in Celsius
+        RH: Union[Raster, np.ndarray] = None,  # relative humidity as a proportion
+        NDVI_minimum: Union[Raster, np.ndarray] = None,  # minimum NDVI
+        NDVI_maximum: Union[Raster, np.ndarray] = None,  # maximum NDVI
+        PAR_albedo: Union[Raster, np.ndarray] = None, # surface albedo in visible wavelengths (initialized to surface albedo if left as None)
+        NIR_albedo: Union[Raster, np.ndarray] = None, # surface albedo in near-infrared wavelengths (initialized to surface albedo if left as None)
+        COT: Union[Raster, np.ndarray] = None,  # cloud optical thickness
+        AOT: Union[Raster, np.ndarray] = None,  # aerosol optical thickness
+        vapor_gccm: Union[Raster, np.ndarray] = None,  # water vapor in g/ccm
+        ozone_cm: Union[Raster, np.ndarray] = None,  # ozone in cm
+        KG_climate: Union[Raster, np.ndarray] = None,  # KG climate
+        canopy_height_meters: Union[Raster, np.ndarray] = None,  # canopy height in meters
+        Ca: Union[Raster, np.ndarray] = None,  # atmospheric CO2 concentration in ppm
+        wind_speed_mps: Union[Raster, np.ndarray] = None,  # wind speed in meters per second
+        SZA_deg: Union[Raster, np.ndarray] = None,  # solar zenith angle in degrees
+        canopy_temperature_C: Union[Raster, np.ndarray] = None, # canopy temperature in Celsius (initialized to surface temperature if left as None)
+        soil_temperature_C: Union[Raster, np.ndarray] = None, # soil temperature in Celsius (initialized to surface temperature if left as None)
+        C4_fraction: Union[Raster, np.ndarray] = None,  # fraction of C4 plants
+        carbon_uptake_efficiency: Union[Raster, np.ndarray] = None,  # intrinsic quantum efficiency for carbon uptake
+        kn: np.ndarray = None,
+        ball_berry_intercept_C3: np.ndarray = None,  # Ball-Berry intercept for C3 plants
+        ball_berry_intercept_C4: Union[np.ndarray, float] = BALL_BERRY_INTERCEPT_C4, # Ball-Berry intercept for C4 plants
+        ball_berry_slope_C3: np.ndarray = None,  # Ball-Berry slope for C3 plants
+        ball_berry_slope_C4: np.ndarray = None,  # Ball-Berry slope for C4 plants
+        peakVCmax_C3_μmolm2s1: np.ndarray = None,  # peak maximum carboxylation rate for C3 plants
+        peakVCmax_C4_μmolm2s1: np.ndarray = None,  # peak maximum carboxylation rate for C4 plants
+        CI: Union[Raster, np.ndarray] = None,
+        C4_fraction_scale_factor: float = C4_FRACTION_SCALE_FACTOR,
+        MODISCI_connection: MODISCI = None,
+        NASADEM_connection: NASADEMConnection = None,
+        resampling: str = RESAMPLING,
+        GEDI_download_directory: str = GEDI_DOWNLOAD_DIRECTORY,
+        offline_mode: bool = False) -> dict:
+    results = {}
+
+    if (day_of_year is None or hour_of_day is None) and time_UTC is not None and geometry is not None:
+        day_of_year = calculate_solar_day_of_year(time_UTC=time_UTC, geometry=geometry)
+        hour_of_day = calculate_solar_hour_of_day(time_UTC=time_UTC, geometry=geometry)
+
+    if time_UTC is None and day_of_year is None and hour_of_day is None:
+        raise ValueError("no time given between time_UTC, day_of_year, and hour_of_day")
+
+    # calculate solar zenith angle if not provided
+    if SZA_deg is None:
+        SZA_deg = calculate_SZA_from_DOY_and_hour(geometry.lat, geometry.lon, day_of_year, hour_of_day)
+
+    if isinstance(SZA_deg, np.ndarray):
+        # If array contains string representations, convert them first
+        if SZA_deg.dtype == object or SZA_deg.dtype.kind in ['U', 'S']:
+            # Handle string arrays by converting each element
+            # This handles cases like '[71.46303285]' or '71.46303285'
+            SZA_deg = np.array([float(str(x).strip('[]')) for x in SZA_deg], dtype=np.float32)
+        else:
+            # cast SZA_deg numpy array to float32
+            SZA_deg = SZA_deg.astype(np.float32)
+
+    print(type(SZA_deg))
+    print(SZA_deg.dtype if isinstance(SZA_deg, np.ndarray) else type(SZA_deg))
+
+    check_distribution(SZA_deg, "SZA_deg")
+    results["SZA_deg"] = SZA_deg
+
+
+    if CI is None and geometry is not None:
+        if offline_mode:
+            raise MissingOfflineParameter("CI not provided in offline mode")
+        
+        if MODISCI_connection is None:
+            MODISCI_connection = MODISCI()
+
+        CI = MODISCI_connection.CI(geometry=geometry, resampling=resampling)
+
+    check_distribution(CI, "CI")
+    results["CI"] = CI
+
+    if elevation_m is None and geometry is not None:
+        if offline_mode:
+            raise MissingOfflineParameter("elevation_m not provided in offline mode")
+        
+        if NASADEM_connection is None:
+            NASADEM_connection = NASADEMConnection()
+
+        elevation_m = NASADEM_connection.elevation_m(geometry=geometry)
+
+    check_distribution(elevation_m, "elevation_m")
+    results["elevation_m"] = elevation_m
+
+
+    # load minimum NDVI if not provided
+    if NDVI_minimum is None and geometry is not None:
+        NDVI_minimum = load_NDVI_minimum(geometry=geometry, resampling=resampling)
+
+    check_distribution(NDVI_minimum, "NDVI_minimum")
+    results["NDVI_minimum"] = NDVI_minimum
+
+    # load maximum NDVI if not provided
+    if NDVI_maximum is None and geometry is not None:
+        NDVI_maximum = load_NDVI_maximum(geometry=geometry, resampling=resampling)
+
+    check_distribution(NDVI_maximum, "NDVI_maximum")
+    results["NDVI_maximum"] = NDVI_maximum
+
+    # load C4 fraction if not provided
+    if C4_fraction is None:
+        C4_fraction = load_C4_fraction(
+            geometry=geometry, 
+            resampling=resampling,
+            scale_factor=C4_fraction_scale_factor
+        )
+
+    check_distribution(C4_fraction, "C4_fraction")
+    results["C4_fraction"] = C4_fraction
+
+    # load carbon uptake efficiency if not provided
+    if carbon_uptake_efficiency is None:
+        carbon_uptake_efficiency = load_carbon_uptake_efficiency(geometry=geometry, resampling=resampling)
+    
+    check_distribution(carbon_uptake_efficiency, "carbon_uptake_efficiency")
+    results["carbon_uptake_efficiency"] = carbon_uptake_efficiency
+
+    # load kn if not provided
+    if kn is None:
+        kn = load_kn(geometry=geometry, resampling=resampling)
+
+    check_distribution(kn, "kn")
+    results["kn"] = kn
+
+    # load peak VC max for C3 plants if not provided
+    if peakVCmax_C3_μmolm2s1 is None:
+        peakVCmax_C3_μmolm2s1 = load_peakVCmax_C3(geometry=geometry, resampling=resampling)
+
+    check_distribution(peakVCmax_C3_μmolm2s1, "peakVCmax_C3_μmolm2s1")
+    results["peakVCmax_C3_μmolm2s1"] = peakVCmax_C3_μmolm2s1
+
+    # load peak VC max for C4 plants if not provided
+    if peakVCmax_C4_μmolm2s1 is None:
+        peakVCmax_C4_μmolm2s1 = load_peakVCmax_C4(geometry=geometry, resampling=resampling)
+
+    check_distribution(peakVCmax_C4_μmolm2s1, "peakVCmax_C4_μmolm2s1")
+    results["peakVCmax_C4_μmolm2s1"] = peakVCmax_C4_μmolm2s1
+
+    # load Ball-Berry slope for C3 plants if not provided
+    if ball_berry_slope_C3 is None:
+        ball_berry_slope_C3 = load_ball_berry_slope_C3(geometry=geometry, resampling=resampling)
+    
+    check_distribution(ball_berry_slope_C3, "ball_berry_slope_C3")
+    results["ball_berry_slope_C3"] = ball_berry_slope_C3
+
+    # load Ball-Berry slope for C4 plants if not provided
+    if ball_berry_slope_C4 is None:
+        ball_berry_slope_C4 = load_ball_berry_slope_C4(geometry=geometry, resampling=resampling)
+
+    check_distribution(ball_berry_slope_C4, "ball_berry_slope_C4")
+    results["ball_berry_slope_C4"] = ball_berry_slope_C4
+
+    # load Ball-Berry intercept for C3 plants if not provided
+    if ball_berry_intercept_C3 is None:
+        ball_berry_intercept_C3 = load_ball_berry_intercept_C3(geometry=geometry, resampling=resampling)
+
+    check_distribution(ball_berry_intercept_C3, "ball_berry_intercept_C3")
+    results["ball_berry_intercept_C3"] = ball_berry_intercept_C3
+
+    # load koppen geiger climate classification if not provided
+    if KG_climate is None:
+        KG_climate = load_koppen_geiger(geometry=geometry)
+
+    check_distribution(np.float32(KG_climate), "KG_climate")
+    results["KG_climate"] = KG_climate
+
+    # load canopy height in meters if not provided
+    if canopy_height_meters is None:
+        canopy_height_meters = load_canopy_height(
+            geometry=geometry, 
+            resampling=resampling,
+            source_directory=GEDI_download_directory
+        )
+
+    # canopy height defaults to zero
+    canopy_height_meters = np.where(np.isnan(canopy_height_meters), 0, canopy_height_meters)
+
+    check_distribution(canopy_height_meters, "canopy_height_meters")
+    results["canopy_height_meters"] = canopy_height_meters
+
+
+    # Retrieve GEOS-5 FP inputs if not provided
+    GEOS5FP_inputs = retrieve_BESS_JPL_GEOS5FP_inputs(
+        time_UTC=time_UTC,
+        geometry=geometry,
+        albedo=albedo,
+        GEOS5FP_connection=GEOS5FP_connection,
+        Ta_C=Ta_C,
+        RH=RH,
+        COT=COT,
+        AOT=AOT,
+        vapor_gccm=vapor_gccm,
+        ozone_cm=ozone_cm,
+        PAR_albedo=PAR_albedo,
+        NIR_albedo=NIR_albedo,
+        Ca=Ca,
+        wind_speed_mps=wind_speed_mps,
+        resampling=resampling,
+        offline_mode=offline_mode
+    )
+    
+    # Merge GEOS-5 FP inputs into results dictionary
+    results.update(GEOS5FP_inputs)
+
+    # canopy temperature defaults to surface temperature
+    if canopy_temperature_C is None:
+        canopy_temperature_C = ST_C
+
+    check_distribution(canopy_temperature_C, "canopy_temperature_C")
+    results["canopy_temperature_C"] = canopy_temperature_C
+
+    # soil temperature defaults to surface temperature
+    if soil_temperature_C is None:
+        soil_temperature_C = ST_C
+
+    check_distribution(soil_temperature_C, "soil_temperature_C")
+    results["soil_temperature_C"] = soil_temperature_C
+
+    return results

BESS_JPL/soil_energy_balance.py

@@ -0,0 +1,35 @@
+import numpy as np
+
+
+def soil_energy_balance(
+        Ts_K: np.ndarray,
+        Ta_K: np.ndarray,
+        G_Wm2: np.ndarray,
+        VPD_Pa: np.ndarray,
+        RH: np.ndarray,
+        gamma: np.ndarray,
+        Cp: np.ndarray,
+        rhoa: np.ndarray,
+        desTa: np.ndarray,
+        Rs: np.ndarray,
+        ASW_soil_Wm2: np.ndarray,
+        ALW_soil_Wm2: np.ndarray,
+        Ls: np.ndarray,
+        epsa: np.ndarray):
+    # Net radiation
+    # Rn = Rnet - Rn_Sun - Rn_Sh
+    sigma = 5.670373e-8  # [W m-2 K-4] (Wiki)
+    Rn_soil_Wm2 = np.clip(ASW_soil_Wm2 + ALW_soil_Wm2 - Ls - 4.0 * epsa * sigma * (Ta_K ** 3) * (Ts_K - Ta_K), 0, None)
+    # G = Rn * 0.35
+
+    # Latent heat
+    LE_soil_Wm2 = desTa / (desTa + gamma) * (Rn_soil_Wm2 - G_Wm2) * (RH ** (VPD_Pa / 1000.0))  # (Ryu et al., 2011)
+    LE_soil_Wm2 = np.clip(LE_soil_Wm2, 0, Rn_soil_Wm2)
+    # Sensible heat
+    H_soil_Wm2 = np.clip(Rn_soil_Wm2 - G_Wm2 - LE_soil_Wm2, 0, Rn_soil_Wm2)
+
+    # Update temperature
+    dT = np.clip(Rs / (rhoa * Cp) * H_soil_Wm2, -20, 20)
+    Ts_K = Ta_K + dT
+
+    return Rn_soil_Wm2, LE_soil_Wm2, Ts_K

BESS_JPL/verify.py

@@ -0,0 +1,127 @@
+def verify() -> bool:
+    """
+    Verifies the correctness of the BESS model implementation by comparing
+    its outputs to a reference dataset.
+
+    This function loads a known input table and the corresponding expected output table.
+    It runs the model on the input data, then compares the resulting outputs to the
+    reference outputs for key variables using strict numerical tolerances. If all
+    outputs match within tolerance, the function returns True. Otherwise, it prints
+    which column failed and returns False.
+
+    Returns:
+        bool: True if all model outputs match the reference outputs within tolerance, False otherwise.
+    """
+    import pandas as pd
+    import numpy as np
+    from .ECOv002_calval_BESS_inputs import load_ECOv002_calval_BESS_inputs
+    from .process_BESS_table import process_BESS_table
+    import os
+
+    # Load input and output tables
+    input_df = load_ECOv002_calval_BESS_inputs()
+    module_dir = os.path.dirname(os.path.abspath(__file__))
+    output_file_path = os.path.join(module_dir, "ECOv002-cal-val-BESS-JPL-outputs.csv")
+    output_df = pd.read_csv(output_file_path)
+
+    # Check that input dataset contains all required inputs for BESS model
+    # This ensures no on-the-fly data retrieval is needed
+    required_inputs = [
+        # Core inputs
+        'ST_C', 'NDVI', 'albedo', 'time_UTC',
+        # Geometry (either geometry or lat/lon)
+        # 'geometry' or both 'lat' and 'lon' required - checked separately
+        'day_of_year', 'hour_of_day',
+        # Meteorological inputs
+        'Ta_C',  # or 'Ta'
+        'RH', 'elevation_m',  # or 'elevation_km'
+        'wind_speed_mps',
+        # Atmospheric inputs from GEOS5FP
+        'COT', 'AOT', 'vapor_gccm', 'ozone_cm',
+        # Note: Ca defaults to 400 ppm if not provided
+        # Radiation inputs
+        'PAR_albedo', 'NIR_albedo',
+        # Vegetation parameters
+        'NDVI_minimum', 'NDVI_maximum', 'C4_fraction',
+        'carbon_uptake_efficiency', 'kn',
+        # Photosynthesis parameters
+        'peakVCmax_C3', 'peakVCmax_C4',
+        'ball_berry_slope_C3', 'ball_berry_slope_C4', 'ball_berry_intercept_C3',
+        # Other required inputs
+        'KG_climate', 'CI', 'canopy_height_meters', 'SZA_deg',
+        # Note: canopy_temperature_C and soil_temperature_C default to ST_C if not provided
+    ]
+    
+    # Check for alternative column names
+    missing_inputs = []
+    for col in required_inputs:
+        # Handle alternative names
+        if col == 'Ta_C' and 'Ta_C' not in input_df and 'Ta' not in input_df:
+            missing_inputs.append('Ta_C (or Ta)')
+        elif col == 'elevation_m' and 'elevation_m' not in input_df and 'elevation_km' not in input_df:
+            missing_inputs.append('elevation_m (or elevation_km)')
+        elif col not in ['Ta_C', 'elevation_m'] and col not in input_df:
+            missing_inputs.append(col)
+    
+    # Check geometry requirement
+    has_geometry = 'geometry' in input_df or ('lat' in input_df and 'lon' in input_df)
+    if not has_geometry:
+        missing_inputs.append('geometry (or lat and lon)')
+    
+    if missing_inputs:
+        print("Input verification failed: Missing required inputs for BESS model.")
+        print("The following inputs are missing from the input dataset:")
+        for inp in missing_inputs:
+            print(f"  - {inp}")
+        print("\nThese inputs must be present to run the model without on-the-fly data retrieval.")
+        return False
+    
+    print("Input verification passed: All required BESS inputs are present.")
+    
+    # Run the model on the input table
+    model_df = process_BESS_table(input_df, offline_mode=True)
+
+    # Columns to compare (model outputs)
+    output_columns = [
+        "G_Wm2",
+        "Rn_Wm2",
+        "LE_Wm2"
+    ]
+
+    # Compare each output column and collect mismatches
+    mismatches = []
+    for col in output_columns:
+        if col not in model_df or col not in output_df:
+            mismatches.append((col, 'missing_column', None))
+            continue
+        model_vals = model_df[col].values
+        ref_vals = output_df[col].values
+        # Use numpy allclose for floating point comparison
+        # Tolerances account for minor platform/version differences (macOS vs Ubuntu, different numpy versions)
+        if not np.allclose(model_vals, ref_vals, rtol=1e-4, atol=1e-7, equal_nan=True):
+            # Find indices where values differ
+            diffs = np.abs(model_vals - ref_vals)
+            max_diff = np.nanmax(diffs)
+            idxs = np.where(~np.isclose(model_vals, ref_vals, rtol=1e-4, atol=1e-7, equal_nan=True))[0]
+            mismatch_info = {
+                'indices': idxs.tolist(),
+                'model_values': model_vals[idxs].tolist(),
+                'ref_values': ref_vals[idxs].tolist(),
+                'diffs': diffs[idxs].tolist(),
+                'max_diff': float(max_diff)
+            }
+            mismatches.append((col, 'value_mismatch', mismatch_info))
+    if mismatches:
+        print("Verification failed. Details:")
+        for col, reason, info in mismatches:
+            if reason == 'missing_column':
+                print(f"  Missing column: {col}")
+            elif reason == 'value_mismatch':
+                print(f"  Mismatch in column: {col}")
+                print(f"    Max difference: {info['max_diff']}")
+                print(f"    Indices off: {info['indices']}")
+                print(f"    Model values: {info['model_values']}")
+                print(f"    Reference values: {info['ref_values']}")
+                print(f"    Differences: {info['diffs']}")
+        return False
+    return True

BESS_JPL/version.py

@@ -0,0 +1,3 @@
+from importlib.metadata import version
+
+__version__ = version("BESS-JPL")

bess_jpl-1.26.0.dist-info/METADATA

@@ -0,0 +1,102 @@
+Metadata-Version: 2.4
+Name: BESS-JPL
+Version: 1.26.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
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: check-distribution>=1.2.0
+Requires-Dist: daylight-evapotranspiration
+Requires-Dist: ECOv002-calval-tables>=1.6.0
+Requires-Dist: FLiESANN>=2.13.0
+Requires-Dist: gedi-canopy-height>=1.1.0
+Requires-Dist: GEOS5FP>=2.19.0
+Requires-Dist: koppengeiger>=1.1.0
+Requires-Dist: MODISCI>=1.4.0
+Requires-Dist: monte-carlo-sensitivity>=1.5.0
+Requires-Dist: NASADEM>=1.3.0
+Requires-Dist: numpy
+Requires-Dist: pytictoc
+Requires-Dist: rasters>=1.9.0
+Requires-Dist: seaborn
+Requires-Dist: solar-apparent-time>=1.5.3
+Provides-Extra: dev
+Requires-Dist: build; extra == "dev"
+Requires-Dist: pytest>=6.0; extra == "dev"
+Requires-Dist: pytest-cov; extra == "dev"
+Requires-Dist: jupyter; extra == "dev"
+Requires-Dist: pytest; extra == "dev"
+Requires-Dist: twine; extra == "dev"
+Dynamic: license-file
+
+# Breathing Earth System Simulator (BESS) Model Python Implementation
+
+[![CI](https://github.com/JPL-Evapotranspiration-Algorithms/breathing-earth-system-simulator/actions/workflows/ci.yml/badge.svg)](https://github.com/JPL-Evapotranspiration-Algorithms/breathing-earth-system-simulator/actions/workflows/ci.yml)
+
+This software package is a Python implementation of the Breathing Earth System Simulator (BESS) model. It was re-implemented in Python by Gregory Halverson at Jet Propulsion Laboratory based on [MATLAB code](https://www.environment.snu.ac.kr/bess-flux) produced by Youngryel Ryu at Seoul University. The BESS model was designed to quantify global gross primary productivity (GPP) and evapotranspiration (ET) using MODIS with a spatial resolution of 1–5 km and a temporal resolution of 8 days. It couples atmospheric and canopy radiative transfer processes with photosynthesis, stomatal conductance, and transpiration models on sunlit and shaded portions of vegetation and soil. An artificial neural network emulator of Hideki Kobayashi's Forest Light Environmental Simulator (FLiES) radiative transfer model is used to estimate incoming solar radiation. This implementation of BESS was designed to process GPP at fine spatial resolution using surface temperature from the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission and normalized difference vegetation index (NDVI) and albedo from the [Spatial Timeseries for Automated high-Resolution multi-Sensor (STARS) data fusion system](https://github.com/STARS-Data-Fusion). The software was developed as part of a research grant by the NASA Research Opportunities in Space and Earth Sciences (ROSES) program. It was designed for use by the ECOSTRESS mission as a precursor for the Surface Biology and Geology (SBG) mission. However, it may also be useful for general remote sensing and GIS projects in Python. This package can be utilized for remote sensing research in Jupyter notebooks and deployed for operations in data processing pipelines. This software is being released according to the SPD-41 open-science requirements of NASA-funded ROSES projects.
+
+Gregory H. Halverson (they/them)<br>
+[gregory.h.halverson@jpl.nasa.gov](mailto:gregory.h.halverson@jpl.nasa.gov)<br>
+Lead developer<br>
+NASA Jet Propulsion Laboratory 329G
+
+Youngryel Ryu (he/him)<br>
+[yryu@snu.ac.kr](mailto:yryu@snu.ac.kr)<br>
+BESS algorithm inventor<br>
+Seoul National University
+
+Hideki Kobayashi (he/him)<br>
+[hkoba@jamstec.go.jp](mailto:hkoba@jamstec.go.jp)<br>
+FLiES algorithm inventor<br>
+Japan Agency for Marine-Earth Science and Technology
+
+Robert Freepartner (he/him)<br>
+[robert.freepartner@jpl.nasa.gov](robert.freepartner@jpl.nasa.gov)<br>
+MATLAB to python translation<br>
+Raytheon
+
+Joshua Fisher (he/him)<br>
+[jbfisher@chapman.edu](mailto:jbfisher@chapman.edu)<br>
+Concept development and project management<br>
+Chapman University
+
+Kerry Cawse-Nicholson (she/her)<br>
+[kerry-anne.cawse-nicholson@jpl.nasa.gov](mailto:kerry-anne.cawse-nicholson@jpl.nasa.gov)<br>
+Project management<br>
+NASA Jet Propulsion Laboratory 329G
+
+Zoe Pierrat (she/her)<br>
+[zoe.a.pierrat@jpl.nasa.gov](mailto:zoe.a.pierrat@jpl.nasa.gov)<br>
+Algorithm maintenance<br>
+NASA Jet Propulsion Laboratory 329G
+
+Claire Villanueva-Weeks (she/her)<br>
+[claire.s.villanueva-weeks@jpl.nasa.gov](mailto:claire.s.villanueva-weeks@jpl.nasa.gov)<br>
+Code maintenance<br>
+NASA Jet Propulsion Laboratory 329G
+
+## Installation
+
+Use the pip package manager to install this package
+
+```
+pip install breathing-earth-system-simulator
+```
+
+## References
+
+The following scientific references provide detailed information about the BESS model and its underlying algorithms:
+
+1. Ryu, Y., Baldocchi, D. D., Kobayashi, H., van Ingen, C., Li, J., Black, T. A., ... & Ueyama, M. (2011). Integration of MODIS land and atmosphere products with a coupled-process model to estimate gross primary productivity and evapotranspiration from 1 km to global scales. *Remote Sensing of Environment, 115*(8), 1865-1874. https://doi.org/10.1016/j.rse.2011.03.009
+
+2. Kobayashi, H., Ryu, Y., Baldocchi, D. D., Welles, J. M., & Norman, J. M. (2012). On the correct estimation of gap fraction: How to remove scattered radiation in gap fraction measurements? *Agricultural and Forest Meteorology, 160*, 14-25. https://doi.org/10.1016/j.agrformet.2012.03.001
+
+3. Fisher, J. B., Lee, B., Purdy, A. J., Halverson, G. H., Dohlen, M. B., & Tu, K. P. (2020). ECOSTRESS: NASA's next generation mission to measure evapotranspiration from the International Space Station. *Water Resources Research, 56*(4), e2019WR026058. https://doi.org/10.1029/2019WR026058
+
+4. Ryu, Y., Jiang, C., Kobayashi, H., & Detto, M. (2018). Modis-derived global land products of shortwave radiation and diffuse and total photosynthetically active radiation at 5 km resolution from 2000. *Remote Sensing of Environment, 204*, 812-825. https://doi.org/10.1016/j.rse.2017.09.021
+
+5. Kobayashi, H., & Iwabuchi, H. (2008). A coupled 1-D atmosphere and canopy radiative transfer model for an atmosphere with a nonlambertian surface. *Journal of Quantitative Spectroscopy and Radiative Transfer, 109*(17-18), 2955-2970. https://doi.org/10.1016/j.jqsrt.2008.07.008

bess_jpl-1.26.0.dist-info/RECORD

@@ -0,0 +1,73 @@
+BESS_JPL/BESS_JPL.py,sha256=l5ubRRpPYlZxhmHjsB-88J5_2exgNLB5lvvVVO10fjw,1720
+BESS_JPL/C3_photosynthesis.py,sha256=VmBYwxCK1URhA2_Etb96JXSqubyvuftjyuhHha9QOXQ,6506
+BESS_JPL/C4_fraction.jpeg,sha256=ECaEYWA8MnNd5z0Yq3beeUCR93KtTDxM8MrH-YYon3Y,2746
+BESS_JPL/C4_fraction.tif,sha256=DFS-J08JCjnf0hi_T_dFridAegVTiU3FGyMBjS9yc1w,25446
+BESS_JPL/C4_fraction.tif.aux.xml,sha256=5Ytn5UtZdGWQUKNmg_kiOU2zRQsm82sY7bq-j3TS1QY,358
+BESS_JPL/C4_photosynthesis.py,sha256=NtT-fV8DQbLXMyePGGDSwYadbsDY0ZW3WFGJ9pxxL40,7035
+BESS_JPL/ECOv002-cal-val-BESS-JPL-GEOS5FP-inputs.csv,sha256=aw7YL3S81qu3SUNp3DHG7Q4HtOdOWxBLycEcKUJIBrA,104716
+BESS_JPL/ECOv002-cal-val-BESS-JPL-inputs.csv,sha256=3oQVYY_dSzIyclORmE181gtbVvWcj8fTFR9E_XW9gjQ,1317067
+BESS_JPL/ECOv002-cal-val-BESS-JPL-outputs.csv,sha256=W5cDy2I3X37z3YaY0ehEeprljIOqfGgEC_5NlpSoUVY,2984084
+BESS_JPL/ECOv002-cal-val-FLiESANN-inputs.csv,sha256=AybkhW0EmPGlkoYe_Ph0aoh7TrA9yd0ozJ6cvDwS4fE,946508
+BESS_JPL/ECOv002-static-tower-BESS-JPL-inputs.csv,sha256=uRrCRPfbnsU_a3Otnnpk47gcBpTrS079EZPrB0KNc8s,33603
+BESS_JPL/ECOv002_calval_BESS_inputs.py,sha256=9XTFJHn9EczVNLr2WDjsJ4b5nHZonkxWpS5yEogo2RY,1109
+BESS_JPL/ECOv002_static_tower_BESS_inputs.py,sha256=jvMGMWVo0UnNsGXebrYDVTgJXroH-MTkLj2zYcCYADk,622
+BESS_JPL/FVC_from_NDVI.py,sha256=LYI1ESsrO8UxlPOr_dxJYDeUsrwjVAwrV_OhL-tcAWc,567
+BESS_JPL/LAI_from_NDVI.py,sha256=uwV1voc2q45hmsSHzIhehA3ZHLivzSeoST-ht9ecSIE,783
+BESS_JPL/NDVI_maximum.jpeg,sha256=6kc0-eE_WOU31bUztdjoWsZ3ZmqHAG-ZwS7lMiZyBGA,763946
+BESS_JPL/NDVI_maximum.tif,sha256=yRgdpB-xoVNMveIUapCrnGAwHX7WkH6pVXu4HECuFfc,19664541
+BESS_JPL/NDVI_minimum.jpeg,sha256=AfpRU8PikwExBJ2IdcFKewg6aTcp3gx3y3EQMzeA7vg,644871
+BESS_JPL/NDVI_minimum.tif,sha256=wPV0CFfQSLfiTAEDlwIAI857pA7aw0iqWhVlOAZFsJg,18762242
+BESS_JPL/__init__.py,sha256=UaK3mWbc5jNoA4kNzIWF9h7CQs64O_chKrpuxNkH48o,95
+BESS_JPL/ball_berry_intercept_C3.jpeg,sha256=c8jlRD3i8cnWXQNRUfyNTuObEpptcG6p6kQd2RGMK70,634200
+BESS_JPL/ball_berry_intercept_C3.tif,sha256=O7ATUIKtXcIjkUI2lloTiHAgvRk4gWyWjYGh5Xc_CRY,1026642
+BESS_JPL/ball_berry_slope_C3.jpeg,sha256=KdFUYG4h1nRNb_eK7w-M2kqChJnhbcGmal9Z2Jos1nI,225597
+BESS_JPL/ball_berry_slope_C3.tif,sha256=PEAZ09AMMbDUXi5R4FsyMQGsanLd4Ptzk-pRHnaa-oU,407646
+BESS_JPL/ball_berry_slope_C4.jpeg,sha256=JJJn3C17duBOlEe98PCyycBT3SS-YBhHJpJZYTq4GMw,1009208
+BESS_JPL/ball_berry_slope_C4.tif,sha256=UedoUOnOjBxQ9pAaf5nXgDa8wvKOHVBD3sVToV12aM8,1012948
+BESS_JPL/calculate_VCmax.py,sha256=dQ6uGzgqrv2rAaUlhxuLauWhMSuUGvb9fgS2pRYugLY,4129
+BESS_JPL/calculate_bulk_aerodynamic_resistance.py,sha256=5j0YJGI_dEnMDud4pv0AyOf03hYxv5cJ6m0nbIJtlWU,5668
+BESS_JPL/calculate_friction_velocity.py,sha256=njTyo9H5tXxqS7-VNGOqf5pCqTV3Y5vaNwYEmF1kxsI,4932
+BESS_JPL/canopy_energy_balance.py,sha256=aNnODLouI33GkVrqmNaviWGzyPvV0DgnbV4iQncKsw4,4945
+BESS_JPL/canopy_longwave_radiation.py,sha256=Hv5TUJ1O_PBIYRmx8kn6Ake8cyKEL9h0cfN4wsC0UgM,5436
+BESS_JPL/canopy_shortwave_radiation.py,sha256=WmzH6birUGkT4O0q7eG1Eo4ed3xi9LVKTydIf9r8J_s,13280
+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=hfTug6k7ch_i2k-HoICVJZTJW310Uc5esK258qaEN-M,16758
+BESS_JPL/colors.py,sha256=vRZoazJwgSa6v3sC678GoyVPRUxVfYg1JMibzI3jHiw,635
+BESS_JPL/constants.py,sha256=LYtSoFHA8zaerW4BOSY8mRcKvqpf3foIX5NX8xwVVFM,688
+BESS_JPL/exceptions.py,sha256=KpYl0RC4FtfmXrMHpRhVb6InwbrALwh3WfZIfbdbpgA,132
+BESS_JPL/generate_BESS_GEOS5FP_inputs.py,sha256=7J5Cv6t6fZZbpETOquTKWf604dvLC1vePE3y8almUK8,1685
+BESS_JPL/generate_BESS_inputs_table.py,sha256=QqYUG98bWXLoTp5TqBfDPI8xXoOu9kEQG04g7oirkBI,9763
+BESS_JPL/generate_input_dataset.py,sha256=4YbPdKvcY69PQVjTDKVb5pYvt7DbMLSzAZcsRIQFY1E,11768
+BESS_JPL/generate_output_dataset.py,sha256=tz7pM3JF2fI1HqKWT3ub1P40VQYzOLoT-qrTV4zeIOI,839
+BESS_JPL/interpolate_C3_C4.py,sha256=adIyk03Yt6SJ_jFCUm6ivK9QpZYw10C3wJKpx1O56Pk,394
+BESS_JPL/kn.jpeg,sha256=wsI40jabT9_iPpQmOxd3GZzfAN7oZ2DubiARmULFUS0,854097
+BESS_JPL/kn.tif,sha256=nZAIAms7LRDVsadxa4z8yRqOudWu3AjJSm6l1T3ywuY,819988
+BESS_JPL/load_C4_fraction.py,sha256=qemgzsyDVKvzvql1p_N7g23L8KYvbLw-qxjWlUC3p_Q,624
+BESS_JPL/load_NDVI_maximum.py,sha256=HMOvChb9CWItkDpGCZtZ9Ode2N4Q_MllmnOgiAt5yzQ,528
+BESS_JPL/load_NDVI_minimum.py,sha256=ed_qmkNCIy-H7Q5POA4mSlUvYUTknKiz31ioXZGpbvU,528
+BESS_JPL/load_ball_berry_intercept_C3.py,sha256=6auXZHej127rhUCFXe7vdEg9C4CosFSHLTcdqKs06Ls,388
+BESS_JPL/load_ball_berry_slope_C3.py,sha256=9nmbdkGnJ7UWXwktpkE8rL1LGA1vi71YZdzpme04P84,380
+BESS_JPL/load_ball_berry_slope_C4.py,sha256=lsKwrkbChEm3sspUgAXjm86Q4VBirwFOZStiO9NuGx8,380
+BESS_JPL/load_carbon_uptake_efficiency.py,sha256=RbWxLbqHnxiHw7XHyFBuWRToc8hDjUR6EUUntiJfe0M,390
+BESS_JPL/load_kn.py,sha256=39RkQkkdmey0IPInziJI0kSiI8alRk5Hz42pzm-qihU,346
+BESS_JPL/load_peakVCmax_C3.py,sha256=A-Wg4PjPgogfpRsAMo6ZNEXM7G3KbnF7ZRV3k_jTtEA,401
+BESS_JPL/load_peakVCmax_C4.py,sha256=iVmJTitkMkB4AjfjanAi3yv1lWO4sSF6J7cjuRmaCwU,401
+BESS_JPL/meteorology.py,sha256=B7tguKu_1H0bWKJI0YQpSZzk9cXS-wc6W4OIe_5igds,14799
+BESS_JPL/model.py,sha256=r4fpvRngARgkQsDmbNVPH7KeC_lHYQiwYheNOI0nanE,25838
+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=U3SMtCEESm7bAu3Vaccdy_MnL2l2Er7LdbZnbryYBkI,13039
+BESS_JPL/process_paw_and_gao_LE.py,sha256=e-hMowIWC-aORGOyJU-omGsFk0Pp8YmTjl5wwQiOXOQ,1760
+BESS_JPL/retrieve_BESS_JPL_GEOS5FP_inputs.py,sha256=mS6IcXELzj2j8fHAAVwregI-FU3g_9l4tolu24x1_l8,10651
+BESS_JPL/retrieve_BESS_inputs.py,sha256=VubZ3ui6ZYnLGpjDjgrcwVRRobeEslCG6fpIPIzzK6A,11924
+BESS_JPL/soil_energy_balance.py,sha256=u5YMywVQPard8DSutBuI37WfyDcCtm64YwxaHZDm-rs,1095
+BESS_JPL/verify.py,sha256=RvGVvaQ68WOfA42KwdZCKOBnydBri4dQfpcYCwDA01M,5628
+BESS_JPL/version.py,sha256=oPRGJspBmrZU7gLv4L_3WPX3aGwee0inKgfzcyvkXcc,74
+bess_jpl-1.26.0.dist-info/licenses/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
+bess_jpl-1.26.0.dist-info/METADATA,sha256=TfruStH4jfszvGX0-wjbRnDKv6BIueGhK2PaThb0UFU,6543
+bess_jpl-1.26.0.dist-info/WHEEL,sha256=qELbo2s1Yzl39ZmrAibXA2jjPLUYfnVhUNTlyF1rq0Y,92
+bess_jpl-1.26.0.dist-info/top_level.txt,sha256=GaKnzt-BBktYn1o-w4Qzh_jHxse4Y3ACOxZQrB2ufhc,9
+bess_jpl-1.26.0.dist-info/RECORD,,

bess_jpl-1.26.0.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (80.10.1)
+Root-Is-Purelib: true
+Tag: py3-none-any
+