From 61c97e33577b8b9d6b40c5d1e7e601e127dc467f Mon Sep 17 00:00:00 2001 From: Lorenzo Figini Date: Fri, 2 Oct 2015 15:09:10 +0000 Subject: [PATCH] corrected bug in equinum_fpol: vacuum value was returned for psin=0 --- src/gray.f | 2 +- src/interp_eqprof.f90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/gray.f b/src/gray.f index cc31738..0144309 100644 --- a/src/gray.f +++ b/src/gray.f @@ -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 diff --git a/src/interp_eqprof.f90 b/src/interp_eqprof.f90 index e681cae..01733c9 100644 --- a/src/interp_eqprof.f90 +++ b/src/interp_eqprof.f90 @@ -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)