[Fortran] 纯文本查看 复制代码
G(17)=45. * 3.14 / 180. - delt1
G(20)=0.85 - TAN(delt1) / ( TAN(beta) + TAN(delt1) )
PI_1=sin(alpha - theta) + f2 * sin(delt2) * cos(alpha - theta)
PI_2=cos(alpha - theta) - f2 * sin(delt2) * sin(alpha - theta) * cos(beta)
A1=PI_2 * COTAN(alpha) + PI_1 * cos(beta)
B1=( f2 * COTAN(alpha) * ( cos(delt2) + sin(delt2) * sin(alpha - theta) * sin(beta) ) - PI_1 * sin(beta) )
C1=( 1 / f1 ) * ( f2 * cos(beta) * ( cos(delt2) + sin(delt2) * sin(alpha - theta) * sin(beta) ) - PI_2 * sin(beta) )
delt1=ASIN( A1 / SQRT( A1**2 + B1**2) ) - Beta - ASIN( C1 / SQRT( A1**2 + B1**2) )
PRINT*,'delt1 = ', delt1
delt2=ATAN( sin(beta) * sin(alpha) / cos(theta) )
PRINT*,'delt2 = ', delt2