fix using interpos for linear interpolation with constant extrapolation

matlab/TCV_IMAS/tcv_get_ids_summary.m

@@ -100,7 +100,7 @@ for i=1:numel(boundary_fieldnames)
     if ~isstruct(ids_summary.boundary.(boundary_fieldnames{i}).value)
       ids_summary.boundary.(boundary_fieldnames{i}).value = boundary.(boundary_fieldnames{i}).data;
       ids_summary.boundary.(boundary_fieldnames{i}).value = interpos(21,boundary.(boundary_fieldnames{i}).t, ...
-          boundary.(boundary_fieldnames{i}).data,ids_summary.time,'linear',NaN);
+          boundary.(boundary_fieldnames{i}).data,ids_summary.time);
       ids_summary.boundary.(boundary_fieldnames{i}).source = ['gdat request: ' boundary_desc.(boundary_fieldnames{i})];
     else
       special_fields{end+1} = boundary_fieldnames{i};
@@ -140,7 +140,7 @@ for i=1:numel(gas_injection_rates_fieldnames)
     if ~isstruct(ids_summary.gas_injection_rates.(gas_injection_rates_fieldnames{i}).value)
       ids_summary.gas_injection_rates.(gas_injection_rates_fieldnames{i}).value = interpos(21, ...
           gas_injection_rates.(gas_injection_rates_fieldnames{i}).t,gas_injection_rates.(gas_injection_rates_fieldnames{i}).data, ...
-          ids_summary.time,'linear',NaN);
+          ids_summary.time);
       ids_summary.gas_injection_rates.(gas_injection_rates_fieldnames{i}).source = ['gdat request: ' gas_injection_rates_desc.(gas_injection_rates_fieldnames{i})];
     else
       special_fields{end+1} = gas_injection_rates_fieldnames{i};
@@ -255,7 +255,7 @@ for i=1:numel(global_quantities_fieldnames)
         ~isempty(global_quantities.(global_quantities_fieldnames{i}).data)
       ids_summary.global_quantities.(global_quantities_fieldnames{i}).value = interpos(21, ...
           global_quantities.(global_quantities_fieldnames{i}).t,global_quantities.(global_quantities_fieldnames{i}).data, ...
-          ids_summary.time,'linear',NaN);
+          ids_summary.time);
       ids_summary.global_quantities.(global_quantities_fieldnames{i}).source = ['gdat request: ' global_quantities_desc.(global_quantities_fieldnames{i})];
     elseif ~isempty(global_quantities.(global_quantities_fieldnames{i}).data)
       special_fields{end+1} = global_quantities_fieldnames{i};
@@ -280,7 +280,7 @@ for i=1:numel(volume_average_fieldnames)
     if ~isstruct(ids_summary.volume_average.(volume_average_fieldnames{i}).value)
       ids_summary.volume_average.(volume_average_fieldnames{i}).value = interpos(21, ...
           volume_average.(volume_average_fieldnames{i}).t,volume_average.(volume_average_fieldnames{i}).data, ...
-          ids_summary.time,'linear',NaN);
+          ids_summary.time);
       ids_summary.volume_average.(volume_average_fieldnames{i}).source = ['gdat request: ' volume_average_desc.(volume_average_fieldnames{i})];
     else
       special_fields{end+1} = volume_average_fieldnames{i};
@@ -298,9 +298,9 @@ line_average_fieldnames = fieldnames(line_average);
 for i=1:numel(line_average_fieldnames)
   if ~any(strcmp(line_average_fieldnames{i},special_fields))
     if ~isstruct(ids_summary.line_average.(line_average_fieldnames{i}).value)
-      ids_summary.line_average.(line_average_fieldnames{i}).value = interp1( ...
+      ids_summary.line_average.(line_average_fieldnames{i}).value = interpos(21, ...
           line_average.(line_average_fieldnames{i}).t,line_average.(line_average_fieldnames{i}).data, ...
-          ids_summary.time,'linear',NaN);
+          ids_summary.time);
       ids_summary.line_average.(line_average_fieldnames{i}).source = ['gdat request: ' line_average_desc.(line_average_fieldnames{i})];
     else
       special_fields{end+1} = line_average_fieldnames{i};