wrfrun 0.1.7__py3-none-any.whl → 0.1.8__py3-none-any.whl

wrfrun/extension/goos_sst/core.py

@@ -0,0 +1,108 @@
+"""
+wrfrun.extension.goos_sst.core
+##############################
+
+Implementation of ``extension.goos_sst``'s core functionality.
+
+
+.. autosummary::
+    :toctree: generated/
+    
+    merge_era5_goos_sst_grib
+"""
+
+from os.path import dirname
+
+import cfgrib as cf
+import numpy as np
+from pandas import to_datetime
+from seafog import goos_sst_find_data, goos_sst_parser
+from xarray import DataArray
+
+from .utils import create_sst_grib
+from wrfrun.utils import logger
+
+
+def merge_era5_goos_sst_grib(surface_grib_path: str, save_path: str, sst_data_save_path: str | None = None, resolution="low"):
+    """
+    This function reads ERA5 skin temperature (SKT) data from the GRIB file,
+    interpolates it to the same grib of NEAR-GOOS sea surface temperature (SST) data,
+    merges them and creates a new GRIB file with the new data.
+
+    :param surface_grib_path: GRIB file which contains SKT data.
+    :type surface_grib_path: str
+    :param save_path: Save path of the new GRIB file.
+    :type save_path: str
+    :param sst_data_save_path: Save path of downloaded NEAR-GOOS SST data.
+                               If None, save data to the parent directory of ``save_path``.
+    :type sst_data_save_path: str | None
+    :param resolution: Resolution of downloaded NEAR-GOO SST data.
+                       Please check ``seafog.goos_sst_find_data`` for more information.
+    :type resolution: str
+    """
+    dataset_list = cf.open_datasets(surface_grib_path)
+
+    dataset = None
+    for _ds in dataset_list:
+        if "skt" in _ds:
+            dataset = _ds
+            break
+
+    if dataset is None:
+        logger.error(f"'skt' data not found in {surface_grib_path}")
+        raise ValueError
+
+    skt: DataArray = dataset["skt"]     # type: ignore
+
+    longitude_start, longitude_end = skt["longitude"][0].data, skt["longitude"][-1].data    # type: ignore
+    latitude_start, latitude_end = skt["latitude"][-1].data, skt["latitude"][0].data    # type: ignore
+
+    data_time = skt["time"].to_numpy()  # type: ignore
+    data_time = to_datetime(data_time).strftime("%Y-%m-%d %H:%M")
+
+    if sst_data_save_path is None:
+        sst_data_save_path = dirname(save_path)
+
+    # read the first time data, so we can interpolate skt data to the same resolution as sst data.
+    _data = goos_sst_find_data(data_time[0], sst_data_save_path, resolution=resolution, show_progress=False)
+    _data = goos_sst_parser(_data, resolution=resolution)
+    _data = _data.loc[latitude_start:latitude_end, longitude_start:longitude_end]
+
+    skt = skt.interp({"longitude": _data["longitude"], "latitude": _data["latitude"]})
+    sst = np.zeros_like(skt)
+    temp_data = np.zeros_like(skt[0])
+
+    temp_data[:] = skt.to_numpy()[0]
+    index = ~np.isnan(_data.to_numpy())
+    temp_data[index] = _data.to_numpy()[index]
+    sst[0] = temp_data
+
+    # loop other time.
+    for time_index, _data_time in enumerate(data_time[1:], start=1):
+        _data = goos_sst_find_data(_data_time, sst_data_save_path, resolution=resolution, show_progress=False)
+        _data = goos_sst_parser(_data, resolution=resolution)
+        _data = _data.loc[latitude_start:latitude_end, longitude_start:longitude_end]
+
+        temp_data[:] = skt.to_numpy()[time_index]
+        index = ~np.isnan(_data.to_numpy())
+        temp_data[index] = _data.to_numpy()[index]
+        sst[time_index] = temp_data
+
+    sst = DataArray(
+        name="sst",
+        data=sst,
+        dims=["time", "latitude", "longitude"],
+        coords={
+            "time": skt["time"],
+            "latitude": skt["latitude"],
+            "longitude": skt["longitude"]
+        }, attrs={
+            "units": "K",
+            "long_name": "Sea surface temperature"
+        }
+    )
+
+    create_sst_grib(sst, save_path=save_path)
+
+
+__all__ = ["merge_era5_goos_sst_grib"]

wrfrun/extension/goos_sst/res/Vtable.ERA_GOOS_SST

