1 |
guez |
98 |
module ord_coordm_m |
2 |
|
|
|
3 |
|
|
implicit none |
4 |
|
|
|
5 |
|
|
contains |
6 |
|
|
|
7 |
|
|
|
8 |
|
|
!*********************************** |
9 |
|
|
|
10 |
|
|
function ord_coordm(xi) |
11 |
|
|
|
12 |
|
|
! From dyn3d/ord_coordm.F, version 1.1.1.1 2004/05/19 12:53:06 |
13 |
|
|
! Author : P. Le Van |
14 |
|
|
|
15 |
|
|
! This procedure converts to degrees, if necessary, and inverts the |
16 |
|
|
! order. |
17 |
|
|
|
18 |
|
|
use nr_util, only: pi |
19 |
|
|
|
20 |
|
|
|
21 |
|
|
REAL, intent(in):: xi(:) ! angle, in rad or degrees |
22 |
|
|
REAL ord_coordm(size(xi)) ! angle, in degrees |
23 |
|
|
|
24 |
|
|
!----------------------------- |
25 |
|
|
|
26 |
|
|
IF (xi(1) < 6.5) THEN |
27 |
|
|
! "xi" is in rad |
28 |
|
|
ord_coordm(:) = xi(size(xi):1:-1) * 180. / pi |
29 |
|
|
else |
30 |
|
|
! "xi" is in degrees |
31 |
|
|
ord_coordm(:) = xi(size(xi):1:-1) |
32 |
|
|
ENDIF |
33 |
|
|
|
34 |
|
|
END function ord_coordm |
35 |
|
|
|
36 |
|
|
end module ord_coordm_m |