-
Notifications
You must be signed in to change notification settings - Fork 14
/
config_ext.f90
50 lines (40 loc) · 1.9 KB
/
config_ext.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
!#######################################################################
MODULE config_ext
!########################################################################
IMPLICIT NONE
CONTAINS
!#######################################################################
FUNCTION Conf_p_avrg(S,x) RESULT(pav)
!========================================================================
! Function returns the value of <p>, i.e. the flux averaged pitch
! <p>=<sqrt(1-x*b)>,
! where p=vpar/v, 0<=x<=1, x = (1-p*p)/b, and b = B(s,theta,phi)/Bmax(s)
!========================================================================
USE const_and_precisions
IMPLICIT NONE
REAL(wp_), INTENT(in) :: S,x
REAL(wp_) :: pav
!========================================================================
! pav = ?
!========================================================================
END FUNCTION Conf_p_avrg
!########################################################################
!########################################################################
!########################################################################
FUNCTION Conf_p_avrg_integr(S,x) RESULT(ipav)
!========================================================================
! Function returns the value of Integral(from x to 1) dx'/p_avrg(s,x') ,
! where p_avrg(s,x) = <sqrt(1-x*B(s,th,ph)/Bmax(s))> 0<=x<=1, 0<=s<=1,
! s is the normalized magnetic flux
!========================================================================
USE const_and_precisions
IMPLICIT NONE
REAL(wp_), INTENT(in) :: S,x
REAL(wp_) :: ipav
!========================================================================
! ipav = ?
!========================================================================
END FUNCTION Conf_p_avrg_integr
!#######################################################################
END MODULE config_ext
!#######################################################################