From dshepher@aoc.nrao.edu Tue Nov 20 17:27:28 2001 Date: Tue, 20 Nov 2001 17:33:19 -0500 (EST) From: Debra Shepherd To: Ed Fomalont Cc: cbrogan@cv3.cv.nrao.edu, cchandle@cv3.cv.nrao.edu, dshepher@cv3.cv.nrao.edu, gtaylor@cv3.cv.nrao.edu, gvanmoor@cv3.cv.nrao.edu, jhibbard@cv3.cv.nrao.edu, jmangum@cv3.cv.nrao.edu, jmcmulli@cv3.cv.nrao.edu, julvesta@cv3.cv.nrao.edu, kdyer@cv3.cv.nrao.edu, mclausse@cv3.cv.nrao.edu, pjewell@cv3.cv.nrao.edu, rzavala@cv3.cv.nrao.edu, smyers@cv3.cv.nrao.edu, wbrisken@cv3.cv.nrao.edu Subject: AIPS++ testing notes 20nov01 Last time I was in a NAUG meeting, many people mentioned that they like to see the testing notes. Here are notes on my first try at GBT OTF continuum data using Joe's new scripts. Summary: Joe filled the data for me and fixed some header problems. I converted the measurement set to an image using glish scripts for GBT continuum processing. Worked like a charm. Ran into one, non-critical bug when I tried to do something that wasn't in a normal data path (and wasn't needed). Still need to scale the data to Janskys -- maps come out in T_ANT. To scale the data I still need a scaling number from GB (incorporating T_CAL, efficiency, effective telescope area, etc). Detailed notes: ================================ gbtcalutils.g is the updated script to reduce GBT continuum data that Joe created. Put this in the directory in which you are working. source /aips++/stable/aipsinit.csh /export/data_1/aips/dss/app/orion.gb AIPS++ Version 1.6 Build 363 glish -l gbtcalutils.g NORMAL: gbtcalutils (gc) is ready for use # Look at the available tools in gc = GBT continuum imaging: - field_names(gc) checkms setdata contcal makeimage covercorr plotsource debug # Make the 1st map: gc.setdata('oria1_DCR'); NORMAL: successful readonly open of default-locked table oria1_DCR/NRAO_GBT_GLISH: 28 columns, 41 rows NORMAL: successful readonly open of default-locked table oria1_DCR/NRAO_GBT_IF: 20 columns, 82 rows NORMAL: successful read/write open of default-locked table oria1_DCR/FEED: 12 columns, 1 rows NORMAL: successful read/write open of default-locked table oria1_DCR/POLARIZATION: 4 columns, 1 rows NORMAL: successful read/write open of default-locked table oria1_DCR/SPECTRAL_WINDOW: 14 columns, 1 rows NORMAL: successful readonly open of default-locked table oria1_DCR/POINTING: 15 columns, 4874 rows # This appears to attach the data to the gc tool # do the basic calibration with "receptor 1" data to get T_ANT # from measured voltages with the noise-diode ON and OFF. gc.contcal(); NORMAL: Starting server imager NORMAL: Server started: /aips++/stable/linux/bin/imager (AIPS++ version: 1.6 (build #363)) NORMAL: Opening MeasurementSet oria1_DCR NORMAL: Successfully closed empty server: imager NORMAL: successful read/write open of default-locked table oria1_DCR: 28 columns, 19516 rows # imager is opened and closed ==> CORRECTED & MODEL_DATA columns are created now. # T_CAL default = 1, receptor default = 1, T_ANT now calculated assuming T_CAL=1 # Select, regrid, and construct an "unnormalized" image gc.makeimage('orion1im'); NORMAL: Starting server imager NORMAL: Server started: /aips++/stable/linux/bin/imager (AIPS++ version: 1.6 (build #363)) NORMAL: Opening MeasurementSet oria1_DCR rm: cannot remove `orion1im': No such file or directory rm: cannot remove `orion1im_weight': No such file or directory NORMAL: Starting imager::setdata NORMAL: Selecting data NORMAL: Performing selection on MeasurementSet NORMAL: Selecting on field and spectral window ids NORMAL: By selection 19516 rows are reduced to 4879 NORMAL: Finished imager::setdata 4.08 real 0.69 user 0.37 system NORMAL: successful readonly open of default-locked table oria1_DCR/NRAO_GBT_GLISH: 28 columns, 41 rows NORMAL: Defining image properties WARN: nx = 66 is not composite; nx = 72 or 64 will be more efficient WARN: ny = 66 is not composite; ny = 72 or 64 will be more efficient NORMAL: Setting processing options NORMAL: Starting imager::weight NORMAL: Weighting MS: IMAGING_WEIGHT column will be changed NORMAL: Natural weighting NORMAL: Sum of weights = 4879 NORMAL: Finished imager::weight 9.38 real 6.02 user 1.12 system NORMAL: Starting imager::makeimage NORMAL: Calculating image (without full skyequation) NORMAL: Making single dish image from corrected data NORMAL: Image is : orion1im NORMAL: Frequency = 8.4351 GHz, synthesized continuum bandwidth = 0.05 GHz NORMAL: Image polarization = Stokes I NORMAL: Frequency = 8.4351 GHz, synthesized continuum bandwidth = 0.05 GHz NORMAL: Image polarization = Stokes I NORMAL: Performing Single Dish gridding with convolution function SF Gridding will use specified common tangent point: 05:35:17.47 -05.23.07.50 J2000 NORMAL: Using prolate spheriodal wave function as the convolution function NORMAL: Transforming I only NORMAL: Finished imager::makeimage 13.06 real 2.06 user 0.52 system NORMAL: Starting imager::makeimage NORMAL: Calculating image (without full skyequation) NORMAL: Making single dish coverage function NORMAL: Image is : orion1im_weight NORMAL: Frequency = 8.4351 GHz, synthesized continuum bandwidth = 0.05 GHz NORMAL: Image polarization = Stokes I NORMAL: Frequency = 8.4351 GHz, synthesized continuum bandwidth = 0.05 GHz NORMAL: Image polarization = Stokes I NORMAL: Using prolate spheriodal wave function as the convolution function NORMAL: Transforming I only NORMAL: Finished imager::makeimage 2.12 real 1.54 user 0.15 system NORMAL: Successfully closed empty server: imager # Created 2 files: the uncorrected image: orion1im/ # & the coverage map: orion1im_weight/ # Create a corrected image by dividing the uncorrected image by the coverage map: gc.covercorr(); NORMAL: Starting server app_image NORMAL: Server started: /aips++/stable/linux/bin/app_image (AIPS++ version: 1.6 (build #363)) NORMAL: Selected bounding box : [1, 1, 1, 1] to [66, 66, 1, 1] (05:35:57.980, -05.33.12.210, I, 8.435100e+09Hz to 05:34:38.209, -05.13.20.955, I, 8.435100e+09Hz) NORMAL: Creating new statistics storage lattice of shape [9] NORMAL: NORMAL: Number points = 4.356000e+03 Sum = 4.355996e+03 NORMAL: Mean = 9.999990e-01 NORMAL: Variance = 1.923394e-01 Std dev = 4.385652e-01 NORMAL: Rms = 1.091922e+00 NORMAL: NORMAL: Minimum value 0.000000e+00 at [1, 1, 1, 1] (05:35:57.980, -05.33.12.210, I, 8.435100e+09Hz) NORMAL: Maximum value 1.426234e+00 at [7, 16, 1, 1] (05:35:50.610, -05.28.37.331, I, 8.435100e+09Hz) NORMAL: NORMAL: Creating image `orion1im_corr' of shape [66, 66, 1, 1] NORMAL: Created mask `mask0' # created a "corrected image" orion1im_corr/ # view the image: include 'image.g' im := image('orion1im_corr'); im.view() # works beautifully # image is what I expect # Try plotsource to fit a source: gc.plotsource('orion1im_corr',nsource=1); rm: cannot remove `tst.tbl': No such file or directory NORMAL: successful creation of default-locked table tst.tbl: 4 columns, 0 rows NORMAL: created table tst.tbl with zero rows NORMAL: Starting image::findsources NORMAL: Selected bounding box : [1, 1, 1, 1] to [66, 66, 1, 1] (05:35:57.980, -05.33.12.210, I, 8.435100e+09Hz to 05:34:38.209, -05.13.20.955, I, 8.435100e+09Hz) NORMAL: Found 1 sources SEVERE: Cannot convert units of brightness to Jy - will assume Jy NORMAL: Finished image::findsources 1.36 real 0.02 user 0.05 system NORMAL: Starting server componentlist NORMAL: Server started: /aips++/stable/linux/bin/componentlist (AIPS++ version: 1.6 (build #363)) NORMAL: defaultentryparser (dep) ready F # Only one source was found (that is all I asked for). OK # there are low-level stripes in the background in the direction of the otf # scanning. # stats function now in the tools option. # rms = 22.16 # max = 84.82 # min = 17.29 # matches Joe's # end of image processing -- looks beautiful. # Trying starting plotsource with nsource=0: gc.plotsource('orion1im_corr',nsource=0); NORMAL: successful creation of default-locked table tst.tbl: 4 columns, 0 rows NORMAL: created table tst.tbl with zero rows NORMAL: Starting image::findsources NORMAL: Selected bounding box : [1, 1, 1, 1] to [66, 66, 1, 1] (05:35:57.980, -05.33.12.210, I, 8.435100e+09Hz to 05:34:38.209, -05.13.20.955, I, 8.435100e+09Hz) LocalExec::SetStatus: abnormal child termination for /aips++/stable/linux/bin/app_image Server 'app_image' has failed unexpectedly! You will need to create the relevant tool again. If that causes unexpected behavior, please restart AIPS++ Please submit a bug-report using bug() if you can reproduce the problem. : no sources found File: ./gbtcalutils.g, Line 320 Stack: .() submitted bug() report end of testing.