seagliderOG1 demo

The purpose of this notebook is to demonstrate the functionality of seagliderOG1 to convert from Seaglider basestation files to OG1 format.

OG1 format is a newly agreed format (since June 2024) for glider data sets from various platforms (e.g., Seaglider, Slocum, Seaexplorer). It lives on github here: (https://github.com/OceanGlidersCommunity/OG-format-user-manual).
OG1 manual: https://oceangliderscommunity.github.io/OG-format-user-manual/OG_Format.html

The test case is to convert sg015 data from the Labrador Sea in September 2004.

The demo is organised to show

Conversion of a single dive cycle (single p*.nc file)
Conversion for a folder of local dive-cycle files (full mission of p*.nc files)
Download from remote server + conversion (directory with full mission of p*.nc files)

Options are provided to only load e.g. 10 files, but note that OG1 format expects a full mission.

[1]:

import pathlib
import sys

script_dir = pathlib.Path().parent.absolute()
parent_dir = script_dir.parents[0]
sys.path.append(str(parent_dir))
sys.path.append(str(parent_dir) + '/seagliderOG1')
print(parent_dir)
print(sys.path)

import xarray as xr
import os
import pooch
from seagliderOG1 import readers, writers, plotters, tools
from seagliderOG1 import convertOG1, vocabularies

/home/runner/work/seagliderOG1/seagliderOG1
['/home/runner/micromamba/envs/TEST/lib/python312.zip', '/home/runner/micromamba/envs/TEST/lib/python3.12', '/home/runner/micromamba/envs/TEST/lib/python3.12/lib-dynload', '', '/home/runner/micromamba/envs/TEST/lib/python3.12/site-packages', '/home/runner/work/seagliderOG1/seagliderOG1', '/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1']

/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1/tools.py:549: SyntaxWarning: invalid escape sequence '\d'
  fill_val = 2 ** (int(re.findall("\d+", str(new_dtype))[0]) - 1) - 1

[ ]:

[2]:

# Specify the path for writing datafiles
data_path = os.path.join(parent_dir, 'data')

Reading basestation files

This has three ways to load a glider dataset.

Load an example dataset using seagliderOG1.fetchers.load_sample_dataset

Alternatively, use your own with e.g. ds = xr.open_dataset('/path/to/yourfile.nc')

Load datasets from a remote directory (URL)

[6]:

# Specify the server where data are located
server = "https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/015/20040924/"

# Load and concatenate all datasets from the server, optionally specifying the range of profiles to load
list_datasets = readers.read_basestation(server, start_profile=500, end_profile=503)

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
Cell In[6], line 5
      2 server = "https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/015/20040924/"
      4 # Load and concatenate all datasets from the server, optionally specifying the range of profiles to load
----> 5 list_datasets = readers.read_basestation(server, start_profile=500, end_profile=503)

AttributeError: module 'seagliderOG1.readers' has no attribute 'read_basestation'

Convert to OG1 format

Process:

For one basestation dataset, split the dataset by dimension (split_ds)
Transform into OG1 format: dataset with dims sg_data_point
- Change the dimension to N_MEASUREMENTS
- Rename variables according to vocabularies.standard_names
- Assign variable attributes according to vocabularies.vocab_attrs. (Note: This could go wrong since it makes assumptions about the input variables. May need additional handling.)
Add missing mandatory variables:
- From split_ds[(gps_info,)], add the LATITUDE_GPS, LONGITUDE_GPS and TIME_GPS (Note: presently TIME_GPS is stripped before saving, but TIME values contain TIME_GPS)
- Create PROFILE_NUMBER and PHASE
- Calculate DEPTH_Z which is positive up
Update attributes for the file.
- Combines creator and contributor from original attributes into contributor
- Adds contributing_institutions based on institution
- Reformats time in time_coverage_* and start_time–> start_date
- Adds date_modified
- Renames comments–>history, site–>summary
- Adds title, platform, platform_vocabulary, featureType, Conventions, rtqc_method* according to OceanGliders format
- Retains naming_authority, institution, project, geospatial_* as OG attributes
- Retains extra attributes: license, keywords, keywords_vocabulary, file_version, acknowledgement, date_created, disclaimer

Future behaviour to be added:

Retain the variables starting with sg_cal and check whether they vary over the mission (shouldn’t)
Add sensors, using information in the split_ds with no dimensions
- Need (from sg_cal_constants: sg_cal plus volmax, vbd_cnts_per_cc, therm_expan, t_*, mass, hd_*, ctcor, cpcor, c_*, abs_compress, a, Tcor, Soc, Pcor, Foffset)
- Maybe also reviewed, magnetic_variation (which will change with position), log_D_FLARE, flight_avg_speed_north and flight_avg_speed_east also with _gsm, depth_avg_curr_north and depth_avg_curr_east also with _gsm, wlbb2f - means sensor sg_cal_mission_title sg_cal_id_str calibcomm_oxygen calibcomm sbe41 means ?? hdm_qc glider

Convert a single (sample) dataset

[7]:

# Loads one dataset (p0150500_20050213.nc)
ds = readers.load_sample_dataset()

ds_OG1 = convertOG1.convert_to_OG1(ds)

# Check the results - uncomment the following lines to either generate a plot or show the variables.
plotters.plot_profile_depth(ds_OG1)

/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:58: FutureWarning: In a future version of xarray decode_timedelta will default to False rather than None. To silence this warning, set decode_timedelta to True, False, or a 'CFTimedeltaCoder' instance.
  return xr.open_dataset(file_path)
/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:58: FutureWarning: In a future version of xarray decode_timedelta will default to False rather than None. To silence this warning, set decode_timedelta to True, False, or a 'CFTimedeltaCoder' instance.
  return xr.open_dataset(file_path)
/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:58: FutureWarning: In a future version of xarray decode_timedelta will default to False rather than None. To silence this warning, set decode_timedelta to True, False, or a 'CFTimedeltaCoder' instance.
  return xr.open_dataset(file_path)
/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:58: FutureWarning: In a future version of xarray decode_timedelta will default to False rather than None. To silence this warning, set decode_timedelta to True, False, or a 'CFTimedeltaCoder' instance.
  return xr.open_dataset(file_path)
/home/runner/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:58: FutureWarning: In a future version of xarray decode_timedelta will default to False rather than None. To silence this warning, set decode_timedelta to True, False, or a 'CFTimedeltaCoder' instance.
  return xr.open_dataset(file_path)

---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
Cell In[7], line 4
      1 # Loads one dataset (p0150500_20050213.nc)
      2 ds = readers.load_sample_dataset()
----> 4 ds_OG1 = convertOG1.convert_to_OG1(ds)
      6 # Check the results - uncomment the following lines to either generate a plot or show the variables.
      7 plotters.plot_profile_depth(ds_OG1)

File ~/work/seagliderOG1/seagliderOG1/seagliderOG1/convertOG1.py:37, in convert_to_OG1(list_of_datasets, contrib_to_append)
     35 for ds1_base in list_of_datasets:
     36     varlist = list(set(varlist + list(ds1_base.variables)))
---> 37     ds_new, attr_warnings, sg_cal, dc_other, dc_log = process_dataset(ds1_base, firstrun)
     38     if ds_new:
     39         processed_datasets.append(ds_new)

File ~/work/seagliderOG1/seagliderOG1/seagliderOG1/convertOG1.py:174, in process_dataset(ds1_base, firstrun)
    172 divenum = ds1_base.attrs['dive_number']
    173 # Split the dataset by unique dimensions
--> 174 split_ds = tools.split_by_unique_dims(ds1)
    175 ds_sgdatapoint = split_ds[('sg_data_point',)]
    176 # Extract the gps_info from the split dataset

NameError: name 'ds1' is not defined

[8]:

# Print to screen a table of the variables and variable attributes
#plotters.show_attributes(ds_OG1)
ds_OG1

---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
Cell In[8], line 3
      1 # Print to screen a table of the variables and variable attributes
      2 #plotters.show_attributes(ds_OG1)
----> 3 ds_OG1

NameError: name 'ds_OG1' is not defined

Convert mission from a local directory of basestation files

For local data in the directory input_dir
Creates a plot of ctd_depth against ctd_time.

[9]:

# Specify the input directory on your local machine
input_dir = "/Users/eddifying/Dropbox/data/sg015-ncei-download"

# Load and concatenate all datasets in the input directory
# Optionally, specify the range of profiles to load (start_profile, end_profile)
list_datasets = readers.read_basestation(input_dir, start_profile=500, end_profile=503)

# Convert the list of datasets to OG1
ds_OG1 = convertOG1.convert_to_OG1(list_datasets)

# Generate a simple plot
plotters.plot_profile_depth(ds_OG1)
#plotters.show_contents(ds_all,'attrs')

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
Cell In[9], line 6
      2 input_dir = "/Users/eddifying/Dropbox/data/sg015-ncei-download"
      4 # Load and concatenate all datasets in the input directory
      5 # Optionally, specify the range of profiles to load (start_profile, end_profile)
----> 6 list_datasets = readers.read_basestation(input_dir, start_profile=500, end_profile=503)
      8 # Convert the list of datasets to OG1
      9 ds_OG1 = convertOG1.convert_to_OG1(list_datasets)

AttributeError: module 'seagliderOG1.readers' has no attribute 'read_basestation'

Convert mission from the NCEI server (with p*nc files)

Data from the sg015 mission in the Labrador Sea (https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.nodc:0111844), dataset identifier gov.noaa.nodc:0111844.

[10]:

# Specify the server where data are located
server = "https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/015/20040924/"

# Load and concatenate all datasets from the server, optionally specifying the range of profiles to load
list_datasets = readers.read_basestation(server, start_profile=500, end_profile=503)

# Convert the list of datasets to OG1
ds_OG1 = convertOG1.convert_to_OG1(list_datasets)

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
Cell In[10], line 5
      2 server = "https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/015/20040924/"
      4 # Load and concatenate all datasets from the server, optionally specifying the range of profiles to load
----> 5 list_datasets = readers.read_basestation(server, start_profile=500, end_profile=503)
      7 # Convert the list of datasets to OG1
      8 ds_OG1 = convertOG1.convert_to_OG1(list_datasets)

AttributeError: module 'seagliderOG1.readers' has no attribute 'read_basestation'

Saving data

Due to problems with writing xarray datasets as netCDF when attributes are not of a specified type (str, Number, np.ndarray, np.number, list, tuple), a function was written save_dataset.

[11]:

# Write the file
# This writer catches errors in data types (DType errors) when using xr.to_netcdf()
# The solution is to convert them to strings, which may be undesired behaviour
output_file = os.path.join(data_path, 'demo_test.nc')
if os.path.exists(output_file):
    os.remove(output_file)

writers.save_dataset(ds_OG1, output_file);

---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
Cell In[11], line 8
      5 if os.path.exists(output_file):
      6     os.remove(output_file)
----> 8 writers.save_dataset(ds_OG1, output_file);

NameError: name 'ds_OG1' is not defined

[12]:

# Load the data saved
ds1 = xr.open_dataset(output_file)

# Generate a simple plot
#plotters.show_contents(ds_all,'attrs')
plotters.plot_depth_colored(ds1, color_by='PROFILE_NUMBER')

---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/file_manager.py:211, in CachingFileManager._acquire_with_cache_info(self, needs_lock)
    210 try:
--> 211     file = self._cache[self._key]
    212 except KeyError:

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/lru_cache.py:56, in LRUCache.__getitem__(self, key)
     55 with self._lock:
---> 56     value = self._cache[key]
     57     self._cache.move_to_end(key)

KeyError: [<class 'netCDF4._netCDF4.Dataset'>, ('/home/runner/work/seagliderOG1/seagliderOG1/data/demo_test.nc',), 'r', (('clobber', True), ('diskless', False), ('format', 'NETCDF4'), ('persist', False)), '2174b3b0-d56a-4819-a507-d7b74ddfa984']

During handling of the above exception, another exception occurred:

FileNotFoundError                         Traceback (most recent call last)
Cell In[12], line 2
      1 # Load the data saved
----> 2 ds1 = xr.open_dataset(output_file)
      4 # Generate a simple plot
      5 #plotters.show_contents(ds_all,'attrs')
      6 plotters.plot_depth_colored(ds1, color_by='PROFILE_NUMBER')

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/api.py:686, in open_dataset(filename_or_obj, engine, chunks, cache, decode_cf, mask_and_scale, decode_times, decode_timedelta, use_cftime, concat_characters, decode_coords, drop_variables, inline_array, chunked_array_type, from_array_kwargs, backend_kwargs, **kwargs)
    674 decoders = _resolve_decoders_kwargs(
    675     decode_cf,
    676     open_backend_dataset_parameters=backend.open_dataset_parameters,
   (...)
    682     decode_coords=decode_coords,
    683 )
    685 overwrite_encoded_chunks = kwargs.pop("overwrite_encoded_chunks", None)
--> 686 backend_ds = backend.open_dataset(
    687     filename_or_obj,
    688     drop_variables=drop_variables,
    689     **decoders,
    690     **kwargs,
    691 )
    692 ds = _dataset_from_backend_dataset(
    693     backend_ds,
    694     filename_or_obj,
   (...)
    704     **kwargs,
    705 )
    706 return ds

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/netCDF4_.py:666, in NetCDF4BackendEntrypoint.open_dataset(self, filename_or_obj, mask_and_scale, decode_times, concat_characters, decode_coords, drop_variables, use_cftime, decode_timedelta, group, mode, format, clobber, diskless, persist, auto_complex, lock, autoclose)
    644 def open_dataset(
    645     self,
    646     filename_or_obj: str | os.PathLike[Any] | ReadBuffer | AbstractDataStore,
   (...)
    663     autoclose=False,
    664 ) -> Dataset:
    665     filename_or_obj = _normalize_path(filename_or_obj)
--> 666     store = NetCDF4DataStore.open(
    667         filename_or_obj,
    668         mode=mode,
    669         format=format,
    670         group=group,
    671         clobber=clobber,
    672         diskless=diskless,
    673         persist=persist,
    674         auto_complex=auto_complex,
    675         lock=lock,
    676         autoclose=autoclose,
    677     )
    679     store_entrypoint = StoreBackendEntrypoint()
    680     with close_on_error(store):

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/netCDF4_.py:452, in NetCDF4DataStore.open(cls, filename, mode, format, group, clobber, diskless, persist, auto_complex, lock, lock_maker, autoclose)
    448     kwargs["auto_complex"] = auto_complex
    449 manager = CachingFileManager(
    450     netCDF4.Dataset, filename, mode=mode, kwargs=kwargs
    451 )
--> 452 return cls(manager, group=group, mode=mode, lock=lock, autoclose=autoclose)

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/netCDF4_.py:393, in NetCDF4DataStore.__init__(self, manager, group, mode, lock, autoclose)
    391 self._group = group
    392 self._mode = mode
--> 393 self.format = self.ds.data_model
    394 self._filename = self.ds.filepath()
    395 self.is_remote = is_remote_uri(self._filename)

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/netCDF4_.py:461, in NetCDF4DataStore.ds(self)
    459 @property
    460 def ds(self):
--> 461     return self._acquire()

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/netCDF4_.py:455, in NetCDF4DataStore._acquire(self, needs_lock)
    454 def _acquire(self, needs_lock=True):
--> 455     with self._manager.acquire_context(needs_lock) as root:
    456         ds = _nc4_require_group(root, self._group, self._mode)
    457     return ds

File ~/micromamba/envs/TEST/lib/python3.12/contextlib.py:137, in _GeneratorContextManager.__enter__(self)
    135 del self.args, self.kwds, self.func
    136 try:
--> 137     return next(self.gen)
    138 except StopIteration:
    139     raise RuntimeError("generator didn't yield") from None

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/file_manager.py:199, in CachingFileManager.acquire_context(self, needs_lock)
    196 @contextlib.contextmanager
    197 def acquire_context(self, needs_lock=True):
    198     """Context manager for acquiring a file."""
--> 199     file, cached = self._acquire_with_cache_info(needs_lock)
    200     try:
    201         yield file

File ~/micromamba/envs/TEST/lib/python3.12/site-packages/xarray/backends/file_manager.py:217, in CachingFileManager._acquire_with_cache_info(self, needs_lock)
    215     kwargs = kwargs.copy()
    216     kwargs["mode"] = self._mode
--> 217 file = self._opener(*self._args, **kwargs)
    218 if self._mode == "w":
    219     # ensure file doesn't get overridden when opened again
    220     self._mode = "a"

File src/netCDF4/_netCDF4.pyx:2521, in netCDF4._netCDF4.Dataset.__init__()

File src/netCDF4/_netCDF4.pyx:2158, in netCDF4._netCDF4._ensure_nc_success()

FileNotFoundError: [Errno 2] No such file or directory: '/home/runner/work/seagliderOG1/seagliderOG1/data/demo_test.nc'

Run multiple missions

[13]:

# Add these to existing attributes - update to your details
contrib_to_append = vocabularies.contrib_to_append
print(contrib_to_append)

{'contributor_name': 'Eleanor Frajka-Williams', 'contributor_email': 'eleanorfrajka@gmail.com', 'contributor_role': 'Data scientist', 'contributor_role_vocabulary': 'http://vocab.nerc.ac.uk/search_nvs/W08', 'contributing_institutions': 'University of Hamburg - Institute of Oceanography', 'contributing_institutions_vocabulary': 'https://edmo.seadatanet.org/report/1156', 'contributing_institutions_role': 'Data scientist', 'contributing_institutions_role_vocabulary': 'http://vocab.nerc.ac.uk/search_nvs/W08'}

[14]:

# Specify a list of servers or local directories
input_locations = [
    # Either Iceland, Faroes or RAPID/MOCHA
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/005/20090829/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/005/20080606/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/005/20081106/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/012/20070831/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/014/20080214/",  # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/014/20080222/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/016/20061112/",  # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/016/20090605/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/016/20071113/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/016/20080607/",  # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/033/20100518/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/033/20100903/", # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/101/20081108/",     # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/101/20061112/",    # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/101/20070609/",   # done
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/102/20061112/",  # done
    # Labrador Sea
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/015/20040924/",
    "https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/014/20040924/",
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/008/20031002/",
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/004/20031002/",
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/016/20050406/",
    # RAPID/MOCHA
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/033/20100729/",
    #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/034/20110128/",
]

for input_loc in input_locations:
    # Example usage
    ds_all = convertOG1.process_and_save_data(input_loc, output_dir=data_path, save=True,  run_quietly=True)

---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
Cell In[14], line 33
      2 input_locations = [
      3     # Either Iceland, Faroes or RAPID/MOCHA
      4     #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/005/20090829/", # done
   (...)
     28     #"https://www.ncei.noaa.gov/data/oceans/glider/seaglider/uw/034/20110128/",
     29 ]
     31 for input_loc in input_locations:
     32     # Example usage
---> 33     ds_all = convertOG1.process_and_save_data(input_loc, output_dir=data_path, save=True,  run_quietly=True)

File ~/work/seagliderOG1/seagliderOG1/seagliderOG1/convertOG1.py:600, in process_and_save_data(input_location, save, output_dir, run_quietly)
    584 """
    585 Processes and saves data from the specified input location.
    586 This function loads and concatenates datasets from the server, converts them to OG1 format,
   (...)
    596 xarray.Dataset: The processed dataset.
    597 """
    599 # Load and concatenate all datasets from the server
--> 600 ds1_base = readers.load_first_basestation_file(input_location)
    602 # Convert the list of datasets to OG1
    603 ds1_og1, varlist = convert_to_OG1(ds1_base)

File ~/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:161, in load_first_basestation_file(source)
    159 filename = file_list[0]
    160 start_profile = _profnum_from_filename(filename)
--> 161 datasets = load_basestation_files(source, start_profile, start_profile)
    162 return datasets[0]

File ~/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:183, in load_basestation_files(source, start_profile, end_profile)
    181 for file in filtered_files:
    182     if source.startswith("http://") or source.startswith("https://"):
--> 183         ds = load_sample_dataset(file)
    184     else:
    185         ds = xr.open_dataset(os.path.join(source, file))

File ~/work/seagliderOG1/seagliderOG1/seagliderOG1/readers.py:61, in load_sample_dataset(dataset_name)
     59 else:
     60     msg = f"Requested sample dataset {dataset_name} not known. Specify one of the following available datasets: {list(data_source_og.registry.keys())}"
---> 61     raise KeyError(msg)

KeyError: "Requested sample dataset p0140001_20040924.nc not known. Specify one of the following available datasets: ['p0330001_20100903.nc', 'p0330002_20100903.nc', 'p0330003_20100903.nc', 'p0330004_20100904.nc', 'p0330005_20100904.nc', 'p0330006_20100904.nc', 'p0330007_20100905.nc', 'p0330008_20100905.nc', 'p0330009_20100905.nc', 'p0330010_20100905.nc', 'p0330011_20100905.nc', 'p0330012_20100905.nc', 'p0330013_20100906.nc', 'p0330014_20100906.nc', 'p0330015_20100906.nc', 'p0330016_20100906.nc', 'p0330017_20100906.nc', 'p0330018_20100906.nc', 'p0330019_20100906.nc']"

[ ]: