full_profile_profile.pyf
520 lines
| 20.9 KiB
| text/plain
|
TextLexer
r1609 | ! -*- 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/ | ||||