[daip] SPFLG trouble and possible fix

Olaf Wucknitz wucknitz at astro.uni-bonn.de
Tue Apr 7 04:55:19 EDT 2009


Hi Eric,

> I think I see a better fix - passing in the list of source numbers found in
> the global read of the data at the beginning.  This involves changing the call
> sequence of GTBWRT and GRIDTB subroutines plus changes in SPFLG, TVFLG, and
> UVIMG.  It is all rather tricky I fear.  Do you have a modest data set that
> would show this problem?

I tried that, but only for SPFLG and subroutines. I don't know if TVFLG 
is also affected. This approach worked fine, but the result was the same 
as with the simpler fix. (in my case)

I boiled down my data set to a reasonable size of 9MB. You can find it 
here:

http://www.astro.uni-bonn.de/~wucknitz/temp/SPFLG_TEST.FITS

Here is the test code that may be used to show the problem:

task 'fitld'
default
datain 'PWD:SPFLG_TEST.FITS'
outd 1
outname 'GW019B'
outclass 'UVAVG'
go; wait


task 'spflg'
default
ind 1
inname 'GW019B'
inclass 'UVAVG'
docat -1
dparm(6)=16
go


Without additional smoothing (AVG=1), it shows the complete time range 
(01/05:05:09 - 01/10:01:25). Once I set AVG to > 1 (via ENTER SMOOTH 
TIME), only 01/05:05:17 - 01/08:03:41 is shown.

By the way: Is it possible to extend the allowed number of time samples 
from 16384 to something larger? It is easy to change the numbers (16384 +- 
a bit) in SPFLG.FOR, but I don't know if other things have to be modified 
too.


Cheers,
Olaf




More information about the Daip mailing list