From 9eb42eb6498f17cb66a2474e4a787a30aabe9868 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Wed, 27 Apr 2016 09:30:02 +0000
Subject: [PATCH] add warning if problem with fir ratio with Thomson, keep ne
 as Nans but can use ne_raw

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@5701 d63d8f72-b253-0410-a779-e742ad2e26cf
---
 crpptbx/TCV/gdat_tcv.m | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index dd01515f..a413e5ec 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -1762,6 +1762,13 @@ if strcmp(data_request_eff(1:2),'ne')
   gdat_data.error_bar_raw = gdat_data.error_bar;
   gdat_data.error_bar = gdat_data.error_bar_raw * diag(tracefirrat_data);
   gdat_data.data_fullpath=[gdat_data.data_fullpath '; fir_thom_rat ; _raw without firrat'];
+  ij=find(~isnan(tracefirrat_data));
+  if isempty(ij)
+    gdat_data.data_fullpath=[gdat_data.data_fullpath ' WARNING use ne_raw since firrat has NaNs thus ne=NaNs; fir_thom_rat ; _raw without firrat'];
+    disp('***********************************************************************')
+    disp('WARNING: ne from Thomson has fir_thom_rat with Nans so only ne_raw can be used');
+    disp('***********************************************************************')
+  end
 end
 z=mdsdata('\diagz::thomson_set_up:vertical_pos');
 gdat_data.dim=[{z};{time}];
-- 
GitLab