[daip] RUNNING RM synthesis tasks (or getting ready to do so)

Eric Perlman eperlman at fit.edu
Wed Jun 8 15:25:47 EDT 2016


Hi,

I am trying to do a rotation measure synthesis analysis on some recent 
VLBA data.  I would like to do it using DOFARS but am having trouble 
getting my data ready to do so.  I am using AIPS version 31-DEC-2015.

To this point, I have done the following:

1) Run SPLIT (with APARM=0) on the data, which has been polarization 
calibrated.

2) Obtained a high-quality self-calibrated Stokes I image, using 
multiple iterations of self-cal.

3) Self-calibrated the data with that Stokes I image.

4) Run IMAGR with no cleaning, FREQID=-1 and FQTOL=-1 to create Stokes 
Q, U, V images.
     This created an image that has a header that looks like this:
AIPS 1: Image=NGC4845   (MA)         Filename=NGC4845     .QIM001. 1
AIPS 1: Telescope=VLBA               Receiver=VLBA
AIPS 1: Observer=BP191               User #= 1806
AIPS 1: Observ. date=08-OCT-2015     Map date=06-JUN-2016
AIPS 1: Minimum=-3.21130268E-03      Maximum= 2.84131849E-03 JY/BEAM
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type    Pixels   Coord value     at Pixel     Coord incr Rotat
AIPS 1: RA---SIN  1024    12 58 01.197     512.00      -0.000300 0.00
AIPS 1: DEC--SIN  1024    01 34 32.420     513.00       0.000300 0.00
AIPS 1: FREQ       128   4.8520000E+09       1.00  2.5000000E+05 0.00
AIPS 1: STOKES       1   2.0000000E+00       1.00  1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Map type=DIRTY               Number of iterations=       0
AIPS 1: Conv size= 0.00296 X 0.00126   Position angle=   0.61
AIPS 1: Rest freq   4852.000         Vel type: OPTICAL wrt YOU
AIPS 1: Alt ref. value  0.00000E+00  wrt pixel    1.00
AIPS 1: Maximum version number of extension files of type CG is   1
AIPS 1: Maximum version number of extension files of type HI is   1
AIPS 1: Keyword = 'WTNOISE '  value =  1.19735482D+00
AIPS 1: Keyword = 'SUMWTIN '  value =  4.81890338D+00
AIPS 1: Keyword = 'OLDRFQ  '  value =  4.85200000D+09
AIPS 1: Keyword = 'CORRELAT'  value = 'DIFX    '
AIPS 1: Keyword = 'MAXABSU '  value = -1.000000E+00

When I try to run FQUBE the following error messages appear:

LOCALH> FQUBE1: Task FQUBE  (release of 31DEC15) begins
LOCALH> FQUBE1: NGC4845     .QIM001.   2 MA 1806 NOT FOUND
LOCALH> FQUBE1: NGC4845     .QIM001.   3 MA 1806 NOT FOUND
LOCALH> FQUBE1: NGC4845     .QIM001.   4 MA 1806 NOT FOUND
LOCALH> FQUBE1: NGC4845     .QIM001.   5 MA 1806 NOT FOUND
LOCALH> FQUBE1: ZCREAT: INVALID BYTE REQUEST =           0
LOCALH> FQUBE1: ZCREAT ERROR     4 FOR NGC4845     .FQUBE .   1
LOCALH> FQUBE1: Purports to die of UNNATURAL causes

If I try to explicitly tell FQUBE to use the FREQ axis in the image 
(with IN2SEQ=INSEQ) I get:

LOCALH> FQUBE1: quasar 31DEC15 TST: Cpu=      0.1  Real=      1 
IO=         2
LOCALH> FQUBE1: Task FQUBE  (release of 31DEC15) begins
LOCALH> FQUBE1: ZCREAT: INVALID BYTE REQUEST =           0
LOCALH> FQUBE1: ZCREAT ERROR     4 FOR NGC4845     .FQUBE .   1
LOCALH> FQUBE1: Purports to die of UNNATURAL causes

Trying to use MCUBE by itself runs successfully but creates:

