[daip] Re: cl2hf

del at cfa.harvard.edu del at cfa.harvard.edu
Wed Oct 3 18:48:02 EDT 2001


Hi Leonia,

Firstly, apologies for the delay.  It's been a VERY busy day here.

When I wrote:
> 
> >DATUTC is defined as REAL, while in fact it should be
> >defined as DOUBLE PRECISION.
> 

The statement is based on the the TYPE output from OGET.
Specifically, SN2HF (and, apparently, CL2HF) gets DATUTC from the
line:

      CALL OGET (ANTAB, 'KEY.DATUTC', TYPE, DIM, DATUTC, CDUMMY, IERR)

If you look at variable TYPE upon subroutine return, it contains
the value 1, which according to OGET documentation refers to a
double precision variable.  I can add that when DATUTC was defined
as real and the value of DATUTC was nonzero, the value of DATUTC
returned from the above subroutine call was garbage (something of
order 10^19), certainly not the value I read when looking at the
AN table via PRTAB.  When I redefined DATUTC as double precision,
the value returned from the subroutine call matched the value I
saw in the AN table.

Perhaps DATUTC is supposed to be REAL but some other location in
AIPS wrongfully defines it as DOUBLE PRECISION.

I maintain that, based on my (albeit limited) understanding, the
value of DATUTC stored in the AN table should be 0 regardless of whether
it is REAL or DOUBLE PRECISION.  DATUTC is not the same thing as
UT1UTC (which I state because in the last few experiments I've looked at,
DATUTC has been set equal to UT1UTC).

If you change AIPS so that DATUTC is defined as REAL throughout, please
let me know so that I can make the change in SN2HF.

Thanks for your prompt attention to this matter.

Dan




More information about the Daip mailing list