ladim 2.0.1__tar.gz → 2.0.2__tar.gz

{ladim-2.0.1 → ladim-2.0.2}/PKG-INFO

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

{ladim-2.0.1 → ladim-2.0.2}/ladim/__init__.py

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

{ladim-2.0.1 → ladim-2.0.2}/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
@@ -242,7 +244,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 +257,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 +359,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 +373,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 +385,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 +394,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 +430,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 +489,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

{ladim-2.0.1 → ladim-2.0.2}/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,11 +36,11 @@ class Grid:
 
     def __init__(self, config):
 
-        logging.info("Initializing zROMS grid object")
+        logger.info("Initializing zROMS grid object")
         try:
             ncid = Dataset(config["gridforce"]["grid_file"])
         except OSError:
-            logging.error(
+            logger.error(
                 "Grid file {} not found".format(config["gridforce"]["grid_file"])
             )
             raise SystemExit(1)
@@ -188,7 +191,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 +202,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
@@ -250,14 +253,14 @@ class Forcing:
                 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 +272,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 +365,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 +402,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-2.0.1 → ladim-2.0.2}/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"""
 
@@ -67,23 +68,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 +93,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 +101,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-2.0.1 → ladim-2.0.2}/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 → ladim-2.0.2}/ladim.egg-info/PKG-INFO

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