##// END OF EJS Templates
schain
RSS
add full_profile_profile.pyf
rflores -
r1609:df12699ff15e
parent child
Show More
Side by Side Unified
Context file:
r1609:df12699ff15e Loading diff...:
Show file before | Show file after | Hide comments
This diff has been collapsed as it changes many lines, (520 lines changed) Show them Hide them
@@ -0,0 +1,520
1 ! -*- f90 -*-
2 ! Note: the context of this file is case sensitive.
3
4 python module fdjac2__user__routines
5 interface fdjac2_user_interface
6 subroutine fcn(m,n,x,wa,iflag) ! in :full_profile_profile:lmdif1.f:fdjac2:unknown_interface
7 integer, optional,check(len(wa)>=m),depend(wa) :: m=len(wa)
8 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
9 real dimension(n) :: x
10 real dimension(m) :: wa
11 integer :: iflag
12 end subroutine fcn
13 end interface fdjac2_user_interface
14 end python module fdjac2__user__routines
15 python module lmdif__user__routines
16 interface lmdif_user_interface
17 subroutine fcn(m,n,x,fvec,iflag) ! in :full_profile_profile:lmdif1.f:lmdif:unknown_interface
18 integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
19 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
20 real dimension(n) :: x
21 real dimension(m) :: fvec
22 integer :: iflag
23 end subroutine fcn
24 end interface lmdif_user_interface
25 end python module lmdif__user__routines
26 python module full_profile_profile ! in
27 interface ! in :full_profile_profile
28 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
29 complex dimension(4,nhts,ibits) :: acf_sum
30 real dimension(nhts,ibits),depend(nhts,ibits) :: acf_err
31 real dimension(nhts),depend(nhts) :: power
32 real dimension(nhts),depend(nhts) :: en
33 real dimension(ibits),depend(ibits) :: alag
34 real dimension(nhts),depend(nhts) :: thb2
35 real dimension(nhts),depend(nhts) :: bfm2
36 real dimension(nacf),depend(nacf) :: ote
37 real dimension(nacf),depend(nacf) :: ete
38 real dimension(nacf),depend(nacf) :: oti
39 real dimension(nacf),depend(nacf) :: eti
40 real dimension(nacf),depend(nacf) :: oph
41 real dimension(nacf),depend(nacf) :: eph
42 real dimension(nacf),depend(nacf) :: ophe
43 real dimension(nacf),depend(nacf) :: ephe
44 real dimension(nhts),depend(nhts) :: range2
45 real :: ut
46 integer, optional,check(shape(acf_sum,1)==nhts),depend(acf_sum) :: nhts=shape(acf_sum,1)
47 integer :: nacf
48 integer, optional,check(shape(acf_sum,2)==ibits),depend(acf_sum) :: ibits=shape(acf_sum,2)
49 real dimension(nhts,ibits),depend(nhts,ibits), :: acf_avg_real
50 real dimension(1) :: status
51 real :: chi2
52 real dimension(91) :: densp
53 real dimension(91) :: tep
54 real dimension(91) :: trp
55 real dimension(91) :: tip
56 real dimension(91) :: hfp
57 real dimension(91) :: hefp
58 real dimension(91) :: altp
59 real :: r0
60 real :: dr
61 real :: wl
62 real dimension(16,91) :: plag
63 real dimension(16,91) :: plag_errors
64 real :: sconst
65 real dimension(91) :: edensp
66 real dimension(91) :: etep
67 real dimension(91) :: etip
68 real dimension(91) :: ehfp
69 real dimension(91) :: ehefp
70 real dimension(85) :: bfld_prof
71 real dimension(85) :: alpha_prof
72 integer :: imode
73 real :: uttime
74 real dimension(34) :: ta
75 real dimension(30,5) :: bcoef
76 real :: te
77 real dimension(10) :: ti
78 real dimension(10) :: fi
79 real :: ven
80 real dimension(10) :: vin
81 real :: alpha
82 real :: dens
83 real :: bfld
84 integer :: nion
85 integer dimension(10) :: wi
86 real :: ak
87 common /chisq/ chi2
88 common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
89 common /data/ plag,plag_errors
90 common /sys/ sconst
91 common /errs/ edensp,etep,etip,ehfp,ehefp
92 common /mag/ bfld_prof,alpha_prof
93 common /mode/ imode
94 common /utime/ uttime
95 common /spline/ ta,bcoef
96 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
97 end subroutine profile
98 subroutine grid ! in :full_profile_profile:full_profile_profile.f
99 real dimension(34) :: ta
100 real dimension(30,5) :: bcoef
101 real dimension(91) :: densp
102 real dimension(91) :: tep
103 real dimension(91) :: trp
104 real dimension(91) :: tip
105 real dimension(91) :: hfp
106 real dimension(91) :: hefp
107 real dimension(91) :: altp
108 real :: r0
109 real :: dr
110 real :: wl
111 real dimension(16,91) :: plag
112 real dimension(16,91) :: plag_errors
113 real :: sconst
114 real dimension(91) :: edensp
115 real dimension(91) :: etep
116 real dimension(91) :: etip
117 real dimension(91) :: ehfp
118 real dimension(91) :: ehefp
119 real :: te
120 real dimension(10) :: ti
121 real dimension(10) :: fi
122 real :: ven
123 real dimension(10) :: vin
124 real :: alpha
125 real :: dens
126 real :: bfld
127 integer :: nion
128 integer dimension(10) :: wi
129 real :: ak
130 real :: chi2
131 real :: uttime
132 common /spline/ ta,bcoef
133 common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
134 common /data/ plag,plag_errors
135 common /sys/ sconst
136 common /errs/ edensp,etep,etip,ehfp,ehefp
137 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
138 common /chisq/ chi2
139 common /utime/ uttime
140 end subroutine grid
141 subroutine propagate(xe) ! in :full_profile_profile:full_profile_profile.f
142 real dimension(150) :: xe
143 real dimension(91) :: densp
144 real dimension(91) :: tep
145 real dimension(91) :: trp
146 real dimension(91) :: tip
147 real dimension(91) :: hfp
148 real dimension(91) :: hefp
149 real dimension(91) :: altp
150 real :: r0
151 real :: dr
152 real :: wl
153 real dimension(16,91) :: plag
154 real dimension(16,91) :: plag_errors
155 real :: sconst
156 real dimension(91) :: edensp
157 real dimension(91) :: etep
158 real dimension(91) :: etip
159 real dimension(91) :: ehfp
160 real dimension(91) :: ehefp
161 real dimension(34) :: ta
162 real dimension(30,5) :: bcoef
163 common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
164 common /data/ plag,plag_errors
165 common /sys/ sconst
166 common /errs/ edensp,etep,etip,ehfp,ehefp
167 common /spline/ ta,bcoef
168 end subroutine propagate
169 subroutine fcn_lpreg(m,n,x,fvec,iflag) ! in :full_profile_profile:full_profile_profile.f
170 integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
171 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
172 real dimension(n) :: x
173 real dimension(m) :: fvec
174 integer :: iflag
175 real dimension(34) :: ta
176 real dimension(30,5) :: bcoef
177 real dimension(91) :: densp
178 real dimension(91) :: tep
179 real dimension(91) :: trp
180 real dimension(91) :: tip
181 real dimension(91) :: hfp
182 real dimension(91) :: hefp
183 real dimension(91) :: altp
184 real :: r0
185 real :: dr
186 real :: wl
187 real dimension(16,91) :: plag
188 real dimension(16,91) :: plag_errors
189 real :: sconst
190 real dimension(91) :: edensp
191 real dimension(91) :: etep
192 real dimension(91) :: etip
193 real dimension(91) :: ehfp
194 real dimension(91) :: ehefp
195 real :: te
196 real dimension(10) :: ti
197 real dimension(10) :: fi
198 real :: ven
199 real dimension(10) :: vin
200 real :: alpha
201 real :: dens
202 real :: bfld
203 integer :: nion
204 integer dimension(10) :: wi
205 real :: ak
206 real :: chi2
207 real :: uttime
208 common /spline/ ta,bcoef
209 common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
210 common /data/ plag,plag_errors
211 common /sys/ sconst
212 common /errs/ edensp,etep,etip,ehfp,ehefp
213 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
214 common /chisq/ chi2
215 common /utime/ uttime
216 end subroutine fcn_lpreg
217 subroutine get_scale(plag2) ! in :full_profile_profile:full_profile_profile.f
218 real dimension(16,75) :: plag2
219 real dimension(91) :: densp
220 real dimension(91) :: tep
221 real dimension(91) :: trp
222 real dimension(91) :: tip
223 real dimension(91) :: hfp
224 real dimension(91) :: hefp
225 real dimension(91) :: altp
226 real :: r0
227 real :: dr
228 real :: wl
229 real dimension(16,91) :: plag
230 real dimension(16,91) :: plag_errors
231 real :: sconst
232 real dimension(91) :: edensp
233 real dimension(91) :: etep
234 real dimension(91) :: etip
235 real dimension(91) :: ehfp
236 real dimension(91) :: ehefp
237 common /fpa/ densp,tep,trp,tip,hfp,hefp,altp,r0,dr,wl
238 common /data/ plag,plag_errors
239 common /sys/ sconst
240 common /errs/ edensp,etep,etip,ehfp,ehefp
241 end subroutine get_scale
242 function enorm(n,x) ! in :full_profile_profile:lmdif1.f
243 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
244 real dimension(n) :: x
245 real :: enorm
246 end function enorm
247 subroutine fdjac2(fcn,m,n,x,fvec,fjac,ldfjac,iflag,epsfcn,wa) ! in :full_profile_profile:lmdif1.f
248 use fdjac2__user__routines
249 external fcn
250 integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
251 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
252 real dimension(n) :: x
253 real dimension(m) :: fvec
254 real dimension(ldfjac,n),depend(n) :: fjac
255 integer, optional,check(shape(fjac,0)==ldfjac),depend(fjac) :: ldfjac=shape(fjac,0)
256 integer :: iflag
257 real :: epsfcn
258 real dimension(m),depend(m) :: wa
259 end subroutine fdjac2
260 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
261 use lmdif__user__routines
262 external fcn
263 integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
264 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
265 real dimension(n) :: x
266 real dimension(m) :: fvec
267 real :: ftol
268 real :: xtol
269 real :: gtol
270 integer :: maxfev
271 real :: epsfcn
272 real dimension(n),depend(n) :: diag
273 integer :: mode
274 real :: factor
275 integer :: nprint
276 integer :: info
277 integer :: nfev
278 real dimension(ldfjac,n),depend(n) :: fjac
279 integer, optional,check(shape(fjac,0)==ldfjac),depend(fjac) :: ldfjac=shape(fjac,0)
280 integer dimension(n),depend(n) :: ipvt
281 real dimension(n),depend(n) :: qtf
282 real dimension(n),depend(n) :: wa1
283 real dimension(n),depend(n) :: wa2
284 real dimension(n),depend(n) :: wa3
285 real dimension(m),depend(m) :: wa4
286 end subroutine lmdif
287 subroutine lmpar(n,r,ldr,ipvt,diag,qtb,delta,par,x,sdiag,wa1,wa2) ! in :full_profile_profile:lmdif1.f
288 integer, optional,check(shape(r,1)==n),depend(r) :: n=shape(r,1)
289 real dimension(ldr,n) :: r
290 integer, optional,check(shape(r,0)==ldr),depend(r) :: ldr=shape(r,0)
291 integer dimension(n),depend(n) :: ipvt
292 real dimension(n),depend(n) :: diag
293 real dimension(n),depend(n) :: qtb
294 real :: delta
295 real :: par
296 real dimension(n),depend(n) :: x
297 real dimension(n),depend(n) :: sdiag
298 real dimension(n),depend(n) :: wa1
299 real dimension(n),depend(n) :: wa2
300 end subroutine lmpar
301 subroutine qrfac(m,n,a,lda,pivot,ipvt,lipvt,rdiag,acnorm,wa) ! in :full_profile_profile:lmdif1.f
302 integer :: m
303 integer, optional,check(shape(a,1)==n),depend(a) :: n=shape(a,1)
304 real dimension(lda,n) :: a
305 integer, optional,check(shape(a,0)==lda),depend(a) :: lda=shape(a,0)
306 logical :: pivot
307 integer dimension(lipvt) :: ipvt
308 integer, optional,check(len(ipvt)>=lipvt),depend(ipvt) :: lipvt=len(ipvt)
309 real dimension(n),depend(n) :: rdiag
310 real dimension(n),depend(n) :: acnorm
311 real dimension(n),depend(n) :: wa
312 end subroutine qrfac
313 subroutine qrsolv(n,r,ldr,ipvt,diag,qtb,x,sdiag,wa) ! in :full_profile_profile:lmdif1.f
314 integer, optional,check(shape(r,1)==n),depend(r) :: n=shape(r,1)
315 real dimension(ldr,n) :: r
316 integer, optional,check(shape(r,0)==ldr),depend(r) :: ldr=shape(r,0)
317 integer dimension(n),depend(n) :: ipvt
318 real dimension(n),depend(n) :: diag
319 real dimension(n),depend(n) :: qtb
320 real dimension(n),depend(n) :: x
321 real dimension(n),depend(n) :: sdiag
322 real dimension(n),depend(n) :: wa
323 end subroutine qrsolv
324 function spmpar(i) ! in :full_profile_profile:lmdif1.f
325 integer :: i
326 real :: spmpar
327 end function spmpar
328 function cdtr1(depth) ! in :full_profile_profile:fitacf.f
329 real :: depth
330 real :: cdtr1
331 end function cdtr1
332 function czte1(zlag,tr) ! in :full_profile_profile:fitacf.f
333 real :: zlag
334 real :: tr
335 real :: czte1
336 end function czte1
337 subroutine fcn(m,n,x,fvec,iflag) ! in :full_profile_profile:fitacf.f
338 integer, optional,check(len(fvec)>=m),depend(fvec) :: m=len(fvec)
339 integer, optional,check(len(x)>=n),depend(x) :: n=len(x)
340 real dimension(n) :: x
341 real dimension(m) :: fvec
342 integer :: iflag
343 real :: te
344 real dimension(10) :: ti
345 real dimension(10) :: fi
346 real :: ven
347 real dimension(10) :: vin
348 real :: alpha
349 real :: dens
350 real :: bfld
351 integer :: nion
352 integer dimension(10) :: wi
353 real :: ak
354 real dimension(100) :: tau
355 real dimension(100) :: rho
356 real dimension(100) :: sigma2
357 real dimension(10) :: params
358 integer dimension(10) :: ifit
359 real :: chisq
360 real dimension(10000) :: ev
361 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
362 common /fitter/ tau,rho,sigma2,params,ifit
363 common /errs/ chisq
364 common /trans/ ev
365 end subroutine fcn
366 function cj_ion(theta,psi) ! in :full_profile_profile:fitacf.f
367 real :: theta
368 real :: psi
369 complex :: cj_ion
370 end function cj_ion
371 function cj_electron(theta,phi,psi,alpha) ! in :full_profile_profile:fitacf.f
372 real :: theta
373 real :: phi
374 real :: psi
375 real :: alpha
376 integer :: imode
377 complex :: cj_electron
378 common /mode/ imode
379 end function cj_electron
380 function y_ion(theta,psi) ! in :full_profile_profile:fitacf.f
381 real :: theta
382 real :: psi
383 complex :: y_ion
384 end function y_ion
385 function y_electron(theta,phi,psi,alpha) ! in :full_profile_profile:fitacf.f
386 real :: theta
387 real :: phi
388 real :: psi
389 real :: alpha
390 complex :: y_electron
391 end function y_electron
392 function spect1(omega) ! in :full_profile_profile:fitacf.f
393 real :: omega
394 real :: te
395 real dimension(10) :: ti
396 real dimension(10) :: fi
397 real :: ven
398 real dimension(10) :: vin
399 real :: alpha
400 real :: dens
401 real :: bfld
402 integer :: nion
403 integer dimension(10) :: wi
404 real :: ak
405 integer :: imode
406 real :: spect1
407 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
408 common /mode/ imode
409 end function spect1
410 subroutine acf2(wl,tau,te1,ti1,fi1,ven1,vin1,wi1,nion1,alpha1,dens1,bfld1,acf) ! in :full_profile_profile:fitacf.f
411 real :: wl
412 real :: tau
413 real :: te1
414 real dimension(nion1) :: ti1
415 real dimension(nion1),depend(nion1) :: fi1
416 real :: ven1
417 real dimension(nion1),depend(nion1) :: vin1
418 integer dimension(nion1),depend(nion1) :: wi1
419 integer, optional,check(len(ti1)>=nion1),depend(ti1) :: nion1=len(ti1)
420 real :: alpha1
421 real :: dens1
422 real :: bfld1
423 real :: acf
424 real :: te
425 real dimension(10) :: ti
426 real dimension(10) :: fi
427 real :: ven
428 real dimension(10) :: vin
429 real :: alpha
430 real :: dens
431 real :: bfld
432 integer :: nion
433 integer dimension(10) :: wi
434 real :: ak
435 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
436 end subroutine acf2
437 subroutine gaussq(tau,acf) ! in :full_profile_profile:fitacf.f
438 real :: tau
439 real :: acf
440 real :: te
441 real dimension(10) :: ti
442 real dimension(10) :: fi
443 real :: ven
444 real dimension(10) :: vin
445 real :: alpha
446 real :: dens
447 real :: bfld
448 integer :: nion
449 integer dimension(10) :: wi
450 real :: ak
451 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
452 end subroutine gaussq
453 subroutine lagp(plag,wl,r0,dr,nl,nrange) ! in :full_profile_profile:lagp.f
454 real dimension(nl,nrange) :: plag
455 real :: wl
456 real :: r0
457 real :: dr
458 integer, optional,check(shape(plag,0)==nl),depend(plag) :: nl=shape(plag,0)
459 integer, optional,check(shape(plag,1)==nrange),depend(plag) :: nrange=shape(plag,1)
460 real :: te
461 real dimension(10) :: ti
462 real dimension(10) :: fi
463 real :: ven
464 real dimension(10) :: vin
465 real :: alpha
466 real :: dens
467 real :: bfld
468 integer :: nion
469 integer dimension(10) :: wi
470 real :: ak
471 real dimension(34) :: ta
472 real dimension(30,5) :: bcoef
473 real dimension(85) :: bfld_prof
474 real dimension(85) :: alpha_prof
475 common /spec/ te,ti,fi,ven,vin,alpha,dens,bfld,nion,wi,ak
476 common /spline/ ta,bcoef
477 common /mag/ bfld_prof,alpha_prof
478 end subroutine lagp
479 function atanh(x) ! in :full_profile_profile:lagp.f
480 real :: x
481 real :: atanh
482 end function atanh
483 subroutine get_spline(alt,dens,te,ti,hf,hef) ! in :full_profile_profile:lagp.f
484 real :: alt
485 real :: dens
486 real :: te
487 real :: ti
488 real :: hf
489 real :: hef
490 real dimension(34) :: ta
491 real dimension(30,5) :: bcoef
492 common /spline/ ta,bcoef
493 end subroutine get_spline
494 function bvalue(t,bcoef,n,k,x,jderiv) ! in :full_profile_profile:lagp.f
495 real dimension(n+k),depend(n,k) :: t
496 real dimension(n) :: bcoef
497 integer, optional,check(len(bcoef)>=n),depend(bcoef) :: n=len(bcoef)
498 integer :: k
499 real :: x
500 integer :: jderiv
501 real :: bvalue
502 end function bvalue
503 subroutine interv(xt,lxt,x,left,mflag) ! in :full_profile_profile:lagp.f
504 real dimension(lxt) :: xt
505 integer, optional,check(len(xt)>=lxt),depend(xt) :: lxt=len(xt)
506 real :: x
507 integer :: left
508 integer :: mflag
509 end subroutine interv
510 function r1mach(i) ! in :full_profile_profile:r1mach.f
511 integer :: i
512 integer :: cray1
513 real :: r1mach
514 common /d8mach/ cray1
515 end function r1mach
516 end interface
517 end python module full_profile_profile
518
519 ! This file was auto-generated with f2py (version:2).
520 ! See http://cens.ioc.ee/projects/f2py2e/
General Comments 0
You need to be logged in to leave comments. Login now