From 6a937eeb008c87b785c10a801e8f111e42da5e7f Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Wed, 16 Jan 2019 09:04:09 +0000
Subject: [PATCH] fix finding sign of q
git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@11288 d63d8f72-b253-0410-a779-e742ad2e26cf
---
crpptbx/AUG/gdat_aug.m | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/crpptbx/AUG/gdat_aug.m b/crpptbx/AUG/gdat_aug.m
index 3f9caaf6..2b37df26 100644
--- a/crpptbx/AUG/gdat_aug.m
+++ b/crpptbx/AUG/gdat_aug.m
@@ -931,9 +931,9 @@ elseif strcmp(mapping_for_aug.method,'switchcase')
Lpf1 = Lpf1_t(it);
% Qpsi and similar data is on (time,radius) with radius being: LCFS..Lpf_points dummy LCFS..SOL part
% change it to (radial,time) and use only Lpf+1 points up to LCFS
- ijok=find(qpsi.value(:,1)); % note: eqr fills in only odd points radially
+ ijok=find(abs(qpsi.value)>0); % note: eqr fills in only odd points radially
% set NaNs to zeroes
- if qpsi.value(ijok(1),1)<0
+ if qpsi.value(ijok(1))<0
gdat_data.qvalue(:,it) = max(qpsi.value(it,Lpf1:-1:1)',-max_qValue);
else
gdat_data.qvalue(:,it) = min(qpsi.value(it,Lpf1:-1:1)',max_qValue);
--
GitLab