Handling of situations where all values for a variable are nans

matlab/TCV_IMAS/tcv_get_ids_summary.m

@@ -252,7 +252,12 @@ for i=1:numel(global_quantities_fieldnames)
         ~isempty(global_quantities.(global_quantities_fieldnames{i}).data)
       % setup mask to get rid of nans
       mask = ~isnan(global_quantities.(global_quantities_fieldnames{i}).data);
-      % interpolate quantity on ids_summary.time
+      % assign a flat 1 to the array if it is all nans
+      if all(mask(:) == 0)
+	fprintf('Warning, %s is filled with NaNs, substituting 1\n', global_quantities_fieldnames{i});
+        global_quantities.(global_quantities_fieldnames{i}).data(:)=1;
+	mask(:)=1;
+      end
       ids_summary.global_quantities.(global_quantities_fieldnames{i}).value = ...
         interpos(21,global_quantities.(global_quantities_fieldnames{i}).t(mask),...
         global_quantities.(global_quantities_fieldnames{i}).data(mask), ...