[daip] Question about DECONV.FOR {External}

Hanno Spreeuw h.spreeuw at esciencecenter.nl
Thu Oct 16 09:16:39 EDT 2025


Hello,

I have been running DECONV.FOR for some time and its output always 
seemed fine.
I guess the current version - 31DEC25/APL/SUB/DECONV.FOR - does not 
differ from the previous ones.

*Yet, I am a bit confused about its output when the convolved source is 
aligned with the restoring beam.*
_Then the returned, i.e. deconvolved, position angle is 90° larger than 
those position angles._

E.g.
/deconv(7.1941910149411346, 6.826481869094818, 86.02929942633546, 
5.699093590220353, 4.181319049159876, 86.7612722265415)/
returns
/(5.396174907684326, 4.390378952026367, 177.1446533203125, 0)
/So rpa=/177.1446533203125/ is 90° larger than one would expect.
Could you please help me to explain this?

Kind regards,
Hanno Spreeuw

PS To reproduce using Python, see the attached source file, which is the 
original DECONV.FOR with a few additions for f2py.
/python -m numpy.f2py -c DECONV.FOR -m deconv_python
ipython
import numpy as np
import deconv_python
# print(deconv_python.deconv.__doc__)
deconv_python.deconv(7.1941910149411346, 6.826481869094818, 
86.02929942633546, 5.699093590220353, 4.181319049159876, 86.7612722265415)
---> (5.396174907684326, 4.390378952026367, 177.1446533203125, 0)/

-- 
handtekening - eScience Center | Enabling digitally enhanced research
	
Hanno Spreeuw
eScience Research Engineer
+31 6 21 537 257
Netherlands eScience Center
www.esciencecenter.nl <https://www.esciencecenter.nl/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20251016/8ae16eb5/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DECONV.FOR
Type: text/x-fortran
Size: 4185 bytes
Desc: not available
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20251016/8ae16eb5/attachment.bin>


More information about the Daip mailing list