feat: Extrapolate TOA/ETA differences beyond TOA range

lib/+artoa/+controller/+show/+satfixEtaToToa/calculateDifference.m

@@ -6,6 +6,17 @@ global artoaWorkspace artoaConfig artoaDataInput;
 
 %% Get required variables
 fnames = fieldnames(artoaWorkspace.filteredSoundsources);
+launchDateVec = artoaWorkspace.float.launchtime;
+if length(launchDateVec) == 5
+    launchDateVec(6) = 0;
+end
+launchDate = artoa.convert.datevec2rd(launchDateVec);
+
+recoveryDateVec = artoaWorkspace.float.recovertime;
+if length(recoveryDateVec) == 5
+    recoveryDateVec(6) = 0;
+end
+recoveryDate = artoa.convert.datevec2rd(recoveryDateVec);
 
 %% Get soundspeed
 floatSoundspeed = artoaWorkspace.editOffsets.offsets.AppliedSoundspeed(1,1);
@@ -35,10 +46,10 @@ for i = 1:length(fnames)
     toas = artoaWorkspace.toaData.toa(selection);
     toaDates = ceil(artoaWorkspace.toaData.toaDate(selection));
     if length(toaDates) > 1
-
         % interpolate
-        toas = interp1(toaDates, toas, min(toaDates):1:max(toaDates))';
-        toaDates = [min(toaDates):1:max(toaDates)]';
+        newToaDates = floor(launchDate):1:ceil(recoveryDate);
+        toas = interp1(toaDates, toas, newToaDates, 'linear', 'extrap')';
+        toaDates = newToaDates';
     end
     dates = intersect(toaDates, satfixEtas.(fnames{i}).satDate);
     indicesToa = ismember(toaDates, dates);