PM-JPL 1.4.0__py3-none-any.whl → 1.5.1__py3-none-any.whl

PMJPL/process_daily_ET_table.py

@@ -0,0 +1,40 @@
+import pandas as pd
+
+from sun_angles import SHA_deg_from_DOY_lat, daylight_from_SHA, sunrise_from_SHA
+from verma_net_radiation import daily_Rn_integration_verma
+from daily_evapotranspiration_upscaling import daily_ET_from_daily_LE
+
+from meteorology_conversion import celcius_to_kelvin
+
+def process_daily_ET_table(input_df: pd.DataFrame) -> pd.DataFrame:
+    hour_of_day = input_df.hour_of_day
+    DOY = input_df.doy
+    lat = input_df.lat
+    LE = input_df.LE
+    Rn = input_df.Rn
+    EF = LE / Rn
+
+    SHA_deg = SHA_deg_from_DOY_lat(DOY=DOY, latitude=lat)
+    sunrise_hour = sunrise_from_SHA(SHA_deg)
+    daylight_hours = daylight_from_SHA(SHA_deg)
+
+    Rn_daylight = daily_Rn_integration_verma(
+        Rn=Rn,
+        hour_of_day=hour_of_day,
+        DOY=DOY,
+        lat=lat,
+        sunrise_hour=sunrise_hour,
+        daylight_hours=daylight_hours
+    )
+
+    LE_daylight = EF * Rn_daylight
+    ET = daily_ET_from_daily_LE(LE_daylight, daylight_hours)
+
+    output_df = input_df.copy()
+    output_df["EF"] = EF
+    output_df["sunrise_hour"] = sunrise_hour
+    output_df["daylight_hours"] = daylight_hours
+    output_df["Rn_daylight"] = Rn_daylight
+    output_df["ET"] = ET
+
+    return output_df

PMJPL/version.py

@@ -0,0 +1,4 @@
+from os.path import join, abspath, dirname
+
+from importlib.metadata import version
+__version__ = version("PM-JPL")

