|
|
! -*- 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/
|
|
|
|
