NASA Ocean Color

ocssw V2022
       SUBROUTINE pout(T,X,LUN)  
 C  VERSION OF 4/13/87
 C  PURPOSE
 C    PRINTOUT ROUTINE FOR TRAJECTORY PROGRAM   
 C  INPUT 
 C    T      = CURRENT TIME (SEC)
 C    X      = 6-D CARTESIAN COORD X,Y,Z,XD,YD,ZD (KM,KM/SEC)
 C    SEE ROUTINES ASAP AND DER FOR EXPLANATION OF COMMON BLOCK
 C    VARIABLES
 C  OUTPUT
 C    NONE  
 C  CALL SUBROUTINES  
 C    EQNOX 
 C  REFERENCES
 C    JPL EM 312/87-153, 20 APRIL 1987
 C  ANALYSIS
 C    J. H. KWOK - JPL  
 C  PROGRAMMER
 C    J. H. KWOK - JPL  
 C  PROGRAM MODIFICATIONS 
 C    NONE  
 C  COMMENTS  
 C    NONE
 C   
 C  MODIFICATION HISTORY
 C   Changed TIME common block name to TIMECMN, to eliminate linker 
 C   linker warnings.  B. A. Franz, GSC, November 14, 1997.
 C
  
       IMPLICIT DOUBLE PRECISION (a-h,o-z)   
       dimension x(6)
       dimension y(18)
       common/option/l,m,ires,isun,imoon,iephem,idrag,idens,isrp,iorb
      1  ,iprint,inode,iplot
       common/timecmn/ti,tf,tr  
       common/pltcon/ge,re,rate,pm,aj2,ellip,ratm
       common/count/node,nrev
       common/case/icase
       DATA hts,dth,rtd,pi,tpi,one/3.6d3,24.d0,57.29577951308232d0,
      1  3.141592653589793d0,6.283185307179586d0,1.d0/
       DATA small/1.d-3/
       thour=(t-ti)/hts
       dday=thour/dth
       iday=idint(dday)
       hour=thour-iday*dth
       CALL eqnox(x,ge,y)
       ha=dmod(pm+(t-tr)*rate,tpi)
       enode=dmod(y(9)-ha+tpi,tpi)
       along=datan2(x(2),x(1))-ha
       along=dmod(along+tpi+tpi,tpi)
       alat=dasin(x(3)/y(16))
       argl=dmod(y(14)-y(9)+tpi,tpi)
       ealt=y(16)-dsqrt(re**2*(one-ellip**2)/(one-(ellip*dcos(alat))**2))
       palt=y(1)*(one-y(7))-re
       alt=y(16)-re
 C
 C *** DOUBLE CHECK APSIS CROSSING, THIS IS NEEDED IN CASE SOME
 C *** PERTURBATION ON NEARLY CIRCULAR ORBIT IS CHANGING MEAN ANOMALY
 C *** BACKWARDS
 C
       IF (iprint.EQ.1) THEN
       IF (icase.EQ.3.OR.icase.EQ.4) THEN
         IF (dabs(y(11)).LT.small.OR.dabs(y(11)-tpi).LT.small) icase=3
         IF (dabs(y(11)-pi).LT.small) icase=4
       ENDIF
       ENDIF
       IF (iprint.EQ.1.OR.inode.EQ.1) THEN
 C
 C *** CHANGE ICASE IF FIRST PRINT OUT IS A NODAL OR APSIS CROSSING
 C
       IF (node.EQ.0.AND.(dabs(argl).LT.small.OR.dabs(argl-tpi).LT.small)
      1) icase=1
       IF (node.EQ.0.AND.dabs(argl-pi).LT.small) icase=2
       IF (nrev.EQ.0.AND.(dabs(y(11)).LT.small.OR.dabs(y(11)-tpi).LT.smal
      1l)) icase=3
       IF (nrev.EQ.0.AND.dabs(y(11)-pi).LT.small) icase=4
       ENDIF
       IF (iprint.EQ.1.OR.inode.EQ.1) THEN
       IF (icase.EQ.0) GO TO 105
       GO TO (101,102,103,104) icase
   101 node=node+1
       WRITE(*,1001)node
       GO TO 105
   102 WRITE(*,1002)node
       GO TO 105
   103 nrev=nrev+1
       WRITE(*,1003)nrev
       GO TO 105
   104 WRITE(*,1004)nrev
   105 CONTINUE
       ENDIF
 C
 C *** START UNIT CONVERSION
 C
       y(6)=y(6)*rtd
       DO 10 i=8,15
    10 y(i)=y(i)*rtd
       y(18)=y(18)/hts
       argl=argl*rtd
       ha=ha*rtd
       enode=enode*rtd
       along=along*rtd
       alat=alat*rtd
 c      WRITE(*,1000)IDAY,HOUR
 c      WRITE(*,2000)X
       WRITE(lun,2001)iday,hour,x
 c      WRITE(*,3000)(Y(I),I=1,6)
 c      WRITE(*,4000)(Y(I),I=7,15),ARGL
 c      WRITE(*,5000)(Y(I),I=16,17),ENODE,ALONG,ALAT,HA,ALT,EALT,PALT
 c     1  ,Y(18)
       IF (iplot.EQ.1) THEN
         WRITE(8,9000)icase,node,nrev,dday,thour,y(1)
      1  ,(y(i),i=7,11),along,alat,enode,ealt,palt,alt
         icase=0
       ENDIF
  1001 FORMAT(/5x,'PRINTING INFORMATION AT ASCENDING NODE # ',i5)
  1002 FORMAT(/5x,'PRINTING INFORMATION AT DESCENDING NODE #',i5)
  1003 FORMAT(/5x,'PRINTING INFORMATION AT PERIAPSIS OF REV #',i5)
  1004 FORMAT(/5x,'PRINTING INFORMATION AT APOAPSIS OF REV # ',i5)
  1000 FORMAT(/5x,i5,' DAYS',f12.6,' HOURS FROM EPOCH')
  2001 FORMAT(i3,f11.7,3f11.4,3f11.7)
  2000 FORMAT(1p,5x,'CARTESIAN COORD X, Y, Z, XD, YD, ZD',/,6d22.14)
  3000 FORMAT(1p,5x,'EQUINOCTIAL ELEMENTS A, H, K, P, Q, MEAN LONG',/
      1,6d22.14)
  4000 FORMAT(1p,5x,'CLASSICAL ELEMENTS E, I, NODE, W, MA, TA, EA, TRUE L
      1ONG, ECC LONG, ARG OF LAT',/,5d22.14/5d22.14) 
  5000 FORMAT(1p,5x,'OTHER PARAMETERS R, V, ENODE, LONG, LAT, HOUR ANGLE/
      1 ALT, ELLIPSOIDAL ALT, PERIAPSIS ALT, PERIOD',/,6d22.14/6d22.14)
  9000 FORMAT(3i5,2f12.6,f14.5,f12.8,7f12.5,3f14.5)
       RETURN
       END
eqnox
subroutine eqnox(X, GM, Y)
Definition: eqnox.f:2
node
Definition: FP_failed_pkt_queue.h:66
pout
subroutine pout(T, X, LUN)
Definition: pout.f:2
#define pi
Definition: vincenty.c:23
#define re
Definition: l1_czcs.c:695
dmod
double dmod(double a, double p)
Description:
Definition: nav.c:23