[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