chore: logging 5 - remove prints from local_modules directory

src/icesat2_tracks/local_modules/JONSWAP_gamma.py

@@ -1,3 +1,4 @@
+import logging
 import os
 from scipy.constants import g
 
@@ -16,6 +17,7 @@
 import matplotlib.pyplot as plt
 import numpy as np
 
+_logger = logging.getLogger(__name__)
 
 # %%
 def normalize_time(time):
@@ -77,7 +79,7 @@ def gamma_time_normlized_amp_shifted(time, gammapar=2, loc=0.2, scale=0.1):
     """
     gamma_mod = gamma(gammapar, loc=loc[0, :], scale=scale).pdf(time[0, :])
     dloc = loc[0, :] - time[0, :][gamma_mod.argmax()]
-    # print(dloc)
+    # _logger.debug(dloc)
     gamma_mod = gamma(gammapar, loc=loc + dloc, scale=scale).pdf(time)
     return gamma_mod / gamma_mod.max()
 
@@ -140,7 +142,7 @@ def pierson_moskowitz_fetch_limit(f, X, U):
 
     alpha = 0.076 * (g * X / U**2) ** (-0.22)  # non-dimentional
     wp = 7.0 * np.pi * (g / U) * (g * X / U**2) ** (-0.33)
-    print("wp=" + str(wp))
+    _logger.debug("wp=%s", wp)
 
     sigma_p = 0.07
     sigma_pp = 0.09
@@ -173,8 +175,8 @@ def JONSWAP_default(f, X, U, gamma=3.3):
 
     alpha = 0.076 * (g * X / U**2) ** (-0.22)  # non-dimentional
     wp = 7.0 * np.pi * (g / U) * (g * X / U**2) ** (-0.33)
-    # print('wp='+str(wp))
-    # print('nu='+str(wp*U/(g)))
+    # _logger.debug('wp='+str(wp))
+    # _logger.debug('nu='+str(wp*U/(g)))
 
     sigma_p = 0.07
     sigma_pp = 0.09
@@ -257,21 +259,21 @@ def gamma_time_JONSWAP_default(
     # intersectf=intersect-intersect/slopet
     intersectF = -intersectT * slope_t
     pfreq = time * slope_t + intersectF
-    # print('intersect F=' + str(intersectF))
+    # _logger.debug('intersect F=' + str(intersectF))
 
     # intersectf=intersect#-intersect/slope
     slopeF = 1 / slope_t
     pfreq_forgamma = f * slopeF + intersectT
 
     # rint(pfreq.shape)
     tt, line = np.meshgrid(time, pfreq_forgamma)
-    # print(line)
-    # print(tt)
+    # _logger.debug(line)
+    # _logger.debug(tt)
     func_t = gamma_time_normlized_amp_shifted(
         tt, gammapar=tgammapar, loc=line, scale=tscale
     )
     # func_t_temp= (tamp*np.exp(- (time-t_center )**4 / tsigma ))
-    # print(func_t.shape)
+    # _logger.debug(func_t.shape)
 
     """ Define X(f_max and U) here """
 
@@ -338,21 +340,21 @@ def gamma_time_JONSWAP_nondim(
     # intersectf=intersect-intersect/slopet
     intersectF = -intersectT * slope_t
     pfreq = time * slope_t + intersectF
-    # print('intersect F=' + str(intersectF))
+    # _logger.debug('intersect F=' + str(intersectF))
 
     # intersectf=intersect#-intersect/slope
     slopeF = 1 / slope_t
     pfreq_forgamma = f * slopeF + intersectT
 
     # rint(pfreq.shape)
     tt, line = np.meshgrid(time, pfreq_forgamma)
-    # print(line)
-    # print(tt)
+    # _logger.debug(line)
+    # _logger.debug(tt)
     func_t = gamma_time_normlized_amp_shifted(
         tt, gammapar=tgammapar, loc=line, scale=tscale
     )
     # func_t_temp= (tamp*np.exp(- (time-t_center )**4 / tsigma ))
-    # print(func_t.shape)
+    # _logger.debug(func_t.shape)
 
     def X(f_max, U10):
         return 3.5**3.0 * g**2.0 / (U10 * f_max**3.0)
@@ -445,7 +447,7 @@ def residual_JONSWAP_default_gamma(
     # tt, tt= np.meshgrid(time, ff)
     model1d = model.reshape(model.shape[0] * model.shape[1])
     model1d[np.isnan(model1d)] = 0
-    # print(model1d)
+    # _logger.debug(model1d)
 
     if data is not None:
         if np.size(data.shape) != 1:
@@ -475,7 +477,7 @@ def residual_JONSWAP_default_gamma(
     if data is None:
         return model1d
     if weight is not None:
-        # print('use weight')
+        # _logger.debug('use weight')
         d = (model1d - data1d) * weight1d
         d[nan_track] = np.nan
         return d
@@ -557,7 +559,7 @@ def residual_JONSWAP_nondim_gamma(
     if data is None:
         return model1d
     if weight is not None:
-        # print('use weight')
+        # _logger.debug('use weight')
         d = (model1d - data1d) * weight1d
         d[nan_track] = np.nan
         return d
@@ -719,7 +721,7 @@ def Jm_regulizer(value_dict, prior):
 
     Jm = list()
     for k, I in prior.items():
-        # print(I)
+        # _logger.debug(I)
         if type(I["m_err"]) is float:
             Jm.append((I["m0"] - vd[k]) / I["m_err"])
         else:
@@ -770,7 +772,7 @@ def reg_func(p0, pi, p_err_unit):
             return (pi / p0 - 1) / ppercent
 
     for k, I in prior.items():
-        # print(I)
+        # _logger.debug(I)
         if type(I["m_err"]) is float:
             Jm.append(reg_func(I["m0"], vd[k], I["m_err"]))
             # Jm.append(    (I['m0']- vd[k] ) / I['m_err']    )
@@ -804,7 +806,7 @@ def reg_func(p0, pi, p_err_unit):
         vd[k] = I * np.random.rand()
 
     Jm = Jm_regulizer(vd, priors)
-    print(Jm)
+    _logger.debug(Jm)
 
 
 # %%
@@ -943,7 +945,7 @@ def cost(value_dict, time, f, data=None, weight=None, prior=False, eps=None):
 #     if prior is not None:
 #         Jm=Jm_regulizer(vdict, prior, alpha=5.0)
 #
-#     #print(Jm)
+#     #_logger.debug(Jm)
 #
 #     if data is None:
 #         return model1d
@@ -955,7 +957,7 @@ def cost(value_dict, time, f, data=None, weight=None, prior=False, eps=None):
 #     if (weight is not None) and (eps is None):
 #         d=(model1d - data1d)*weight
 #         d[nan_track]=np.nan
-#         #print(np.concatenate(( d,Jm )) )
+#         #_logger.debug(np.concatenate(( d,Jm )) )
 #         return np.concatenate(( d, Jm ))
 #
 #     if (eps is None) and (weight is None):

src/icesat2_tracks/local_modules/m_colormanager_ph3.py

@@ -19,7 +19,7 @@ def ase_to_json(path):
     A=swatch.parse(path)
 
     for i in A[0]['swatches']:
-        _logger.debug(i['name'] + '  ' + str(i['data']['values']))
+        _logger.debug("%s %s", i['name'], i['data']['values'])
 
     return A
 
@@ -32,7 +32,7 @@ def json_save(name, path, data, return_name=False):
     full_name= (os.path.join(full_name_root+ '.json'))
     with open(full_name, 'w') as outfile:
         json.dump(data, outfile)
-    _logger.debug(f'save at: {full_name}')
+    _logger.debug('save at: %s', full_name)
     if return_name:
         return full_name_root
     else:
@@ -45,7 +45,7 @@ def json_load(name, path):
 
     with open(full_name, 'r') as ifile:
         data=json.load(ifile)
-    _logger.debug('loaded from: {full_name}')
+    _logger.debug('loaded from: %s', full_name)
     return data
 
 
@@ -55,7 +55,7 @@ class color:
         def __init__(self, path=None, name=None):
             self.white=(1,1,1)
             if (path is not None) & (name is not None):
-                _logger.debug('color theme: '+name)
+                _logger.debug('color theme: %s', name)
                 try:
                     theme=json_load(name, path)
                     for k, v in theme.items():
@@ -177,13 +177,13 @@ def colormaps(self,n, gamma=None):
 
         def show(self):
             for key in self.__dict__.keys():
-                _logger.debug(key)
+                _logger.debug("%s",key)
 
             _logger.debug('  rels dict:')
             for key in self.rels.keys():
-                _logger.debug('  '+key)
+                _logger.debug('  %s', key)
 
-            #print(self.__dict__)
+            #_logger.debug(self.__dict__)
         def plot(self):
             dd=self.__dict__.copy()
             dd_colors=dd['rels']