@@ -0,0 +1,7 @@
+GRIB | Level| Level| Level| metgrid  |  metgrid | metgrid                                  |
+Code | Code |   1  |   2  | Name     |  Units   | Description                              |
+-----+------+------+------+----------+----------+------------------------------------------+
+  34 |  1   |   0  |      | SST      | K        | Sea-Surface Temperature                  |
+-----+------+------+------+----------+----------+------------------------------------------+
+
+# Vtable setting for GRIB file created by wrfrun.extension.sst

wrfrun/extension/goos_sst/res/__init__.py

@@ -0,0 +1,9 @@
+from os.path import abspath, dirname
+
+
+_PATH = abspath(dirname(__file__))
+
+VTABLE_ERA_GOOS_SST = f"{_PATH}/Vtable.ERA_GOOS_SST"
+
+
+__all__ = ["VTABLE_ERA_GOOS_SST"]

wrfrun/extension/goos_sst/utils.py

@@ -0,0 +1,81 @@
+from cfgrib.xarray_to_grib import to_grib
+from numpy.dtypes import DateTime64DType
+from pandas import to_datetime
+from xarray import DataArray, Dataset
+
+from wrfrun.utils import logger
+
+
+def create_sst_grib(data: DataArray, save_path: str):
+    """
+    Create a dataset and save it to a GRIB file.
+
+    Args:
+        data: DataArray data, it should at least contain three dimensions: ``["time", "latitude", "longitude"]``.
+        save_path: GRIB file path.
+
+    """
+    # check the data's dimensions.
+    for _dim in ["time", "longitude", "latitude"]:
+        if _dim not in data.dims:
+            logger.error(f"Essential dimension '{_dim}' not found in data")
+            raise KeyError
+
+    if not isinstance(data["time"].dtype, DateTime64DType):     # type: ignore
+        data = data.assign_coords(time=to_datetime(data["time"].data))
+
+    data = data.sortby(data["latitude"], ascending=False)
+
+    longitude_start, longitude_stop = data["longitude"][0].data, data["longitude"][-1].data
+    latitude_start, latitude_stop = data["latitude"][0].data, data["latitude"][-1].data
+    delta_longitude = data["longitude"][1].data - data["longitude"][0].data
+    delta_latitude = data["latitude"][0].data - data["latitude"][1].data
+    longitude_length = data["longitude"].size
+    latitude_length = data["latitude"].size
+    points_number = data.size
+
+    data = data.assign_attrs(**{
+        "units": "K",
+        "long_name": "Sea surface temperature",
+        "standard_name": "Sea surface temperature",
+        # The following keys and values will be used in GRIB.
+        "GRIB_paramId": 34,
+        "GRIB_shortName": "sst",
+        "GRIB_units": "K",
+        "GRIB_name": "Sea surface temperature",
+        "GRIB_stepUnits": 1,
+        "GRIB_stepType": "instant",
+        "GRIB_gridType": "regular_ll",
+        "GRIB_iDirectionIncrementInDegrees": delta_longitude,
+        "GRIB_iScanNegatively": 0,
+        "GRIB_jDirectionIncrementInDegrees": delta_latitude,
+        "GRIB_jScanPositively": 0,
+        "GRIB_latitudeOfFirstGridPointInDegrees": latitude_start,
+        "GRIB_latitudeOfLastGridPointInDegrees": latitude_stop,
+        "GRIB_longitudeOfFirstGridPointInDegrees": longitude_start,
+        "GRIB_longitudeOfLastGridPointInDegrees": longitude_stop,
+        "GRIB_Ny": latitude_length,
+        "GRIB_Nx": longitude_length,
+        "GRIB_typeOfLevel": "surface",
+        # The following keys and values can't be found at ECMWF websites.
+        "GRIB_cfName": "unknown",
+        "GRIB_cfVarName": "sst",
+        "GRIB_dataType": "an",  # Analysis data, defined at https://codes.ecmwf.int/grib/format/mars/type/
+        "GRIB_gridDefinitionDescription": "Latitude/Longitude Grid",
+        # "GRIB_missingValue": -9999,
+        "GRIB_numberOfPoints": points_number,
+        "GRIB_totalNumber": 0,
+        "GRIB_uvRelativeToGird": 0
+    })
+
+    to_grib(
+        Dataset({
+            "sst": data
+        }, attrs={
+            "GRIB_centre": "ecmf",
+            "GRIB_edition": 1,
+        }), save_path
+    )
+
+
+__all__ = ["create_sst_grib"]

wrfrun/model/plot.py

