File: getdphi.f90

package info (click to toggle)
wsjtx 2.7.0%2Brepack-1
  • links: PTS, VCS
  • area: main
  • in suites: trixie
  • size: 70,440 kB
  • sloc: cpp: 75,379; f90: 46,460; python: 27,241; ansic: 13,367; fortran: 2,382; makefile: 197; sh: 133
file content (18 lines) | stat: -rwxr-xr-x 266 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
subroutine getdphi(qphi)

  real qphi(12)

  s=0.
  c=0.
  do i=1,12
     th=i*30/57.2957795
     s=s+qphi(i)*sin(th)
     c=c+qphi(i)*cos(th)
  enddo

  dphi=57.2957795*atan2(s,c)
  write(*,1010) nint(dphi)
1010 format('!Best-fit Dphi =',i4,' deg')

  return
  end