- Timestamp:
- 09/05/08 11:49:42 (16 years ago)
- Branches:
- master
- Children:
- fafe4d6
- Parents:
- 2ebb8b6
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
difang.f
r2ebb8b6 rbd2278d 1 c**************************************************************2 c 3 cThis file contains the subroutines: difang,addang4 c 5 cCopyright 2003 Frank Eisenmenger, U.H.E. Hansmann,6 cShura Hayryan, Chin-Ku7 cCopyright 2007 Frank Eisenmenger, U.H.E. Hansmann,8 cJan H. Meinke, Sandipan Mohanty9 c 10 c**************************************************************1 ! ************************************************************** 2 ! 3 ! This file contains the subroutines: difang,addang 4 ! 5 ! Copyright 2003 Frank Eisenmenger, U.H.E. Hansmann, 6 ! Shura Hayryan, Chin-Ku 7 ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann, 8 ! Jan H. Meinke, Sandipan Mohanty 9 ! 10 ! ************************************************************** 11 11 12 12 real*8 function difang(a1,a2) 13 13 14 c......................................................15 cPURPOSE: difang = a2 - a1 with: -pi < difang <= pi16 c17 cINPUT: a1,a2-two angles [rad.]18 c 19 cCALLS: none20 c 21 c......................................................14 ! ...................................................... 15 ! PURPOSE: difang = a2 - a1 with: -pi < difang <= pi 16 ! 17 ! INPUT: a1,a2-two angles [rad.] 18 ! 19 ! CALLS: none 20 ! 21 ! ...................................................... 22 22 23 23 implicit real*8 (a-h,o-z) 24 24 25 25 parameter (pi=3.141592653589793d0, 26 #pi2=2.d0*pi)26 & pi2=2.d0*pi) 27 27 28 28 d=mod((a2-a1),pi2) … … 35 35 return 36 36 end 37 c*********************************37 ! ********************************* 38 38 real*8 function addang(a1,a2) 39 39 40 c......................................................41 cPURPOSE: addang = a1 + a2 with: -pi < addang <= pi42 c43 cINPUT: a1,a2-two angles [rad.]44 c 45 cCALLS: none46 c 47 c......................................................40 ! ...................................................... 41 ! PURPOSE: addang = a1 + a2 with: -pi < addang <= pi 42 ! 43 ! INPUT: a1,a2-two angles [rad.] 44 ! 45 ! CALLS: none 46 ! 47 ! ...................................................... 48 48 49 49 implicit real*8 (a-h,o-z) 50 50 51 51 parameter (pi=3.141592653589793d0, 52 #pi2=2.d0*pi)52 & pi2=2.d0*pi) 53 53 54 54 d=mod((a1+a2),pi2)
Note:
See TracChangeset
for help on using the changeset viewer.