search

ratt-ru / meqtrees

A library for implementing radio astronomical Measurement Equations
http://meqtrees.net
10 stars 2 forks source link

(r8829) pybeams_fits (turbo-sim.py) doesn't like my (hitherto working) beam FITS files #851

Closed gijzelaerr closed 10 years ago

gijzelaerr commented 10 years ago
at 2013-01-07 20:38:51 Ian Heywood reported:

(r8829) pybeams_fits (turbo-sim.py) doesn't like my (hitherto working) beam FITS files

gijzelaerr commented 10 years ago

Original comment thread migrated from bugzilla

at 2013-01-07 20:38:51 Ian Heywood replied:

I've been scratching my head long enough to consider this more of a bug report in good-faith than a cry for help: I can't get pybeams_fits to pick up any beam FITS images. It's throwing this error, which we've seen before:

Traceback (most recent call last): File "/home/ianh/trunk/Cattery/Siamese/OMS/InterpolatedBeams.py", line 512, in get_result l,dl = unite_shapes(l,dl); File "/home/ianh/trunk/Cattery/Siamese/OMS/InterpolatedBeams.py", line 37, in unite_shapes if a.shape == b.shape: AttributeError: 'int' object has no attribute 'shape'

This happens on both my beast machine (with and without complete explicit paths) as you can see here:

http://www-astro.physics.ox.ac.uk/~ianh/missingbeams1.jpg

and on my laptop, as you can see here:

http://www-astro.physics.ox.ac.uk/~ianh/missingbeams2.jpg

The former is using FITS images that I've generated myself from the cassbeam output, but I don't think this is the problem as the latter is using the X-band beams that we used for the WHDF paper. These were definitely working at CALIM time as I used this very setup to do the source suppression simulations.

Apologies if I'm doing something stupid, and thanks for your time.

Cheers.

at 2013-01-07 21:00:47 Oleg Smirnov replied:

Ah hang on, have you got only a single source in your LSM?

It's definitely a bug, but I have a feeling it only affects the single-source case.

at 2013-01-07 21:07:02 Ian Heywood replied:

In the two examples I screenshotted I was using the gridded sky module with tens of sources. I just did a check and the error persists using Tigger format sky models with both single and multiple sources. It was also definitely working with a single-source Tigger LSM for the source suppression sims as I was just applying the beam effects to the single bright confusing source.

Cheers.

at 2013-01-07 21:12:04 Oleg Smirnov replied:

Hmmm, weird... ok, next question, is sky rotation on?

at 2013-01-07 21:16:50 Ian Heywood replied:

Yes, sky rotation is on.

at 2013-01-08 10:34:57 Oleg Smirnov replied:

All right -- can you either set it up on jake or elwood, or upload a tarball with your MS, your fits files and your .tdl.conf?

at 2013-01-09 14:05:35 Ian Heywood replied:

I've set up both the L-band and X-band sims on Jake.

/home/ianh/Data/VLA_Beams/L-band/ /home/ianh/Data/VLA_Beams/X-band/

If you enter either of these folders and source run_sim.sh then the bug should spew forth to the terminal window.

I've also left a 1665 MHz cassbeam sim and the script I used to generate the L-band FITS images in:

/home/ianh/Data/VLA_Beams/

although it seems to hang on Jake when it comes to generating the plots (a matplotlib error maybe, certainly not critical -- I can modify it so it just generates the FITS files easily enough).

Cheers.

at 2013-01-10 17:49:32 Oleg Smirnov replied:

OK, fixed it in r8832.

One of these how-did-this-ever-work scenarios. Perhaps that's because it should only have affected sims without pointing error, and we always ran it with?