@@ -5,8 +5,8 @@ from shutil import copyfile, move
 from wrfrun.core import WRFRUNConfig
 from wrfrun.res import EXT_NCL_PLOT_SCRIPT
 from wrfrun.utils import call_subprocess, check_path, logger
-from .wrf.namelist import prepare_wps_namelist
 from .base import NamelistName
+from .wrf.namelist import prepare_wps_namelist
 
 
 def plot_domain_area():

wrfrun/model/wrf/core.py

@@ -3,7 +3,7 @@ from os.path import abspath, basename, exists
 from shutil import copyfile, move, rmtree
 from typing import Optional
 
-from wrfrun.core import ExecutableBase, FileConfigDict, InputFileError, NamelistError, WRFRUNConfig, WRFRUNExecDB
+from wrfrun.core import ExecutableBase, FileConfigDict, InputFileError, NamelistIDError, WRFRUNConfig, WRFRUNExecDB
 from wrfrun.utils import logger
 from ._metgrid import reconcile_namelist_metgrid
 from ._ndown import process_after_ndown
@@ -626,7 +626,7 @@ class DFI(ExecutableBase):
 
         if not WRFRUNConfig.register_custom_namelist_id("dfi"):
             logger.error("Can't register namelist for DFI.")
-            raise NamelistError("Can't register namelist for DFI.")
+            raise NamelistIDError("Can't register namelist for DFI.")
         WRFRUNConfig.update_namelist(self.custom_config["namelist"], "dfi")
 
     def before_exec(self):
@@ -659,7 +659,7 @@ class DFI(ExecutableBase):
 
             if not WRFRUNConfig.register_custom_namelist_id("dfi"):
                 logger.error("Can't register namelist for DFI.")
-                raise NamelistError("Can't register namelist for DFI.")
+                raise NamelistIDError("Can't register namelist for DFI.")
 
             prepare_dfi_namelist()
 

wrfrun/model/wrf/namelist.py

@@ -79,18 +79,14 @@ def prepare_wps_namelist():
             "dy": wrf_config["domain"]["dy"],
             "ref_lat": wrf_config["domain"]["ref_lat"],
             "ref_lon": wrf_config["domain"]["ref_lon"],
+            "map_proj": wrf_config["domain"]["map_proj"],
+            "truelat1": wrf_config["domain"]["truelat1"],
+            "truelat2": wrf_config["domain"]["truelat2"],
             "stand_lon": wrf_config["domain"]["stand_lon"],
             "geog_data_path": wrf_config["geog_data_path"]
         }
     }
 
-    # # use loop to process config of map_proj
-    for key in wrf_config["domain"]["map_proj"]:
-        if key == "name":
-            update_value["geogrid"]["map_proj"] = wrf_config["domain"]["map_proj"][key]
-        else:
-            update_value["geogrid"][key] = wrf_config["domain"]["map_proj"][key]
-
     # # update namelist
     WRFRUNConfig.update_namelist(update_value, "wps")
 

wrfrun/model/wrf/scheme.py

@@ -308,4 +308,88 @@ class SchemeSurfaceLayer:
         return integer_label_map[key]
 
 
-__all__ = ["SchemeCumulus", "SchemeLandSurfaceModel", "SchemeLongWave", "SchemePBL", "SchemeShortWave", "SchemeSurfaceLayer"]
+@dataclass()
+class SchemeMicrophysics:
+    """
+    Microphysics schemes.
+    """
+    OFF = 0
+    KESSLER = 1
+    PURDUE_LIN = 2
+    WSM_3_CLASS_ICE = 3
+    WSM_5_CLASS = 4
+    FERRIER_ETA = 5
+    WSM_6_CLASS_GRAUPEL = 6
+    GODDARD_4_ICE = 7
+    THOMPSON_GRAUPEL = 8
+    MILBRANDT_YAU_2_MOMENT = 9
+    MORRISON_2_MOMENT = 10
+    CAM_V51_5_CLASS = 11
+    SBU_YLIN_5_CLASS = 13
+    WDM_5_CLASS = 14
+    HIGH_FERRIER_ADVECTION = 15
+    WDM_6_CLASS = 16
+    NSSL_2_MOMENT_4_ICE = 18
+    WSM7 = 24
+    WDM7 = 26
+    AEROSOL_AWARE_THOMPSON = 28
+    HUJI = 30
+    THOMPSON_HAIL_GRAUPEL_AEROSOL = 38
+    MORRISON_WITH_CESM_AEROSOL = 40
+    P3_1_ICE_1_MOMENT_QCLOUD = 50
+    P3_1_ICE_2_MOMENT_QCLOUD = 51
+    P3_2_ICE_2_MOMENT_QCLOUD = 52
+    P3_1_ICE_3_MOMENT_ICE_2_MOMENT_QCLOUD = 53
+    JENSEN_ISHMAEL = 55
+    NTU = 56
+
+    @classmethod
+    def get_scheme_id(cls, key: str = "lin"):
+        """Get corresponding integer label for microphysics scheme
+
+        Args:
+            key (str, optional): Name of microphysics scheme. Defaults to "PURDUE_LIN".
+        """
+        # Here is the map of scheme name and integer label in WRF
+        # Reference link: https://www2.mmm.ucar.edu/wrf/users/wrf_users_guide/build/html/namelist_variables.html
+        integer_label_map = {
+            "off": cls.OFF,
+            "kessler": cls.KESSLER,
+            "lin": cls.PURDUE_LIN,
+            "wsm3": cls.WSM_3_CLASS_ICE,
+            "wsm5": cls.WSM_5_CLASS,
+            "ferrier": cls.FERRIER_ETA,
+            "wsm6": cls.WDM_6_CLASS,
+            "goddard4": cls.GODDARD_4_ICE,
+            "thompson": cls.THOMPSON_GRAUPEL,
+            "milbrandt": cls.MILBRANDT_YAU_2_MOMENT,
+            "morrison2": cls.MORRISON_2_MOMENT,
+            "cam": cls.CAM_V51_5_CLASS,
+            "sbu": cls.SBU_YLIN_5_CLASS,
+            "wdm5": cls.WDM_5_CLASS,
+            "high_ferrier": cls.HIGH_FERRIER_ADVECTION,
+            "wdm6": cls.WDM_6_CLASS,
+            "nssl": cls.NSSL_2_MOMENT_4_ICE,
+            "wsm7": cls.WSM7,
+            "wdm7": cls.WDM7,
+            "thompson_aerosol": cls.AEROSOL_AWARE_THOMPSON,
+            "cesm_morrison": cls.MORRISON_WITH_CESM_AEROSOL,
+            "p311": cls.P3_1_ICE_1_MOMENT_QCLOUD,
+            "p312": cls.P3_1_ICE_2_MOMENT_QCLOUD,
+            "p322": cls.P3_2_ICE_2_MOMENT_QCLOUD,
+            "p3132": cls.P3_1_ICE_3_MOMENT_ICE_2_MOMENT_QCLOUD,
+            "jensen": cls.JENSEN_ISHMAEL,
+            "ntu": cls.NTU,
+        }
+
+        # check if key is in map
+        if key not in integer_label_map:
+            logger.error(
+                f"Key error: {key}. Valid key: {list(integer_label_map.keys())}"
+            )
+            raise KeyError
+
+        return integer_label_map[key]
+
+
+__all__ = ["SchemeCumulus", "SchemeLandSurfaceModel", "SchemeLongWave", "SchemePBL", "SchemeShortWave", "SchemeSurfaceLayer", "SchemeMicrophysics"]

wrfrun/plot/wps.py

@@ -12,29 +12,29 @@ from xarray import DataArray
 from wrfrun.utils import logger
 
 # Here defines the colormap for landuse and soiltype
