4 REAL(rprec) :: chi(nthet), mu0
5 INTEGER :: ipsi, jthet, j, i
11 chi(jthet)=2.*pi*(jthet-1)/(nthet-1)
19 arcrad(ipsi,jthet)=arcrad(ipsi-1,jthet)+
20 & sqrt((xs(ipsi,jthet)-xs(ipsi-1,jthet))**2
21 & +(zs(ipsi,jthet)-zs(ipsi-1,jthet))**2)
28 jacob(i,j)=arcsur(i,nthet)/(2.*pi*bps(i,j))
33 vprime(ipsi)=2*pi*trap(nthet,chi,jacob(ipsi,1:nthet))
36 volume(ipsi)=trap(ipsi,psiv(1:ipsi),vprime(1:ipsi))
38 vprime(1)=axisv(psiv,vprime,
SIZE(psiv))
41 tflx(ipsi)=trap(ipsi,psiv(1:ipsi),qsfin(1:ipsi))