From e75df9fe5458f647dded8a9fe89d1fdd9afd39b1 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Tue, 20 Aug 2024 15:38:18 +0200
Subject: [PATCH] fix using interpos for linear interpolation with constant
extrapolation
---
matlab/TCV_IMAS/tcv_get_ids_summary.m | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv_get_ids_summary.m b/matlab/TCV_IMAS/tcv_get_ids_summary.m
index 5bc8827a..7fc129ec 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_summary.m
+++ b/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};
--
GitLab