refactor: use print() instead of sys.stdout.write()

flopy/mbase.py

@@ -5,7 +5,6 @@
 
 """
 import abc
-import sys
 import os
 import shutil
 import threading
@@ -1203,8 +1202,8 @@ def change_model_ws(self, new_pth=None, reset_external=False):
         # --reset the model workspace
         old_pth = self._model_ws
         self._model_ws = new_pth
-        line = f"\nchanging model workspace...\n   {new_pth}\n"
-        sys.stdout.write(line)
+        if self.verbose:
+            print(f"\nchanging model workspace...\n   {new_pth}")
         # reset the paths for each package
         for pp in self.packagelist:
             pp.fn_path = os.path.join(self.model_ws, pp.file_name[0])

flopy/modflow/mfbas.py

@@ -9,7 +9,6 @@
 """
 
 import re
-import sys
 import numpy as np
 from ..pakbase import Package
 from ..utils import Util3d
@@ -324,7 +323,7 @@ def load(cls, f, model, ext_unit_dict=None, check=True, **kwargs):
         """
 
         if model.verbose:
-            sys.stdout.write("loading bas6 package file...\n")
+            print("loading bas6 package file...")
 
         # parse keywords
         if "nlay" in kwargs:

flopy/modflow/mfbcf.py

@@ -1,5 +1,3 @@
-import sys
-
 import numpy as np
 
 from ..pakbase import Package
@@ -349,7 +347,7 @@ def load(cls, f, model, ext_unit_dict=None):
         """
 
         if model.verbose:
