diff --git a/schainf/Ffiles/full_profile_profile.pyf b/schainf/Ffiles/full_profile_profile.pyf
new file mode 100644
index 0000000..ed0d093
--- /dev/null
+++ b/schainf/Ffiles/full_profile_profile.pyf
@@ -0,0 +1,520 @@
+!    -*- f90 -*-
+! Note: the context of this file is case sensitive.
+
+python module fdjac2__user__routines 
+    interface fdjac2_user_interface 
+        subroutine fcn(m,n,x,wa,iflag) ! in :full_profile_profile:lmdif1.f:fdjac2:unknown_interface
+            integer, optional,check(len(wa)>=m),depend(wa) :: m=len(wa)
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real dimension(m) :: wa
+            integer :: iflag
+        end subroutine fcn
+    end interface fdjac2_user_interface
+end python module fdjac2__user__routines
+python module lmdif__user__routines 
+    interface lmdif_user_interface 
+        subroutine fcn(m,n,x,fvec,iflag) ! in :full_profile_profile:lmdif1.f:lmdif:unknown_interface
+            integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real dimension(m) :: fvec
+            integer :: iflag
+        end subroutine fcn
+    end interface lmdif_user_interface
+end python module lmdif__user__routines
+python module full_profile_profile ! in 
+    interface  ! in :full_profile_profile
+        subroutine profile(acf_sum,acf_err,power,en,alag,thb2,bfm2,ote,ete,oti,eti,oph,eph,ophe,ephe,range2,ut,nhts,nacf,ibits,acf_avg_real,status) ! in :full_profile_profile:full_profile_profile.f
+            complex dimension(4,nhts,ibits) :: acf_sum
+            real dimension(nhts,ibits),depend(nhts,ibits) :: acf_err
+            real dimension(nhts),depend(nhts) :: power
+            real dimension(nhts),depend(nhts) :: en
+            real dimension(ibits),depend(ibits) :: alag
+            real dimension(nhts),depend(nhts) :: thb2
+            real dimension(nhts),depend(nhts) :: bfm2
+            real dimension(nacf),depend(nacf) :: ote
+            real dimension(nacf),depend(nacf) :: ete
+            real dimension(nacf),depend(nacf) :: oti
+            real dimension(nacf),depend(nacf) :: eti
+            real dimension(nacf),depend(nacf) :: oph
+            real dimension(nacf),depend(nacf) :: eph
+            real dimension(nacf),depend(nacf) :: ophe
+            real dimension(nacf),depend(nacf) :: ephe
+            real dimension(nhts),depend(nhts) :: range2
+            real :: ut
+            integer, optional,check(shape(acf_sum,1)==nhts),depend(acf_sum) :: nhts=shape(acf_sum,1)
+            integer :: nacf
+            integer, optional,check(shape(acf_sum,2)==ibits),depend(acf_sum) :: ibits=shape(acf_sum,2)
+            real dimension(nhts,ibits),depend(nhts,ibits), :: acf_avg_real	    
+            real dimension(1) :: status
+            real :: chi2
+            real dimension(91) :: densp
+            real dimension(91) :: tep
+            real dimension(91) :: trp
+            real dimension(91) :: tip
+            real dimension(91) :: hfp
+            real dimension(91) :: hefp
+            real dimension(91) :: altp
+            real :: r0
+            real :: dr
+            real :: wl
+            real dimension(16,91) :: plag
+            real dimension(16,91) :: plag_errors
+            real :: sconst
+            real dimension(91) :: edensp
+            real dimension(91) :: etep
+            real dimension(91) :: etip
+            real dimension(91) :: ehfp
+            real dimension(91) :: ehefp
+            real dimension(85) :: bfld_prof
+            real dimension(85) :: alpha_prof
+            integer :: imode
+            real :: uttime
+            real dimension(34) :: ta
+            real dimension(30,5) :: bcoef
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            common /chisq/ chi2
+            common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
+            common /data/ plag,plag_errors
+            common /sys/ sconst
+            common /errs/ edensp,etep,etip,ehfp,ehefp
+            common /mag/ bfld_prof,alpha_prof
+            common /mode/ imode
+            common /utime/ uttime
+            common /spline/ ta,bcoef
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+        end subroutine profile
+        subroutine grid ! in :full_profile_profile:full_profile_profile.f
+            real dimension(34) :: ta
+            real dimension(30,5) :: bcoef
+            real dimension(91) :: densp
+            real dimension(91) :: tep
+            real dimension(91) :: trp
+            real dimension(91) :: tip
+            real dimension(91) :: hfp
+            real dimension(91) :: hefp
+            real dimension(91) :: altp
+            real :: r0
+            real :: dr
+            real :: wl
+            real dimension(16,91) :: plag
+            real dimension(16,91) :: plag_errors
+            real :: sconst
+            real dimension(91) :: edensp
+            real dimension(91) :: etep
+            real dimension(91) :: etip
+            real dimension(91) :: ehfp
+            real dimension(91) :: ehefp
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            real :: chi2
+            real :: uttime
+            common /spline/ ta,bcoef
+            common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
+            common /data/ plag,plag_errors
+            common /sys/ sconst
+            common /errs/ edensp,etep,etip,ehfp,ehefp
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+            common /chisq/ chi2
+            common /utime/ uttime
+        end subroutine grid
+        subroutine propagate(xe) ! in :full_profile_profile:full_profile_profile.f
+            real dimension(150) :: xe
+            real dimension(91) :: densp
+            real dimension(91) :: tep
+            real dimension(91) :: trp
+            real dimension(91) :: tip
+            real dimension(91) :: hfp
+            real dimension(91) :: hefp
+            real dimension(91) :: altp
+            real :: r0
+            real :: dr
+            real :: wl
+            real dimension(16,91) :: plag
+            real dimension(16,91) :: plag_errors
+            real :: sconst
+            real dimension(91) :: edensp
+            real dimension(91) :: etep
+            real dimension(91) :: etip
+            real dimension(91) :: ehfp
+            real dimension(91) :: ehefp
+            real dimension(34) :: ta
+            real dimension(30,5) :: bcoef
+            common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
+            common /data/ plag,plag_errors
+            common /sys/ sconst
+            common /errs/ edensp,etep,etip,ehfp,ehefp
+            common /spline/ ta,bcoef
+        end subroutine propagate
+        subroutine fcn_lpreg(m,n,x,fvec,iflag) ! in :full_profile_profile:full_profile_profile.f
+            integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real dimension(m) :: fvec
+            integer :: iflag
+            real dimension(34) :: ta
+            real dimension(30,5) :: bcoef
+            real dimension(91) :: densp
+            real dimension(91) :: tep
+            real dimension(91) :: trp
+            real dimension(91) :: tip
+            real dimension(91) :: hfp
+            real dimension(91) :: hefp
+            real dimension(91) :: altp
+            real :: r0
+            real :: dr
+            real :: wl
+            real dimension(16,91) :: plag
+            real dimension(16,91) :: plag_errors
+            real :: sconst
+            real dimension(91) :: edensp
+            real dimension(91) :: etep
+            real dimension(91) :: etip
+            real dimension(91) :: ehfp
+            real dimension(91) :: ehefp
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            real :: chi2
+            real :: uttime
+            common /spline/ ta,bcoef
+            common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
+            common /data/ plag,plag_errors
+            common /sys/ sconst
+            common /errs/ edensp,etep,etip,ehfp,ehefp
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+            common /chisq/ chi2
+            common /utime/ uttime
+        end subroutine fcn_lpreg
+        subroutine get_scale(plag2) ! in :full_profile_profile:full_profile_profile.f
+            real dimension(16,75) :: plag2
+            real dimension(91) :: densp
+            real dimension(91) :: tep
+            real dimension(91) :: trp
+            real dimension(91) :: tip
+            real dimension(91) :: hfp
+            real dimension(91) :: hefp
+            real dimension(91) :: altp
+            real :: r0
+            real :: dr
+            real :: wl
+            real dimension(16,91) :: plag
+            real dimension(16,91) :: plag_errors
+            real :: sconst
+            real dimension(91) :: edensp
+            real dimension(91) :: etep
+            real dimension(91) :: etip
+            real dimension(91) :: ehfp
+            real dimension(91) :: ehefp
+            common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
+            common /data/ plag,plag_errors
+            common /sys/ sconst
+            common /errs/ edensp,etep,etip,ehfp,ehefp
+        end subroutine get_scale
+        function enorm(n,x) ! in :full_profile_profile:lmdif1.f
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real :: enorm
+        end function enorm
+        subroutine fdjac2(fcn,m,n,x,fvec,fjac,ldfjac,iflag,epsfcn,wa) ! in :full_profile_profile:lmdif1.f
+            use fdjac2__user__routines
+            external fcn
+            integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real dimension(m) :: fvec
+            real dimension(ldfjac,n),depend(n) :: fjac
+            integer, optional,check(shape(fjac,0)==ldfjac),depend(fjac) :: ldfjac=shape(fjac,0)
+            integer :: iflag
+            real :: epsfcn
+            real dimension(m),depend(m) :: wa
+        end subroutine fdjac2
+        subroutine lmdif(fcn,m,n,x,fvec,ftol,xtol,gtol,maxfev,epsfcn,diag,mode,factor,nprint,info,nfev,fjac,ldfjac,ipvt,qtf,wa1,wa2,wa3,wa4) ! in :full_profile_profile:lmdif1.f
+            use lmdif__user__routines
+            external fcn
+            integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real dimension(m) :: fvec
+            real :: ftol
+            real :: xtol
+            real :: gtol
+            integer :: maxfev
+            real :: epsfcn
+            real dimension(n),depend(n) :: diag
+            integer :: mode
+            real :: factor
+            integer :: nprint
+            integer :: info
+            integer :: nfev
+            real dimension(ldfjac,n),depend(n) :: fjac
+            integer, optional,check(shape(fjac,0)==ldfjac),depend(fjac) :: ldfjac=shape(fjac,0)
+            integer dimension(n),depend(n) :: ipvt
+            real dimension(n),depend(n) :: qtf
+            real dimension(n),depend(n) :: wa1
+            real dimension(n),depend(n) :: wa2
+            real dimension(n),depend(n) :: wa3
+            real dimension(m),depend(m) :: wa4
+        end subroutine lmdif
+        subroutine lmpar(n,r,ldr,ipvt,diag,qtb,delta,par,x,sdiag,wa1,wa2) ! in :full_profile_profile:lmdif1.f
+            integer, optional,check(shape(r,1)==n),depend(r) :: n=shape(r,1)
+            real dimension(ldr,n) :: r
+            integer, optional,check(shape(r,0)==ldr),depend(r) :: ldr=shape(r,0)
+            integer dimension(n),depend(n) :: ipvt
+            real dimension(n),depend(n) :: diag
+            real dimension(n),depend(n) :: qtb
+            real :: delta
+            real :: par
+            real dimension(n),depend(n) :: x
+            real dimension(n),depend(n) :: sdiag
+            real dimension(n),depend(n) :: wa1
+            real dimension(n),depend(n) :: wa2
+        end subroutine lmpar
+        subroutine qrfac(m,n,a,lda,pivot,ipvt,lipvt,rdiag,acnorm,wa) ! in :full_profile_profile:lmdif1.f
+            integer :: m
+            integer, optional,check(shape(a,1)==n),depend(a) :: n=shape(a,1)
+            real dimension(lda,n) :: a
+            integer, optional,check(shape(a,0)==lda),depend(a) :: lda=shape(a,0)
+            logical :: pivot
+            integer dimension(lipvt) :: ipvt
+            integer, optional,check(len(ipvt)>=lipvt),depend(ipvt) :: lipvt=len(ipvt)
+            real dimension(n),depend(n) :: rdiag
+            real dimension(n),depend(n) :: acnorm
+            real dimension(n),depend(n) :: wa
+        end subroutine qrfac
+        subroutine qrsolv(n,r,ldr,ipvt,diag,qtb,x,sdiag,wa) ! in :full_profile_profile:lmdif1.f
+            integer, optional,check(shape(r,1)==n),depend(r) :: n=shape(r,1)
+            real dimension(ldr,n) :: r
+            integer, optional,check(shape(r,0)==ldr),depend(r) :: ldr=shape(r,0)
+            integer dimension(n),depend(n) :: ipvt
+            real dimension(n),depend(n) :: diag
+            real dimension(n),depend(n) :: qtb
+            real dimension(n),depend(n) :: x
+            real dimension(n),depend(n) :: sdiag
+            real dimension(n),depend(n) :: wa
+        end subroutine qrsolv
+        function spmpar(i) ! in :full_profile_profile:lmdif1.f
+            integer :: i
+            real :: spmpar
+        end function spmpar
+        function cdtr1(depth) ! in :full_profile_profile:fitacf.f
+            real :: depth
+            real :: cdtr1
+        end function cdtr1
+        function czte1(zlag,tr) ! in :full_profile_profile:fitacf.f
+            real :: zlag
+            real :: tr
+            real :: czte1
+        end function czte1
+        subroutine fcn(m,n,x,fvec,iflag) ! in :full_profile_profile:fitacf.f
+            integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
+            integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
+            real dimension(n) :: x
+            real dimension(m) :: fvec
+            integer :: iflag
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            real dimension(100) :: tau
+            real dimension(100) :: rho
+            real dimension(100) :: sigma2
+            real dimension(10) :: params
+            integer dimension(10) :: ifit
+            real :: chisq
+            real dimension(10000) :: ev
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+            common /fitter/ tau,rho,sigma2,params,ifit
+            common /errs/ chisq
+            common /trans/ ev
+        end subroutine fcn
+        function cj_ion(theta,psi) ! in :full_profile_profile:fitacf.f
+            real :: theta
+            real :: psi
+            complex :: cj_ion
+        end function cj_ion
+        function cj_electron(theta,phi,psi,alpha) ! in :full_profile_profile:fitacf.f
+            real :: theta
+            real :: phi
+            real :: psi
+            real :: alpha
+            integer :: imode
+            complex :: cj_electron
+            common /mode/ imode
+        end function cj_electron
+        function y_ion(theta,psi) ! in :full_profile_profile:fitacf.f
+            real :: theta
+            real :: psi
+            complex :: y_ion
+        end function y_ion
+        function y_electron(theta,phi,psi,alpha) ! in :full_profile_profile:fitacf.f
+            real :: theta
+            real :: phi
+            real :: psi
+            real :: alpha
+            complex :: y_electron
+        end function y_electron
+        function spect1(omega) ! in :full_profile_profile:fitacf.f
+            real :: omega
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            integer :: imode
+            real :: spect1
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+            common /mode/ imode
+        end function spect1
+        subroutine acf2(wl,tau,te1,ti1,fi1,ven1,vin1,wi1,nion1,alpha1,dens1,bfld1,acf) ! in :full_profile_profile:fitacf.f
+            real :: wl
+            real :: tau
+            real :: te1
+            real dimension(nion1) :: ti1
+            real dimension(nion1),depend(nion1) :: fi1
+            real :: ven1
+            real dimension(nion1),depend(nion1) :: vin1
+            integer dimension(nion1),depend(nion1) :: wi1
+            integer, optional,check(len(ti1)>=nion1),depend(ti1) :: nion1=len(ti1)
+            real :: alpha1
+            real :: dens1
+            real :: bfld1
+            real :: acf
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+        end subroutine acf2
+        subroutine gaussq(tau,acf) ! in :full_profile_profile:fitacf.f
+            real :: tau
+            real :: acf
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+        end subroutine gaussq
+        subroutine lagp(plag,wl,r0,dr,nl,nrange) ! in :full_profile_profile:lagp.f
+            real dimension(nl,nrange) :: plag
+            real :: wl
+            real :: r0
+            real :: dr
+            integer, optional,check(shape(plag,0)==nl),depend(plag) :: nl=shape(plag,0)
+            integer, optional,check(shape(plag,1)==nrange),depend(plag) :: nrange=shape(plag,1)
+            real :: te
+            real dimension(10) :: ti
+            real dimension(10) :: fi
+            real :: ven
+            real dimension(10) :: vin
+            real :: alpha
+            real :: dens
+            real :: bfld
+            integer :: nion
+            integer dimension(10) :: wi
+            real :: ak
+            real dimension(34) :: ta
+            real dimension(30,5) :: bcoef
+            real dimension(85) :: bfld_prof
+            real dimension(85) :: alpha_prof
+            common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
+            common /spline/ ta,bcoef
+            common /mag/ bfld_prof,alpha_prof
+        end subroutine lagp
+        function atanh(x) ! in :full_profile_profile:lagp.f
+            real :: x
+            real :: atanh
+        end function atanh
+        subroutine get_spline(alt,dens,te,ti,hf,hef) ! in :full_profile_profile:lagp.f
+            real :: alt
+            real :: dens
+            real :: te
+            real :: ti
+            real :: hf
+            real :: hef
+            real dimension(34) :: ta
+            real dimension(30,5) :: bcoef
+            common /spline/ ta,bcoef
+        end subroutine get_spline
+        function bvalue(t,bcoef,n,k,x,jderiv) ! in :full_profile_profile:lagp.f
+            real dimension(n+k),depend(n,k) :: t
+            real dimension(n) :: bcoef
+            integer, optional,check(len(bcoef)>=n),depend(bcoef) :: n=len(bcoef)
+            integer :: k
+            real :: x
+            integer :: jderiv
+            real :: bvalue
+        end function bvalue
+        subroutine interv(xt,lxt,x,left,mflag) ! in :full_profile_profile:lagp.f
+            real dimension(lxt) :: xt
+            integer, optional,check(len(xt)>=lxt),depend(xt) :: lxt=len(xt)
+            real :: x
+            integer :: left
+            integer :: mflag
+        end subroutine interv
+        function r1mach(i) ! in :full_profile_profile:r1mach.f
+            integer :: i
+            integer :: cray1
+            real :: r1mach
+            common /d8mach/ cray1
+        end function r1mach
+    end interface 
+end python module full_profile_profile
+
+! This file was auto-generated with f2py (version:2).
+! See http://cens.ioc.ee/projects/f2py2e/