LOCALH> MCUBE1: Task MCUBE  (release of 31DEC15) begins
LOCALH> MCUBE1: Create NGC4845     .MCUBE .   1 (MA)  on disk  2 cno   79
LOCALH> MCUBE1: Map plane does not fit into output map.
LOCALH> MCUBE1: Skipping sequence #   1
LOCALH> MCUBE1: Blanking empty output plane    1
LOCALH> MCUBE1: Blanking empty output plane    2
LOCALH> MCUBE1: Blanking empty output plane    3
LOCALH> MCUBE1: Blanking empty output plane    4
LOCALH> MCUBE1: Blanking empty output plane    5
LOCALH> MCUBE1: Blanking empty output plane    6
LOCALH> MCUBE1: Blanking empty output plane    7
LOCALH> MCUBE1: Blanking empty output plane    8
LOCALH> MCUBE1: Blanking empty output plane    9
LOCALH> MCUBE1: Blanking empty output plane   10
LOCALH> MCUBE1: Blanking empty output plane   11
LOCALH> MCUBE1: Blanking empty output plane   12
LOCALH> MCUBE1: Blanking empty output plane   13
LOCALH> MCUBE1: Blanking empty output plane   14
LOCALH> MCUBE1: Blanking empty output plane   15
LOCALH> MCUBE1: Blanking empty output plane   16
LOCALH> MCUBE1: Blanking empty output plane   17
LOCALH> MCUBE1: Blanking empty output plane   18
LOCALH> MCUBE1: Blanking empty output plane   19
LOCALH> MCUBE1: Blanking empty output plane   20
LOCALH> MCUBE1: Blanking empty output plane   21
LOCALH> MCUBE1: Blanking empty output plane   22
LOCALH> MCUBE1: Blanking empty output plane   23
LOCALH> MCUBE1: Blanking empty output plane   24
LOCALH> MCUBE1: Blanking empty output plane   25
LOCALH> MCUBE1: Blanking empty output plane   26
LOCALH> MCUBE1: Blanking empty output plane   27
LOCALH> MCUBE1: Blanking empty output plane   28
LOCALH> MCUBE1: Blanking empty output plane   29
LOCALH> MCUBE1: Blanking empty output plane   30
LOCALH> MCUBE1: Blanking empty output plane   31
LOCALH> MCUBE1: Blanking empty output plane   32
LOCALH> MCUBE1: Blanking empty output plane   33
LOCALH> MCUBE1: Blanking empty output plane   34
LOCALH> MCUBE1: Blanking empty output plane   35
LOCALH> MCUBE1: Blanking empty output plane   36
LOCALH> MCUBE1: Blanking empty output plane   37
LOCALH> MCUBE1: Blanking empty output plane   38
LOCALH> MCUBE1: Blanking empty output plane   39
LOCALH> MCUBE1: Blanking empty output plane   40
LOCALH> MCUBE1: Blanking empty output plane   41
LOCALH> MCUBE1: Blanking empty output plane   42
LOCALH> MCUBE1: Blanking empty output plane   43
LOCALH> MCUBE1: Blanking empty output plane   44
LOCALH> MCUBE1: Blanking empty output plane   45
LOCALH> MCUBE1: Blanking empty output plane   46
LOCALH> MCUBE1: Blanking empty output plane   47
LOCALH> MCUBE1: Blanking empty output plane   48
LOCALH> MCUBE1: Blanking empty output plane   49
LOCALH> MCUBE1: Blanking empty output plane   50
LOCALH> MCUBE1: Blanking empty output plane   51
LOCALH> MCUBE1: Blanking empty output plane   52
LOCALH> MCUBE1: Blanking empty output plane   53
LOCALH> MCUBE1: Blanking empty output plane   54
LOCALH> MCUBE1: Blanking empty output plane   55
LOCALH> MCUBE1: Blanking empty output plane   56
LOCALH> MCUBE1: Blanking empty output plane   57
LOCALH> MCUBE1: Blanking empty output plane   58
LOCALH> MCUBE1: Blanking empty output plane   59
LOCALH> MCUBE1: Blanking empty output plane   60
LOCALH> MCUBE1: Blanking empty output plane   61
LOCALH> MCUBE1: Blanking empty output plane   62
LOCALH> MCUBE1: Blanking empty output plane   63
LOCALH> MCUBE1: Blanking empty output plane   64
LOCALH> MCUBE1: Blanking empty output plane   65
LOCALH> MCUBE1: Blanking empty output plane   66
LOCALH> MCUBE1: Blanking empty output plane   67
LOCALH> MCUBE1: Blanking empty output plane   68
LOCALH> MCUBE1: Blanking empty output plane   69
LOCALH> MCUBE1: Blanking empty output plane   70
LOCALH> MCUBE1: Blanking empty output plane   71
LOCALH> MCUBE1: Blanking empty output plane   72
LOCALH> MCUBE1: Blanking empty output plane   73
LOCALH> MCUBE1: Blanking empty output plane   74
LOCALH> MCUBE1: Blanking empty output plane   75
LOCALH> MCUBE1: Blanking empty output plane   76
LOCALH> MCUBE1: Blanking empty output plane   77
LOCALH> MCUBE1: Blanking empty output plane   78
LOCALH> MCUBE1: Blanking empty output plane   79
LOCALH> MCUBE1: Blanking empty output plane   80
LOCALH> MCUBE1: Blanking empty output plane   81
LOCALH> MCUBE1: Blanking empty output plane   82
LOCALH> MCUBE1: Blanking empty output plane   83
LOCALH> MCUBE1: Blanking empty output plane   84
LOCALH> MCUBE1: Blanking empty output plane   85
LOCALH> MCUBE1: Blanking empty output plane   86
LOCALH> MCUBE1: Blanking empty output plane   87
LOCALH> MCUBE1: Blanking empty output plane   88
LOCALH> MCUBE1: Blanking empty output plane   89
LOCALH> MCUBE1: Blanking empty output plane   90
LOCALH> MCUBE1: Blanking empty output plane   91
LOCALH> MCUBE1: Blanking empty output plane   92
LOCALH> MCUBE1: Blanking empty output plane   93
LOCALH> MCUBE1: Blanking empty output plane   94
LOCALH> MCUBE1: Blanking empty output plane   95
LOCALH> MCUBE1: Blanking empty output plane   96
LOCALH> MCUBE1: Blanking empty output plane   97
LOCALH> MCUBE1: Blanking empty output plane   98
LOCALH> MCUBE1: Blanking empty output plane   99
LOCALH> MCUBE1: Blanking empty output plane  100
LOCALH> MCUBE1: Blanking empty output plane  101
LOCALH> MCUBE1: Blanking empty output plane  102
LOCALH> MCUBE1: Blanking empty output plane  103
LOCALH> MCUBE1: Blanking empty output plane  104
LOCALH> MCUBE1: Blanking empty output plane  105
LOCALH> MCUBE1: Blanking empty output plane  106
LOCALH> MCUBE1: Blanking empty output plane  107
LOCALH> MCUBE1: Blanking empty output plane  108
LOCALH> MCUBE1: Blanking empty output plane  109
LOCALH> MCUBE1: Blanking empty output plane  110
LOCALH> MCUBE1: Blanking empty output plane  111
LOCALH> MCUBE1: Blanking empty output plane  112
LOCALH> MCUBE1: Blanking empty output plane  113
LOCALH> MCUBE1: Blanking empty output plane  114
LOCALH> MCUBE1: Blanking empty output plane  115
LOCALH> MCUBE1: Blanking empty output plane  116
LOCALH> MCUBE1: Blanking empty output plane  117
LOCALH> MCUBE1: Blanking empty output plane  118
LOCALH> MCUBE1: Blanking empty output plane  119
LOCALH> MCUBE1: Blanking empty output plane  120
LOCALH> MCUBE1: Blanking empty output plane  121
LOCALH> MCUBE1: Blanking empty output plane  122
LOCALH> MCUBE1: Blanking empty output plane  123
LOCALH> MCUBE1: Blanking empty output plane  124
LOCALH> MCUBE1: Blanking empty output plane  125
LOCALH> MCUBE1: Blanking empty output plane  126
LOCALH> MCUBE1: Blanking empty output plane  127
LOCALH> MCUBE1: Blanking empty output plane  128
LOCALH> MCUBE1: Appears to have ended successfully

Even with DOCONCAT=2 it seems to not create an image that has a FQID 
axis. Can you help me?

Thanks,

Eric


-- 
Eric S. Perlman, Professor 		eperlman at fit.edu
Physics and Space Sciences Department	Phone:  +1 321 674 7741
Florida Institute of Technology	     	Fax:    +1 321 674 7482
150 West University Boulevard		http://my.fit.edu/~eperlman
Melbourne, FL  32901			http://astro.fit.edu/



More information about the Daip mailing list