-            sys.stdout.write("loading bcf package file...\n")
+            print("loading bcf package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:

flopy/modflow/mfchd.py

@@ -8,7 +8,6 @@
 
 """
 
-import sys
 import numpy as np
 from ..pakbase import Package
 from ..utils import MfList
@@ -266,7 +265,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading chd package file...\n")
+            print("loading chd package file...")
 
         return Package.load(
             f,

flopy/modflow/mfde4.py

@@ -7,7 +7,6 @@
 <http://water.usgs.gov/ogw/modflow-nwt/MODFLOW-NWT-Guide/de4.htm>`_.
 
 """
-import sys
 from ..pakbase import Package
 
 
@@ -256,7 +255,7 @@ def load(cls, f, model, ext_unit_dict=None):
         """
 
         if model.verbose:
-            sys.stdout.write("loading de4 package file...\n")
+            print("loading de4 package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:

flopy/modflow/mfdis.py

@@ -7,8 +7,6 @@
 <http://water.usgs.gov/ogw/modflow/MODFLOW-2005-Guide/index.html?dis.htm>`_.
 
 """
-
-import sys
 import warnings
 
 import numpy as np
@@ -763,7 +761,7 @@ def load(cls, f, model, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading dis package file...\n")
+            print("loading dis package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:

flopy/modflow/mfdisu.py

@@ -3,9 +3,6 @@
 the ModflowDisU class as `flopy.modflow.ModflowDisU`.
 
 """
-
-import sys
-import warnings
 import numpy as np
 from ..pakbase import Package
 from ..utils import Util2d, Util3d, read1d
@@ -552,7 +549,7 @@ def load(cls, f, model, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading disu package file...\n")
+            print("loading disu package file...")
 
         if model.version != "mfusg":
             print(

flopy/modflow/mfdrn.py

@@ -7,7 +7,6 @@
 <http://water.usgs.gov/ogw/modflow/MODFLOW-2005-Guide/index.html?drn.htm>`_.
 
 """
-import sys
 import numpy as np
 from ..pakbase import Package
 from ..utils.util_list import MfList
@@ -323,7 +322,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading drn package file...\n")
+            print("loading drn package file...")
 
         return Package.load(
             f,

flopy/modflow/mfdrt.py

@@ -7,7 +7,6 @@
 <http://water.usgs.gov/ogw/modflow/MODFLOW-2005-Guide/index.html?drt.htm>`_.
 
 """
-import sys
 import numpy as np
 from ..pakbase import Package
 from ..utils.util_list import MfList
@@ -309,7 +308,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading drt package file...\n")
+            print("loading drt package file...")
 
         return Package.load(
             f,

flopy/modflow/mfevt.py

@@ -7,8 +7,6 @@
 <http://water.usgs.gov/ogw/modflow/MODFLOW-2005-Guide/index.html?evt.htm>`_.
 
 """
-import sys
-
 import numpy as np
 from .mfparbc import ModflowParBc as mfparbc
 from ..utils import Transient2d, Util2d
@@ -267,7 +265,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
 
         """
         if model.verbose:
-            sys.stdout.write("loading evt package file...\n")
+            print("loading evt package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:

flopy/modflow/mffhb.py

@@ -7,8 +7,6 @@
 <http://water.usgs.gov/ogw/modflow/MODFLOW-2005-Guide/index.html?fhb.htm>`_.
 
 """
-import sys
-
 import numpy as np
 
 from ..pakbase import Package
@@ -455,7 +453,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
 
         """
         if model.verbose:
-            sys.stdout.write("loading fhb package file...\n")
+            print("loading fhb package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:
@@ -477,7 +475,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
 
         # dataset 1
         if model.verbose:
-            sys.stdout.write("loading fhb dataset 1\n")
+            print("loading fhb dataset 1")
         raw = line.strip().split()
         nbdtim = int(raw[0])
         nflw = int(raw[1])
@@ -493,10 +491,8 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
         flow_aux = []
         if nfhbx1 > 0:
             if model.verbose:
-                sys.stdout.write("loading fhb dataset 2\n")
-            sys.stdout.write(
-                "dataset 2 will not be preserved in the created hfb object.\n"
-            )
+                print("loading fhb dataset 2")
+            print("dataset 2 will not be preserved in the created hfb object.")
             for idx in range(nfhbx1):
                 line = f.readline()
                 raw = line.strip().split()
@@ -510,10 +506,8 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
         head_aux = []
         if nfhbx2 > 0:
             if model.verbose:
-                sys.stdout.write("loading fhb dataset 3\n")
-            sys.stdout.write(
-                "dataset 3 will not be preserved in the created hfb object.\n"
-            )
+                print("loading fhb dataset 3")
+            print("dataset 3 will not be preserved in the created hfb object.")
             for idx in range(nfhbx2):
                 line = f.readline()
                 raw = line.strip().split()
@@ -525,7 +519,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
 
         # Dataset 4a IFHBUN CNSTM IFHBPT
         if model.verbose:
-            sys.stdout.write("loading fhb dataset 4a\n")
+            print("loading fhb dataset 4a")
         line = f.readline()
         raw = line.strip().split()
         ifhbun = int(raw[0])
@@ -539,7 +533,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
 
         # Dataset 4b
         if model.verbose:
-            sys.stdout.write("loading fhb dataset 4b\n")
+            print("loading fhb dataset 4b")
         line = f.readline()
         raw = line.strip().split()
         bdtime = []
@@ -553,7 +547,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
         ds6 = None
         if nflw > 0:
             if model.verbose:
-                sys.stdout.write("loading fhb dataset 5a\n")
+                print("loading fhb dataset 5a")
             # Dataset 5a IFHBUN CNSTM IFHBPT
             line = f.readline()
             raw = line.strip().split()
@@ -567,7 +561,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
             ifhbpt = max(ifhbpt, int(raw[2]))
 
             if model.verbose:
-                sys.stdout.write("loading fhb dataset 5b\n")
+                print("loading fhb dataset 5b")
             dtype = ModflowFhb.get_default_dtype(
                 nbdtim=nbdtim, head=False, structured=model.structured
             )
@@ -598,12 +592,10 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
                     dtype.append((name, np.float32))
                 for naux in range(nfhbx1):
                     if model.verbose:
-                        sys.stdout.write(
-                            f"loading fhb dataset 6a - aux {naux + 1}\n"
-                        )
-                    sys.stdout.write(
+                        print(f"loading fhb dataset 6a - aux {naux + 1}")
+                    print(
                         "dataset 6a will not be preserved in "
-                        "the created hfb object.\n"
+                        "the created hfb object."
                     )
                     # Dataset 6a IFHBUN CNSTM IFHBPT
                     line = f.readline()
@@ -619,12 +611,10 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
                     ifhbpt = max(ifhbpt, int(raw[2]))
 
                     if model.verbose:
-                        sys.stdout.write(
-                            f"loading fhb dataset 6b - aux {naux + 1}\n"
-                        )
-                    sys.stdout.write(
+                        print(f"loading fhb dataset 6b - aux {naux + 1}")
+                    print(
                         "dataset 6b will not be preserved in "
-                        "the created hfb object.\n"
+                        "the created hfb object."
                     )
                     current = np.recarray(nflw, dtype=dtype)
                     for n in range(nflw):
@@ -640,7 +630,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
         ds8 = None
         if nhed > 0:
             if model.verbose:
-                sys.stdout.write("loading fhb dataset 7a\n")
+                print("loading fhb dataset 7a")
             # Dataset 7a IFHBUN CNSTM IFHBPT
             line = f.readline()
             raw = line.strip().split()
@@ -654,7 +644,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
             ifhbpt = max(ifhbpt, int(raw[2]))
 
             if model.verbose:
-                sys.stdout.write("loading fhb dataset 7b\n")
+                print("loading fhb dataset 7b")
             dtype = ModflowFhb.get_default_dtype(
                 nbdtim=nbdtim, head=True, structured=model.structured
             )
@@ -685,12 +675,10 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
                     dtype.append((name, np.float32))
                 for naux in range(nfhbx1):
                     if model.verbose:
-                        sys.stdout.write(
-                            f"loading fhb dataset 8a - aux {naux + 1}\n"
-                        )
-                    sys.stdout.write(
+                        print(f"loading fhb dataset 8a - aux {naux + 1}")
+                    print(
                         "dataset 8a will not be preserved in "
-                        "the created hfb object.\n"
+                        "the created hfb object."
                     )
                     # Dataset 6a IFHBUN CNSTM IFHBPT
                     line = f.readline()
@@ -707,10 +695,8 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
                     ifhbpt = max(ifhbpt, ifhbpt6)
 
                     if model.verbose:
-                        sys.stdout.write(
-                            f"loading fhb dataset 8b - aux {naux + 1}\n"
-                        )
-                    sys.stdout.write(
+                        print(f"loading fhb dataset 8b - aux {naux + 1}")
+                    print(
                         "dataset 8b will not be preserved in "
                         "the created hfb object."
                     )

flopy/modflow/mfflwob.py

@@ -1,5 +1,4 @@
 import os
-import sys
 import numpy as np
 from ..pakbase import Package
 from ..utils import parsenamefile
@@ -414,7 +413,7 @@ def load(cls, f, model, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading flwob package file...\n")
+            print("loading flwob package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:

flopy/modflow/mfgage.py

@@ -8,7 +8,6 @@
 
 """
 import os
-import sys
 
 import numpy as np
 
@@ -318,7 +317,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None):
         """
 
         if model.verbose:
-            sys.stdout.write("loading gage package file...\n")
+            print("loading gage package file...")
 
         openfile = not hasattr(f, "read")
         if openfile:

flopy/modflow/mfghb.py

@@ -7,7 +7,6 @@
 <http://water.usgs.gov/ogw/modflow/MODFLOW-2005-Guide/index.html?ghb.htm>`_.
 
 """
-import sys
 import numpy as np
 from ..pakbase import Package
 from ..utils import MfList
@@ -301,7 +300,7 @@ def load(cls, f, model, nper=None, ext_unit_dict=None, check=True):
         """
 
         if model.verbose:
-            sys.stdout.write("loading ghb package file...\n")
+            print("loading ghb package file...")
 
         return Package.load(
             f,