[daip] visibility weight problem in output of DIFUV

Eric Greisen egreisen at nrao.edu
Thu Nov 12 15:18:24 EST 2009


Huib Intema wrote:
> Dear DAIP,
> 
> I use DIFUV to subtract a model UV database from a real UV database. The 
> model UV database is created with UVSUB with OPCODE = 'MODL', using the 
> real UV database as input, which results in model visibilities having 
> the same weights as in the input. Before (using a mid-2008 copy of 
> AIPS), I needed to run WTMOD on the output of DIFUV with APARM(1)=0.5 to 
> compensate for the unwanted summation of (quadratic?) weights by DIFUV.
> 
> I found that DIFUV in the current 2009 AIPS causes the output visibility 
> (quadratic?) weights to be halved rather than doubled. This I could fix 
> by using WTMOD with APARM(1)=2.0. However, when the input visibility 
> weights are both zero the output weight is NaN. This looks like a zero 
> division error while adding the reciprokes of the input visibility 
> weights rather than the weights themselves.

Sorry - this slipped through the cracks of EVLA needs.  I have fixed the 
computation to be more careful (it would abort on Suns and perhaps 
others).  Tomorrow's MNJ should fix it.

Eric Greisen




More information about the Daip mailing list