1
! @(#)precesopto.prg 19.1 (ES0-DMD) 02/25/03 14:27:39
2
!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
3
!.IDENTIFICATION: precesopto.prg
4
!.PURPOSE: Precess coordinates in table created by CREATE/OPTABLE
6
!.USE: @s precess [inp_tab] [equinox1]
7
!.AUTHOR: Alessandra Gemmo, department of Astronomy, Padova
8
!.VERSION: 910826 AG Creation
9
!------------------------------------------------------------------------------
10
DEFINE/PARAM P1 opto1.tbl T "Enter input table [Default=opto1.tbl]: "
11
DEFINE/PARAM P2 {EQUINOX(2)} N "Enter new equinox of coordinates: "
12
DEFINE/LOCAL INDX/I/1/1 0
13
DEFINE/LOCAL TABEQUI/D/1/1 1950.
15
COMPUTE/KEYW INDX = M$INDEX(P1,".tbl")
17
WRITE/KEYW INPUTFI/C/1/50 'P1'.tbl
19
WRITE/KEYW INPUTFI/C/1/50 'P1'
24
COPY/DK 'INPUTFI' TABEQUI/D/1/1 TABEQUI/D/1/1
25
IF 'TABEQUI(1)' .EQ. 'EQUINOX(2)' THEN
26
WRITE/OUT "*** FATAL: Table has already been precessed !!"
30
IF 'P2' .LT. 1800. THEN
31
WRITE/OUT "*** FATAL: Value of EQUINOX not allowed !!"
33
ELSEIF 'P2' .GT. 2100. THEN
34
WRITE/OUT "*** FATAL: Value of EQUINOX not allowed !!"
39
IF 'P2' .EQ. {EQUINOX(2)} THEN
40
WRITE/DESCR 'INPUTFI' TABEQUI/D/1/1 'EPOCH(1)'
42
WRITE/DESCR 'INPUTFI' TABEQUI/D/1/1 'P2'
43
WRITE/KEYW EQUINOX/D/2/1 'P2'
48
NAME/COLU 'INPUTFI' :RA :OLDRA
49
NAME/COLU 'INPUTFI' :DEC :OLDDEC
50
NAME/COLU 'INPUTFI' :PRA :RA
51
NAME/COLU 'INPUTFI' :PDEC :DEC