From 667f67b8057ad98c1ff2765f30c74d3f026feedb Mon Sep 17 00:00:00 2001
From: Olivier Sauter <olivier.sauter@epfl.ch>
Date: Wed, 19 Mar 2014 17:17:34 +0000
Subject: [PATCH] fix getting plasma edge from 0.999 psiedge to make sure gets
 a closed surface

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@4336 d63d8f72-b253-0410-a779-e742ad2e26cf
---
 crpptbx/AUG/geteqdskAUG.m | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/crpptbx/AUG/geteqdskAUG.m b/crpptbx/AUG/geteqdskAUG.m
index 948e52cb..87b8f41a 100644
--- a/crpptbx/AUG/geteqdskAUG.m
+++ b/crpptbx/AUG/geteqdskAUG.m
@@ -75,7 +75,8 @@ eqdsk.zaxis = zmag.data(itrmag);
 figure
 contour(eqdsk.rmesh,eqdsk.zmesh,eqdsk.psi',100)
 hold
-[hh1 hh2]=contour(eqdsk.rmesh,eqdsk.zmesh,eqdsk.psi',[eqdsk.psiedge eqdsk.psiedge],'k');
+psiedge_eff = 0.001*eqdsk.psiaxis + 0.999*eqdsk.psiedge;
+[hh1 hh2]=contour(eqdsk.rmesh,eqdsk.zmesh,eqdsk.psi',[psiedge_eff psiedge_eff],'k');
 axis equal
 ij=1;
 ij_prev = 0;
-- 
GitLab