STIC-JPL 1.1.0__py3-none-any.whl

STIC_JPL/soil_moisture_initialization.py

@@ -0,0 +1,114 @@
+from typing import Union, Tuple
+import numpy as np
+
+import rasters as rt
+from rasters import Raster
+
+from .vegetation_conversion.vegetation_conversion import FVC_from_NDVI
+
+from .constants import GAMMA_HPA
+from .root_zone_initialization import calculate_root_zone_moisture
+
+def initialize_soil_moisture(
+        delta_hPa: Union[Raster, np.ndarray],  # Rate of change of saturation vapor pressure with temperature (hPa/°C)
+        ST_C: Union[Raster, np.ndarray],  # Surface temperature (°C)
+        Ta_C: Union[Raster, np.ndarray],  # Air temperature (°C)
+        Td_C: Union[Raster, np.ndarray],  # Dewpoint temperature (°C)
+        dTS: Union[Raster, np.ndarray],  # Difference between surface and air temperature (°C)
+        Rn_Wm2: Union[Raster, np.ndarray],  # Net radiation (W/m²)
+        LWnet_Wm2: Union[Raster, np.ndarray],  # Net longwave radiation (W/m²)
+        FVC: Union[Raster, np.ndarray],  # Fractional vegetation cover (unitless)
+        VPD_hPa: Union[Raster, np.ndarray],  # Vapor pressure deficit (hPa)
+        SVP_hPa: Union[Raster, np.ndarray],  # Saturation vapor pressure (hPa)
+        Ea_hPa: Union[Raster, np.ndarray],  # Actual vapor pressure (hPa)
+        Estar_hPa: Union[Raster, np.ndarray],  # Saturation vapor pressure at surface temperature (hPa)
+        gamma_hPa: Union[Raster, np.ndarray, float] = GAMMA_HPA  # Psychrometric constant (hPa/°C)
+    ) -> Tuple[Union[Raster, np.ndarray]]:
+    """
+    This function estimates the soil moisture availability (SM) based on thermal IR and meteorological information. 
+    The estimated SM is treated as initial SM, which will be later estimated through iteration in the actual ET estimation loop to establish feedback between SM and biophysical states.
+
+    Parameters:
+    delta (np.ndarray): Rate of change of saturation vapor pressure with temperature (kPa/°C)
+    ST_C (np.ndarray): Surface temperature (°C)
+    Ta_C (np.ndarray): Air temperature (°C)
+    Td_C (np.ndarray): Dewpoint temperature (°C)
+    dTS (np.ndarray): Difference between surface and air temperature (°C)
+    Rn (np.ndarray): Net radiation (W/m²)
+    LWnet (np.ndarray): Net longwave radiation (W/m²)
+    NDVI (np.ndarray): Normalized Difference Vegetation Index (unitless)
+    VPD_hPa (np.ndarray): Vapor pressure deficit (hPa)
+    SVP_hPa (np.ndarray): Saturation vapor pressure (hPa)
+    Ea_hPa (np.ndarray): Actual vapor pressure (hPa)
+    Estar (np.ndarray): Saturation vapor pressure at surface temperature (hPa)
+
+    Returns:
+    SM (np.ndarray): soil moisture (m³/m³)
+    Mrz (np.ndarray): The rootzone moisture (m³/m³)
+    Ms (np.ndarray): The surface moisture (m³/m³)
+    Ep_PT (np.ndarray): The potential evaporation (Priestley-Taylor eqn.) (mm/day)
+    Ds (np.ndarray): The vapor pressure deficit at surface (hPa)
+    s1 (np.ndarray): The slope of SVP at dewpoint temperature (hPa/K)
+    s3 (np.ndarray): The slope of SVP at surface temperature (hPa/K)
+    Tsd_C (np.ndarray): The surface dewpoint temperature (°C)
+    """
+    # Compute the surface dewpoint temperature
+    s11 = (45.03 + 3.014 * Td_C + 0.05345 * Td_C ** 2 + 0.00224 * Td_C ** 3) * 1e-2  # slope of SVP at TD (hpa/K)
+    s22 = (Estar_hPa - Ea_hPa) / (ST_C - Td_C)
+    s33 = (45.03 + 3.014 * ST_C + 0.05345 * ST_C ** 2 + 0.00224 * ST_C ** 3) * 1e-2  # slope of SVP at TS (hpa/K)
+    s44 = (SVP_hPa - Ea_hPa) / (Ta_C - Td_C)
+    Tsd_C = (Estar_hPa - Ea_hPa - s33 * ST_C + s11 * Td_C) / (s11 - s33)  # Surface dewpoint temperature (degC)
+
+    # Calculate the surface moisture or wetness
+    Msurf = (s11 / s22) * ((Tsd_C - Td_C) / (ST_C - Td_C))  # Surface wetness
+    Msurf = rt.clip(Msurf, 0.0001, 0.9999)
+
+    # Calculate the surface vapor pressure and deficit
+    esurf = Ea_hPa + Msurf * (Estar_hPa - Ea_hPa)
+    Dsurf = esurf - Ea_hPa
+
+    # Separate the soil and canopy wetness to form a composite surface moisture
+    Ms = Msurf
+    Mcan = FVC * Msurf
+    Msoil = (1 - FVC) * Msurf
+
+    TdewIndex = (ST_C - Tsd_C) / (Ta_C - Td_C)  # % TdewIndex > 1 signifies super dry condition
+    Ep_PT = (1.26 * delta_hPa * Rn_Wm2) / (delta_hPa + gamma_hPa)  # Potential evaporation (Priestley-Taylor eqn.)
+
+    # Adjust surface wetness based on certain conditions
+    Ms = rt.where((FVC <= 0.25) & (TdewIndex < 1), Msoil, Ms)
+    Mcan = rt.where((FVC <= 0.25) & (TdewIndex < 1), 0, Mcan)
+    Ms = rt.where((FVC <= 0.25) & (Ta_C > 10) & (Td_C < 0) & (LWnet_Wm2 < -125), Msoil, Ms)
+    Mcan = rt.where((FVC <= 0.25) & (Ta_C > 10) & (Td_C < 0) & (LWnet_Wm2 < -125), 0, Mcan)
+    
+    # Calculate the root-zone moisture
+    SMrz = calculate_root_zone_moisture(
+        delta_hPa=delta_hPa, 
+        ST_C=ST_C, 
+        Ta_C=Ta_C, 
+        Td_C=Td_C, 
+        s11=s11, 
+        s33=s33, 
+        s44=s44, 
+        Tsd_C=Tsd_C,
+        gamma_hPa=gamma_hPa
+    )
+
+    # Combine soil moisture to account for hysteresis and initial estimation of surface vapor pressure
+    SM = Ms
+    SM = rt.where((Ep_PT > Rn_Wm2) & (dTS > 0), SMrz, SM)
+    SM = rt.where((Ep_PT > Rn_Wm2) & (FVC <= 0.25), SMrz, SM)
+    SM = rt.where((Ep_PT > Rn_Wm2) & (Dsurf > VPD_hPa), SMrz, SM)
+    SM = rt.where((FVC <= 0.25) & (dTS > 0) & (Ta_C > 10) & (Td_C < 0) & (LWnet_Wm2 < -125), SMrz, SM)
+    SM = rt.where((FVC <= 0.25) & (dTS > 0) & (Ta_C > 10) & (Td_C < 0) & (Dsurf > VPD_hPa), SMrz, SM)
+    SM = rt.where((Ep_PT < Rn_Wm2) & (FVC <= 0.25) & (Dsurf > VPD_hPa), SMrz, SM)
+    
+    es = Ea_hPa + SM * (Estar_hPa - Ea_hPa)
+    
+    # vapor pressure deficit at surface
+    Ds = (Estar_hPa - es)  
+    
+    s1 = s11
+    s3 = s33
+
+    return SM, SMrz, Ms, Ep_PT, Ds, s1, s3, Tsd_C

