diff --git a/pyglider/seaexplorer.py b/pyglider/seaexplorer.py
index 6451af7..3a6a422 100644
--- a/pyglider/seaexplorer.py
+++ b/pyglider/seaexplorer.py
@@ -302,8 +302,7 @@ def _remove_fill_values(df, fill_value=9999):
 
 
 def raw_to_timeseries(indir, outdir, deploymentyaml, kind='raw',
-                      profile_filt_time=100, profile_min_time=300,
-                      fnamesuffix=''):
+                      profile_filt_time=100, profile_min_time=300, interpolate=False, fnamesuffix=''):
     """
     A little different than above, for the 4-file version of the data set.
     """
@@ -339,7 +338,11 @@ def raw_to_timeseries(indir, outdir, deploymentyaml, kind='raw',
     indctd = np.where(~np.isnan(vals))[0]
     ds['time'] = (('time'), sensor.select('time').to_numpy()[indctd, 0], attr)
     thenames = list(ncvar.keys())
-    for i in ['time', 'timebase', 'keep_variables']:
+    # Check yaml to see if interpolate has been set to True
+    if "interpolate" in thenames:
+        if ncvar["interpolate"]:
+            interpolate = True
+    for i in ['time', 'timebase', 'keep_variables', 'interpolate']:
         if i in thenames:
             thenames.remove(i)
     for name in thenames:
@@ -371,7 +374,21 @@ def raw_to_timeseries(indir, outdir, deploymentyaml, kind='raw',
                     )[:-1, :]
                     val2 = sensor_sub_grouped.select(sensorname).to_numpy()[:, 0]
                     val = _interp_gli_to_pld(sensor_sub_grouped, sensor, val2, indctd)
-                val = val[indctd]
+                if interpolate and not np.isnan(val).all():
+                    time_original = sensor.select('time').to_numpy()[:, 0]
+                    time_var = time_original[np.where(~np.isnan(val))[0]]
+                    var_non_nan = val[np.where(~np.isnan(val))[0]]
+                    time_timebase = sensor.select('time').to_numpy()[indctd, 0]
+                    if val.dtype == "<M8[us]":
+                        # for datetime, must convert to numerical, interpolate, then convert back
+                        us_since_1970 = (var_non_nan - np.datetime64("1970-01-01")).astype(int)
+                        val_int = np.interp(time_timebase.astype(float), time_var.astype(float), us_since_1970)
+                        val_us = val_int.astype("timedelta64[us]")
+                        val = np.datetime64("1970-01-01") + val_us
+                    else:
+                        val = np.interp(time_timebase.astype(float), time_var.astype(float), var_non_nan)
+                else:
+                    val = val[indctd]
 
                 ncvar['method'] = 'linear fill'
             else:
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.10.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.10.gz
new file mode 100644
index 0000000..3815411
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.10.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.11.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.11.gz
new file mode 100644
index 0000000..de0b6cd
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.11.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.12.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.12.gz
new file mode 100644
index 0000000..0e3ab44
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.12.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.13.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.13.gz
new file mode 100644
index 0000000..c7b6720
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.13.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.14.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.14.gz
new file mode 100644
index 0000000..72eb16a
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.gli.sub.14.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.10.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.10.gz
new file mode 100644
index 0000000..f3b0f61
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.10.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.11.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.11.gz
new file mode 100644
index 0000000..26d145a
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.11.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.12.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.12.gz
new file mode 100644
index 0000000..a4dbbbf
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.12.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.13.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.13.gz
new file mode 100644
index 0000000..a03e090
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.13.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.14.gz b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.14.gz
new file mode 100644
index 0000000..c2de77a
Binary files /dev/null and b/tests/example-data/example-seaexplorer-raw/delayed_raw/sea046.45.pld1.raw.14.gz differ
diff --git a/tests/example-data/example-seaexplorer-raw/deployment.yml b/tests/example-data/example-seaexplorer-raw/deployment.yml
new file mode 100644
index 0000000..ed82130
--- /dev/null
+++ b/tests/example-data/example-seaexplorer-raw/deployment.yml
@@ -0,0 +1,338 @@
+metadata:
+  # https://github.com/ioos/ioosngdac/wiki/NGDAC-NetCDF-File-Format-Version-2
+  acknowledgement: Funding from Fisheries and Oceans Canada, Canadian Foundation
+                   for Innovation, BC Knowledge Development Fund
+  comment:        "Calvert Island deployment"
+  contributor_name: James Pegg, Jody Klymak, Tetjana Ross, Hayley Dosser
+  contributor_role: Lead Technician, Principal Investigator, Co-PI, Ocean Glider Scientist
+  creator_email: jklymak@uvic.ca
+  creator_name:  Jody Klymak
+  creator_url:   http://cproof.uvic.ca
+  # date_created etc: added automatically
+  # numbers must be enclosed by strings so YAML keeps as strings
+  deployment_id: '45'
+  deployment_name: 'dfo-bb046-20200908'
+  deployment_start: '2020-09-08'
+  deployment_end: '2020-09-30'
+  format_version: IOOS_Glider_NetCDF_v2.0.nc
+  # id filled automatically...
+  glider_name: dfo-bb
+  glider_serial: '046'
+  glider_model: SeaExplorer
+  glider_instrument_name: seaexplorer
+  glider_wmo: "4803918"
+  institution: C-PROOF
+  keywords: "AUVS, Autonomous Underwater Vehicles, Oceans, Ocean Pressure,
+             Water Pressure, Oceans, Ocean Temperature, Water Temperature,
+             Oceans, Salinity/Density, Conductivity, Oceans,
+             Salinity/Density, Density, Oceans, Salinity/Density, Salinity"
+  keywords_vocabulary: GCMD Science Keywords
+  license: "This data may be redistributed and used without restriction or
+            warranty"
+  metadata_link: "https://cproof.uvic.ca"
+  Metadata_Conventions: CF-1.6, Unidata Dataset Discovery v1.0
+  naming_authority: "ca.uvic.cproof"
+  platform_type:    "SeaExplorer Glider"
+  processing_level: "Data provided as is with no expressed or implied
+                     assurance of quality assurance or quality control."
+  project: Calvert Island
+  project_url: http://cproof.uvic.ca
+  publisher_email: jklymak@uvic.ca
+  publisher_name:  Jody Klymak
+  publisher_url:   http://cproof.uvic.ca
+  references:     cproof toolbox URL
+  # https://www.nodc.noaa.gov/General/NODC-Archive/seanamelist.txt
+  sea_name:   BC Coastal Waters
+  source:     Observational data from a profiling glider.
+  standard_name_vocabulary: CF STandard Name Table v49
+  summary: Transect from Calvert Island to beyond shelf break in Queen Charlotte Sound.
+  transmission_system: IRIDIUM
+  wmo_id: "4803918"
+
+
+glider_devices:
+  pressure:
+    make: Micron
+    model: Pressure
+    serial: ' '
+  ctd:
+    make: Seabird
+    model: GPCTD
+    serial: '0278'
+    long_name: Seabird GPCTD
+    make_model: Seabird GPCTD
+    factory_calibrated: "Yes"
+    calibration_date: "20181102"
+    calibration_report: " "
+    comment:   " "
+  optics:
+    make: Wetlabs
+    model: FLBBCDSLC
+    serial: '4741'
+    calibration_date: "20180807"
+  oxygen:
+    make: AROD_FT
+    model: Optode4831
+    serial: '0022'
+    calibration_date: "20191207"
+
+# map between glider variables and netcdf variables.  This shouldn't
+# change too much.
+netcdf_variables:
+  # Time and Place:
+  timebase:
+    source: GPCTD_TEMPERATURE
+  time:
+    source:        time
+    long_name:     Time
+    standard_name: time
+    calendar:      gregorian
+    units:         seconds since 1970-01-01T00:00:00Z
+    axis:          T
+    observation_type: "measured"
+    coordinates:   time depth latitude longitude
+
+  latitude:
+    source:       NAV_LATITUDE
+    long_name:    latitude
+    standard_name: latitude
+    units:        degrees_north
+    axis:         Y
+    coordinates:   time depth latitude longitude
+    conversion:   nmea2deg
+    comment:     "Estimated between surface fixes"
+    observation_type: measured
+    platform:     platform
+    reference:    WGS84
+    valid_max:    "90.0"
+    valid_min:    "-90.0"
+    coordinate_reference_frame:  urn:ogc:crs:EPSG::4326
+
+  longitude:
+    source:       NAV_LONGITUDE
+    long_name:    longitude
+    standard_name: longitude
+    units:        degrees_east
+    axis:         X
+    coordinates:  time depth latitude longitude
+    conversion:   nmea2deg
+    comment:     "Estimated between surface fixes"
+    observation_type: measured
+    platform:     platform
+    reference:    WGS84
+    valid_max:    "180.0"
+    valid_min:    "-180.0"
+    coordinate_reference_frame:  urn:ogc:crs:EPSG::4326
+
+  heading:
+    source:       Heading
+    long_name:    glider heading angle
+    standard_name: platform_orientation
+    units:        rad
+    coordinates:  time depth latitude longitude
+
+  pitch:
+    source:       Pitch
+    long_name:    glider pitch angle
+    standard_name: platform_pitch_angle
+    units:        rad
+    coordinates:  time depth latitude longitude
+
+  roll:
+    source:       Roll
+    long_name:    glider roll angle
+    standard_name: platform_roll_angle
+    units:        rad
+    coordinates:  time depth latitude longitude
+
+  # data parameters
+  conductivity:
+    source:       GPCTD_CONDUCTIVITY
+    long_name:    water conductivity
+    standard_name: sea_water_electrical_conductivity
+    units:        S m-1
+    coordinates:  time depth latitude longitude
+    instrument:    instrument_ctd
+    valid_min:    "0."
+    valid_max:    "10."
+    observation_type: "measured"
+    accuracy:      "0.0003"
+    precision:     "0.0001"
+    resolution:     "0.00002"
+
+  temperature:
+    source:       GPCTD_TEMPERATURE
+    long_name:    water temperature
+    standard_name: sea_water_temperature
+    units:        Celsius
+    coordinates:  time depth latitude longitude
+    instrument:   instrument_ctd
+    valid_min:    "-5"
+    valid_max:    "50"
+    observation_type: "measured"
+    accuracy:      "0.002"
+    precision:     "0.001"
+    resolution:    "0.0002"
+
+  pressure:
+    source:       GPCTD_PRESSURE
+    long_name:    water pressure
+    standard_name:  sea_water_pressure
+    units:        dbar
+    coordinates:  time depth latitude longitude
+    valid_min:    "0"
+    valid_max:    "2000"
+    positive:      "down"
+    reference_datum:  "sea-surface"
+    instrument:     "instrument_ctd"
+    observation_type: "measured"
+    accuracy:         "1"
+    precision:        "2"
+    resolution:       "0.02"
+    comment:          "ctd pressure sensor"
+
+# optics:
+  chlorophyll:
+    source:       FLBBCD_CHL_SCALED
+    long_name:    chlorophyll
+    standard_name: concentration_of_chlorophyll_in_sea_water
+    units:        mg m-3
+    coordinates:  time depth latitude longitude
+
+  cdom:
+    source:  FLBBCD_CDOM_SCALED
+    long_name:    CDOM
+    units:        ppb
+    coordinates:  time depth latitude longitude
+
+  backscatter_700:
+    source:       FLBBCD_BB_700_SCALED
+    long_name:    700 nm wavelength backscatter
+    units:         "1"
+    coordinates:  time depth latitude longitude
+
+# Oxygen
+  oxygen_concentration:
+    source:       AROD_FT_DO
+    long_name:    oxygen concentration
+    standard_name: mole_concentration_of_dissolved_molecular_oxygen_in_sea_water
+    units:        umol l-1
+    coordinates:   time depth latitude longitude
+
+  temperature_oxygen:
+    source:  AROD_FT_TEMP
+    long_name:    oxygen sensor temperature
+    standard_name: temperature_of_sensor_for_oxygen_in_sea_water
+    units:        Celsius
+    coordinates:   time depth latitude longitude
+
+profile_variables:
+  # variables for extract_L1timeseries_profiles processing step...
+  profile_id:
+    comment: Sequential profile number within the trajectory.  This value is unique in each file that is part of a single trajectory/deployment.
+    long_name: 'Profile ID'
+    valid_max: '2147483647'
+    valid_min: '1'
+
+  profile_time:
+    comment:           Timestamp corresponding to the mid-point of the profile
+    long_name:         Profile Center Time
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     time
+    units:             seconds since 1970-01-01T00:00:00Z
+
+  profile_time_start:
+    comment:           Timestamp corresponding to the start of the profile
+    long_name:         Profile Start Time
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     time
+    units:             seconds since 1970-01-01T00:00:00Z
+
+  profile_time_end:
+    comment:           Timestamp corresponding to the end of the profile
+    long_name:         Profile End Time
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     time
+    units:             seconds since 1970-01-01T00:00:00Z
+
+  profile_lat:
+    comment:           Value is interpolated to provide an estimate of the latitude at the mid-point of the profile
+    long_name:         Profile Center Latitude
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     latitude
+    units:             degrees_north
+    valid_max:         "90.0"
+    valid_min:         "-90.0"
+
+  profile_lon:
+    comment:           Value is interpolated to provide an estimate of the latitude at the mid-point of the profile
+    long_name:         Profile Center Longitude
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     longitude
+    units:             degrees_east
+    valid_max:         "180.0"
+    valid_min:         "-180.0"
+
+  u:
+    comment:  The depth-averaged current is an estimate of the net current measured while the glider is underwater.  The value is calculated over the entire underwater segment, which may consist of 1 or more dives.
+    long_name:         Depth-Averaged Eastward Sea Water Velocity
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     eastward_sea_water_velocity
+    units:             m s-1
+    valid_max:         "10.0"
+    valid_min:         "-10.0"
+
+  v:
+    comment:  The depth-averaged current is an estimate of the net current measured while the glider is underwater.  The value is calculated over the entire underwater segment, which may consist of 1 or more dives.
+    long_name:         Depth-Averaged Northward Sea Water Velocity
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     northward_sea_water_velocity
+    units:             m s-1
+    valid_max:         "10.0"
+    valid_min:         "-10.0"
+
+  lon_uv:
+    comment:           Not computed
+    long_name:         Longitude
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     longitude
+    units:             degrees_east
+    valid_max:         "180.0"
+    valid_min:         "-180.0"
+
+  lat_uv:
+    comment:           Not computed
+    long_name:         Latitude
+    observation_type:  calculated
+    platform:          platform
+    standard_name:     latitude
+    units:             degrees_north
+    valid_max:         "90.0"
+    valid_min:         "-90.0"
+
+  time_uv:
+    comment:       Not computed
+    long_name:     Time
+    standard_name: time
+    calendar:      gregorian
+    units:         seconds since 1970-01-01T00:00:00Z
+    observation_type: calculated
+
+  instrument_ctd:
+    comment:    pumped CTD
+    calibration_date:     "2018-11-02"
+    calibration_report:   
+    factory_calibrated:  "yes"
+    long_name:           Seabird Glider Payload CTD
+    make_model:          Seabird GPCTD
+    platform:            platform
+    serial_number:       "0278"
+    type:                platform
diff --git a/tests/example-data/example-seaexplorer-raw/process_delayed.py b/tests/example-data/example-seaexplorer-raw/process_delayed.py
new file mode 100644
index 0000000..479dbff
--- /dev/null
+++ b/tests/example-data/example-seaexplorer-raw/process_delayed.py
@@ -0,0 +1,25 @@
+import logging
+import os
+import pyglider.seaexplorer as seaexplorer
+
+logging.basicConfig(level='INFO')
+
+rawdir  = './delayed_raw/'
+rawncdir     = './realtime_rawnc/'
+deploymentyaml = './deployment.yml'
+l0tsdir    = './L0-timeseries/'
+profiledir = './L0-profiles/'
+griddir    = './L0-gridfiles/'
+
+
+if __name__ == '__main__':
+    # clean last processing
+    os.system('rm ' + rawncdir + '* ' + l0tsdir + '* ' + profiledir + '* ' +
+              griddir + '* ')
+    # Convert from csv to parquet
+    seaexplorer.raw_to_rawnc(rawdir, rawncdir, deploymentyaml)
+    # merge individual parquet files
+    seaexplorer.merge_parquet(rawncdir, rawncdir, deploymentyaml, kind='raw')
+
+    # Make level-1 timeseries netcdf file from the raw files
+    outname = seaexplorer.raw_to_timeseries(rawncdir, l0tsdir, deploymentyaml, kind='raw')
diff --git a/tests/results/example-seaexplorer-raw/L0-timeseries/dfo-bb046-20200908.nc b/tests/results/example-seaexplorer-raw/L0-timeseries/dfo-bb046-20200908.nc
new file mode 100644
index 0000000..55d56cd
Binary files /dev/null and b/tests/results/example-seaexplorer-raw/L0-timeseries/dfo-bb046-20200908.nc differ
diff --git a/tests/test_pyglider.py b/tests/test_pyglider.py
index 1313252..e5663ca 100644
--- a/tests/test_pyglider.py
+++ b/tests/test_pyglider.py
@@ -2,6 +2,7 @@
 from pathlib import Path
 import pytest
 import numpy as np
+import yaml
 
 import pyglider.seaexplorer as seaexplorer
 import pyglider.slocum as slocum
@@ -57,6 +58,98 @@ def test_example_seaexplorer_metadata():
     assert output.attrs == test_data.attrs
 
 
+# Test that interpolation over nans does not change the output with nrt data
+with open(deploymentyaml) as fin:
+    deployment = yaml.safe_load(fin)
+interp_yaml = str(example_dir / 'example-seaexplorer/deploymentRealtimeInterp.yml')
+deployment['netcdf_variables']["interpolate"] = True
+with open(interp_yaml, "w") as fout:
+    yaml.dump(deployment, fout)
+l0tsdir_interp = str(example_dir / 'example-seaexplorer/L0-timeseries-test-interp/') + '/'
+
+outname_interp = seaexplorer.raw_to_L0timeseries(rawncdir, l0tsdir_interp, interp_yaml, kind='sub')
+output_interp = xr.open_dataset(outname_interp)
+
+@pytest.mark.parametrize("var", variables)
+def test_example_seaexplorer_interp_nrt(var):
+    assert output_interp[var].attrs == test_data[var].attrs
+    if var not in ['time']:
+        np.testing.assert_allclose(output_interp[var].values, test_data[var].values, rtol=1e-5)
+    else:
+        dt0 = output_interp[var].values - np.datetime64('2000-01-01')
+        dt1 = test_data[var].values - np.datetime64('2000-01-01')
+        assert np.allclose(
+            np.array(dt0, dtype='float64'),
+            np.array(dt1, dtype='float64'))
+
+
+# Test raw (full resolution) seaexplorer data.
+
+rawdir = str(example_dir / 'example-seaexplorer-raw/delayed_raw/') + '/'
+rawncdir = str(example_dir / 'example-seaexplorer-raw/delayed_rawnc/') + '/'
+deploymentyaml_raw = str(example_dir / 'example-seaexplorer-raw/deployment.yml')
+l0tsdir = str(example_dir / 'example-seaexplorer-raw/L0-timeseries-test/') + '/'
+seaexplorer.raw_to_rawnc(rawdir, rawncdir, deploymentyaml_raw)
+seaexplorer.merge_parquet(rawncdir, rawncdir, deploymentyaml_raw, kind='raw')
+outname_raw = seaexplorer.raw_to_L0timeseries(rawncdir, l0tsdir,
+                                          deploymentyaml_raw, kind='raw')
+output_raw = xr.open_dataset(outname_raw)
+# Open test data file
+test_data_raw = xr.open_dataset(
+    library_dir /
+    'tests/results/example-seaexplorer-raw/L0-timeseries/dfo-bb046-20200908.nc')
+
+@pytest.mark.parametrize("var", variables)
+def test_example_seaexplorer_raw(var):
+    # Test that each variable and its coordinates match
+    assert output_raw[var].attrs == test_data_raw[var].attrs
+    if var not in ['time']:
+        np.testing.assert_allclose(output_raw[var].values, test_data_raw[var].values, rtol=1e-5)
+    else:
+        dt0 = output_raw[var].values - np.datetime64('2000-01-01')
+        dt1 = test_data_raw[var].values - np.datetime64('2000-01-01')
+        assert np.allclose(
+            np.array(dt0, dtype='float64'),
+            np.array(dt1, dtype='float64'))
+
+
+def test_example_seaexplorer_metadata_raw():
+    # Test that attributes match. Have to remove creation and issue dates first
+    output_raw.attrs.pop('date_created')
+    output_raw.attrs.pop('date_issued')
+    test_data_raw.attrs.pop('date_created')
+    test_data_raw.attrs.pop('date_issued')
+    assert output_raw.attrs == test_data_raw.attrs
+
+
+# Test nan interpolation on raw data.
+
+# Test that interpolation over nans in raw data results in a greater or equal number of non-nan values
+with open(deploymentyaml_raw) as fin:
+    deployment_raw = yaml.safe_load(fin)
+interp_yaml = str(example_dir / 'example-seaexplorer-raw/deploymentDelayedInterp.yml')
+deployment_raw['netcdf_variables']["interpolate"] = True
+with open(interp_yaml, "w") as fout:
+    yaml.dump(deployment_raw, fout)
+l0tsdir_interp_raw = str(example_dir / 'example-seaexplorer-raw/L0-timeseries-test-interp/') + '/'
+
+outname_interp_raw = seaexplorer.raw_to_L0timeseries(rawncdir, l0tsdir_interp_raw, interp_yaml, kind='raw')
+output_interp_raw = xr.open_dataset(outname_interp_raw)
+
+@pytest.mark.parametrize("var", variables)
+def test_example_seaexplorer_interp_raw(var):
+    assert output_interp_raw[var].attrs == test_data_raw[var].attrs
+    if var not in ['time']:
+        assert np.count_nonzero(~np.isnan(output_interp_raw[var].values)) >= \
+               np.count_nonzero(~np.isnan(test_data_raw[var].values))
+    else:
+        dt0 = output_interp_raw[var].values - np.datetime64('2000-01-01')
+        dt1 = test_data_raw[var].values - np.datetime64('2000-01-01')
+        assert np.allclose(
+            np.array(dt0, dtype='float64'),
+            np.array(dt1, dtype='float64'))
+
+
 # Create an L0 timeseries from slocum data and test that the resulting netcdf is
 # identical to the test data
 cacdir = str(example_dir / 'example-slocum/cac/') + '/'