From tcornwel@cv3.cv.nrao.edu Wed Apr 19 14:06:44 2000 Date: Wed, 19 Apr 2000 14:03:54 -0600 From: Tim Cornwell To: Steven T. Myers , Tim Cornwell , Athol Kemball Cc: Greg Taylor Subject: RE: aips++ test log 14Apr00 - 18Apr00 [ The following text is in the "iso-8859-1" character set. ] [ Your display is set for the "US-ASCII" character set. ] [ Some characters may be displayed incorrectly. ] Steve, Some comments and replies. Cheers, Tim > >>>ENHANCEMENT: VISPLOT AS IT STANDS IS BARELY USEABLE. DEVELOPMENT > OF A MORE COMPREHENSIVE AND INTUITIVE VISIBILITY VISUALIZATION > TOOL WOULD BE PREFERABLE We agreed and it's planned to happen. visplot is an interim solution. > > 3. Now try imager. Use script > > imgr:=imager('J0739+016.ms') # make imager object > > ERROR: unhandled mount type > (eventually ended with 1711 errors!) > > I knew this was going to be a problem, as Tim Pearson put the mount > type in the FITS header as > > imgr.summary() > > ================================================================== > ============= > MeasurementSet Name: J0739+016.ms MS Version 1 > ================================================================== > ============= > Observer: Project: > Array: CBI > Fields: 1 > ID Name Right Ascension Declination Type > 1 J0739+016 07:39:18.03 +01.37.04.62 J2000 > Data records: 210100 Total integration time = 1089.25 seconds > Observed from 01-Mar-2000/23:58:22 to 02-Mar-2000/00:16:31 > Spectral Windows: 10 > Ref.Freq RestFreq #Chans Resolution TotalBW Correlations > 35500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 34500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 33500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 32500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 31500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 30500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 29500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 28500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 27500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > 26500 MHz continuum 1 1e+06 kHz 1e+06 kHz LL > Feeds = 12: printing first row only > Antenna Spectral Window # Receptors Polarizations > 1 -1 2 [ R, L] > Antennas: 13 > 1-13: RX0 RX1 RX2 RX3 RX4 RX5 RX6 RX7 RX8 RX9 RX10 RX11 RX12 > ------------------------------------------------------------------ > ------------- > > Tables: > MAIN 210100 rows > ANTENNA 13 rows > ARRAY 1 row > FEED 12 rows > FIELD 1 row > OBSERVATION 1 row > OBSLOG 25 rows > SOURCE (see FIELD) > SPECTRAL_WINDOW 10 rows > SYSCAL > WEATHER > ================================================================== > ============= > > General: > MeasurementSet is J0739+016.ms > Beam fit is not valid > Image: parameters not yet set (use setimage) > Data selection settings: (use setdata to change) > All data selected > Options settings: (use setoptions to change) > Gridding cache has 4194304 complex pixels, in tiles of 16 > pixels on a side > Gridding convolution function is Spheroidal wave function > No primary beam correction will be made > Image plane padding : 1 > > > Seems to have gone ahead and made the imager object anyway! > Continue... > > There seem to be 10 spectral windows (the 10 CBI channels) each of > 1 channel 1 GHz wide. > > imgr.setimage(cellx='1arcmin', celly='1arcmin', nx=128, ny=128, > stokes='LL', spwid=[1,2,3,4,5,6,7,8,9,10], fieldid=1) > > ERROR: Illegal Stokes string LL > > Manual sez "Allowed: 'I'|'IV'|'IQU'|'IQUV' " > WHAT, NO RR OR LL???? You can make an RR and LL complex dirty image using imager.image setting the compleximage. Set Stokes = 'IV'. It then keeps the compleximage around. > > I guess "I" will do what I want, however one of the 13 CBI > receivers is cross-polarized (to R) versus the 12 others (in L) > in order to do cmbr polarization. Have to watch that this gets > treated right. To do this right, we'll need an accurate description of the mount type plus the polarization behavior. With that it should be quite straightforward to add. > > I wonder how you get rid of the frame. Does "frame" have a done > function? Well, after reading under the Glish/Tk in the Glish > manual about the frame, I dont see anything. Do it manually... To get rid of a frame, set it to false i.e. myframe := F; This is documented in the Glish User Manual. > 7. 2000.04.18 : Try clean again, make ms from scratch > > include 'synthesis.g' > dowait:=T > m:=fitstoms(msfile='J0739+016.ms',fitsfile='J0739+016.uvf'); > m.close() > imgr:=imager('J0739+016.ms') > > As before, gave a whole lot of errors "unhandled mount type". > Noted strange behavior that it made 3 small SEVERE ERROR windows > in addition to the one it was printing in, but the extras went > away when the errors stopped being generated. > > Also gave errors (reported in the glish window) > error, self.vsb is not an agent > WheneverStmt ref count > error, self.text is not an agent > WheneverStmt ref count > > >>>BUG: IS THIS BEHAVIOR EXPECTED? That's a bug that I think has been fixed in the last couple of days. > > imgr.restore(model='J0739+016.clean', image='J0739+016.clean.restored', > bmaj=bmaj, bmin=bmin, bpa=bpa) > > THIS GIVES ERROR > > : named arg "bmaj" does not match any formal in call to > function (val model = , val complist = , val image = , val > residual = , val async = !dowait) { > { > { > { > { > (self.restoreRec.model := model); > (self.restoreRec.complist := complist);} > > (self.restoreRec.image := image);} > > (self.restoreRec.residual := residual);} > > (returnval := defaultservers.run(self.agent, self.restoreRec, async));} > > return returnval} > > File: imager.g, Line 239 > Stack: .() > > EVEN THOUGH > > - print bmaj,bmin,bpa > [value=256.004608, unit=arcsec] [value=239.322708, unit=arcsec] > [value=27.4288387, unit=deg] > > Lets try > > imgr.restore(model='J0739+016.clean', image='J0739+016.clean.restored', > bmaj='247arcsec', bmin='247arcsec') > > >>>BUG: STILL GIVES SAME ERROR, EVEN THOUGH THE DOCUMENTATION SAYS TO > DO IT THIS WAY. WHAT IS WRONG? The documentation is out of date! You should use setbeam to set the beam and then restore to do the restore. I'll fix that. > BOOM. LOOKS LIKE THERE HAVE BEEN SOME CHANGES NOT PROPAGATED > THROUGH CORRECTLY EVEN THOUGH THIS IS THE "STABLE" VERSION. Yes. It's a bug in the documentation. > > Fire up the Tool Manager for cl object. Aha, the dophase input > has been changed to phaseonly. > > >>>BUG: THERE SHOULD BE A MORE ELEGANT WAY THAT A FUNCTION CAN > LET YOU KNOW YOU GAVE AN UNKNOWN INPUT PARAMETER!!!! Yes!