1 SUBROUTINE pout(T,X,LUN)
30 IMPLICIT DOUBLE PRECISION (a-h,o-z)
33 common/option/l,m,ires,isun,imoon,iephem,idrag,idens,isrp,iorb
35 common/timecmn/ti,tf,tr
36 common/pltcon/ge,
re,rate,pm,aj2,ellip,ratm
37 common/count/
node,nrev
39 DATA hts,dth,rtd,
pi,tpi,one/3.6d3,24.d0,57.29577951308232d0,
40 1 3.141592653589793d0,6.283185307179586d0,1.d0/
47 ha=dmod(pm+(t-tr)*rate,tpi)
48 enode=dmod(y(9)-ha+tpi,tpi)
49 along=datan2(x(2),x(1))-ha
50 along=dmod(along+tpi+tpi,tpi)
51 alat=dasin(x(3)/y(16))
52 argl=dmod(y(14)-y(9)+tpi,tpi)
53 ealt=y(16)-dsqrt(
re**2*(one-ellip**2)/(one-(ellip*dcos(alat))**2))
54 palt=y(1)*(one-y(7))-
re
62 IF (icase.EQ.3.OR.icase.EQ.4)
THEN
63 IF (dabs(y(11)).LT.small.OR.dabs(y(11)-tpi).LT.small) icase=3
64 IF (dabs(y(11)-
pi).LT.small) icase=4
67 IF (iprint.EQ.1.OR.inode.EQ.1)
THEN
71 IF (
node.EQ.0.AND.(dabs(argl).LT.small.OR.dabs(argl-tpi).LT.small)
73 IF (
node.EQ.0.AND.dabs(argl-
pi).LT.small) icase=2
74 IF (nrev.EQ.0.AND.(dabs(y(11)).LT.small.OR.dabs(y(11)-tpi).LT.smal
76 IF (nrev.EQ.0.AND.dabs(y(11)-
pi).LT.small) icase=4
78 IF (iprint.EQ.1.OR.inode.EQ.1)
THEN
79 IF (icase.EQ.0)
GO TO 105
80 GO TO (101,102,103,104) icase
106 WRITE(lun,2001)iday,hour,x
112 WRITE(8,9000)icase,
node,nrev,dday,thour,y(1)
113 1 ,(y(i),i=7,11),along,alat,enode,ealt,palt,alt
116 1001
FORMAT(/5x,
'PRINTING INFORMATION AT ASCENDING NODE # ',i5)
117 1002
FORMAT(/5x,
'PRINTING INFORMATION AT DESCENDING NODE #',i5)
118 1003
FORMAT(/5x,
'PRINTING INFORMATION AT PERIAPSIS OF REV #',i5)
119 1004
FORMAT(/5x,
'PRINTING INFORMATION AT APOAPSIS OF REV # ',i5)
120 1000
FORMAT(/5x,i5,
' DAYS',f12.6,
' HOURS FROM EPOCH')
121 2001
FORMAT(i3,f11.7,3f11.4,3f11.7)
122 2000
FORMAT(1p,5x,
'CARTESIAN COORD X, Y, Z, XD, YD, ZD',/,6d22.14)
123 3000
FORMAT(1p,5x,
'EQUINOCTIAL ELEMENTS A, H, K, P, Q, MEAN LONG',/
125 4000
FORMAT(1p,5x,
'CLASSICAL ELEMENTS E, I, NODE, W, MA, TA, EA, TRUE L
126 1ONG, ECC LONG, ARG OF LAT',/,5d22.14/5d22.14)
127 5000
FORMAT(1p,5x,
'OTHER PARAMETERS R, V, ENODE, LONG, LAT, HOUR ANGLE/
128 1 ALT, ELLIPSOIDAL ALT, PERIAPSIS ALT, PERIOD',/,6d22.14/6d22.14)
129 9000
FORMAT(3i5,2f12.6,f14.5,f12.8,7f12.5,3f14.5)