ladim 2.0.1__py3-none-any.whl → 2.0.3__py3-none-any.whl

ladim/__init__.py

@@ -1,3 +1,3 @@
-__version__ = '2.0.1'
+__version__ = '2.0.3'
 
 from .main import main, run

ladim/config.py

@@ -28,98 +28,118 @@ def configure(module_conf):
 
 def _versioned_configure(config_dict):
     if config_dict['version'] == 1:
-        config_dict = _convert_1_to_2(config_dict)
+        config_dict = convert_1_to_2(config_dict)
 
     return config_dict
 
 
-def _convert_1_to_2(c):
+def dict_get(d, items, default=None):
+    if isinstance(items, str):
+        items = [items]
+
+    for item in items:
+        try:
+            return dict_get_single(d, item)
+        except KeyError:
+            pass
+
+    return default
+
+
+def dict_get_single(d, item):
+    tokens = str(item).split(sep='.')
+    sub_dict = d
+    for t in tokens:
+        if t in sub_dict:
+            sub_dict = sub_dict[t]
+        else:
+            raise KeyError
+
+    return sub_dict
+
+
+def convert_1_to_2(c):
+
+    out = {}
+
     # Read timedelta
-    dt_value, dt_unit = c['numerics']['dt']
-    dt_sec = np.timedelta64(dt_value, dt_unit).astype('timedelta64[s]').astype('int64')
+    dt_sec = None
+    if 'numerics' in c:
+        if 'dt' in c['numerics']:
+            dt_value, dt_unit = c['numerics']['dt']
+            dt_sec = np.timedelta64(dt_value, dt_unit).astype('timedelta64[s]').astype('int64')
+
+    out['version'] = 2
+
+    out['solver'] = {}
+    out['solver']['start'] = dict_get(c, 'time_control.start_time')
+    out['solver']['stop'] = dict_get(c, 'time_control.stop_time')
+    out['solver']['step'] = dt_sec
+    out['solver']['seed'] = dict_get(c, 'numerics.seed')
+    out['solver']['order'] = ['release', 'forcing', 'output', 'tracker', 'ibm', 'state']
+
+    out['grid'] = {}
+    out['grid']['file'] = dict_get(c, [
+        'files.grid_file', 'gridforce.grid_file',
+        'files.input_file', 'gridforce.input_file'])
+    out['grid']['legacy_module'] = dict_get(c, 'gridforce.module', '') + '.Grid'
+    out['grid']['start_time'] = np.datetime64(dict_get(c, 'time_control.start_time', '1970'), 's')
+
+    out['forcing'] = {}
+    out['forcing']['file'] = dict_get(c, ['gridforce.input_file', 'files.input_file'])
+    out['forcing']['legacy_module'] = dict_get(c, 'gridforce.module', '') + '.Forcing'
+    out['forcing']['start_time'] = np.datetime64(dict_get(c, 'time_control.start_time', '1970'), 's')
+    out['forcing']['stop_time'] = np.datetime64(dict_get(c, 'time_control.stop_time', '1970'), 's')
+    out['forcing']['dt'] = dt_sec
+    out['forcing']['ibm_forcing'] = dict_get(c, 'gridforce.ibm_forcing', [])
 
-    # Read output variables
-    outvars = dict()
-    outvar_names = c['output_variables'].get('particle', []) + c['output_variables'].get('instance', [])
+    out['output'] = {}
+    out['output']['file'] = dict_get(c, 'files.output_file')
+    out['output']['frequency'] = dict_get(c, 'output_variables.outper')
+    out['output']['variables'] = {}
+
+    # Convert output variable format spec
+    outvar_names = dict_get(c, 'output_variables.particle', []).copy()
+    outvar_names += dict_get(c, 'output_variables.instance', [])
     for v in outvar_names:
