-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathIons_param_Sphe.f90
46 lines (46 loc) · 1.85 KB
/
Ions_param_Sphe.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
!***********************************************************************
! *
SUBROUTINE Ions_param_Sphe &
(ion_number,q_ion,R_ion,GTheta_ion,Phy_ion)
! *
! Written by G. Gaigalas and D. Kato *
! NIFS, Japan, September 2013 *
! The last modification made by G. Gaigalas June 2018 *
! *
!***********************************************************************
!-----------------------------------------------
! M o d u l e s
!-----------------------------------------------
USE vast_kind_param, ONLY: DOUBLE
USE CONS_C, ONLY: ONE
USE CrysPAR_C, ONLY: QQ, Vector_L, X1, X2, X3
IMPLICIT NONE
!-----------------------------------------------
! D u m m y A r g u m e n t s
!-----------------------------------------------
INTEGER, INTENT(IN) :: ion_number
REAL(DOUBLE), INTENT(OUT) :: q_ion, R_ion, GTheta_ion, Phy_ion
!-----------------------------------------------
! L o c a l V a r i a b l e s
!-----------------------------------------------
REAL(DOUBLE) :: PI
!-----------------------------------------------
!
PI = 4.0D00*ATAN (ONE)
!
q_ion = QQ(ion_number)
R_ion = Vector_L*X1(ion_number)
GTheta_ion = X2(ion_number)
Phy_ion = X3(ion_number)
!
! Pervedame angstremus i atominius vienetus
!
R_ion = R_ion/0.52917706
!
! Pervedame laipsnius i radianus
!
GTheta_ion = GTheta_ion*PI / 180.D00
Phy_ion = Phy_ion*PI / 180.D00
!
RETURN
END SUBROUTINE Ions_param_Sphe