From 860143c5d8ba08cc47f04d28091f7b2ddbc1bba5 Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Thu, 3 Nov 2016 16:21:00 +0000
Subject: [PATCH] get Ipi from IDG in case of IDE equilibrium
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@6698 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/AUG/gdat_aug.m | 10 +++++++++-
1 file changed, 9 insertions(+), 1 deletion(-)
diff --git a/crpptbx/AUG/gdat_aug.m b/crpptbx/AUG/gdat_aug.m
index 86538d89..71e6e268 100644
--- a/crpptbx/AUG/gdat_aug.m
+++ b/crpptbx/AUG/gdat_aug.m
@@ -969,7 +969,15 @@ elseif strcmp(mapping_for_aug.method,'switchcase')
[equil_Zcoil,e]=rdaAUG_eff(shot,DIAG,'Zcl',exp_name_eff);
gdat_data.Zcoils=equil_Zcoil.value;
%
- [IpiPSI,e]=rdaAUG_eff(shot,DIAG,'IpiPSI',exp_name_eff);
+ if strcmp(DIAG,'IDE')
+ [IpiPSI,e]=rdaAUG_eff(shot,'IDG','Itor',exp_name_eff);
+ if length(IpiPSI.value)~=NTIME
+ disp('Problem with nb time points of IDG/Itor with respect to IDE NTIME, check with O. Sauter')
+ return
+ end
+ else
+ [IpiPSI,e]=rdaAUG_eff(shot,DIAG,'IpiPSI',exp_name_eff);
+ end
gdat_data.Ip = IpiPSI.value(1:NTIME);
%
gdat_data.data = gdat_data.qvalue; % put q in data
--
GitLab