############## JHIBBARD AIPS++ TESTING LOG ################
##############   FOR MONTH OF APRIL, 2002  ################

Spent most of the end of March, begining of April testing various view
tools, including imageprofileplotter, and sending in enhancement
requests. For this, I did *not* use cubes made in aips++, but used my
NGC4038CD datacubes that I had on disk. I would still like to go back
and make sure that the solutions to the gaussian fitting to line
profiles agrees with what I got from karma+IDL, but later.

George just put a new VLA reduction cookbook into the system, so I am
reading through that.  First off, it occurs to me that perhaps I can
get by my previous problems by simply copying channels 3-122 from my
ms to a new ms, and maybe then I can calibrate it. So lets try that
first.

Put these lines into cal4038.g:

################################

BCHAN			:= 3;  
ECHAN			:= 122;
NCHAN                   := ECHAN-BCHAN+1;

include 'ms.g'

myms:=fitstoms(fitsfile=FITSFILENAME, msfile=spaste(BASENAME,'.tempms'));
myms.selectchannel(nchan=NCHAN, 
     start=BCHAN, 
     width=1, 
     inc=1);

include 'table.g';
ok:=tablecopy(tablename   =spaste(BASENAME,'.tempms'), 
              newtablename=spaste(BASENAME,'.ms'), 
              deep=T) 
myms.done();

################################

Didn't work. Left a message w/George; resubmited as enhancement
to "AOCso03259: would like an ms.copy function". 

4/25/02: 
Read data into AIPS, copied to dataset only using channels 3-122. 
File is called /export/data_1/jhibbard/DATA/n4038/AH582/AH582.960603.GDLINE'.
Edited cal4038.g to read this in. No channel flagging, but
keep other flag commands. Things to check:

1) Does msplot let me see calibrated data?

   No. Tried to browse the table via dc->view, but can't see corrected
   rows. So this is not linked to the fact that I deleted channels
   before. Here is the error message: 

log window:
successful readonly open of default-locked table /export/data_1/jhibbard/aips++/n4038test.ms: <fail>: shape: Table::lock cannot sync; another process changed the number of columns
	File:	tableserver.g, Line 71
	Stack:	             .(), tableserver.g line 645
		             .(), tableserver.g line 652
		             .(), tableserver.g line 155
		        spaste(), tableserver.g line 156
		     tablenote(), tableserver.g line 158
		... columns, 77213 rows

If I get out of aips++ and back in, **THEN** I can see the corrected
and model table columns, and plot them in msplot.

I update this information to AOCso03457

2) Does ch0 map look right (specline_mapch0.g)?

Yep. 

Specifically, this is what I got before (2002-01-03-hibbard.txt)
(note: AIPS was selfcal'ed):

Old numbers:
#JEH>              aips++        AIPS
#JEH> noise:        2.4 mJy/b     0.63 mJy/b
#JEH> n4038 peak:  99.8 mJy/b    122.6 mJy/b
#JEH> n4038 flux: 511.9 mJy      589.4 mJy
#JEH> sw peak:     61.0 mJy/b     66.3 mJy/b
#JEH> sw flux:     61.5 mJy       70.8 mJy
#JEH> se peak:    105.6 mJy/b    113.3 mJy/b
#JEH> se flux:    115.6 mJy      114.9 mJy
#JEH> east peak:   30.7 mJy/b     32.3 mJy/b
#JEH> east flux:   27.4 mJy       31.6 mJy

New numbers:
#JEH> noise:        1.9 mJy/b     0.63 mJy/b
#JEH> n4038 peak:  98.9 mJy/b    122.6 mJy/b
#JEH> n4038 flux: 496.0 mJy      589.4 mJy
#JEH> sw peak:     64.5 mJy/b     66.3 mJy/b
#JEH> sw flux:     70.6 mJy       70.8 mJy
#JEH> se peak:    110.5 mJy/b    113.3 mJy/b
#JEH> se flux:    118.1 mJy      114.9 mJy
#JEH> east peak:   26.1 mJy/b     32.3 mJy/b
#JEH> east flux:   26.6 mJy       31.6 mJy


Restoring Beam   : 40.7747 arcsec, 19.7435 arcsec, 81.0597 deg 
          AIPS was 51.1744 arcsec, 41.2073 arcsec, 15.0826 deg
          Diff in beam areas is AIPS:aips++ = 2.62:1
Flux can differ by up to 15%. Should check this out, but
perhaps I need to selfcal in aips++ first. Or go back to
AIPS and do reduction w/o selfcal.


3) Does cube look flat (specline_mapcal.g, specline_mapcube.g)?

   mapcal: yep, looks pretty flat.
   mapcube: had to change channel ranges - channels 119-120 are
  bad. Should look into that. After redoing, looked good, but need
  to subtract continuum, so ran specline_subcont.g.
  Basically ok, but clearly some problems. Why is beam so small?
  Stripes. But at least basic calibration is working. Now lets
  go back to the data and see if we can find bad data.

  First, I'll update all the specline_*.g scripts to use the
  new dataset and flagging commands.