pypromice 1.6.0__py3-none-any.whl → 1.7.0__py3-none-any.whl

pypromice/core/qc/github_data_issues.py

@@ -159,7 +159,7 @@ def adjustData(ds, adj_dir, var_list=[], skip_var=[]):
         adj_info.loc[adj_info.t0.isnull()|(adj_info.t0==''), "t0"] = None
 		
         # if "*" is in the variable name then we interpret it as regex
-        selec =  adj_info['variable'].str.contains('\*') & (adj_info['variable'] != "*")
+        selec = adj_info['variable'].str.contains(r'\*') & (adj_info['variable'] != "*")
         for ind in adj_info.loc[selec, :].index:
             line_template = adj_info.loc[ind, :].copy()
             regex = adj_info.loc[ind, 'variable']
@@ -209,23 +209,11 @@ def adjustData(ds, adj_dir, var_list=[], skip_var=[]):
 
                 if func == "add":
                     ds_out[var].loc[index_slice] = ds_out[var].loc[index_slice].values + val
-                    # flagging adjusted values
-                    # if var + "_adj_flag" not in ds_out.columns:
-                    #     ds_out[var + "_adj_flag"] = 0
-                    # msk = ds_out[var].loc[index_slice])].notnull()
-                    # ind = ds_out[var].loc[index_slice])].loc[msk].time
-                    # ds_out.loc[ind, var + "_adj_flag"] = 1
 
                 if func == "multiply":
                     ds_out[var].loc[index_slice] = ds_out[var].loc[index_slice].values * val
                     if "DW" in var:
                         ds_out[var].loc[index_slice] = ds_out[var].loc[index_slice] % 360
-                    # flagging adjusted values
-                    # if var + "_adj_flag" not in ds_out.columns:
-                    #     ds_out[var + "_adj_flag"] = 0
-                    # msk = ds_out[var].loc[index_slice].notnull()
-                    # ind = ds_out[var].loc[index_slice].loc[msk].time
-                    # ds_out.loc[ind, var + "_adj_flag"] = 1
 
                 if func == "min_filter":
                     tmp = ds_out[var].loc[index_slice].values
@@ -277,6 +265,27 @@ def adjustData(ds, adj_dir, var_list=[], skip_var=[]):
                     ds_out[var2].loc[index_slice] = val_var
                     ds_out[var].loc[index_slice] = val_var2
 
+                if "delete_when_same_as_" in func:
+                    var2 = func.replace('delete_when_same_as_','')
+                    tmp = ds_out[var].loc[index_slice]
+                    msk = np.abs(tmp - ds_out[var2].loc[index_slice]) < val
+                    tmp = tmp.where(~msk)
+                    # remove isolated singletons and pairs surrounded by NaNs
+                    m1 = tmp.notnull() & tmp.shift(time=1).isnull() & tmp.shift(time=-1).isnull()
+                    
+                    m2_first  = (tmp.notnull()
+                                 & tmp.shift(time=1).isnull()     # left is NaN
+                                 & tmp.shift(time=-1).notnull()   # right is value
+                                 & tmp.shift(time=-2).isnull())   # right+1 is NaN
+                    
+                    m2_second = (tmp.notnull()
+                                 & tmp.shift(time=-1).isnull()    # right is NaN
+                                 & tmp.shift(time=1).notnull()    # left is value
+                                 & tmp.shift(time=2).isnull())    # left-1 is NaN
+                    
+                    tmp = tmp.where(~(m1 | m2_first | m2_second))
+                    ds_out[var].loc[index_slice] = tmp.values
+
                 if func == "rotate":
                     ds_out[var].loc[index_slice] = (ds_out[var].loc[index_slice].values + val) % 360
 

pypromice/core/resampling.py

