Skip to content
Snippets Groups Projects
Commit 667f67b8 authored by Olivier Sauter's avatar Olivier Sauter
Browse files

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
parent 47452573
No related branches found
No related tags found
No related merge requests found
......@@ -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;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment