ADAS Subroutine a8optm
subroutine a8optm(itype,xa,oa,n,s,b0,bp0,f10,f20,f30,xk0,ifail)
c-----------------------------------------------------------------------
c
c ********************* fortran77 subroutine a8optm *******************
c
c purpose: to find the best approximate form parameters for neutral
c atoms by varing the matching position.
c
c calling program: adas108.for
c
c input:
c (i*4) itype = type of transition (1=dipole,2=non-dipole
c non-spin change, 3=spin change, 4=null)
c (r*8) xa() = x-parameters for cross-section
c (r*8) oa = collisions strengths for transition
c (i*4) n = no of collision strengths
c (r*8) s = line strength if dipole transition
c (i*4) ifail = failure code on entry (ifail=0 two point
c fit, ifail=-1 one point fit)
c output:
c (r*8) b0 = threshold form parameter
c (r*8) bp0 = matching parameter
c (r*8) f10 = threshold form parameter
c (r*8) f20 = asymptotic form parameter
c (r*8) f30 = asymptotic form parameter
c (r*8) xk0 = optimum matching x-value
c (i*4) ifail = failure code on exit
c (ifail=0 successful two point fit
c ifail=1 converted to one point fit)
c
c
c routines:
c a8slvf adas solves for asymptotic parms f2 and f3
c a8slv2 adas solves for the parms f1,f2,f3,b
c
c author: Hugh Summers, University of Strathclyde ext.4196
c
c
c version 1.1 date: 19/07/99
c modified: Hugh Summers
c - first release
c
c-----------------------------------------------------------------------
INTEGER IFAIL, ITYPE, N
REAL*8 B0, BP0, F10, F20
REAL*8 F30, OA(N), S, XA(N)
REAL*8 XK0