@@ -0,0 +1,142 @@
+import datetime
+import numpy as np
+import pandas as pd
+
+DEFAULT_COMPLETENESS_THRESHOLDS = {
+    "default": 0.8,
+    "albedo": 0.2,
+    "p_u": 0.5,
+    "p_l": 0.5,
+    "z_boom_u": 0.1,
+    "z_boom_l": 0.1,
+    "z_boom_cor_u": 0.1,
+    "z_boom_cor_l": 0.1,
+    "z_stake": 0.1,
+    "z_stake_cor": 0.1,
+    "z_surf_combined": 0.1,
+    "t_i_1": 0.1,
+    "t_i_2": 0.1,
+    "t_i_3": 0.1,
+    "t_i_4": 0.1,
+    "t_i_5": 0.1,
+    "t_i_6": 0.1,
+    "t_i_7": 0.1,
+    "t_i_8": 0.1,
+    "t_i_9": 0.1,
+    "t_i_10": 0.1,
+    "t_i_11": 0.1,
+    "gps_lat": 0.1,
+    "gps_lon": 0.1,
+    "gps_alt": 0.1,
+    "batt_v": 0.1,
+}
+
+ALLOWED_TIME_STAMP_DURATIONS = (
+    datetime.timedelta(minutes=10),
+    datetime.timedelta(minutes=30),
+    datetime.timedelta(hours=1),
+    datetime.timedelta(hours=6),
+    datetime.timedelta(days=1),
+)
+
+
+def classify_timestamp_durations(
+        index: pd.DatetimeIndex,
+) -> pd.TimedeltaIndex:
+    """
+    Classifies the durations between consecutive timestamps in a given DatetimeIndex.
+
+    The function computes the time differences between consecutive timestamps and
+    checks if these differences belong to a predefined set of allowed durations.
+    It performs backward filling to handle missing values
+
+    Parameters
+    ----------
+    index : pd.DatetimeIndex
+        A pandas DatetimeIndex containing the timestamps to classify.
+
+    Returns
+    -------
+    pd.TimedeltaIndex
+        A TimedeltaIndex containing the classified durations for the corresponding
+        timestamps in the input index.
+    """
+    return pd.TimedeltaIndex(
+        index.to_series()
+        .diff()
+        .where(lambda d: d.isin(ALLOWED_TIME_STAMP_DURATIONS))
+        .bfill()
+    )
+
+
+def get_completeness_mask(
+    data_frame: pd.DataFrame,
+    resample_offset: str,
+    completeness_thresholds: dict[str, float] = DEFAULT_COMPLETENESS_THRESHOLDS,
+    *,
+    atol: float = 1e-9,
+) -> pd.DataFrame:
+    """
+    Returns a completeness mask for the given DataFrame based on the specified
+    resampling offset, completeness threshold, and tolerance for over-completeness.
+
+    This function evaluates the completeness of timestamped data, ensuring that
+    records match the expected durations defined by the `resample_offset`. It
+    computes whether each resampled group of data satisfies the completeness
+    constraints defined by the `completeness_threshold` and `atol`.
+
+    Parameters
+    ----------
+    data_frame : pd.DataFrame
+        Input data containing a DatetimeIndex and associated values. The index must
+        be a DatetimeIndex as the function relies on timestamp durations for
+        computations.
+    resample_offset : str
+        Offset string defining resampling frequency. Examples include 'MS' (month
+        start) or other Pandas-compatible offset strings.
+    completeness_threshold : float, optional
+        Dictionary containing the variable-specific minimum completeness ratio
+        required to consider a time period as valid. Must contain a key 'default'
+        used for variables not explicitly listed.
+        Defaults to the dictionary `DEFAULT_COMPLETENESS_THRESHOLD`.
+    atol : float, optional
+        Absolute tolerance for over-completeness. Specifies an allowable margin by
+        which completeness can exceed 1. Defaults to 1e-9.
+
+    Returns
+    -------
+    pd.DataFrame
+        A DataFrame containing Boolean values, where True indicates that the data
+        for the corresponding time period satisfies the completeness constraints,
+        while False indicates the data is either under-complete or over-complete.
+    """
+    if resample_offset in ['MS', 'ME']:
+        offset_timedelta = datetime.timedelta(days=30)
+        # Increase tolerance for overcomplete values in monthly resampling
+        # to handle months with 31 days.
+        atol = 1/30 + atol
+    else:
+        offset_timedelta = pd.to_timedelta(resample_offset)
+
+    index = data_frame.index
+    assert isinstance(index, pd.DatetimeIndex)
+
+    timestamp_durations = classify_timestamp_durations(index)
+    timestamp_coverage = timestamp_durations / np.array(offset_timedelta)
+    data_frame_is_valid = data_frame.notna()
+
+    completeness = (
+        data_frame_is_valid
+        .mul(timestamp_coverage, axis=0)
+        .resample(resample_offset).sum()
+    )
+
+    thresholds = pd.Series(
+        {col: completeness_thresholds.get(col, completeness_thresholds["default"])
+         for col in data_frame.columns}
+    )
+
+    is_under_complete = completeness.lt(thresholds, axis=1)
+    is_over_complete = completeness.gt(1 + atol)
+    completeness_mask = ~(is_under_complete | is_over_complete)
+    return completeness_mask

