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

corrected bug in equinum_fpol: vacuum value was returned for psin=0

Browse Source
This commit is contained in:
Lorenzo Figini 2015-10-02 15:09:10 +00:00
parent dcd3badba9
commit 61c97e3357
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/gray.f
View File

@ -660,7 +660,7 @@ c
use coreprofiles, only : psdbnd,read_profiles,tene_scal,set_prfspl
use interp_eqprof, only : rmxm,zbmin,zbmax,rcen,psia,rmaxis,
. zmaxis,rbbbs,zbbbs,read_eqdsk,change_cocos,eq_scal,set_eqspl,
. btrcen,psinr=>psiq,qpsi=>q
. btrcen,psinr=>psiq,qpsi=>q,equinum_fpol
use reflections, only : rlim,zlim,nlim,alloc_lim
use beamdata, only : nrayr,nrayth,nstep
implicit none

2
src/interp_eqprof.f90
View File

@ -440,7 +440,7 @@ contains
real(wp_), dimension(1) :: rrs,ffspl
real(wp_), dimension(nsf) :: wrkfd
!
if(psinv.le.1.0_wp_.and.psinv.gt.0.0_wp_) then
if(psinv.le.1.0_wp_.and.psinv.ge.0.0_wp_) then
rrs(1)=psinv
call splev(tfp,nsf,cfp,3,rrs,ffspl,1,ier)
fpolv=ffspl(1)