corrected check on exit from plasma

This commit is contained in:
Lorenzo Figini 2012-08-29 16:34:47 +00:00
parent 1b037fa88a
commit c78d88c7a8

View File

@ -257,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.1.0d0) iov(j,k)=2 if(iov(j,k).eq.1.and.psinv.ge.psdbnd) 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