pypromice/core/variables/station_boom_height.py

@@ -1,4 +1,4 @@
-__all__ = ["adjust", "adjust_and_include_uncorrected_values"]
+__all__ = ["adjust", "include_uncorrected_values"]
 
 import numpy as np
 import xarray as xr
@@ -25,25 +25,51 @@ def adjust(z_boom: xr.DataArray,
     return z_boom * ((air_temperature + T_0)/T_0)**0.5
 
 
-def adjust_and_include_uncorrected_values(z_boom: xr.DataArray,
-           air_temperature: xr.DataArray
+def include_uncorrected_values(
+    z_boom: xr.DataArray,
+    z_boom_cor: xr.DataArray,
+    air_temperature_other_level: xr.DataArray = None,
+    t_rad: xr.DataArray = None,
+    T_0: float = 273.15,
 ) -> xr.DataArray:
-    """Adjust sonic ranger readings for sensitivity to air temperature,
-    and retain uncorrected values where air temperature measurements
-    are not available.
+    """
+    Adjust sonic ranger readings for sensitivity to air temperature and
+    retain uncorrected values where temperature measurements are unavailable.
 
     Parameters
     ----------
     z_boom : xr.DataArray
-        Station boom height from sonic ranger
-    air_temperature : xr.DataArray
-        Air temperature
+        Uncorrected station boom height from sonic ranger
+    z_boom_cor : xr.DataArray
+        Boom height corrected with air_temperature_1
+    air_temperature_other_level : xr.DataArray, optional
+        Secondary air temperature
+    t_rad : xr.DataArray, optional
+        Radiative temperature
+    T_0 : float, optional
+        Reference temperature in Kelvin (default 273.15)
 
     Returns
     -------
     xr.DataArray
-        Adjusted station boom height
+        Corrected boom height with fallback where needed
     """
-    return xr.where(air_temperature.notnull(),
-                    z_boom * ((air_temperature + T_0)/T_0)**0.5,
-                    z_boom)
+    if air_temperature_other_level is None:
+        air_temperature_other_level = xr.full_like(z_boom, np.nan)
+    if t_rad is None:
+        t_rad = xr.full_like(z_boom, np.nan)
+    else:
+        t_rad = t_rad.clip(max=0)
+    
+    # first gap filling using values corrected with air_temperature_other_level
+    z_boom_cor_w_ta2 = z_boom * ((air_temperature_other_level + T_0) / T_0) ** 0.5
+    z_boom_cor=z_boom_cor.fillna(z_boom_cor_w_ta2)
+
+    # second gap filling using values corrected with t_rad
+    z_boom_cor_w_t_rad = z_boom * ((t_rad + T_0) / T_0) ** 0.5
+    z_boom_cor=z_boom_cor.fillna(z_boom_cor_w_t_rad)
+
+    # third gap filling using uncorrected values
+    z_boom_cor=z_boom_cor.fillna(z_boom)
+
+    return z_boom_cor

pypromice/pipeline/L2toL3.py

@@ -172,16 +172,28 @@ def process_surface_height(ds, data_adjustments_dir, station_config={}):
     ds['z_surf_1'] = ('time', ds['z_boom_u'].data * np.nan)
     ds['z_surf_2'] = ('time', ds['z_boom_u'].data * np.nan)
 
-    z_boom_best_u = station_boom_height.adjust_and_include_uncorrected_values(ds["z_boom_u"], ds["t_u"])
+    z_boom_best_u = station_boom_height.include_uncorrected_values(
+                                ds["z_boom_u"],
+                                ds["z_boom_cor_u"],
+                                ds["t_u"],
+                                ds["t_l"] if "t_l" in ds.data_vars else None,
+                                ds["t_rad"] if "t_rad" in ds.data_vars else None)
+
+
+
+    if 'z_stake' in ds.data_vars and ds.z_stake.notnull().any():
+        # Calculate stake boom height correction with uncorrected values where needed
+        z_stake_best = station_boom_height.include_uncorrected_values(
+                                    ds["z_stake"],
+                                    ds["z_stake_cor"],
+                                    ds["t_u"],
+                                    ds["t_l"] if "t_l" in ds.data_vars else None,
+                                    ds["t_rad"] if "t_rad" in ds.data_vars else None)
 
     if ds.attrs['site_type'] == 'ablation':
         # Calculate surface heights for ablation sites
         ds['z_surf_1'] = 2.6 - z_boom_best_u
         if ds.z_stake.notnull().any():
-
-            # Calculate stake boom height correction with uncorrected values where needed
-            z_stake_best = station_boom_height.adjust_and_include_uncorrected_values(ds["z_stake"], ds["t_u"])
-
             first_valid_index = ds.time.where((z_stake_best + z_boom_best_u).notnull(), drop=True).data[0]
             ds['z_surf_2'] = ds.z_surf_1.sel(time=first_valid_index) + z_stake_best.sel(time=first_valid_index) - z_stake_best
 
@@ -195,14 +207,18 @@ def process_surface_height(ds, data_adjustments_dir, station_config={}):
         ds['z_surf_1'] = z_boom_best_u.sel(time=first_valid_index) - z_boom_best_u
 
         if 'z_stake' in ds.data_vars and ds.z_stake.notnull().any():
-            z_stake_best = station_boom_height.adjust_and_include_uncorrected_values(ds["z_stake"], ds["t_u"])
             first_valid_index = ds.time.where(z_stake_best.notnull(), drop=True).data[0]
             ds['z_surf_2'] = z_stake_best.sel(time=first_valid_index) - z_stake_best
 
         if 'z_boom_l' in ds.data_vars:
 
             # Calculate lower boom height correction with uncorrected values where needed
-            z_boom_best_l = station_boom_height.adjust_and_include_uncorrected_values(ds["z_boom_l"], ds["t_l"])
+            z_boom_best_l = station_boom_height.include_uncorrected_values(
+                                        ds["z_boom_l"],
+                                        ds["z_boom_cor_l"],
+                                        ds["t_l"],
+                                        ds["t_u"] if "t_u" in ds.data_vars else None,
+                                        ds["t_rad"] if "t_rad" in ds.data_vars else None)
 
             # need a combine first because KAN_U switches from having a z_stake_best
             # to having a z_boom_best_l

pypromice/pipeline/get_l2tol3.py

@@ -92,7 +92,7 @@ def get_l2tol3(config_folder: Path|str, inpath, outpath, variables, metadata, da
     if outpath is not None:
         prepare_and_write(l3, outpath, v, m, '60min')
         prepare_and_write(l3, outpath, v, m, '1D')
-        prepare_and_write(l3, outpath, v, m, 'ME')
+        prepare_and_write(l3, outpath, v, m, 'MS')
     return l3
 
 def main():

pypromice/pipeline/resample.py

@@ -1,18 +1,15 @@
 #!/usr/bin/env python3
 # -*- coding: utf-8 -*-
-"""
-Created on Mon Jun 10 10:58:39 2024
-
-@author: pho
-"""
 import logging
 import numpy as np
 import pandas as pd
 import xarray as xr
+
+from pypromice.core.resampling import get_completeness_mask, DEFAULT_COMPLETENESS_THRESHOLDS
 from pypromice.core.variables.wind import calculate_directional_wind_speed
 logger = logging.getLogger(__name__)
 
-def resample_dataset(ds_h, t):
+def resample_dataset(ds_h, t, completeness_thresholds=DEFAULT_COMPLETENESS_THRESHOLDS):
     '''Resample L2 AWS data, e.g. hourly to daily average. This uses pandas
     DataFrame resampling at the moment as a work-around to the xarray Dataset
     resampling. As stated, xarray resampling is a lengthy process that takes
@@ -27,6 +24,11 @@ def resample_dataset(ds_h, t):
     t : str
         Resample factor( "60min", "1D" or "MS"), same variable definition as in
         pandas.DataFrame.resample()
+    completeness_thresholds : Dict
+        A dict with, for each variable, the lower limit of completness of an
+        hourly/daily/monthly aggregate (nr of samples in aggregate / expected
+        nr of samples). Aggregates below that limit are replaced by NaNs.
+        Must include a "default" value used for variables not listed explicitly.
 
     Returns
     -------
@@ -62,6 +64,15 @@ def resample_dataset(ds_h, t):
         if var in df_h.columns:
             df_resampled[var] = df_h[var].resample(t).sum()
 
+    # Apply completeness filter based on the the data frame time index
+    completeness_mask = get_completeness_mask(
+        data_frame=df_h,
+        resample_offset=t,
+        completeness_thresholds=completeness_thresholds,
+    )
+
+    df_resampled[~completeness_mask] = np.nan
+
     # taking the 10 min data and using it as instantaneous values:
     is_10_minutes_timestamp = (ds_h.time.diff(dim='time') / np.timedelta64(1, 's') == 600)
     if (t == '60min') and is_10_minutes_timestamp.any():

pypromice/resources/variable_aliases_GC-Net.csv

@@ -27,8 +27,8 @@ dlhf_u,LHF
 dlhf_l,
 dshf_u,SHF
 dshf_l,
-z_boom_u,HW2
-z_boom_l,HW1
+z_boom_cor_u,HW2
+z_boom_cor_l,HW1
 precip_u,
 precip_u_cor,
 precip_l,

pypromice/resources/variables.csv

@@ -34,13 +34,13 @@ dlhf_l,surface_downward_latent_heat_flux,Latent heat flux (lower boom),W m-2,mod
 dshf_u,surface_downward_sensible_heat_flux,Sensible heat flux (upper boom),W m-2,modelResult,time,FALSE,L3 or later,,,,all,0,0,1,4
 dshf_l,surface_downward_sensible_heat_flux,Sensible heat flux (lower boom),W m-2,modelResult,time,FALSE,L3 or later,,,,two-boom,0,0,1,4
 z_boom_u,distance_to_surface_from_boom,Upper boom height,m,physicalMeasurement,time,TRUE,,0.3,10,z_boom_cor_u,all,1,1,1,4
-z_boom_cor_u,distance_to_surface_from_boom_corrected,Upper boom height - corrected,m,modelResult,time,TRUE,,0.3,10,"",all,1,1,1,4
+z_boom_cor_u,distance_to_surface_from_boom_corrected,Upper boom height - corrected,m,modelResult,time,TRUE,,0.3,10,z_boom_u,all,1,1,1,4
 z_boom_q_u,distance_to_surface_from_boom_quality,Upper boom height (quality),-,qualityInformation,time,TRUE,L0 or L2,,,,all,1,1,0,4
 z_boom_l,distance_to_surface_from_boom,Lower boom height,m,physicalMeasurement,time,TRUE,,0.3,5,z_boom_cor_l,two-boom,1,1,1,4
-z_boom_cor_l,distance_to_surface_from_boom_corrected,Lower boom height - corrected,m,modelResult,time,TRUE,,0.3,5,"",two-boom,1,1,1,4
+z_boom_cor_l,distance_to_surface_from_boom_corrected,Lower boom height - corrected,m,modelResult,time,TRUE,,0.3,5,z_boom_l,two-boom,1,1,1,4
 z_boom_q_l,distance_to_surface_from_boom_quality,Lower boom height (quality),-,qualityInformation,time,TRUE,L0 or L2,,,,two-boom,1,1,0,4
 z_stake,distance_to_surface_from_stake_assembly,Stake height,m,physicalMeasurement,time,TRUE,,0.3,8,z_stake_cor,one-boom,1,1,1,4
-z_stake_cor,distance_to_surface_from_stake_assembly_corrected,Stake height - corrected,m,physicalMeasurement,time,TRUE,,0.3,8,,one-boom,1,1,1,4
+z_stake_cor,distance_to_surface_from_stake_assembly_corrected,Stake height - corrected,m,physicalMeasurement,time,TRUE,,0.3,8,z_stake,one-boom,1,1,1,4
 z_stake_q,distance_to_surface_from_stake_assembly_quality,Stake height (quality),-,qualityInformation,time,TRUE,L0 or L2,,,,one-boom,1,1,0,4
 z_pt,depth_of_pressure_transducer_in_ice,Depth of pressure transducer in ice,m,physicalMeasurement,time,FALSE,,0,30,z_pt_cor,one-boom,1,1,1,4
 z_pt_cor,depth_of_pressure_transducer_in_ice_corrected,Depth of pressure transducer in ice - corrected,m,modelResult,time,FALSE,L2 or later,0,30,,one-boom,0,1,1,4

{pypromice-1.6.0.dist-info → pypromice-1.7.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pypromice
-Version: 1.6.0
+Version: 1.7.0
 Summary: PROMICE/GC-Net data processing toolbox
 Home-page: https://github.com/GEUS-Glaciology-and-Climate/pypromice
 Author: GEUS Glaciology and Climate

{pypromice-1.6.0.dist-info → pypromice-1.7.0.dist-info}/RECORD

@@ -1,8 +1,9 @@
 pypromice/__init__.py,sha256=X2LaniNJv4iVYqRtmf2jyGMQLaj59bIgkhlWhT6LCgQ,74
 pypromice/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 pypromice/core/dependency_graph.py,sha256=bqoXasC8pg5ipjBd6rqDhfHwIq11t2_cFlNT72ncw4w,3135
+pypromice/core/resampling.py,sha256=KyYk6HWFqJo60CPiX-gBd0uwLx5iPrYJBMQ3Sqb9yFg,4645
 pypromice/core/qc/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-pypromice/core/qc/github_data_issues.py,sha256=gbbF62oMMWbXiLnsrs60vXbwfAqSUP113plhidVeqCU,13353
+pypromice/core/qc/github_data_issues.py,sha256=mTnJArv89p8vSzJdWFRkA4kBAYnFcR-I-Xz4QzlrifA,13857
 pypromice/core/qc/persistence.py,sha256=Y9CmAAPHNqEjdsZW4LEAVm7cnMEBK-zwEAD_UBDASVw,6466
 pypromice/core/qc/value_clipping.py,sha256=KGLLN54-QeD4TQ-Dd-NeooYG3kdEC6SPhz4LT2ZsRi4,1533
 pypromice/core/qc/percentiles/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -16,7 +17,7 @@ pypromice/core/variables/humidity.py,sha256=TfQfGbQ0gyj_A9zRChGei47oZYMRAR13eIzj
 pypromice/core/variables/precipitation.py,sha256=iH7m1sJDwbstZuAutO2HCRVo3RqaSuDG33YERBtvQag,3467
 pypromice/core/variables/pressure_transducer_depth.py,sha256=SadI_8oRKmMI2dnsYfYDJD7AnE6-i1bv0MOmO5xms8E,2418
 pypromice/core/variables/radiation.py,sha256=11iGl2cq5WGJxnPwTAy3H4Y5bU-QXTaMaEife1BPNA0,14318
-pypromice/core/variables/station_boom_height.py,sha256=5bomhyXLSlBhcL9EaMInkcriNeatuvnUqu98hX-0n0s,1317
+pypromice/core/variables/station_boom_height.py,sha256=5uGalC-uVEdJFMPQ-cuayP5kpFu8mvou10w38X7kd8Y,2281
 pypromice/core/variables/station_pose.py,sha256=3aVb8hywsFBCDorgHPW_nNdf_pyyil6ziJsP3UBhje4,13112
 pypromice/core/variables/wind.py,sha256=-dpaBOtqGyPDIU4O1HbbWRzlLNRC2a50OdnZhIaWdeI,1701
 pypromice/io/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -37,28 +38,28 @@ pypromice/io/ingest/l0_repository.py,sha256=jmlQMSETj9aBA1iOC9sPl4-ia9UfsRHTWiR0
 pypromice/io/ingest/toa5.py,sha256=GiDCUMJpvbTRs0PzU1wz6uPkYkih9XjDVPcGbl__RtU,2857
 pypromice/pipeline/L0toL1.py,sha256=g_pH9sVPJMfYKdFgJpT7WW5rqJntJPKsxpw9N_nKs-o,13816
 pypromice/pipeline/L1toL2.py,sha256=lNSC3QOQf-ZZpFlCW6PFeRIkWtGWlDMBXWzNBU2yKAc,10158
-pypromice/pipeline/L2toL3.py,sha256=zDJggJMojKtK1VPkWoQBEgUNBbm02b2-xR-ajW2E3tg,60846
+pypromice/pipeline/L2toL3.py,sha256=yEvjIaheTdT136HEMJXmqAuReYgnx0pF6VhG23VFj00,61650
 pypromice/pipeline/__init__.py,sha256=yqYenngaSNFtpnAsfc953142P84ocq3ykvlsCIbsw3g,151
 pypromice/pipeline/aws.py,sha256=5Xv7XVf5PvcSAnEu_jPpjDWV4U0p_UvdpOzCtmH0ICU,5092
 pypromice/pipeline/get_l2.py,sha256=VV4iV3P34HFSOi3jr0IQrNXjBMgCx5GpF0uwNbT84Ck,3072
-pypromice/pipeline/get_l2tol3.py,sha256=RLyyV2iCqCXInc91sMgxltTmAD8hAt6NU3-meIR9kac,4497
+pypromice/pipeline/get_l2tol3.py,sha256=KIsScgatZbccz7ypZHBGiibbresJ9HWq9Sv4Ys2LcHs,4497
 pypromice/pipeline/join_l2.py,sha256=AQL9wVqRFbdB6bevuOArg9knZCv-jDVwypMETjuDqBM,4561
 pypromice/pipeline/join_l3.py,sha256=Wd5slBCo1XgzExOjPu9IRLWnk7NJMTwDcckhithvw0E,20316
-pypromice/pipeline/resample.py,sha256=SMq9Rh-Uy6DbtQdDv062t3n-AChRHx_J3AR914sMkLc,8099
+pypromice/pipeline/resample.py,sha256=tq6GvnmbQJsFvdJxZF2DrwTTmeNRVig_LHZQboKWreU,8818
 pypromice/pipeline/utilities.py,sha256=m-BaHWMKBGzTHa09w-49yqWRulXex5TTScg7IZu8fSY,1248
 pypromice/resources/__init__.py,sha256=MpKmvV11R4tcqmyvJpXQt-_It3oRI0WEIQNbMST--4w,907
 pypromice/resources/file_attributes.csv,sha256=ISKR-Ax12CT9tQD38ByNyvWPLLpMszpWXwviPyhuUaI,7018
-pypromice/resources/variable_aliases_GC-Net.csv,sha256=ts10e0AI1Fz_TtEj8qdRFIZBTkE8pt_VtJGp5hSBUBI,993
-pypromice/resources/variables.csv,sha256=gUmA37r_k9JEP5Tvog3gUC6lROwz-w9DfX0FW3RA-FQ,14166
+pypromice/resources/variable_aliases_GC-Net.csv,sha256=mIAnBy2XvnGGy4_HODHlbZCvk0jq3v9gAncMwyDmyqI,1001
+pypromice/resources/variables.csv,sha256=hcd0WcY9vaG6MhqnCUNjbFxXNg7p5d_ZIyaebLKy0s8,14185
 pypromice/tx/__init__.py,sha256=-62bhHWJGfzFh5JwHcLqRj2jcGzmqzYOLWByhO706YY,30
 pypromice/tx/get_l0tx.py,sha256=b34-96KGshTyTN2tBFaAIBl7oZZzbRB_JR7sXtDNfXA,6957
 pypromice/tx/get_msg.py,sha256=OGS60OHjy4Wf8JExTfOdK-9xhjFdjhuChxoTSPe_MjI,3417
 pypromice/tx/payload_formats.csv,sha256=A46-XcYdpe9-gzmADylP2UVizLi_UphF-BPT5u3Lyn8,7903
 pypromice/tx/payload_types.csv,sha256=C1-xCmHytAqqAzgzPwBLWqabzWu6s6tKAd8AjVd935s,457
 pypromice/tx/tx.py,sha256=asbgXVI5vurKM-WVACTfpKRt-70wtzVvSbvjvYufajI,34416
-pypromice-1.6.0.dist-info/licenses/LICENSE.txt,sha256=gXf5dRMhNSbfLPYYTY_5hsZ1r7UU1OaKQEAQUhuIBkM,18092
-pypromice-1.6.0.dist-info/METADATA,sha256=9Y4jNegT_VKzN0dQdqFSvOWcLSONdsva-tNfvG5qiB0,4958
-pypromice-1.6.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-pypromice-1.6.0.dist-info/entry_points.txt,sha256=aU2SG5A0wI2n4HE9fotG-y5yeIqJ0G1-8UrEgr1GUQk,535
-pypromice-1.6.0.dist-info/top_level.txt,sha256=cBdfwgSbWDQq3a07nKRjrfmLC7jdaYXs98GG58HpTks,10
-pypromice-1.6.0.dist-info/RECORD,,
+pypromice-1.7.0.dist-info/licenses/LICENSE.txt,sha256=gXf5dRMhNSbfLPYYTY_5hsZ1r7UU1OaKQEAQUhuIBkM,18092
+pypromice-1.7.0.dist-info/METADATA,sha256=sU7zLwdi_aCLLUUUagNsQgIFVEwae5acab3osEHJDGM,4958
+pypromice-1.7.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+pypromice-1.7.0.dist-info/entry_points.txt,sha256=aU2SG5A0wI2n4HE9fotG-y5yeIqJ0G1-8UrEgr1GUQk,535
+pypromice-1.7.0.dist-info/top_level.txt,sha256=cBdfwgSbWDQq3a07nKRjrfmLC7jdaYXs98GG58HpTks,10
+pypromice-1.7.0.dist-info/RECORD,,