istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

added check (z only) to avoid absorption outside the plasma boundary

Browse Source
This commit is contained in:
Lorenzo Figini 2012-08-29 16:22:02 +00:00
parent 6c4682fb56
commit 1b037fa88a
1 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
src/gray.f
View File

@ -209,6 +209,7 @@ c
common/index_rt/index_rt common/index_rt/index_rt
common/ipass/ipass common/ipass/ipass
common/rwallm/rwallm common/rwallm/rwallm
common/bound/zbmin,zbmax
pabstot=0.0d0 pabstot=0.0d0
currtot=0.0d0 currtot=0.0d0
@ -236,10 +237,12 @@ c ierr=0
psjki(j,k,i)=psinv psjki(j,k,i)=psinv
rrm=sqrt(xv(1)**2+xv(2)**2)/100.d0 rrm=sqrt(xv(1)**2+xv(2)**2)/100.d0
zzm=xv(3)/100.d0
if(j.eq.1) rrm11=rrm if(j.eq.1) rrm11=rrm
if (iwarm.gt.0.and.i.gt.1) then if (iwarm.gt.0.and.i.gt.1) then
if(psinv.ge.0.and.psinv.le.1.0d0) then if(psinv.ge.0.and.psinv.le.1.0d0.and.
. zzm.ge.zbmin.and.zzm.le.zbmax) then
call pabs_curr(i,j,k) call pabs_curr(i,j,k)
iiv(j,k)=i iiv(j,k)=i
else else
@ -254,7 +257,7 @@ c ierr=0
if(i.gt.1.and.psinv.ge.0.and.psinv.lt.psdbnd) if(i.gt.1.and.psinv.ge.0.and.psinv.lt.psdbnd)
. iov(j,k)=1 . iov(j,k)=1
if(iov(j,k).eq.1.and.psinv.ge.psdbnd) iov(j,k)=2 if(iov(j,k).eq.1.and.psinv.ge.1.0d0) iov(j,k)=2
c iov=0 initially, iov=1 first entrance in plasma, c iov=0 initially, iov=1 first entrance in plasma,
c iov=2 first exit from plasma, iov=3 after 2nd entrance into plasma c iov=2 first exit from plasma, iov=3 after 2nd entrance into plasma