-        outvars[v] = c['output_variables'][v].copy()
-    for v in c['output_variables'].get('particle', []):
-        outvars[v]['kind'] = 'initial'
-    if 'release_time' in outvars and 'units' in outvars['release_time']:
-        outvars['release_time']['units'] = 'seconds since 1970-01-01'
+        out['output']['variables'][v] = c['output_variables'][v].copy()
+        if v == 'release_time' and 'units' in c['output_variables'][v]:
+            out['output']['variables'][v]['units'] = 'seconds since 1970-01-01'
+    for v in dict_get(c, 'output_variables.particle', []):
+        out['output']['variables'][v]['kind'] = 'initial'
+
+    out['tracker'] = {}
+    out['tracker']['method'] = dict_get(c, 'numerics.advection')
+    out['tracker']['diffusion'] = dict_get(c, 'numerics.diffusion')
 
     # Read release config
-    relconf = dict(
-        file=c['files']['particle_release_file'],
-        frequency=c['particle_release'].get('release_frequency', [0, 's']),
-        colnames=c['particle_release']['variables'],
-        formats={
-            c['particle_release'][v]: v
-            for v in c['particle_release']['variables']
-            if v in c['particle_release'].keys()
-        },
-    )
-    if c['particle_release'].get('release_type', '') != 'continuous':
-        del relconf['frequency']
-    ibmvars = c.get('state', dict()).get('ibm_variables', [])
-    ibmvars += c.get('ibm', dict()).get('variables', [])
-    relconf['defaults'] = {
+    out['release'] = {}
+    out['release']['file'] = dict_get(c, 'files.particle_release_file')
+    out['release']['colnames'] = dict_get(c, 'particle_release.variables', [])
+    if dict_get(c, 'particle_release.release_type', '') == 'continuous':
+        out['release']['frequency'] = dict_get(c, 'particle_release.release_frequency', [0, 's'])
+    out['release']['formats'] = {
+        c.get('particle_release', {})[v]: v
+        for v in dict_get(c, 'particle_release.variables', [])
+        if v in c.get('particle_release', {}).keys()
+    }
+    out['release']['defaults'] = {
         k: np.float64(0)
-        for k in ibmvars
-        if k not in relconf['colnames']
+        for k in dict_get(c, 'state.ibm_variables', []) + dict_get(c, 'ibm.variables', [])
+        if k not in out['release']['colnames']
     }
 
-    # Read ibm config
-    ibmconf_legacy = c.get('ibm', dict()).copy()
-    if 'module' in ibmconf_legacy:
-        ibmconf_legacy['ibm_module'] = ibmconf_legacy.pop('module')
-    ibmconf = dict()
-    if 'ibm_module' in ibmconf_legacy:
-        ibmconf['module'] = 'ladim.ibms.LegacyIBM'
-        ibmconf['legacy_module'] = ibmconf_legacy['ibm_module']
-        ibmconf['conf'] = dict(
-            dt=dt_sec,
-            output_instance=c.get('output_variables', {}).get('instance', []),
-            nc_attributes={k: v for k, v in outvars.items()}
-        )
-        ibmconf['conf']['ibm'] = {
+    out['ibm'] = {}
+    if 'ibm' in c:
+        out['ibm']['module'] = 'ladim.ibms.LegacyIBM'
+        out['ibm']['legacy_module'] = dict_get(c, ['ibm.ibm_module', 'ibm.module'])
+        out['ibm']['conf'] = {}
+        out['ibm']['conf']['dt'] = dt_sec
+        out['ibm']['conf']['output_instance'] = dict_get(c, 'output_variables.instance', [])
+        out['ibm']['conf']['nc_attributes'] = {
             k: v
-            for k, v in ibmconf_legacy.items()
-            if k != 'ibm_module'
+            for k, v in out['output']['variables'].items()
         }
+        out['ibm']['conf']['ibm'] = {k: v for k, v in c['ibm'].items() if k != 'ibm_module'}
 
-    config_dict = dict(
-        version=2,
-        solver=dict(
-            start=c['time_control']['start_time'],
-            stop=c['time_control']['stop_time'],
-            step=dt_sec,
-            seed=c['numerics'].get('seed', None),
-            order=['release', 'forcing', 'output', 'tracker', 'ibm', 'state'],
-        ),
-        grid=dict(
-            file=c['gridforce']['input_file'],
-            legacy_module=c['gridforce']['module'] + '.Grid',
-            start_time=np.datetime64(c['time_control']['start_time'], 's'),
-        ),
-        forcing=dict(
-            file=c['gridforce']['input_file'],
-            legacy_module=c['gridforce']['module'] + '.Forcing',
-            start_time=np.datetime64(c['time_control']['start_time'], 's'),
-            stop_time=np.datetime64(c['time_control']['stop_time'], 's'),
-            dt=dt_sec,
-            ibm_forcing=c['gridforce'].get('ibm_forcing', []),
-        ),
-        release=relconf,
-        output=dict(
-            file=c['files']['output_file'],
-            frequency=c['output_variables']['outper'],
-            variables=outvars,
-        ),
-        tracker=dict(
-            method=c['numerics']['advection'],
-            diffusion=c['numerics']['diffusion'],
-        ),
-        ibm=ibmconf,
-    )
-    return config_dict
+    return out

ladim/gridforce/ROMS.py

@@ -15,10 +15,12 @@ import glob
 import logging
 import numpy as np
 from netCDF4 import Dataset, num2date
-
 from ladim.sample import sample2D, bilin_inv
 
 
+logger = logging.getLogger(__name__)
+
+
 class Grid:
     """Simple ROMS grid object
 
@@ -34,7 +36,7 @@ class Grid:
 
     def __init__(self, config):
 
-        logging.info("Initializing ROMS-type grid object")
+        logger.info("Initializing ROMS-type grid object")
 
         # Grid file
         if "grid_file" in config["gridforce"]:
@@ -44,14 +46,14 @@ class Grid:
             files.sort()
             grid_file = files[0]
         else:
-            logging.error("No grid file specified")
+            logger.error("No grid file specified")
             raise SystemExit(1)
 
         try:
             ncid = Dataset(grid_file)
             ncid.set_auto_mask(False)
         except OSError:
-            logging.error("Could not open grid file " + grid_file)
+            logger.error("Could not open grid file " + grid_file)
             raise SystemExit(1)
 
         # Subgrid, only considers internal grid cells
@@ -180,6 +182,8 @@ class Grid:
         """Return the depth of grid cells"""
         I = X.round().astype(int) - self.i0
         J = Y.round().astype(int) - self.j0
+        I = np.minimum(np.maximum(I, 0), self.H.shape[1] - 1)
+        J = np.minimum(np.maximum(J, 0), self.H.shape[0] - 1)
         return self.H[J, I]
 
     def lonlat(self, X, Y, method="bilinear"):
@@ -189,6 +193,8 @@ class Grid:
         # else: containing grid cell, less accurate
         I = X.round().astype("int") - self.i0
         J = Y.round().astype("int") - self.j0
+        I = np.minimum(np.maximum(I, 0), self.lon.shape[1] - 1)
+        J = np.minimum(np.maximum(J, 0), self.lon.shape[0] - 1)
         return self.lon[J, I], self.lat[J, I]
 
     def ingrid(self, X: np.ndarray, Y: np.ndarray) -> np.ndarray:
@@ -204,6 +210,11 @@ class Grid:
         """Returns True for points on land"""
         I = X.round().astype(int) - self.i0
         J = Y.round().astype(int) - self.j0
+
+        # Constrain to valid indices
+        I = np.minimum(np.maximum(I, 0), self.M.shape[-1] - 1)
+        J = np.minimum(np.maximum(J, 0), self.M.shape[-2] - 1)
+
         return self.M[J, I] < 1
 
     # Error if point outside
@@ -242,7 +253,7 @@ class Forcing:
 
     def __init__(self, config, grid):
 
-        logging.info("Initiating forcing")
+        logger.info("Initiating forcing")
 
         self._grid = grid  # Get the grid object, make private?
         # self.config = config["gridforce"]
@@ -255,9 +266,9 @@ class Forcing:
         files = self.find_files(config["gridforce"])
         numfiles = len(files)
         if numfiles == 0:
-            logging.error("No input file: {}".format(config["gridforce"]["input_file"]))
+            logger.error("No input file: {}".format(config["gridforce"]["input_file"]))
             raise SystemExit(3)
-        logging.info("Number of forcing files = {}".format(numfiles))
+        logger.info("Number of forcing files = {}".format(numfiles))
 
         # ---------------------------
         # Overview of all the files
@@ -357,6 +368,7 @@ class Forcing:
         all_frames = []  # All time frames
         num_frames = {}  # Number of time frames in each file
         for fname in files:
+            logger.info(f'Open forcing file {fname}')
             with Dataset(fname) as nc:
                 new_times = nc.variables["ocean_time"][:]
                 num_frames[fname] = len(new_times)
@@ -370,11 +382,11 @@ class Forcing:
         if np.any(I):
             i = I.nonzero()[0][0] + 1   # Index of first out-of-order frame
             oooframe = str(all_frames[i]).split('.')[0]  # Remove microseconds
-            logging.info(f"Time frame {i} = {oooframe} out of order")
-            logging.critical("Forcing time frames not strictly sorted")
+            logger.info(f"Time frame {i} = {oooframe} out of order")
+            logger.critical("Forcing time frames not strictly sorted")
             raise SystemExit(4)
 
-        logging.info(f"Number of available forcing times = {len(all_frames)}")
+        logger.info(f"Number of available forcing times = {len(all_frames)}")
         return all_frames, num_frames
 
     @staticmethod
@@ -382,8 +394,8 @@ class Forcing:
 
         time0 = all_frames[0]
         time1 = all_frames[-1]
-        logging.info(f"First forcing time = {time0}")
-        logging.info(f"Last forcing time = {time1}")
+        logger.info(f"First forcing time = {time0}")
+        logger.info(f"Last forcing time = {time1}")
         start_time = np.datetime64(config["start_time"])
         dt = np.timedelta64(int(config["dt"]), "s")
 
@@ -391,10 +403,10 @@ class Forcing:
         # ------------------------------------------------------
 
         if time0 > start_time:
-            logging.error("No forcing at start time")
+            logger.error("No forcing at start time")
             raise SystemExit(3)
         if time1 < config["stop_time"]:
-            logging.error("No forcing at stop time")
+            logger.error("No forcing at stop time")
             raise SystemExit(3)
 
         # Make a list steps of the forcing time steps
@@ -427,7 +439,7 @@ class Forcing:
         interpolate_velocity_in_time = True
         interpolate_ibm_forcing_in_time = False
 
-        logging.debug("Updating forcing, time step = {}".format(t))
+        logger.debug("Updating forcing, time step = {}".format(t))
         if t in self.steps:  # No time interpolation
             self.U = self.Unew
             self.V = self.Vnew
@@ -486,7 +498,7 @@ class Forcing:
 
         # Handle file opening/closing
         # Always read velocity before other fields
-        logging.info("Reading velocity for time step = {}".format(n))
+        logger.info("Reading velocity for time step = {}".format(n))
 
         # If finished a file or first read (self._nc == "")
         if not self._nc:  # First read
@@ -799,6 +811,11 @@ def sample3D(F, X, Y, K, A, method="bilinear"):
     # else:  method == 'nearest'
     I = X.round().astype("int")
     J = Y.round().astype("int")
+
+    # Constrain to valid indices
+    I = np.minimum(np.maximum(I, 0), F.shape[-1] - 1)
+    J = np.minimum(np.maximum(J, 0), F.shape[-2] - 1)
+
     return F[K, J, I]
 
 

ladim/gridforce/zROMS.py

@@ -18,6 +18,9 @@ from netCDF4 import Dataset, num2date
 from ladim.sample import sample2D
 
 
+logger = logging.getLogger(__name__)
+
+
 class Grid:
     """Simple ROMS grid object
 
@@ -33,13 +36,24 @@ class Grid:
 
     def __init__(self, config):
 
-        logging.info("Initializing zROMS grid object")
+        logger.info("Initializing zROMS grid object")
+
+        # Grid file
+        if "grid_file" in config["gridforce"]:
+            grid_file = config["gridforce"]["grid_file"]
+        elif "input_file" in config["gridforce"]:
+            files = glob.glob(config["gridforce"]["input_file"])
+            files.sort()
+            grid_file = files[0]
+        else:
+            logger.error("No grid file specified")
+            raise SystemExit(1)
+
         try:
-            ncid = Dataset(config["gridforce"]["grid_file"])
+            ncid = Dataset(grid_file)
+            ncid.set_auto_mask(False)
         except OSError:
-            logging.error(
-                "Grid file {} not found".format(config["gridforce"]["grid_file"])
-            )
+            logger.error("Could not open grid file " + grid_file)
             raise SystemExit(1)
 
         # Subgrid, only considers internal grid cells
@@ -188,7 +202,7 @@ class Forcing:
 
     def __init__(self, config, grid):
 
-        logging.info("Initiating forcing")
+        logger.info("Initiating forcing")
 
         self._grid = grid  # Get the grid object, make private?
 
@@ -199,9 +213,9 @@ class Forcing:
         files.sort()
         numfiles = len(files)
         if numfiles == 0:
-            logging.error("No input file: {}".format(config["gridforce"]["input_file"]))
+            logger.error("No input file: {}".format(config["gridforce"]["input_file"]))
             raise SystemExit(3)
-        logging.info("Number of forcing files = {}".format(numfiles))
+        logger.info("Number of forcing files = {}".format(numfiles))
 
         # ----------------------------------------
         # Open first file for some general info
@@ -244,20 +258,20 @@ class Forcing:
         num_frames = []  # Available time frames in each file
         # change_times = []     # Times for change of file
         for fname in files:
-            print(fname)
+            logging.info(f'Load {fname}')
             with Dataset(fname) as nc:
                 # new_times = nc.variables['ocean_time'][:]
                 new_times = nc.variables["time"][:]
                 times.extend(new_times)
                 num_frames.append(len(new_times))
-        logging.info("Number of available forcing times = {:d}".format(len(times)))
+        logger.info("Number of available forcing times = {:d}".format(len(times)))
 
         # Find first/last forcing times
         # -----------------------------
         time0 = num2date(times[0], time_units)
         time1 = num2date(times[-1], time_units)
-        logging.info("time0 = {}".format(str(time0)))
-        logging.info("time1 = {}".format(str(time1)))
+        logger.info("time0 = {}".format(str(time0)))
+        logger.info("time1 = {}".format(str(time1)))
         # print(time0)
         # print(time1)
         start_time = np.datetime64(config["start_time"])
@@ -269,10 +283,10 @@ class Forcing:
         # Use logging module for this
 
         if time0 > start_time:
-            logging.error("No forcing at start time")
+            logger.error("No forcing at start time")
             raise SystemExit(3)
         if time1 < config["stop_time"]:
-            logging.error("No forcing at stop time")
+            logger.error("No forcing at stop time")
             raise SystemExit(3)
 
         # Make a list steps of the forcing time steps
@@ -362,7 +376,7 @@ class Forcing:
         interpolate_velocity_in_time = True
         interpolate_ibm_forcing_in_time = False
 
-        logging.debug("Updating forcing, time step = {}".format(t))
+        logger.debug("Updating forcing, time step = {}".format(t))
         if t in self.steps:  # No time interpolation
             self.U = self.Unew
             self.V = self.Vnew
@@ -399,7 +413,7 @@ class Forcing:
 
         # Handle file opening/closing
         # Always read velocity before other fields
-        logging.debug("Reading velocity for time step = {}".format(n))
+        logger.debug("Reading velocity for time step = {}".format(n))
         first = True
         if first:  # Open file initiallt
             self._nc = Dataset(self._files[self.file_idx[n]])

ladim/main.py

@@ -12,13 +12,14 @@ Lagrangian Advection and Diffusion Model
 # ---------------------------------
 
 import logging
-
 import ladim
-
 from .config import configure
 from .model import Model
 
 
+logger = logging.getLogger(__name__)
+
+
 def main(config_stream, loglevel=logging.INFO):
     """Main function for LADiM"""
 
@@ -46,7 +47,9 @@ def run():
 
     logging.basicConfig(
         level=logging.INFO,
-        format='%(levelname)s:%(module)s - %(message)s')
+        format='%(asctime)s  %(levelname)s:%(module)s - %(message)s',
+        datefmt='%Y-%m-%d %H:%M:%S',
+    )
 
     # ====================
     # Parse command line
@@ -67,23 +70,23 @@ def run():
 
     args = parser.parse_args()
 
-    logging.info(" ================================================")
-    logging.info(" === Lagrangian Advection and Diffusion Model ===")
-    logging.info(" ================================================\n")
+    logger.info(" ================================================")
+    logger.info(" === Lagrangian Advection and Diffusion Model ===")
+    logger.info(" ================================================\n")
 
-    logging.info(f"ladim path: {ladim.__file__.strip('__init.py__')}")
-    logging.info(f"ladim version:  {ladim.__version__}\n")
-    logging.info(f"python version:  {sys.version.split()[0]}\n")
+    logger.info(f"ladim path: {ladim.__file__.strip('__init.py__')}")
+    logger.info(f"ladim version:  {ladim.__version__}\n")
+    logger.info(f"python version:  {sys.version.split()[0]}\n")
 
-    logging.info(f"  Configuration file: {args.config_file}")
-    logging.info(f"  loglevel = {logging.getLevelName(args.loglevel)}")
+    logger.info(f"  Configuration file: {args.config_file}")
+    logger.info(f"  loglevel = {logging.getLevelName(args.loglevel)}")
 
     # =============
     # Sanity check
     # =============
 
     if not Path(args.config_file).exists():
-        logging.critical(f'Configuration file {args.config_file} not found')
+        logger.critical(f'Configuration file {args.config_file} not found')
         raise SystemExit(1)
 
     # ===================
@@ -92,7 +95,7 @@ def run():
 
     # Start message
     now = datetime.datetime.now().replace(microsecond=0)
-    logging.info(f'LADiM simulation starting, wall time={now}')
+    logger.info(f'LADiM simulation starting, wall time={now}')
 
     fp = open(args.config_file, encoding='utf8')
     ladim.main(config_stream=fp, loglevel=args.loglevel)
@@ -100,4 +103,4 @@ def run():
     # Reset logging and print final message
     logging.getLogger().setLevel(logging.INFO)
     now = datetime.datetime.now().replace(microsecond=0)
-    logging.info(f'LADiM simulation finished, wall time={now}')
+    logger.info(f'LADiM simulation finished, wall time={now}')

ladim/model.py

@@ -1,4 +1,7 @@
 import importlib
+import importlib.util
+import sys
+from pathlib import Path
 
 from typing import TYPE_CHECKING
 if TYPE_CHECKING:
@@ -101,7 +104,26 @@ class Model:
 
 def load_class(name):
     pkg, cls = name.rsplit(sep='.', maxsplit=1)
-    return getattr(importlib.import_module(pkg), cls)
+
+    # Check if "pkg" is an existing file
+    spec = None
+    module_name = None
+    file_name = pkg + '.py'
+    if Path(file_name).exists():
+        # This can return None if there were import errors
+        module_name = pkg
+        spec = importlib.util.spec_from_file_location(module_name, file_name)
+
+    # If pkg can not be interpreted as a file, use regular import
+    if spec is None:
+        return getattr(importlib.import_module(pkg), cls)
+
+    # File import
+    else:
+        module = importlib.util.module_from_spec(spec)
+        sys.modules[module_name] = module
+        spec.loader.exec_module(module)
+        return getattr(module, cls)
 
 
 class Module:

ladim/state.py

@@ -122,6 +122,14 @@ class DynamicState(State):
             raise AttributeError(f'Attribute not defined: {item}')
         return self[item]
 
+    def __setattr__(self, item, value):
+        if item in list(self.__dict__.keys()) + ['_data', '_model', '_num_released', '_varnames']:
+            super().__setattr__(item, value)
+        elif item in self._data:
+            self._data[item] = value
+        else:
+            raise AttributeError(f"Attribute not defined: '{item}'")
+
     def __contains__(self, item):
         return item in self._data
 

{ladim-2.0.1.dist-info → ladim-2.0.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: ladim
-Version: 2.0.1
+Version: 2.0.3
 Summary: Lagrangian Advection and Diffusion Model
 Home-page: https://github.com/pnsaevik/ladim
 Author: Bjørn Ådlandsvik

{ladim-2.0.1.dist-info → ladim-2.0.3.dist-info}/RECORD

@@ -1,21 +1,21 @@
-ladim/__init__.py,sha256=gItkgENsZYTMR2SxeAgoMMvJ9uTjFAJA29jZoVWvM2k,51
+ladim/__init__.py,sha256=Vn8lFql5FOG5tGQ1TDaKybIZLoYlLHVXvLcitN1QO04,51
 ladim/__main__.py,sha256=8f07EMfxQllDZSgpak5ECyYHnfQFy8LaHl2xdC-aO9c,23
-ladim/config.py,sha256=6b0ikBrBnq_sSElgO2YcVtJHjQU0nyZsVyPV3q9fy5I,4233
+ladim/config.py,sha256=AAmQflh_u8T-0Qudij_cGpuc_WHOKouoLTBtMZw_Bh0,5075
 ladim/forcing.py,sha256=f4PpSwyilSScXeNyorTWLMgVTiat9htSLkCwAkRlJVM,3048
 ladim/grid.py,sha256=m6bQrGJ3cux7rqC8pbRXD86cOI-VQKF-XjP9m1jCIcY,2221
-ladim/main.py,sha256=wO91-nLd1gvF3V20XK5qRvvOIV4xoTOKiWFcbwV2oag,2812
-ladim/model.py,sha256=jpjq_ZSh7ULpwi3_RqDb-p5SG8WcdgCPaBkSpnNWblU,3137
+ladim/main.py,sha256=6_blu3PYnDXaYdPxfZoukWsjN0o9vh7O8_-W2-aguAI,2894
+ladim/model.py,sha256=iXClvieChhipCSZ-dDrmnjqwS4cuM53VpJv7oaJyQ88,3794
 ladim/output.py,sha256=Rz7iujvS7Z3LoABiJduQqyb3zPswNqhhFsywr3MLsBY,8373
 ladim/release.py,sha256=1j__9Gj0BD0CqVCM2KLZhio1Ia-hz1gbUIhTsa0J3Rg,8451
 ladim/sample.py,sha256=n8wRGd_VsW_qyQe1ZoTpmfZcdcwB929vsM8PoKG6JTs,8292
 ladim/solver.py,sha256=sZvYgOxzJ-EItI-IB2y8_z8Tf-SJAQSrmydlhDRa7ZQ,755
-ladim/state.py,sha256=5ICIiujsV3KOAUYagGLK7YdmhcItgJmtntZeR11nIpw,3781
+ladim/state.py,sha256=4XNIIx5sGjlqkZ6bg-dGbqzp8ujFNkHHFL2D9qCQA2w,4119
 ladim/tracker.py,sha256=VVX6T5CqiU6nGSCgLlSCC8w0UYhW273OGFE7ApPjdyI,5091
 ladim/utilities.py,sha256=r7-zShqJhh0cBctDUmtfw-GBOk1eTTYR4S72b0ouiSQ,994
-ladim/gridforce/ROMS.py,sha256=yrMr1GiDgr7VG9V630nHFdeajPY4WHLu653ZRwAk1aw,26888
+ladim/gridforce/ROMS.py,sha256=DF5CSR2iJsPWAmTOrqYavARVtqokQbtAWLSCDWb_9f0,27525
 ladim/gridforce/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 ladim/gridforce/analytical.py,sha256=qI-1LJdjmnwwanzOdrsDZqwGgo73bT75CB7pMaxbHKM,1094
-ladim/gridforce/zROMS.py,sha256=MVA6PQuY1ukvs2E20sWY4kr5-QieeQHTrA5ruxCqbUM,22826
+ladim/gridforce/zROMS.py,sha256=217lDAHY0H7vbez6iLgIHQxK6MrNaZlAtWo6MteKVxU,23225
 ladim/ibms/__init__.py,sha256=GOG75jZDmNEiLr8brxrKqIlqVj-pNR7pnPP8FUKE6hU,565
 ladim/ibms/light.py,sha256=POltHmKkX8-q3t9wXyfcseCKEq9Bq-kX1WEJYsr1lNQ,2737
 ladim/plugins/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -24,9 +24,9 @@ postladim/cellcount.py,sha256=nCFu9iJmprubn4YmPB4W0VO02GfEb90Iif7D49w1Kss,2054
 postladim/kde_plot.py,sha256=GvMWzT6VxIeXKh1cnqaGzR-4jGG_WIHGMLPpRMXIpo4,1628
 postladim/particlefile.py,sha256=0aif9wYUJ-VrpQKeCef8wB5VCiBB-gWY6sxNCUYviTA,4889
 postladim/variable.py,sha256=-2aihoppYMMmpSpCqaF31XvpinTMaH3Y01-USDIkbBc,6587
-ladim-2.0.1.dist-info/LICENSE,sha256=BgtXyjNr6Ly9nQ7ZLXKpV3r5kWRLnh5MiN0dxp0Bvfc,1085
-ladim-2.0.1.dist-info/METADATA,sha256=-183UhMhVbBOfOcjonoP1P8l2q500OF8vexlV1SDAEg,1841
-ladim-2.0.1.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-ladim-2.0.1.dist-info/entry_points.txt,sha256=JDlNJo87GJaOkH0-BpAzTPLCrZcuPSdSlHNQ4XmnoRg,41
-ladim-2.0.1.dist-info/top_level.txt,sha256=TK8Gl7d6MsrAQvqKG4b6YJCbB4UL46Se3SzsI-sJAuc,16
-ladim-2.0.1.dist-info/RECORD,,
+ladim-2.0.3.dist-info/LICENSE,sha256=BgtXyjNr6Ly9nQ7ZLXKpV3r5kWRLnh5MiN0dxp0Bvfc,1085
+ladim-2.0.3.dist-info/METADATA,sha256=QjbEu8Z5DysRolgZ7RuuvAT4CZ0S_8bWPP49ByNIQIE,1841
+ladim-2.0.3.dist-info/WHEEL,sha256=GV9aMThwP_4oNCtvEC2ec3qUYutgWeAzklro_0m4WJQ,91
+ladim-2.0.3.dist-info/entry_points.txt,sha256=JDlNJo87GJaOkH0-BpAzTPLCrZcuPSdSlHNQ4XmnoRg,41
+ladim-2.0.3.dist-info/top_level.txt,sha256=TK8Gl7d6MsrAQvqKG4b6YJCbB4UL46Se3SzsI-sJAuc,16
+ladim-2.0.3.dist-info/RECORD,,

{ladim-2.0.1.dist-info → ladim-2.0.3.dist-info}/WHEEL

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