STIC_JPL/soil_moisture_iteration.py

@@ -0,0 +1,131 @@
+from typing import Union
+import numpy as np
+
+import rasters as rt
+
+from rasters import Raster
+
+from .constants import GAMMA_HPA
+from .root_zone_iteration import calculate_rootzone_moisture
+
+def iterate_soil_moisture(
+        delta_hPa: Union[Raster, np.ndarray], # Rate of change of saturation vapor pressure with temperature (hPa/°C)
+        s1: Union[Raster, np.ndarray], # The slope of SVP at surface temperature (hPa/K)
+        s3: Union[Raster, np.ndarray], # The slope of SVP at dewpoint temperature (hPa/K)
+        ST_C: Union[Raster, np.ndarray], # surface temperature (°C)
+        Ta_C: Union[Raster, np.ndarray], # air temperature (°C)
+        dTS_C: Union[Raster, np.ndarray], # difference between surface and air temperature (°C)
+        Td_C: Union[Raster, np.ndarray], # dewpoint temperature (°C)
+        Tsd_C: Union[Raster, np.ndarray], # surface dewpoint temperature (°C)
+        Rg_Wm2: Union[Raster, np.ndarray], # Incoming solar radiation (W/m^2)
+        Rn_Wm2: Union[Raster, np.ndarray], # Net radiation (W/m^2)
+        LWnet_Wm2: Union[Raster, np.ndarray], # Net longwave radiation (W/m^2)
+        FVC: Union[Raster, np.ndarray], # Fractional vegetation cover (unitless)
+        VPD_hPa: Union[Raster, np.ndarray], # Vapor pressure deficit (hPa)
+        D0_hPa: Union[Raster, np.ndarray], # Vapor pressure deficit at source (hPa)
+        SVP_hPa: Union[Raster, np.ndarray], # Saturation vapor pressure (hPa)
+        Ea_hPa: Union[Raster, np.ndarray], # Actual vapor pressure (hPa)
+        T0: Union[Raster, np.ndarray], # Temperature at source (°C)
+        gamma_hPa: Union[Raster, np.ndarray, float] = GAMMA_HPA # Psychrometric constant (hPa/°C)
+    ) -> Union[Raster, np.ndarray]:
+    """
+    Estimates the soil moisture availability (SM) (or wetness) (value 0 to 1) based on thermal IR and meteorological
+    information. However, this M will be treated as initial M, which will be later on estimated through iteration in the
+    actual ET estimation loop to establish feedback between M and biophysical states.
+
+    Args:
+        delta (np.ndarray): Rate of change of saturation vapor pressure with temperature (kPa/°C).
+        s1 (np.ndarray): The slope of saturation vapor pressure at surface temperature (hPa/K).
+        s3 (np.ndarray): The slope of saturation vapor pressure at dewpoint temperature (hPa/K).
+        ST_C (np.ndarray): Surface temperature in degrees Celsius.
+        Ta_C (np.ndarray): Air temperature in degrees Celsius.
+        dTS (np.ndarray): Difference between surface and air temperature in degrees Celsius.
+        Td_C (np.ndarray): Dewpoint temperature in degrees Celsius.
+        Tsd_C (np.ndarray): Surface dewpoint temperature in degrees Celsius.
+        Rg (np.ndarray): Incoming solar radiation in W/m^2.
+        Rn (np.ndarray): Net radiation in W/m^2.
+        LWnet (np.ndarray): Net longwave radiation in W/m^2.
+        FVC (np.ndarray): Fractional vegetation cover (unitless).
+        VPD_hPa (np.ndarray): Vapor pressure deficit in hPa.
+        D0 (np.ndarray): Vapor pressure deficit at source in hPa.
+        SVP_hPa (np.ndarray): Saturation vapor pressure in hPa.
+        Ea_hPa (np.ndarray): Actual vapor pressure in hPa.
+        T0 (np.ndarray): Temperature at source in degrees Celsius.
+
+    Returns:
+        np.ndarray: Soil moisture availability (SM) (or wetness) (value 0 to 1).
+    """
+    # calculate surface wetness (Msurf)
+    kTSTD = (T0 - Td_C) / (ST_C - Td_C)
+    Msurf = (s1 / s3) * ((Tsd_C - Td_C) / (kTSTD * (ST_C - Td_C)))  # Surface wetness
+    Msurf = rt.where((Rn_Wm2 < 0) & (dTS_C < 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    Msurf = rt.where((Rn_Wm2 < 0) & (dTS_C > 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    Msurf = rt.where((Rn_Wm2 > 0) & (dTS_C < 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    Msurf = rt.where((Rn_Wm2 > 0) & (dTS_C > 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    
+    # solar radiation is only used to correct negative values of surface wetness
+    Msurf = rt.where((Rg_Wm2 > 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    Msurf = rt.where((Rg_Wm2 < 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    
+    Msurf = rt.where((Td_C < 0) & (Msurf < 0), np.abs(Msurf), Msurf)
+    Msurf = rt.clip(Msurf, 0.0001, 1)
+
+    # Separating soil and canopy wetness to form a composite surface moisture
+    Ms = Msurf
+    Mcan = FVC * Msurf
+    Msoil = (1 - FVC) * Msurf
+
+    TdewIndex = (ST_C - Tsd_C) / (Ta_C - Td_C)
+
+    # Potential evaporation (Priestley-Taylor eqn.)
+    Ep_PT = (1.26 * delta_hPa * Rn_Wm2) / (delta_hPa + gamma_hPa)  
+
+    # calculate surface wetness (Ms)
+    Ms = rt.where((Ep_PT > Rn_Wm2) & (FVC <= 0.25) & (TdewIndex < 1), Msoil, Ms)
+    Ms = rt.where((FVC <= 0.25) & (Ta_C > 10) & (Td_C < 0) & (LWnet_Wm2 < -125), Msoil, Ms)
+    Ms = rt.where((Rn_Wm2 > Ep_PT) & (FVC <= 0.25) & (TdewIndex < 1) & (Td_C <= 0), Msoil, Ms)
+    Ms = rt.where((Rn_Wm2 > Ep_PT) & (FVC <= 0.25) & (TdewIndex < 1), Msoil, Ms)
+    Ms = rt.where((D0_hPa > VPD_hPa) & (FVC <= 0.25) & (TdewIndex < 1), Msoil, Ms)
+
+    # calculate canopy wetness (Mcan)
+    Mcan = rt.where((Ep_PT > Rn_Wm2) & (FVC <= 0.25) & (TdewIndex < 1), 0, Mcan)
+    Mcan = rt.where((FVC <= 0.25) & (Ta_C > 10) & (Td_C < 0) & (LWnet_Wm2 < -125), 0, Mcan)
+    Mcan = rt.where((Rn_Wm2 > Ep_PT) & (FVC <= 0.25) & (TdewIndex < 1) & (Td_C <= 0), 0, Mcan)
+    Mcan = rt.where((Rn_Wm2 > Ep_PT) & (FVC <= 0.25) & (TdewIndex < 1), 0, Mcan)
+    Mcan = rt.where((D0_hPa > VPD_hPa) & (FVC <= 0.25) & (TdewIndex < 1), 0, Mcan)
+
+    # calculate rootzone moisture (Mrz)
+    Mrz = calculate_rootzone_moisture(
+        delta_hPa=delta_hPa, 
+        s1_hPa=s1, 
+        s3_hPa=s3, 
+        ST_C=ST_C, 
+        Ta_C=Ta_C, 
+        dTS_C=dTS_C, 
+        Td_C=Td_C, 
+        Tsd_C=Tsd_C, 
+        Rg_Wm2=Rg_Wm2, 
+        Rn_Wm2=Rn_Wm2, 
+        LWnet_Wm2=LWnet_Wm2, 
+        FVC=FVC, 
+        VPD_hPa=VPD_hPa, 
+        D0_hPa=D0_hPa, 
+        SVP_hPa=SVP_hPa, 
+        Ea_hPa=Ea_hPa,
+        gamma_hPa=gamma_hPa
+    )
+
+    TdewIndex = (ST_C - Tsd_C) / (Ta_C - Td_C)
+
+    # Potential evaporation (Priestley-Taylor eqn.)
+    Ep_PT = (1.26 * delta_hPa * Rn_Wm2) / (delta_hPa + gamma_hPa)  
+
+    # COMBINE M to account for Hysteresis and initial estimation of surface vapor pressure
+    SM = Msurf
+    SM = rt.where((Ep_PT > Rn_Wm2) & (dTS_C > 0) & (FVC <= 0.25) & (D0_hPa > VPD_hPa) & (TdewIndex < 1), Mrz, SM)
+    SM = rt.where((FVC <= 0.25) & (dTS_C > 0) & (Ta_C > 10) & (Td_C < 0) & (LWnet_Wm2 < -125) & (D0_hPa > VPD_hPa), Mrz, SM)
+
+    return SM
+
+
+

STIC_JPL/timer/__init__.py

@@ -0,0 +1 @@
+from .timer import *

STIC_JPL/timer/timer.py

@@ -0,0 +1,77 @@
+"""
+This is a minimalistic performance timer.
+"""
+import time
+
+__author__ = "Gregory Halverson"
+
+DEFAULT_FORMAT = "0.2f"
+
+class Timer(object):
+    """
+    This is a minimalistic performance timer.
+    """
+
+    def __init__(self):
+        self._start_time = None
+        self._end_time = None
+        self.start()
+
+    def __enter__(self, *args, **kwargs):
+        self.start()
+        return self
+
+    def __exit__(self, *args, **kwargs):
+        self.end()
+
+    def __repr__(self):
+        # print("Timer.__repr__")
+        return self.__format__(format_string=DEFAULT_FORMAT)
+
+    def __str__(self):
+        # print("Timer.__str__")
+        return self.__repr__()
+
+    def __format__(self, format_string=DEFAULT_FORMAT):
+        if format_string is None or format_string == "":
+            format_string = DEFAULT_FORMAT
+
+        return format(self.duration, format_string)
+
+    @property
+    def now(self):
+        # return datetime.now()
+        return time.perf_counter()
+
+    def start(self):
+        self._start_time = self.now
+
+        return self.start_time
+
+    @property
+    def start_time(self):
+        return self._start_time
+
+    def end(self):
+        self._end_time = self.now
+
+        return self.end_time
+
+    @property
+    def end_time(self):
+        return self._end_time
+
+    @property
+    def duration(self):
+        if self.start_time is None:
+            raise Exception("timer never started")
+
+        if self.end_time is None:
+            end_time = self.now
+        else:
+            end_time = self.end_time
+
+        duration = end_time - self.start_time
+
+        return duration
+

STIC_JPL/vegetation_conversion/__init__.py

@@ -0,0 +1 @@
+from .vegetation_conversion import *

STIC_JPL/vegetation_conversion/vegetation_conversion.py

@@ -0,0 +1,47 @@
+from typing import Union
+import numpy as np
+import rasters as rt
+from rasters import Raster
+
+KPAR = 0.5
+MIN_FIPAR = 0.0
+MAX_FIPAR = 1.0
+MIN_LAI = 0.0
+MAX_LAI = 10.0
+
+def FVC_from_NDVI(NDVI: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
+    """
+    Convert Normalized Difference Vegetation Index (NDVI) to Fractional Vegetation Cover (FVC).
+
+    Parameters:
+        NDVI (Union[Raster, np.ndarray]): Input NDVI data.
+
+    Returns:
+        Union[Raster, np.ndarray]: Converted FVC data.
+    """
+    NDVIv = 0.52  # +- 0.03
+    NDVIs = 0.04  # +- 0.03
+    FVC = rt.clip((NDVI - NDVIs) / (NDVIv - NDVIs), 0.0, 1.0)
+
+    return FVC
+
+def LAI_from_NDVI(
+        NDVI: Union[Raster, np.ndarray],
+        min_fIPAR: float = MIN_FIPAR,
+        max_fIPAR: float = MAX_FIPAR,
+        min_LAI: float = MIN_LAI,
+        max_LAI: float = MAX_LAI) -> Union[Raster, np.ndarray]:
+    """
+    Convert Normalized Difference Vegetation Index (NDVI) to Leaf Area Index (LAI).
+
+    Parameters:
+        NDVI (Union[Raster, np.ndarray]): Input NDVI data.
+
+    Returns:
+        Union[Raster, np.ndarray]: Converted LAI data.
+    """
+    fIPAR = rt.clip(NDVI - 0.05, min_fIPAR, max_fIPAR)
+    fIPAR = np.where(fIPAR == 0, np.nan, fIPAR)
+    LAI = rt.clip(-np.log(1 - fIPAR) * (1 / KPAR), min_LAI, max_LAI)
+
+    return LAI

STIC_JPL/version.txt

@@ -0,0 +1 @@
+1.0.4

stic_jpl-1.1.0.dist-info/METADATA

@@ -0,0 +1,90 @@
+Metadata-Version: 2.4
+Name: STIC-JPL
+Version: 1.1.0
+Summary: Surface Temperature Initiated Closure (STIC) Evapotranspiration Model Python Implementation
+Author-email: Gregory Halverson <gregory.h.halverson@jpl.nasa.gov>, Kaniska Mallick <kaniska.mallick@list.lu>, Madeleine Pascolini-Campbell <madeleine.a.pascolini-campbell@jpl.nasa.gov>, "Claire S. Villanueva-Weeks" <claire.s.villanueva-weeks@jpl.gov>
+Project-URL: Homepage, https://github.com/JPL-Evapotranspiration-Algorithms/STIC-JPL
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+Requires-Dist: colored-logging
+Requires-Dist: ECOv002-CMR
+Requires-Dist: ECOv002-granules
+Requires-Dist: GEOS5FP>=1.1.1
+Requires-Dist: numpy
+Requires-Dist: pandas
+Requires-Dist: rasters
+Requires-Dist: solar-apparent-time
+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"
+
+# Surface Temperature Initiated Closure (STIC) Evapotranspiration Model Python Implementation
+
+[![CI](https://github.com/JPL-Evapotranspiration-Algorithms/STIC/actions/workflows/ci.yml/badge.svg)](https://github.com/JPL-Evapotranspiration-Algorithms/STIC/actions/workflows/ci.yml)
+
+This repository contains the python implementation for the Surface Temperature Initiated Closure (STIC) evapotranspiration model, used by the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) and Surface Biology and Geology (SBG) missions.
+
+This software package is a Python implementation of the Surface Temperature Initiated Closure (STIC) version 1.3 model designed to implement LST to solve the aerodynamic temperature, which is critical for ET modeling. The original STIC model produced by Kaniska Mallick (Mallick et al. 2015; 2018; 2022) was re-implemented from MATLAB code to Python by Gregory Halverson and Madeleine Pascolini-Campbell. The software was developed under a research grant by the NASA Research Opportunities in Space and Earth Sciences (ROSES) program. It is intended for use by the Hyperspectral Thermal Emission Spectrometer (HyTES), MODIS/ASTER (MASTER) Airborne Simulator, Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission and as a precursor for the Surface Biology and Geology (SBG) mission. 
+
+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 Hyperspectral Thermal Emission Spectrometer (HyTES), MODIS/ASTER (MASTER) Airborne Simulator, Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission as a precursor for the Surface Biology and Geology (SBG) mission. 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.
+ 
+The 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
+
+Kaniska Mallick (he/him)<br>
+[kaniska.mallick@list.lu](mailto:kaniska.mallick@list.lu)<br>
+Algorithm inventor<br>
+Luxembourg Institute of Science and Technology
+
+Tian Hu (he/him)<br>
+[tian.hu@list.lu](mailto:tian.hu@list.lu)<br>
+Algorithm developer<br>
+Luxembourg Institute of Science and Technology
+
+Madeleine Pascolini-Campbell (she/her)<br>
+[madeleine.a.pascolini-campbell@jpl.nasa.gov](mailto:madeleine.a.pascolini-campbell@jpl.nasa.gov)<br>
+Algorithm developer<br>
+NASA Jet Propulsion Laboratory 329F
+
+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 the `STIC` PyPi package.
+
+```
+pip install STIC
+```
+
+## Usage
+
+Import the `STIC` function from the `STIC` package.
+
+```
+from STIC import STIC
+```
+
+See the [ECOSTRESS example](ECOSTRESS%20Example.ipynb) for usage.
+
+## References
+ 
+Mallick, K., Boegh, E., Trebs, I., Alfieri, J. G., Kustas, W. P., Prueger, J. H., ... & Jarvis, A. J. (2015). Reintroducing radiometric surface temperature into the P enman‐M onteith formulation. Water Resources Research, 51(8), 6214-6243. https://doi.org/10.1002/2014WR016106 
+
+Mallick, K., Toivonen, E., Trebs, I., Boegh, E., Cleverly, J., Eamus, D., ... & Garcia, M. (2018). Bridging Thermal Infrared Sensing and Physically‐Based Evapotranspiration Modeling: From Theoretical Implementation to Validation Across an Aridity Gradient in Australian Ecosystems. Water Resources Research, 54(5), 3409-3435. https://doi.org/10.1029/2017WR021357
+
+Mallick, K., Baldocchi, D., Jarvis, A., Hu, T., Trebs, I., Sulis, M., ... & Kustas, W. P. (2022). Insights Into the Aerodynamic Versus Radiometric Surface Temperature Debate in Thermal‐Based Evaporation Modeling. Geophysical Research Letters, 49(15), e2021GL097568. https://doi.org/10.1029/2021GL097568
+
+

stic_jpl-1.1.0.dist-info/RECORD

@@ -0,0 +1,28 @@
+STIC_JPL/STIC_JPL.py,sha256=a0EpRTX6Sd75RPIqAzWckIuj9-KwbCE8T1aSEDpLM-c,17434
+STIC_JPL/__init__.py,sha256=FmTp-Ir0Wbae0SbZgxlNIxq48TDrWCev17Gv1hWskxU,216
+STIC_JPL/canopy_air_stream.py,sha256=UYp3l7mt0XH4SqvCGZhRzCIL-v7-Rvxxfm-C81lheSU,1637
+STIC_JPL/closure.py,sha256=AzWoLnpoZ1MlQAKq0CgpmCMfDGVsp8-hYZ50ry33cic,4031
+STIC_JPL/constants.py,sha256=xIVTsdc3_Rf3EPSgSk1QPjWiojFGn4lf1Swffov3oTM,595
+STIC_JPL/diagnostic.py,sha256=RXjtZbCEJMiTyRuiMGUvekDgO_GcrhXVXvorR0WlAsk,2613
+STIC_JPL/initialize_with_solar.py,sha256=6DUUHdVfeKvDk1pf_bOteSW8qd-n9pU4dowO1sMYOpw,3825
+STIC_JPL/initialize_without_solar.py,sha256=M_7QfP2pr2w5sQO7UOW1GPqxpmQkjmsQ526BgYVVALk,4675
+STIC_JPL/iterate_with_solar.py,sha256=WKySvuIJ4HzcDni0UYafmk5lXex2lJfde82kGnGUDI8,6368
+STIC_JPL/iterate_without_solar.py,sha256=09RERpN6tdZr6ORDhMwEFXwx8941vyn1HMJQ-Q_bY6s,6140
+STIC_JPL/net_radiation.py,sha256=Uwudsazul8V-x5t8KQLi3wpYi-wjMwT__eqZCRV2bIw,1187
+STIC_JPL/root_zone_initialization.py,sha256=3JVKNDt3ebIiGVfuhBawjKs-BicNwkbfMprkzSxd4Cg,1581
+STIC_JPL/root_zone_iteration.py,sha256=1XOMFE3-TdJZzUU5vouflUO4NKaZwAW6s1KJrNez3Es,3787
+STIC_JPL/soil_moisture_initialization.py,sha256=YF0GcHkpH5amABNPpxD8f1lLPRyIAkG2dYJTfEuIX4g,5599
+STIC_JPL/soil_moisture_iteration.py,sha256=QJXOPMxxwIIskpx9zLkXUPfuhWgFPUBcRnGhZo2UjAw,6493
+STIC_JPL/version.txt,sha256=rqWtvvV0eFJzxPOcG3aHz4AZk-DLa_Z4GkFonk7bsgw,5
+STIC_JPL/meteorology_conversion/__init__.py,sha256=dLF40imqJCC_l8_QMetbsj1YzkXZuFt9tQ2fsd7Wcdw,38
+STIC_JPL/meteorology_conversion/meteorology_conversion.py,sha256=8vzx71FLLNhR5v6EQBL8OSgAmupQcgdLX2c9tL9TPxA,4202
+STIC_JPL/soil_heat_flux/__init__.py,sha256=bN0dHUfqR9ob68M1TqdjQTK3kmaL5S4P08yK0Thub6k,74
+STIC_JPL/soil_heat_flux/calculate_SEBAL_soil_heat_flux.py,sha256=PZylUGET5-Uj6BGfhOeqS4tp2pGE__gdc5X1nidLjxw,1516
+STIC_JPL/timer/__init__.py,sha256=I_MQKp_aamBLUzZv0psEbRgs6GZLOJd4mmJ7bli0Ikc,21
+STIC_JPL/timer/timer.py,sha256=tn5e3NQmsh55Jp9Fstjf-8KJW4F8UIJs-d_ZLooFYE8,1610
+STIC_JPL/vegetation_conversion/__init__.py,sha256=8gnz0yK_euCV0TtVGDrSGaVfx4g0EJCG2J68ppuA5oc,37
+STIC_JPL/vegetation_conversion/vegetation_conversion.py,sha256=-K3FLV6pXkpxtRxQvimqqkXAgkn9mUUlHQ8FIscq1Zg,1306
+stic_jpl-1.1.0.dist-info/METADATA,sha256=sDO9vRDdU6d-qTxHfk-0ZoCLX5ZyT1o1F086M0hbqWA,5348
+stic_jpl-1.1.0.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+stic_jpl-1.1.0.dist-info/top_level.txt,sha256=9NkchxttzACJcGcAaWzMaZarzX40OXQ216hERNA9LIo,9
+stic_jpl-1.1.0.dist-info/RECORD,,

stic_jpl-1.1.0.dist-info/WHEEL

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

stic_jpl-1.1.0.dist-info/top_level.txt

@@ -0,0 +1 @@
+STIC_JPL