-LANDUSE_CMAP = ListedColormap([
-    [0, .4, 0],         # 1 Evergreen Needleleaf Forest
-    [0, .4, .2],        # 2 Evergreen Broadleaf Forest
-    [.2, .8, .2],       # 3 Deciduous Needleleaf Forest
-    [.2, .8, .4],       # 4 Deciduous Broadleaf Forest
-    [.2, .6, .2],       # 5 Mixed Forests
-    [.3, .7, 0],        # 6 Closed Shrublands
-    [.82, .41, .12],    # 7 Open Shurblands
-    [.74, .71, .41],    # 8 Woody Savannas
-    [1, .84, .0],       # 9 Savannas
-    [0, 1, 0],          # 10 Grasslands
-    [0, 1, 1],          # 11 Permanant Wetlands
-    [1, 1, 0],          # 12 Croplands
-    [1, 0, 0],          # 13 Urban and Built-up
-    [.7, .9, .3],       # 14 Cropland/Natual Vegation Mosaic
-    [1, 1, 1],          # 15 Snow and Ice
-    [.914, .914, .7],   # 16 Barren or Sparsely Vegetated
-    [.5, .7, 1],        # 17 Water (like oceans)
-    [.86, .08, .23],    # 18 Wooded Tundra
-    [.97, .5, .31],     # 19 Mixed Tundra
-    [.91, .59, .48],    # 20 Barren Tundra
-    [0, 0, .88],        # 21 Lake
-])
+LANDUSE_CMAP = ListedColormap((
+    [0, .4, 0],         # 1: Evergreen Needleleaf Forest
+    [0, .4, .2],        # 2: Evergreen Broadleaf Forest
+    [.2, .8, .2],       # 3: Deciduous Needleleaf Forest
+    [.2, .8, .4],       # 4: Deciduous Broadleaf Forest
+    [.2, .6, .2],       # 5: Mixed Forests
+    [.3, .7, 0],        # 6: Closed Shrublands
+    [.82, .41, .12],    # 7: Open Shurblands
+    [.74, .71, .41],    # 8: Woody Savannas
+    [1, .84, .0],       # 9: Savannas
+    [0, 1, 0],          # 10: Grasslands
+    [0, 1, 1],          # 11: Permanent Wetlands
+    [1, 1, 0],          # 12: Croplands
+    [1, 0, 0],          # 13: Urban and Built-up
+    [.7, .9, .3],       # 14: Cropland/Natual Vegation Mosaic
+    [1, 1, 1],          # 15: Snow and Ice
+    [.914, .914, .7],   # 16: Barren or Sparsely Vegetated
+    [.5, .7, 1],        # 17: Water (like oceans)
+    [.86, .08, .23],    # 18: Wooded Tundra
+    [.97, .5, .31],     # 19: Mixed Tundra
+    [.91, .59, .48],    # 20: Barren Tundra
+    [0, 0, .88],        # 21: Lake
+))
 
 LANDUSE_LABELS = [
     'Evergreen Needleleaf Forest',
@@ -61,25 +61,25 @@ LANDUSE_LABELS = [
 ]
 
 
-# Here defines the colormap and labels for soil type
-SOILTYPE_CMAP = ListedColormap([
-    [0, .4, 0],         # 1 Sand
-    [0, .4, .2],        # 2 Loamy Sand
-    [.2, .8, .2],       # 3 Sandy Loam
-    [.2, .8, .4],       # 4 Silt Loam
-    [.2, .6, .2],       # 5 Silt
-    [.3, .7, 0],        # 6 Loam
-    [.82, .41, .12],    # 7 Sandy Clay Loam
-    [.74, .71, .41],    # 8 Silty Clay Loam
-    [1, .84, .0],       # 9 Clay Loam
-    [0, 1, 0],          # 10 Sandy Clay
-    [0, 1, 1],          # 11 Silty Clay
-    [1, 1, 0],          # 12 Clay
-    [1, 0, 0],          # 13 Organic Material
-    [.7, .9, .3],       # 14 Water
-    [1, 1, 1],          # 15 Bedrock
-    [.914, .914, .7]    # 16 Other
-])
+# Here defines the colormap and labels for soil types
+SOILTYPE_CMAP = ListedColormap((
+    [0, .4, 0],         # 1: Sand
+    [0, .4, .2],        # 2: Loamy Sand
+    [.2, .8, .2],       # 3: Sandy Loam
+    [.2, .8, .4],       # 4: Silt Loam
+    [.2, .6, .2],       # 5: Silt
+    [.3, .7, 0],        # 6: Loam
+    [.82, .41, .12],    # 7: Sandy Clay Loam
+    [.74, .71, .41],    # 8: Silty Clay Loam
+    [1, .84, .0],       # 9: Clay Loam
+    [0, 1, 0],          # 10: Sandy Clay
+    [0, 1, 1],          # 11: Silty Clay
+    [1, 1, 0],          # 12: Clay
+    [1, 0, 0],          # 13: Organic Material
+    [.7, .9, .3],       # 14: Water
+    [1, 1, 1],          # 15: Bedrock
+    [.914, .914, .7]    # 16: Other
+))
 
 SOILTYPE_LABELS = [
     "Sand",
@@ -102,13 +102,13 @@ SOILTYPE_LABELS = [
 
 
 def get_cmap_ticks(name: str) -> tuple:
-    """Get corresponding colormap, labels and ticks.
-
-    Args:
-        name (str): Field name.
+    """
+    Get corresponding colormap, labels and ticks.
 
-    Returns:
-        tuple: (colormap, labels, ticks)
+    :param name: File name.
+    :type name: str
+    :return: (colormap, labels, ticks)
+    :rtype: tuple
     """
     # name map
     cmap_ticks_map = {
@@ -141,15 +141,23 @@ def get_cmap_ticks(name: str) -> tuple:
 
 
 def draw_geogrid(data_path: str, field: str, fig: Figure, nrow: int, ncol: int, index: int) -> Tuple[GeoAxes, QuadMesh]:
-    """Draw specific data in geo_em*
-
-    Args:
-        data_path (str): data path.
-        field (str): Which field you want to draw.
-        fig (Figure): Matplotlib figure object.
-        nrow (int): Row number of axes.
-        ncol (int): Column number of axes.
-        index (int): Index of axes.
+    """
+    Draw specific data in geogrid outputs.
+
+    :param data_path: Data path.
+    :type data_path: str
+    :param field: Which field you want to draw.
+    :type field: str
+    :param fig: Matplotlib figure object.
+    :type fig: Figure
+    :param nrow: Row number of axes.
+    :type nrow: int
+    :param ncol: Column number of axes.
+    :type ncol: int
+    :param index: Index of axes.
+    :type index: int
+    :return: (GeoAxes, QuadMesh)
+    :rtype: tuple
     """
     # take out data
     data: DataArray = getvar(Dataset(data_path), field)     # type: ignore
@@ -159,7 +167,7 @@ def draw_geogrid(data_path: str, field: str, fig: Figure, nrow: int, ncol: int,
         data = data.argmax(dim="soil_cat", keep_attrs=True)     # type: ignore
         # get data attrs
         data_attrs = data.attrs
-        data = data + 1
+        data = data + 1     # type: ignore
         data.attrs = data_attrs
 
     # get latitude and longitude

wrfrun/res/config.toml.template

@@ -109,7 +109,9 @@ dx = 9000
 dy = 9000
 
 # Projection.
-map_proj = { name = 'lambert', truelat1 = 34.0, truelat2 = 40.0}
+map_proj = 'lambert'
+truelat1 = 34.0
+truelat2 = 40.0
 
 # Central point of the first area.
 ref_lat = 37.0

{wrfrun-0.1.7.dist-info → wrfrun-0.1.8.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: wrfrun
-Version: 0.1.7
+Version: 0.1.8
 Summary: wrfrun is a comprehensive toolkit for managing and using WRF
 Keywords: WRF
 Author-Email: Syize <syizeliu@gmail.com>
@@ -10,7 +10,7 @@ Project-URL: homepage, https://github.com/Syize/wrfrun
 Project-URL: repository, https://github.com/Syize/wrfrun
 Project-URL: documentation, https://wrfrun.syize.cn
 Project-URL: Bug Tracker, https://github.com/Syize/wrfrun/issues
-Requires-Python: <3.13
+Requires-Python: >=3.10
 Requires-Dist: numpy<2.0.0
 Requires-Dist: xarray
 Requires-Dist: netCDF4

{wrfrun-0.1.7.dist-info → wrfrun-0.1.8.dist-info}/RECORD

@@ -1,38 +1,43 @@
-wrfrun-0.1.7.dist-info/METADATA,sha256=nrPBmS0UxcjIiCN9wLyL798-r3OKpirBJUPlWrzobw8,2912
-wrfrun-0.1.7.dist-info/WHEEL,sha256=5J4neoE7k6LMgx4Fz1FHgBiO3YevhJGtNQ3muDrdLQM,75
+wrfrun-0.1.8.dist-info/METADATA,sha256=_Y0huf7Q8hAsLFsgIlPZYF8FDmZ_DKqsGMxcBaewcpQ,2913
+wrfrun-0.1.8.dist-info/WHEEL,sha256=5J4neoE7k6LMgx4Fz1FHgBiO3YevhJGtNQ3muDrdLQM,75
 wrfrun/__init__.py,sha256=l2zGk2qDWa3EXIlwvzqBYZHOU9273TN75Rmvx_7PAYM,58
-wrfrun/data.py,sha256=AjkDSZqkLROnWLw9qDOBUctoQmJ3R1bkfjnPzLMjPfs,14751
+wrfrun/data.py,sha256=W6IhopHpp0lhvVo-beuO6NsbfsGDmN5MNf8Bqkt7xsA,14845
 wrfrun/pbs.py,sha256=gzRHbYqpU-c8lynu39TUKNvmPntvrPwfQVxK337yvjw,2454
 wrfrun/run.py,sha256=0mvuI98gLjq1gP1eaLUZiU1LJMCjtL668aKZTwu3n2g,9114
 wrfrun/utils.py,sha256=v3ZRu3Rd8-pZjo7I0PrJtgBN1zIBAt_nVPs0BP_DFwY,8118
 wrfrun/workspace.py,sha256=5tWIIFxwAE4hkdg_wJALGDc3P0bfwb_IRTIyPIjN-tE,3105
-wrfrun/core/__init__.py,sha256=bf9Fp6VS5i9xB10lazYBOkJapF2jnQpyiB6qucVjp0I,107
-wrfrun/core/base.py,sha256=S9Mqr3LOpos0CP0TsFwtCQl0tLIObIhjX7VxjS0LRDY,23172
-wrfrun/core/config.py,sha256=v7zdhAthlMmjv_BWINxkmgNHLvKZdzQhgM5B-Nlsl7M,23795
-wrfrun/core/error.py,sha256=fCWVC4-GTyDH2obzu9-D_bzF8_uzqZoUi7Ffy0B5qoc,1324
-wrfrun/core/replay.py,sha256=oyV9QppZG99GbAksQqoRxeBQ3jjXU-7YqkpxHRyfBvc,3229
-wrfrun/core/server.py,sha256=KK_yj4J-5F_s-rKO5FK7e8WojsnuZMbhWtMceWYbaBY,6524
-wrfrun/extension/__init__.py,sha256=alIDGBnxWH4JvP-UW-7N99seBBi0r1GV1h8f1ERFBec,21
+wrfrun/core/__init__.py,sha256=zsitfcmbO_Qf_Uh0FzZ4rceZAo3WEYev32jSjgRjtkw,1325
+wrfrun/core/base.py,sha256=zix413Nc0YrgIR6oakVhTN405KUmZ7OKlJTgL0fNUVc,30307
+wrfrun/core/config.py,sha256=g021WEkRZpvsLyzVHlahL9WtAP7zKG9polW3B84eT3k,29252
+wrfrun/core/error.py,sha256=etnImryyqDUggfZOKK1ynOm71GRSFJVGGOwhuiZahTI,2383
+wrfrun/core/replay.py,sha256=nqVmR3Hb4zO6PnROssv3PKw7gvvdi-UqIsjX-NV1WLI,4782
+wrfrun/core/server.py,sha256=RjzFVlRZg60rLMy4Gy5UhZWPUOQP307rT4ttpvbcecQ,8640
+wrfrun/extension/__init__.py,sha256=YDYCS0KD3e-uLvf1K08ktNfP2sUSKpbgkHrvXj21EPQ,913
 wrfrun/extension/utils.py,sha256=OBIAlAXakbh1QhrF2BFlYAFV062ZfJjX0mYJkmcU_RI,2751
+wrfrun/extension/goos_sst/__init__.py,sha256=M8LWGLXqHyPUWijl5925jFmwJ1ih3kHfT2jTxZzwMBA,2126
+wrfrun/extension/goos_sst/core.py,sha256=D0po6PjBqcEareVQO1VjdIUrFz9lJrnAZ5q2sax78Tw,3794
+wrfrun/extension/goos_sst/utils.py,sha256=h0OBCwxS67kG6YNclZLvyPYnT7uNZd3Uq9f6Vg2EaMc,3052
+wrfrun/extension/goos_sst/res/__init__.py,sha256=M-qN-WI0EIvhYo2hwnnIK_fXL3TbR12KzSC9fvjgDaY,164
+wrfrun/extension/goos_sst/res/Vtable.ERA_GOOS_SST,sha256=aX8RSTWL8GVdxf2TgvIoQx5xDVuLi0F5_Eg2e8Gaa-E,529
 wrfrun/extension/littler/__init__.py,sha256=pG-YPVG0gJIJ6s4xcAz9S_CnUxpUcz33wl9eNUSgxGk,20
 wrfrun/extension/littler/utils.py,sha256=TsN3GiAFszPWF7RNM_hYdARRY_OLL8mYm87UyZupUtU,21615
 wrfrun/model/__init__.py,sha256=lQNNnh3evXZytLAjdnhlIEFWerP1L0XAgDAUQhqIRyE,119
 wrfrun/model/base.py,sha256=WnVNu12ICJgThfZR7kK6_XLCQWVYFsHrdyAz8wl_YRc,218
-wrfrun/model/plot.py,sha256=ATG5nl5kV8YfoeNFIPLCUnnIQIrfIFNGPb9d8raM6S0,1582
+wrfrun/model/plot.py,sha256=Hh-GxQPDzpwAEAq-rFh0VcIWzgMDxZVT0c40SIRa6_o,1582
 wrfrun/model/utils.py,sha256=Pyj9D9NMKw3wQcZcM7pJaLG09Vov447AxFKoVeZZU2s,1155
 wrfrun/model/wrf/__init__.py,sha256=T3S1BD9QDoW4uq871QW4_8-2BvCRWYEecz71EcdbMaI,136
-wrfrun/model/wrf/core.py,sha256=FEzRXUKK5dENxev3wR_GON5n1Zxq4Hj-Bs_BNP7P_ok,32134
+wrfrun/model/wrf/core.py,sha256=srOD7KFF5TgPAcWTw0PLHFjkGjDzb8PYbxuPx5B4bQE,32140
 wrfrun/model/wrf/exec_wrap.py,sha256=gIsGZDAjzRzxt24eo8t-laIK4UBiyKf_ysfAZLlE1og,3918
 wrfrun/model/wrf/geodata.py,sha256=X9OUOm0__5NI7sl2z1F9_ur5wZ4P0HkpBNcRZQYKof0,9740
 wrfrun/model/wrf/_metgrid.py,sha256=i2qILM_1LkdnBZqXXL3swVhMRvhFvq4bOcRPGI6kSzg,2248
-wrfrun/model/wrf/namelist.py,sha256=t_JuWYj4oNNVwzTMKo6BNlv2LIq4t3JC7ZgVkaKfNbA,14923
+wrfrun/model/wrf/namelist.py,sha256=T9b5H6pCljhKUhoR9qwTbdMF04Hrdp3mHwCP-VuhFWE,14791
 wrfrun/model/wrf/_ndown.py,sha256=ZcWttGXiT5FzGCAohmbeNU6bi8GUOUleiTK0ivn2BaI,1548
-wrfrun/model/wrf/scheme.py,sha256=F7HxoaNkkPpjkse5t9sGS3CJDvn_yBdHoaCG7WLdbqY,8524
+wrfrun/model/wrf/scheme.py,sha256=8y85Dbu-GajwjHr3Spw8_tN21xAloD-SnmxJd11nXKE,11254
 wrfrun/model/wrf/vtable.py,sha256=0V4fo0Limnx4oJ2NByq_eHPzj8L4mnHNlrKL34Buw1A,2432
 wrfrun/plot/__init__.py,sha256=9Kn0IgkX10sHEqHJwk7mZV-dP14XMNcvXN8znO89FIw,19
-wrfrun/plot/wps.py,sha256=AQ1B6llNUrPU3PR8XKNLSzs_H-v60a0Wep-8Mcto6NM,5319
+wrfrun/plot/wps.py,sha256=pvkxbh5760AAM4KaPteMzhFniZZFtkYF31xhzSBa7lo,5552
 wrfrun/res/__init__.py,sha256=YBcVaM2EZZJkDDvyKdk-WZmcFQoTFENZKCwcIQzEf3I,1133
-wrfrun/res/config.toml.template,sha256=aVXCSzhbV3P-MQCZOk2g-y0auRmqgr3jW920W3wks9g,4386
+wrfrun/res/config.toml.template,sha256=LTCEwXSFF2QXnwGRoB0Jw5Qs_vXwdb6Zd7fQtNAz9Uw,4374
 wrfrun/res/run.sh.template,sha256=ZN7k0a9ONOtZUQ67JiTva-0DNVJG5O7BdOXQ8-LCscs,247
 wrfrun/res/extension/plotgrids.ncl,sha256=B0mvH1H1j_w7EEana9HmU1XJZtG0w9IccQ54zXpbQG0,7546
 wrfrun/res/namelist/namelist.input.da_wrfvar.template,sha256=Cwc-XPu_spJeQte4duWrulPBOLRMEBtn0mIn0pgMmKY,4912
@@ -43,4 +48,4 @@ wrfrun/res/namelist/namelist.wps.template,sha256=HlA7-SHs4C-cKRb3h6D0Kl1Y-5VSJ1L
 wrfrun/res/namelist/parame.in.template,sha256=vR8JSix20FAKGqj6jK8QuEAeWkGvg8_iptdVlzjPX6o,259
 wrfrun/res/job_scheduler/pbs.template,sha256=m7eyvZkzQRpjs43CvSeZ1biHYNyEJJkgXt3IU5ZhCmg,148
 wrfrun/res/job_scheduler/slurm.template,sha256=BrKL8DcHfpFfv0y0z7Qe_cqjwfzMyF-3V6vM72ehr9Y,312
-wrfrun-0.1.7.dist-info/RECORD,,
+wrfrun-0.1.8.dist-info/RECORD,,