{pm_jpl-1.4.0.dist-info → pm_jpl-1.5.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PM-JPL
-Version: 1.4.0
+Version: 1.5.1
 Summary: JPL implementation of the MOD16 evapotranspiration algorithm for high resolution instantaneous remote sensing imagery
 Author-email: Gregory Halverson <gregory.h.halverson@jpl.nasa.gov>, Qiaozhen Mu <qiaozhen@ntsg.umt.edu>, Maosheng Zhao <zhao@ntsg.umt.edu>, "Steven W. Running" <swr@ntsg.umt.edu>, "Claire S. Villanueva-Weeks" <claire.s.villanueva-weeks@jpl.gov>
 Project-URL: Homepage, https://github.com/JPL-Evapotranspiration-Algorithms/PM-JPL
@@ -9,18 +9,24 @@ Classifier: Operating System :: OS Independent
 Requires-Python: >=3.10
 Description-Content-Type: text/markdown
 License-File: LICENSE
+Requires-Dist: carlson-fractional-vegetation-cover
+Requires-Dist: carlson-leaf-area-index
 Requires-Dist: check-distribution
+Requires-Dist: daily-evapotranspiration-upscaling
 Requires-Dist: ECOv002-CMR>=1.0.5
 Requires-Dist: ECOv002-granules>=1.0.3
 Requires-Dist: ECOv003-granules
 Requires-Dist: GEOS5FP>=1.1.1
 Requires-Dist: MCD12C1_2019_v006
+Requires-Dist: meteorology-conversion
 Requires-Dist: NASADEM
 Requires-Dist: numpy
 Requires-Dist: pandas
+Requires-Dist: priestley-taylor
+Requires-Dist: PTJPL>=1.5.1
 Requires-Dist: rasters>=1.4.6
 Requires-Dist: SEBAL-soil-heat-flux
-Requires-Dist: sun-angles
+Requires-Dist: sun-angles>=1.3.0
 Requires-Dist: verma-net-radiation
 Provides-Extra: dev
 Requires-Dist: build; extra == "dev"

pm_jpl-1.5.1.dist-info/RECORD

@@ -0,0 +1,24 @@
+PMJPL/PMJPL.py,sha256=Hi2zYpOF0Fk2VUbZuQ31uYDYFx1W92t5zxyGcy3E-h8,16079
+PMJPL/VPD_factor.py,sha256=ezHMnnMPGNBt6ZTyiwsTEyMr8fwre6JDsHQea68P7T0,970
+PMJPL/__init__.py,sha256=jke4us8XYaXqU4fL97WYb7Ot1kdD4hV_2NyqFchUh6Q,91
+PMJPL/calculate_gamma.py,sha256=cB7pnbhJM4nArjmtg-gvGR6aj1NvXMOlG2RNZ8oLPFs,5907
+PMJPL/canopy_aerodynamic_resistance.py,sha256=eU5x8riEROyynsMiWDHR_sDCeDexsWDlhuU69tUvC2s,8490
+PMJPL/canopy_conductance.py,sha256=ow7Egm_2wu2Zu5WxkESaUABf7UQA0ldmyPe5Wf4FQ2Y,4814
+PMJPL/constants.py,sha256=T5Da1IzPGM6_E6tBIm2s3CXtwUYP0XXYCHyiPuaLknw,786
+PMJPL/correctance_factor.py,sha256=bVsu23oeTJoeMe-SIsH6cdFyATsLhgwvQUz79yf6kcQ,3181
+PMJPL/interception.py,sha256=2kDqMefH_kzMO9mLVY9BPFSoJE4K0pkQCXwrSi2CR9c,1605
+PMJPL/mod16.csv,sha256=1aVfqLOAJwYuZPaScyqbk4yacObPgp8Q-IPzJm0yReY,847
+PMJPL/parameters.py,sha256=kM26BUDMi-N03Gi2PqQi48kYMKakCB3ReiRMUGrvRQU,1974
+PMJPL/potential_soil_evaporation.py,sha256=O8bHGM4ARqSXJDHu6IWneSEy5qScenotGR3uawijGWg,1829
+PMJPL/process_daily_ET_table.py,sha256=dbKp8LaTXnWaj8gOLueX_AvQpvVVb9OWfz5RqC-kYrc,1197
+PMJPL/soil_moisture_constraint.py,sha256=URUGsE1_p2P1dVQ8qKRcKb93hahYWtAmYqoassRu-PI,667
+PMJPL/tmin_factor.py,sha256=UEeET10ErR2rUr1ViRN8UZmoCACxeeiFAReZXr9Kts4,1595
+PMJPL/transpiration.py,sha256=lhRm2A0yxrlJxDw2he-k3PhE7G-sBsVvAMKo5fmh7wU,1998
+PMJPL/version.py,sha256=7e7k9a2pQYZR4Fyhl6UDY8NQ1lFi5ChR0Zk9Mpl45BU,115
+PMJPL/wet_canopy_resistance.py,sha256=Eape5EBJuoP-IRvD0SWHJKCJWhzm03vcG72OZVYqt_Q,879
+PMJPL/wet_soil_evaporation.py,sha256=mw9JG1oXfVdSWnVRLP2SPd53mSRK96b8k3E1zf0ci7A,1440
+pm_jpl-1.5.1.dist-info/licenses/LICENSE,sha256=xx0jnfkXJvxRnG63LTGOxlggYnIysveWIZ6H3PNdCrQ,11357
+pm_jpl-1.5.1.dist-info/METADATA,sha256=-QtcUy-InQc4tDHH58GFGte0PAvrz5hlnoIbVu0EdI4,4391
+pm_jpl-1.5.1.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+pm_jpl-1.5.1.dist-info/top_level.txt,sha256=YaAFwdYHUxfUW08hiuFquUEdDGrsYn1duuMMjJKh0Ws,6
+pm_jpl-1.5.1.dist-info/RECORD,,

{pm_jpl-1.4.0.dist-info → pm_jpl-1.5.1.dist-info}/WHEEL

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

PMJPL/downscaling/__init__.py

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

PMJPL/downscaling/downscaling.py

@@ -1,271 +0,0 @@
-from datetime import datetime
-from dateutil import parser
-import numpy as np
-from rasters import Raster, RasterGeometry, RasterGrid
-import rasters as rt
-
-DEFAULT_UPSAMPLING = "average"
-DEFAULT_DOWNSAMPLING = "linear"
-
-def bias_correct(
-        coarse_image: Raster,
-        fine_image: Raster,
-        upsampling: str = "average",
-        downsampling: str = "linear",
-        return_bias: bool = False):
-    fine_geometry = fine_image.geometry
-    coarse_geometry = coarse_image.geometry
-    upsampled = fine_image.to_geometry(coarse_geometry, resampling=upsampling)
-    bias_coarse = upsampled - coarse_image
-    bias_fine = bias_coarse.to_geometry(fine_geometry, resampling=downsampling)
-    bias_corrected_fine = fine_image - bias_fine
-
-    if return_bias:
-        return bias_corrected_fine, bias_fine
-    else:
-        return bias_corrected_fine
-
-def linear_downscale(
-        coarse_image: Raster,
-        fine_image: Raster,
-        upsampling: str = "average",
-        downsampling: str = "linear",
-        use_gap_filling: bool = False,
-        apply_scale: bool = True,
-        apply_bias: bool = True,
-        return_scale_and_bias: bool = False) -> Raster:
-    if upsampling is None:
-        upsampling = DEFAULT_UPSAMPLING
-
-    if downsampling is None:
-        downsampling = DEFAULT_DOWNSAMPLING
-
-    coarse_geometry = coarse_image.geometry
-    fine_geometry = fine_image.geometry
-    upsampled = fine_image.to_geometry(coarse_geometry, resampling=upsampling)
-
-    if apply_scale:
-        scale_coarse = coarse_image / upsampled
-        scale_coarse = rt.where(coarse_image == 0, 0, scale_coarse)
-        scale_coarse = rt.where(upsampled == 0, 0, scale_coarse)
-        scale_fine = scale_coarse.to_geometry(fine_geometry, resampling=downsampling)
-        scale_corrected_fine = fine_image * scale_fine
-        fine_image = scale_corrected_fine
-    else:
-        scale_fine = fine_image * 0 + 1
-
-    if apply_bias:
-        upsampled = fine_image.to_geometry(coarse_geometry, resampling=upsampling)
-        bias_coarse = upsampled - coarse_image
-        bias_fine = bias_coarse.to_geometry(fine_geometry, resampling=downsampling)
-        bias_corrected_fine = fine_image - bias_fine
-        fine_image = bias_corrected_fine
-    else:
-        bias_fine = fine_image * 0
-
-    if use_gap_filling:
-        gap_fill = coarse_image.to_geometry(fine_geometry, resampling=downsampling)
-        fine_image = fine_image.fill(gap_fill)
-
-    if return_scale_and_bias:
-        fine_image["scale"] = scale_fine
-        fine_image["bias"] = bias_fine
-
-    return fine_image
-
-def NDVI_to_FVC(NDVI: Raster) -> Raster:
-    NDVIv = 0.52  # +- 0.03
-    NDVIs = 0.04  # +- 0.03
-    FVC = rt.clip((NDVI - NDVIs) / (NDVIv - NDVIs), 0, 1)
-
-    return FVC
-
-def downscale_air_temperature(
-        time_UTC: datetime,
-        Ta_K_coarse: Raster,
-        ST_K: Raster,
-        water: Raster = None,
-        fine_geometry: RasterGeometry = None,
-        coarse_geometry: RasterGeometry = None,
-        resampling: str = None,
-        upsampling: str = None,
-        downsampling: str = None,
-        apply_scale: bool = True,
-        apply_bias: bool = True,
-        return_scale_and_bias: bool = False) -> Raster:
-    """
-    near-surface air temperature (Ta) in Kelvin
-    :param time_UTC: date/time in UTC
-    :param geometry: optional target geometry
-    :param resampling: optional sampling method for resampling to target geometry
-    :return: raster of Ta
-    """
-
-    if isinstance(time_UTC, str):
-        time_UTC = parser.parse(time_UTC)
-
-    if fine_geometry is None:
-        fine_geometry = ST_K.geometry
-
-    if coarse_geometry is None:
-        coarse_geometry = Ta_K_coarse.geometry
-
-    ST_K_water = None
-
-    if water is not None:
-        ST_K_water = rt.where(water, ST_K, np.nan)
-        ST_K = rt.where(water, np.nan, ST_K)
-
-    scale = None
-    bias = None
-
-    Ta_K = linear_downscale(
-        coarse_image=Ta_K_coarse,
-        fine_image=ST_K,
-        upsampling=upsampling,
-        downsampling=downsampling,
-        apply_scale=apply_scale,
-        apply_bias=apply_bias,
-        return_scale_and_bias=return_scale_and_bias
-    )
-
-    if water is not None:
-        Ta_K_water = linear_downscale(
-            coarse_image=Ta_K_coarse,
-            fine_image=ST_K_water,
-            upsampling=upsampling,
-            downsampling=downsampling,
-            apply_scale=apply_scale,
-            apply_bias=apply_bias,
-            return_scale_and_bias=False
-        )
-
-        Ta_K = rt.where(water, Ta_K_water, Ta_K)
-
-    Ta_K.filenames = Ta_K_coarse.filenames
-
-    return Ta_K
-
-
-def downscale_soil_moisture(
-        time_UTC: datetime,
-        fine_geometry: RasterGrid,
-        coarse_geometry: RasterGrid,
-        SM_coarse: Raster,
-        SM_resampled: Raster,
-        ST_fine: Raster,
-        NDVI_fine: Raster,
-        water: Raster,
-        fvlim=0.5,
-        a=0.5,
-        smoothing="linear") -> Raster:
-    fine = fine_geometry
-    ST_fine = ST_fine.mask(~water)
-    NDVI_fine = NDVI_fine.mask(~water)
-    FVC_fine = NDVI_to_FVC(NDVI_fine)
-    soil_fine = FVC_fine < fvlim
-    Tmin_coarse = ST_fine.to_geometry(coarse_geometry, resampling="min")
-    Tmax_coarse = ST_fine.to_geometry(coarse_geometry, resampling="max")
-    Ts_fine = ST_fine.mask(soil_fine)
-    Tsmin_coarse = Ts_fine.to_geometry(coarse_geometry, resampling="min").fill(Tmin_coarse)
-    Tsmax_coarse = Ts_fine.to_geometry(coarse_geometry, resampling="max").fill(Tmax_coarse)
-    ST_coarse = ST_fine.to_geometry(coarse_geometry, resampling="average")
-    SEE_coarse = (Tsmax_coarse - ST_coarse) / rt.clip(Tsmax_coarse - Tsmin_coarse, 1, None)
-    SM_SEE_proportion = (SM_coarse / SEE_coarse).to_geometry(fine, resampling=smoothing)
-    Tsmax_fine = Tsmax_coarse.to_geometry(fine_geometry, resampling=smoothing)
-    Tsrange_fine = (Tsmax_coarse - Tsmin_coarse).to_geometry(fine, resampling=smoothing)
-    SEE_fine = (Tsmax_fine - ST_fine) / rt.clip(Tsrange_fine, 1, None)
-
-    SEE_mean = SEE_coarse.to_geometry(fine, resampling=smoothing)
-    SM_fine = rt.clip(SM_resampled + a * SM_SEE_proportion * (SEE_fine - SEE_mean), 0, 1)
-    SM_fine = SM_fine.mask(~water)
-
-    return SM_fine
-
-
-def downscale_vapor_pressure_deficit(
-        time_UTC: datetime,
-        VPD_Pa_coarse: Raster,
-        ST_K: Raster,
-        fine_geometry: RasterGeometry = None,
-        coarse_geometry: RasterGeometry = None,
-        resampling: str = None,
-        upsampling: str = None,
-        downsampling: str = None,
-        return_scale_and_bias: bool = False) -> Raster:
-    if upsampling is None:
-        upsampling = "average"
-
-    if downsampling is None:
-        downsampling = "linear"
-
-    if fine_geometry is None:
-        fine_geometry = ST_K.geometry
-
-    if coarse_geometry is None:
-        coarse_geometry = VPD_Pa_coarse.geometry
-
-    return linear_downscale(
-        coarse_image=VPD_Pa_coarse,
-        fine_image=ST_K,
-        upsampling=upsampling,
-        downsampling=downsampling,
-        return_scale_and_bias=return_scale_and_bias
-    )
-
-
-def downscale_relative_humidity(
-        time_UTC: datetime,
-        RH_coarse: Raster,
-        SM,
-        ST_K,
-        VPD_kPa,
-        water: Raster = None,
-        fine_geometry: RasterGeometry = None,
-        coarse_geometry: RasterGeometry = None,
-        resampling: str = None,
-        upsampling: str = None,
-        downsampling: str = None) -> Raster:
-    if upsampling is None:
-        upsampling = "average"
-
-    if downsampling is None:
-        downsampling = "linear"
-
-    if fine_geometry is None:
-        fine_geometry = SM.geometry
-
-    if coarse_geometry is None:
-        coarse_geometry = RH_coarse.geometry
-
-    bias_fine = None
-
-    RH_estimate_fine = SM ** (1 / VPD_kPa)
-
-    RH = bias_correct(
-        coarse_image=RH_coarse,
-        fine_image=RH_estimate_fine,
-        upsampling=upsampling,
-        downsampling=downsampling,
-        return_bias=False
-    )
-
-    if water is not None:
-        ST_K_water = rt.where(water, ST_K, np.nan)
-        RH_coarse_complement = 1 - RH_coarse
-        RH_complement_water = linear_downscale(
-            coarse_image=RH_coarse_complement,
-            fine_image=ST_K_water,
-            upsampling=upsampling,
-            downsampling=downsampling,
-            apply_bias=True,
-            return_scale_and_bias=False
-        )
-
-        RH_water = 1 - RH_complement_water
-        RH = rt.where(water, RH_water, RH)
-
-    RH = rt.clip(RH, 0, 1)
-
-    return RH
-

PMJPL/downscaling/linear_downscale.py

@@ -1,71 +0,0 @@
-from rasters import Raster
-import rasters as rt
-
-DEFAULT_UPSAMPLING = "average"
-DEFAULT_DOWNSAMPLING = "linear"
-
-def bias_correct(
-        coarse_image: Raster,
-        fine_image: Raster,
-        upsampling: str = "average",
-        downsampling: str = "linear",
-        return_bias: bool = False):
-    fine_geometry = fine_image.geometry
-    coarse_geometry = coarse_image.geometry
-    upsampled = fine_image.to_geometry(coarse_geometry, resampling=upsampling)
-    bias_coarse = upsampled - coarse_image
-    bias_fine = bias_coarse.to_geometry(fine_geometry, resampling=downsampling)
-    bias_corrected_fine = fine_image - bias_fine
-
-    if return_bias:
-        return bias_corrected_fine, bias_fine
-    else:
-        return bias_corrected_fine
-
-def linear_downscale(
-        coarse_image: Raster,
-        fine_image: Raster,
-        upsampling: str = "average",
-        downsampling: str = "cubic",
-        use_gap_filling: bool = False,
-        apply_scale: bool = True,
-        apply_bias: bool = True,
-        return_scale_and_bias: bool = False) -> Raster:
-    if upsampling is None:
-        upsampling = DEFAULT_UPSAMPLING
-
-    if downsampling is None:
-        downsampling = DEFAULT_DOWNSAMPLING
-
-    coarse_geometry = coarse_image.geometry
-    fine_geometry = fine_image.geometry
-    upsampled = fine_image.to_geometry(coarse_geometry, resampling=upsampling)
-
-    if apply_scale:
-        scale_coarse = coarse_image / upsampled
-        scale_coarse = rt.where(coarse_image == 0, 0, scale_coarse)
-        scale_coarse = rt.where(upsampled == 0, 0, scale_coarse)
-        scale_fine = scale_coarse.to_geometry(fine_geometry, resampling=downsampling)
-        scale_corrected_fine = fine_image * scale_fine
-        fine_image = scale_corrected_fine
-    else:
-        scale_fine = fine_image * 0 + 1
-
-    if apply_bias:
-        upsampled = fine_image.to_geometry(coarse_geometry, resampling=upsampling)
-        bias_coarse = upsampled - coarse_image
-        bias_fine = bias_coarse.to_geometry(fine_geometry, resampling=downsampling)
-        bias_corrected_fine = fine_image - bias_fine
-        fine_image = bias_corrected_fine
-    else:
-        bias_fine = fine_image * 0
-
-    if use_gap_filling:
-        gap_fill = coarse_image.to_geometry(fine_geometry, resampling=downsampling)
-        fine_image = fine_image.fill(gap_fill)
-
-    if return_scale_and_bias:
-        fine_image["scale"] = scale_fine
-        fine_image["bias"] = bias_fine
-
-    return fine_image

PMJPL/evapotranspiration_conversion/__init__.py

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

PMJPL/evapotranspiration_conversion/evapotranspiration_conversion.py

@@ -1,82 +0,0 @@
-from typing import Union
-from sun_angles import SHA_deg_from_DOY_lat, daylight_from_SHA, sunrise_from_SHA
-
-import rasters as rt
-from rasters import Raster
-import numpy as np
-import pandas as pd
-
-from verma_net_radiation import daily_Rn_integration_verma
-
-from ..meteorology_conversion.meteorology_conversion import celcius_to_kelvin
-
-
-# latent heat of vaporization for water at 20 Celsius in Joules per kilogram
-LAMBDA_JKG_WATER_20C = 2450000.0
-
-def lambda_Jkg_from_Ta_K(Ta_K: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    # Calculate the latent heat of vaporization (J kg-1)
-    return (2.501 - 0.002361 * (Ta_K - 273.15)) * 1e6
-
-def lambda_Jkg_from_Ta_C(Ta_C: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    Ta_K = celcius_to_kelvin(Ta_C)
-    lambda_Jkg = lambda_Jkg_from_Ta_K(Ta_K)
-
-    return lambda_Jkg
-
-def daily_ET_from_daily_LE(
-        LE_daylight: Union[Raster, np.ndarray], 
-        daylight_hours: Union[Raster, np.ndarray],
-        lambda_Jkg: float = LAMBDA_JKG_WATER_20C) -> Union[Raster, np.ndarray]:
-    """
-    Calculate daily evapotranspiration (ET) from daily latent heat flux (LE).
-
-    Parameters:
-        LE_daily (Union[Raster, np.ndarray]): Daily latent heat flux.
-        daylight_hours (Union[Raster, np.ndarray]): Length of day in hours.
-        latent_vaporization (float, optional): Latent heat of vaporization. Defaults to LATENT_VAPORIZATION.
-
-    Returns:
-        Union[Raster, np.ndarray]: Daily evapotranspiration in kilograms.
-    """
-    # convert length of day in hours to seconds
-    daylight_seconds = daylight_hours * 3600.0
-
-    # factor seconds out of watts to get joules and divide by latent heat of vaporization to get kilograms
-    ET_daily_kg = rt.clip(LE_daylight * daylight_seconds / LAMBDA_JKG_WATER_20C, 0.0, None)
-
-    return ET_daily_kg
-
-def process_daily_ET_table(input_df: pd.DataFrame) -> pd.DataFrame:
-    hour_of_day = input_df.hour_of_day
-    DOY = input_df.doy
-    lat = input_df.lat
-    LE = input_df.LE
-    Rn = input_df.Rn
-    EF = LE / Rn
-
-    SHA_deg = SHA_deg_from_DOY_lat(DOY=DOY, latitude=lat)
-    sunrise_hour = sunrise_from_SHA(SHA_deg)
-    daylight_hours = daylight_from_SHA(SHA_deg)
-
-    Rn_daylight = daily_Rn_integration_verma(
-        Rn=Rn,
-        hour_of_day=hour_of_day,
-        DOY=DOY,
-        lat=lat,
-        sunrise_hour=sunrise_hour,
-        daylight_hours=daylight_hours
-    )
-
-    LE_daylight = EF * Rn_daylight
-    ET = daily_ET_from_daily_LE(LE_daylight, daylight_hours)
-
-    output_df = input_df.copy()
-    output_df["EF"] = EF
-    output_df["sunrise_hour"] = sunrise_hour
-    output_df["daylight_hours"] = daylight_hours
-    output_df["Rn_daylight"] = Rn_daylight
-    output_df["ET"] = ET
-
-    return output_df
-    

PMJPL/fwet.py

@@ -1,21 +0,0 @@
-import numpy as np
-from typing import Union
-from rasters import Raster
-
-from .constants import RH_THRESHOLD, MIN_FWET
-
-def calculate_fwet(
-        RH: Union[Raster, np.ndarray], 
-        RH_threshold: float = RH_THRESHOLD, 
-        min_fwet: float = MIN_FWET) -> Union[Raster, np.ndarray]:
-    """
-    calculates relative surface wetness
-    :param RH: relative humdity from 0.0 to 1.0
-    :return: relative surface wetness from 0.0 to 1.0
-    """
-    fwet = np.float32(np.clip(RH ** 4.0, min_fwet, None))
-
-    if RH_threshold is not None:
-        fwet = np.where(RH < RH_threshold, min_fwet, fwet)
-
-    return fwet

PMJPL/meteorology_conversion/__init__.py

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

PMJPL/meteorology_conversion/meteorology_conversion.py

@@ -1,123 +0,0 @@
-from typing import Union
-import numpy as np
-import rasters as rt
-from rasters import Raster
-
-# gas constant for dry air in joules per kilogram per kelvin
-RD = 286.9
-
-# gas constant for moist air in joules per kilogram per kelvin
-RW = 461.5
-
-# specific heat of water vapor in joules per kilogram per kelvin
-CPW = 1846.0
-
-# specific heat of dry air in joules per kilogram per kelvin
-CPD = 1005.0
-
-def kelvin_to_celsius(T_K: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    convert temperature in kelvin to celsius.
-    :param T_K: temperature in kelvin
-    :return: temperature in celsius
-    """
-    return T_K - 273.15
-
-def celcius_to_kelvin(T_C: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    convert temperature in celsius to kelvin.
-    :param T_C: temperature in celsius
-    :return: temperature in kelvin
-    """
-    return T_C + 273.15
-
-def calculate_specific_humidity(
-        Ea_Pa: Union[Raster, np.ndarray], 
-        Ps_Pa: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    Calculate the specific humidity of air as a ratio of kilograms of water to kilograms of air.
-    
-    Args:
-        Ea_Pa (Union[Raster, np.ndarray]): Actual water vapor pressure in Pascal.
-        surface_pressure_Pa (Union[Raster, np.ndarray]): Surface pressure in Pascal.
-    
-    Returns:
-        Union[Raster, np.ndarray]: Specific humidity in kilograms of water per kilograms of air.
-    """
-    return ((0.622 * Ea_Pa) / (Ps_Pa - (0.387 * Ea_Pa)))
-
-def calculate_specific_heat(specific_humidity: Union[Raster, np.ndarray]):
-    # calculate specific heat capacity of the air (Cp)
-    # in joules per kilogram per kelvin
-    # from specific heat of water vapor (CPW)
-    # and specific heat of dry air (CPD)
-    Cp_Jkg = specific_humidity * CPW + (1 - specific_humidity) * CPD
-
-    return Cp_Jkg
-
-def calculate_air_density(
-        surface_pressure_Pa: Union[Raster, np.ndarray], 
-        Ta_K: Union[Raster, np.ndarray], 
-        specific_humidity: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    Calculate air density.
-
-    Parameters:
-    surface_pressure_Pa (Union[Raster, np.ndarray]): Surface pressure in Pascal.
-    Ta_K (Union[Raster, np.ndarray]): Air temperature in Kelvin.
-    specific_humidity (Union[Raster, np.ndarray]): Specific humidity.
-
-    Returns:
-    Union[Raster, np.ndarray]: Air density in kilograms per cubic meter.
-    """
-    # numerator: Pa(N / m ^ 2 = kg * m / s ^ 2); denominator: J / kg / K * K)
-    rhoD = surface_pressure_Pa / (RD * Ta_K)
-
-    # calculate air density (rho) in kilograms per cubic meter
-    rho = rhoD * ((1.0 + specific_humidity) / (1.0 + specific_humidity * (RW / RD)))
-
-    return rho
-
-def SVP_kPa_from_Ta_C(Ta_C: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    Calculate the saturation vapor pressure in kiloPascal (kPa) from air temperature in Celsius.
-
-    Parameters:
-    Ta_C (Union[Raster, np.ndarray]): Air temperature in Celsius.
-
-    Returns:
-    Union[Raster, np.ndarray]: Saturation vapor pressure in kPa.
-
-    """
-    SVP_kPa = np.clip(0.611 * np.exp((Ta_C * 17.27) / (Ta_C + 237.7)), 1, None)
-
-    return SVP_kPa
-
-def SVP_Pa_from_Ta_C(Ta_C: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    Calculate the saturation vapor pressure in Pascal (Pa) from the air temperature in Celsius (Ta_C).
-
-    Parameters:
-        Ta_C (Union[Raster, np.ndarray]): Air temperature in Celsius.
-
-    Returns:
-        Union[Raster, np.ndarray]: Saturation vapor pressure in Pascal (Pa).
-    """
-    return SVP_kPa_from_Ta_C(Ta_C) * 1000
-
-def calculate_surface_pressure(elevation_m: Union[Raster, np.ndarray], Ta_C: Union[Raster, np.ndarray]) -> Union[Raster, np.ndarray]:
-    """
-    Calculate surface pressure using elevation and air temperature.
-
-    Parameters:
-        elevation_m (Union[Raster, np.ndarray]): Elevation in meters.
-        Ta_K (Union[Raster, np.ndarray]): Air temperature in Kelvin.
-
-    Returns:
-        Union[Raster, np.ndarray]: Surface pressure in Pascal (Pa).
-    """
-    Ta_K = kelvin_to_celsius(Ta_C)
-    Ps_Pa = 101325.0 * (1.0 - 0.0065 * elevation_m / Ta_K) ** (9.807 / (0.0065 * 287.0))  # [Pa]
-
-    return Ps_Pa
-