WSJT-X/getdphi.f90
Joe Taylor c8e94ea0ee Finished(?) inserting mutex lockouts.
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@1299 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
2009-07-26 16:31:11 +00:00

21 lines
309 B
Fortran

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)
call cs_lock('getdphi')
write(11,1010) nint(dphi)
1010 format('Best-fit Dphi =',i4,' deg')
call cs_unlock
return
end