Shadow of Io Program

SUBROUTINE SHADOW ( TIMEVO, SHADOW )
DOUBLE PRECISION TIMEVO	        ...
INTEGER JUPITE
PARAMETER (JUPITE = 599)        ...
DOUBLE PRECISION RADJUP ( 3 )   ...
CALL BODVAR ( JUPITE, 'RADII', DMYO, RADJUP ) 
TJUPIT = SENT ( JUPITE, VOYGR2, TIMEVO ) 
CALL FINDPV ( JUPITE, TJUPIT, PJUPIT, DMY20 ) 
CALL BODMAT ( JUPITE, TJUPIT, MJUPIT ) 
TIO = SENT ( IO, JUPITE, TJUPIT ) 
CALL FINDPV ( IO, TIO, PIO, DMY30 ) 
TSUN = SENT ( SUN, 10, TIO ) 
CALL FINDPV ( SUN, TSUN, PSUN, DMY40 ) 
CALL VSUB ( PIO, PSUN, DPSPI ) 
CALL VSUB ( PSUN, PJUPIT, DPJPS ) 
CALL MXV ( MJUPIT, DPSPI, XDPSPI ) 
CALL MXV ( MJUPIT, DPJPS, XDPJPS )
CALL SURFPT ( XDPJPS, XDPSPI, RADJUP ( 1 ),
              RADJUP, RADJUP ( 3 ), P, DMY90 )
CALL VSUB ( P, PJUPIT, DPJUPP )
CALL MTXV ( MJUPIT, DPJUPP, SHADOW )
END