Sunspot Band6 Calibration for CASA 6.5.4: Difference between revisions

From CASA Guides
Jump to navigationJump to search
 
(51 intermediate revisions by the same user not shown)
Line 3: Line 3:


==Getting Started==
==Getting Started==
Details of this ALMA observations are provided at [[Sunspot_Band6]]. This portion of the guide will cover the calibration of the raw visibility data.  
Details of these ALMA observations are provided at [[Sunspot_Band6]]. This portion of the guide will cover the calibration of the raw visibility data.
 
<font color="red">
WARNING: The command in the [[#Tsys+Tant calibration of the Sun]] section of this guide takes a long time to run. Depending on your machine, it may be between a couple hours and one day.
</font>


===Download and Unpack the Data===
===Download and Unpack the Data===
Line 49: Line 53:
<source lang="python">
<source lang="python">
#In CASA
#In CASA
asdm ='uid___A002_Xae00c5_X2a8d'  
asdm = 'uid___A002_Xae00c5_X2a8d'  
mso = asdm + '.ms'          #MS original
mso = asdm + '.ms'          #MS original
mss = asdm + '.ms.split'    #MS split
mss = asdm + '.ms.split'    #MS split (science spws)
msc = asdm + '.ms.split.cal' #MS calibrated
msc = mss  + '.cal'         #MS calibrated
</source>
</source>
The raw data have been provided to you in the ASDM format. It is the native format of the data produced by the ALMA observatory.
The raw data have been provided to you in the ASDM format. It is the native format of the data produced by the ALMA observatory.
Line 62: Line 66:
</source>
</source>


To fix the bug in the SYSCAL table times, the following commands are executed.
To check if a bug fix is needed in the SYSCAL table times, the following commands are executed.


<source lang="python">
<source lang="python">
Line 81: Line 85:
The output will be sent to the CASA {{logger_6.5.4}}, or saved in a text file.
The output will be sent to the CASA {{logger_6.5.4}}, or saved in a text file.


* The first section shows information about the project and observing time:
The first section shows information about the project and observing time:


<pre style="background-color: #fffacd;">
<pre style="background-color: #fffacd;">
Line 92: Line 96:
====Scan List====
====Scan List====


* The second section describes details of each scan. This is often easier to read in a text viewer with line wrapping turned off. We see that our Fields have the following associated intents:
The second section describes details of each scan. This is often easier to read in a text viewer with line wrapping turned off. We see that our Fields have the following associated intents:
** Field 0 = Sun: observed for the measurement of zero-signal level [scan 1] and the calibration of ''atmosphere'' [scans 10,13,16]. Scientific observations of the Sun ('''Target''') [scans 11,14,17] are listed many times, which represent a mosaic.
* Field 0 = Sun: observed for the measurement of zero-signal level [scan 1] and the calibration of ''atmosphere'' [scans 10,13,16]. Scientific observations of the Sun ('''Target''') [scans 11,14,17] are listed many times, which represent a mosaic.
** Field 1 = J1924-2914: observed for the calibrations of ''pointing'' [scan 2], ''sideband ratio'' [scan 3], ''atmosphere'' [scan 4], and '''bandpass''' [scan 5]
* Field 1 = J1924-2914: observed for the calibrations of ''pointing'' [scan 2], ''sideband ratio'' [scan 3], ''atmosphere'' [scan 4], and '''bandpass''' [scan 5]
** Field 2 = nrao530: observed for the calibration of ''pointing'' [scan 6], ''atmosphere'' [scan 7], '''flux''' [scan 8], and '''phase''' [scans 9,12,15,18]
* Field 2 = nrao530: observed for the calibration of ''pointing'' [scan 6], ''atmosphere'' [scan 7], '''flux''' [scan 8], and '''phase''' [scans 9,12,15,18]
 
Subscans are NOT necessarily listed in chronological order, but rather ordered by Field ID. For example, for scan 14, Fields 3~9 are observed AFTER Fields 80~150, but listed first. See [[Sunspot_Band6_Imaging_for_CASA_6.5.4#Flag_the_surplus_subscans_and_baselines]].


<pre style="background-color: #fffacd;">
<pre style="background-color: #fffacd;">
Line 114: Line 120:
               19:39:23.9 - 19:49:36.0    11      4 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11      4 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11      5 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11      5 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11      6 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
 
              19:39:23.9 - 19:49:36.0    11     7 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
*** Scan 11, Fields 6~76 trimmed for brevity ***
              19:39:23.9 - 19:49:36.0    11      8 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
 
              19:39:23.9 - 19:49:36.0    11      9 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    10 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    11 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    12 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    13 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    14 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    15 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    16 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    17 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    18 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    19 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    20 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    21 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    22 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    23 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    24 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    25 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    26 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    27 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    28 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    29 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    30 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    31 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    32 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    33 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    34 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    35 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    36 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    37 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    38 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    39 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    40 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    41 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    42 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    43 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    44 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    45 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    46 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    47 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    48 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    49 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    50 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    51 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    52 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    53 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    54 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    55 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    56 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    57 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    58 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    59 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    60 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    61 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    62 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    63 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    64 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    65 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    66 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    67 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    68 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    69 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    70 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    71 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    72 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    73 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    74 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    75 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11    76 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11    77 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11    77 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11    78 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:39:23.9 - 19:49:36.0    11    78 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
Line 201: Line 139:
               19:53:26.0 - 20:03:37.0    14    81 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    81 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    82 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    82 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    83 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
 
              19:53:26.0 - 20:03:37.0    14     84 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
*** Scan 14, Fields 83~147 trimmed for brevity ***
              19:53:26.0 - 20:03:37.0    14    85 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
 
              19:53:26.0 - 20:03:37.0    14    86 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    87 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    88 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    89 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    90 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    91 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    92 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    93 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    94 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    95 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    96 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    97 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    98 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    99 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    100 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    101 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    102 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    103 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    104 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    105 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    106 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    107 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    108 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    109 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    110 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    111 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    112 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    113 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    114 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    115 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    116 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    117 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    118 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    119 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    120 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    121 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    122 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    123 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    124 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    125 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    126 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    127 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    128 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    129 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    130 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    131 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    132 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    133 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    134 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    135 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    136 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    137 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    138 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    139 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    140 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    141 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    142 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    143 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    144 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    145 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    146 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    147 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    148 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    148 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    149 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
               19:53:26.0 - 20:03:37.0    14    149 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
Line 281: Line 157:
====Field List====
====Field List====


* The third section shows the information of the observing targets. This shows that three Sources with 151 Fields were observed: The Sun, J1924-2914, and nrao530 (J1733-1304). The Sun as one Source has many Fields - this is how a mosaic is represented:
The third section shows the information of the observing targets. This shows that three Sources with 151 Fields were observed: The Sun, J1924-2914, and nrao530 (J1733-1304). The Sun as one Source has many Fields - this is how a mosaic is represented:


<pre style="background-color: #fffacd;">
<pre style="background-color: #fffacd;">
Line 292: Line 168:
   4    none Sun                17:44:04.746002 -23.20.37.96880 ICRS    0          36053
   4    none Sun                17:44:04.746002 -23.20.37.96880 ICRS    0          36053
   5    none Sun                17:44:05.707104 -23.20.37.97466 ICRS    0          36053
   5    none Sun                17:44:05.707104 -23.20.37.97466 ICRS    0          36053
  6    none Sun                17:44:06.668207 -23.20.37.98018 ICRS    0          36053
 
  7    none Sun                 17:44:07.629310 -23.20.37.98535 ICRS    0          36053
*** Sun Fields 6~147 trimmed for brevity ***
  8    none Sun                17:44:08.590412 -23.20.37.99016 ICRS    0          36022
 
  9    none Sun                17:44:09.551515 -23.20.37.99463 ICRS    0          36022
  10  none Sun                17:44:10.512617 -23.20.37.99875 ICRS    0          36022
  11  none Sun                17:44:11.473720 -23.20.38.00251 ICRS    0          36053
  12  none Sun                17:44:12.434823 -23.20.38.00593 ICRS    0          36053
  13  none Sun                17:44:13.395926 -23.20.38.00900 ICRS    0          18011
  14  none Sun                17:44:03.555941 -23.20.26.81043 ICRS    0          18011
  15  none Sun                17:44:04.517022 -23.20.26.81682 ICRS    0          18011
  16  none Sun                17:44:05.478102 -23.20.26.82286 ICRS    0          18011
  17  none Sun                17:44:06.439183 -23.20.26.82855 ICRS    0          18042
  18  none Sun                17:44:07.400264 -23.20.26.83389 ICRS    0          18011
  19  none Sun                17:44:08.361345 -23.20.26.83888 ICRS    0          18042
  20  none Sun                17:44:09.322425 -23.20.26.84352 ICRS    0          18011
  21  none Sun                17:44:10.283506 -23.20.26.84781 ICRS    0          18011
  22  none Sun                17:44:11.244587 -23.20.26.85175 ICRS    0          18011
  23  none Sun                17:44:12.205668 -23.20.26.85534 ICRS    0          18011
  24  none Sun                17:44:13.166748 -23.20.26.85858 ICRS    0          18042
  25  none Sun                17:44:14.127829 -23.20.26.86147 ICRS    0          18011
  26  none Sun                17:44:04.288064 -23.20.15.66473 ICRS    0          18011
  27  none Sun                17:44:05.249123 -23.20.15.67095 ICRS    0          18011
  28  none Sun                17:44:06.210181 -23.20.15.67681 ICRS    0          18011
  29  none Sun                17:44:07.171240 -23.20.15.68232 ICRS    0          18042
  30  none Sun                17:44:08.132299 -23.20.15.68749 ICRS    0          18011
  31  none Sun                17:44:09.093358 -23.20.15.69230 ICRS    0          18042
  32  none Sun                17:44:10.054417 -23.20.15.69676 ICRS    0          18011
  33  none Sun                17:44:11.015475 -23.20.15.70087 ICRS    0          18011
  34  none Sun                17:44:11.976534 -23.20.15.70464 ICRS    0          18011
  35  none Sun                17:44:12.937593 -23.20.15.70805 ICRS    0          18011
  36  none Sun                17:44:13.898652 -23.20.15.71111 ICRS    0          18042
  37  none Sun                17:44:04.059128 -23.20.04.51254 ICRS    0          18011
  38  none Sun                17:44:05.020164 -23.20.04.51893 ICRS    0          18011
  39  none Sun                17:44:05.981201 -23.20.04.52496 ICRS    0          18011
  40  none Sun                17:44:06.942238 -23.20.04.53065 ICRS    0          18011
  41  none Sun                17:44:07.903275 -23.20.04.53598 ICRS    0          18042
  42  none Sun                17:44:08.864312 -23.20.04.54097 ICRS    0          18011
  43  none Sun                17:44:09.825349 -23.20.04.54560 ICRS    0          18042
  44  none Sun                17:44:10.786386 -23.20.04.54989 ICRS    0          18011
  45  none Sun                17:44:11.747423 -23.20.04.55382 ICRS    0          18011
  46  none Sun                17:44:12.708460 -23.20.04.55741 ICRS    0          18011
  47  none Sun                17:44:13.669497 -23.20.04.56064 ICRS    0          18011
  48  none Sun                17:44:14.630534 -23.20.04.56353 ICRS    0          18042
  49  none Sun                17:44:04.791228 -23.19.53.36679 ICRS    0          18011
  50  none Sun                17:44:05.752243 -23.19.53.37300 ICRS    0          18011
  51  none Sun                17:44:06.713258 -23.19.53.37886 ICRS    0          18011
  52  none Sun                17:44:07.674273 -23.19.53.38437 ICRS    0          18011
  53  none Sun                17:44:08.635288 -23.19.53.38953 ICRS    0          18042
  54  none Sun                17:44:09.596303 -23.19.53.39433 ICRS    0          18011
  55  none Sun                17:44:10.557318 -23.19.53.39879 ICRS    0          18042
  56  none Sun                17:44:11.518333 -23.19.53.40290 ICRS    0          18011
  57  none Sun                17:44:12.479348 -23.19.53.40666 ICRS    0          18011
  58  none Sun                17:44:13.440364 -23.19.53.41007 ICRS    0          18011
  59  none Sun                17:44:14.401379 -23.19.53.41312 ICRS    0          18011
  60  none Sun                17:44:04.562314 -23.19.42.21455 ICRS    0          18042
  61  none Sun                17:44:05.523307 -23.19.42.22093 ICRS    0          18011
  62  none Sun                17:44:06.484300 -23.19.42.22696 ICRS    0          18011
  63  none Sun                17:44:07.445293 -23.19.42.23264 ICRS    0          18011
  64  none Sun                17:44:08.406286 -23.19.42.23797 ICRS    0          18011
  65  none Sun                17:44:09.367280 -23.19.42.24295 ICRS    0          18042
  66  none Sun                17:44:10.328273 -23.19.42.24758 ICRS    0          18011
  67  none Sun                17:44:11.289266 -23.19.42.25187 ICRS    0          18042
  68  none Sun                17:44:12.250259 -23.19.42.25580 ICRS    0          18011
  69  none Sun                17:44:13.211252 -23.19.42.25938 ICRS    0          18011
  70  none Sun                17:44:14.172245 -23.19.42.26261 ICRS    0          18011
  71  none Sun                17:44:15.133239 -23.19.42.26549 ICRS    0          18011
  72  none Sun                17:44:05.294393 -23.19.31.06875 ICRS    0          18042
  73  none Sun                17:44:06.255364 -23.19.31.07495 ICRS    0          18011
  74  none Sun                17:44:07.216335 -23.19.31.08081 ICRS    0          18011
  75  none Sun                17:44:08.177306 -23.19.31.08631 ICRS    0          18011
  76  none Sun                17:44:09.138278 -23.19.31.09146 ICRS    0          18011
  77  none Sun                17:44:11.060220 -23.19.31.10072 ICRS    0          18011
  78  none Sun                17:44:12.021191 -23.19.31.10482 ICRS    0          18042
  79  none Sun                17:44:12.982163 -23.19.31.10858 ICRS    0          18011
  80  none Sun                17:44:14.662815 -23.19.31.26419 ICRS    0          18042
  81  none Sun                17:44:15.623787 -23.19.31.26724 ICRS    0          18011
  82  none Sun                17:44:05.785179 -23.19.20.06866 ICRS    0          18042
  83  none Sun                17:44:06.746129 -23.19.20.07503 ICRS    0          18011
  84  none Sun                17:44:07.707078 -23.19.20.08105 ICRS    0          18011
  85  none Sun                17:44:08.668028 -23.19.20.08672 ICRS    0          18011
  86  none Sun                17:44:09.628977 -23.19.20.09204 ICRS    0          18011
  87  none Sun                17:44:10.589927 -23.19.20.09701 ICRS    0          18042
  88  none Sun                17:44:11.550877 -23.19.20.10163 ICRS    0          18011
  89  none Sun                17:44:12.511827 -23.19.20.10590 ICRS    0          18011
  90  none Sun                17:44:13.472776 -23.19.20.10982 ICRS    0          18011
  91  none Sun                17:44:14.433726 -23.19.20.11339 ICRS    0          18011
  92  none Sun                17:44:15.394676 -23.19.20.11661 ICRS    0          18042
  93  none Sun                17:44:16.355625 -23.19.20.11948 ICRS    0          18011
  94  none Sun                17:44:06.517236 -23.19.08.92273 ICRS    0          18042
  95  none Sun                17:44:07.478164 -23.19.08.92892 ICRS    0          18011
  96  none Sun                17:44:08.439092 -23.19.08.93476 ICRS    0          18011
  97  none Sun                17:44:09.400019 -23.19.08.94026 ICRS    0          18011
  98  none Sun                17:44:10.360947 -23.19.08.94540 ICRS    0          18011
  99  none Sun                17:44:11.321875 -23.19.08.95019 ICRS    0          18042
  100  none Sun                17:44:12.282803 -23.19.08.95464 ICRS    0          18011
  101  none Sun                17:44:13.243731 -23.19.08.95873 ICRS    0          18011
  102  none Sun                17:44:14.204658 -23.19.08.96247 ICRS    0          18011
  103  none Sun                17:44:15.165586 -23.19.08.96587 ICRS    0          18011
  104  none Sun                17:44:16.126514 -23.19.08.96891 ICRS    0          18042
  105  none Sun                17:44:06.288366 -23.18.57.77032 ICRS    0          18011
  106  none Sun                17:44:07.249271 -23.18.57.77669 ICRS    0          18042
  107  none Sun                17:44:08.210177 -23.18.57.78270 ICRS    0          18011
  108  none Sun                17:44:09.171083 -23.18.57.78837 ICRS    0          18011
  109  none Sun                17:44:10.131989 -23.18.57.79368 ICRS    0          18011
  110  none Sun                17:44:11.092895 -23.18.57.79865 ICRS    0          18011
  111  none Sun                17:44:12.053801 -23.18.57.80326 ICRS    0          18042
  112  none Sun                17:44:13.014707 -23.18.57.80753 ICRS    0          18011
  113  none Sun                17:44:13.975613 -23.18.57.81145 ICRS    0          18011
  114  none Sun                17:44:14.936519 -23.18.57.81501 ICRS    0          18011
  115  none Sun                17:44:15.897425 -23.18.57.81823 ICRS    0          18011
  116  none Sun                17:44:16.858331 -23.18.57.82110 ICRS    0          18042
  117  none Sun                17:44:07.020401 -23.18.46.62434 ICRS    0          18011
  118  none Sun                17:44:07.981285 -23.18.46.63053 ICRS    0          18042
  119  none Sun                17:44:08.942169 -23.18.46.63637 ICRS    0          18011
  120  none Sun                17:44:09.903053 -23.18.46.64185 ICRS    0          18011
  121  none Sun                17:44:10.863937 -23.18.46.64699 ICRS    0          18011
  122  none Sun                17:44:11.824821 -23.18.46.65178 ICRS    0          18011
  123  none Sun                17:44:12.785705 -23.18.46.65622 ICRS    0          18042
  124  none Sun                17:44:13.746589 -23.18.46.66031 ICRS    0          18011
  125  none Sun                17:44:14.707473 -23.18.46.66405 ICRS    0          18011
  126  none Sun                17:44:15.668357 -23.18.46.66744 ICRS    0          18011
  127  none Sun                17:44:16.629241 -23.18.46.67048 ICRS    0          18011
  128  none Sun                17:44:06.791552 -23.18.35.47189 ICRS    0          18042
  129  none Sun                17:44:07.752414 -23.18.35.47825 ICRS    0          18011
  130  none Sun                17:44:08.713276 -23.18.35.48426 ICRS    0          18042
  131  none Sun                17:44:09.674116 -23.18.35.48972 ICRS    0          18011
  132  none Sun                17:44:10.634936 -23.18.35.49465 ICRS    0          18011
  133  none Sun                17:44:11.595756 -23.18.35.49924 ICRS    0          18011
  134  none Sun                17:44:12.556576 -23.18.35.50347 ICRS    0          18011
  135  none Sun                17:44:13.517396 -23.18.35.50736 ICRS    0          18042
  136  none Sun                17:44:14.478216 -23.18.35.51089 ICRS    0          18011
  137  none Sun                17:44:15.439036 -23.18.35.51408 ICRS    0          18011
  138  none Sun                17:44:16.399856 -23.18.35.51691 ICRS    0          18011
  139  none Sun                17:44:17.360676 -23.18.35.51940 ICRS    0          18011
  140  none Sun                17:44:07.523164 -23.18.24.32227 ICRS    0          18042
  141  none Sun                17:44:08.483962 -23.18.24.32807 ICRS    0          18011
  142  none Sun                17:44:09.444760 -23.18.24.33353 ICRS    0          18042
  143  none Sun                17:44:10.405558 -23.18.24.33864 ICRS    0          18011
  144  none Sun                17:44:11.366356 -23.18.24.34339 ICRS    0          18011
  145  none Sun                17:44:12.327154 -23.18.24.34780 ICRS    0          18011
  146  none Sun                17:44:13.287952 -23.18.24.35186 ICRS    0          18011
  147  none Sun                17:44:14.248750 -23.18.24.35557 ICRS    0          18042
   148  none Sun                17:44:15.209549 -23.18.24.35893 ICRS    0          18011
   148  none Sun                17:44:15.209549 -23.18.24.35893 ICRS    0          18011
   149  none Sun                17:44:16.170347 -23.18.24.36194 ICRS    0          18011
   149  none Sun                17:44:16.170347 -23.18.24.36194 ICRS    0          18011
Line 441: Line 178:
====Spectral Window List====
====Spectral Window List====


* The fourth section shows the information of the spws in the dataset. From the scan list, we see that science (Target) scans are done with the spw 0~12. Of these, the spws with 128 channels are 5,7,9,11, which are used for image synthesis. '''Therefore, we will apply calibrations to spws 5,7,9,11 only.''' The data of spws 0,1,2,3 are the data from the square-law detectors of the basebands, which will be used for creating Tsys+Tant tables, and are archived as the auto-correlation data in the dataset.
The fourth section shows the information of the spws in the dataset. From the scan list, we see that science (Target) scans are done with the spw 0~12. Of these, the spws with 128 channels are 5,7,9,11, which are used for image synthesis. '''Therefore, we will apply calibrations to spws 5,7,9,11 only.''' The data of spws 0,1,2,3 are the data from the square-law detectors of the basebands, which will be used for creating Tsys+Tant tables, and are archived as the auto-correlation data in the dataset.


<pre style="background-color: #fffacd;">
<pre style="background-color: #fffacd;">
Line 485: Line 222:
</pre>
</pre>


* The fifth sections lists the spectral windows (spws) used for each Source. We exclude this for brevity.
The fifth sections lists the spectral windows (spws) used for each Source. We exclude this for brevity.


====Antenna List====
====Antenna List====


* The sixth section shows antenna information. Thirty-one antennas were used for the dataset. Note that numbering in python always begins with "0", so the antennas have IDs 0-30. To see what the antenna configuration looked like at the time of this observation, we use the task {{plotants_6.5.4}} (Figure 1):
The sixth section shows antenna information. Thirty-one antennas were used for the dataset. Note that numbering in python always begins with "0", so the antennas have IDs 0-30. To see what the antenna configuration looked like at the time of this observation, we use the task {{plotants_6.5.4}} (Figure 1):


<pre style="background-color: #fffacd;">
<pre style="background-color: #fffacd;">
Line 564: Line 301:
===Tsys calibration of the calibrators===
===Tsys calibration of the calibrators===


[[File:sunspot_tsys.png|thumb|right|'''Fig. 2.''' One example of the Tsys plots.]]
[[File:sunspot_tsys.png|400px|thumb|right|'''Fig. 2.''' One example of the Tsys plots.]]


The System Temperature (Tsys) calibration gives a first-order correction for the atmospheric opacity as a function of time and frequency, and associates weights with each visibility that persist through imaging. The MS dataset contains Tsys measurements; the task {{gencal_6.5.4}} is used to generate a calibration table. Since the dataset obtained with TDM, the data in the channels near the both edges of the spectrum window (~10 channels) are flagged with {{flagdata_6.5.4}}. The plots for checking are created by the ''checkCalTable'' subroutine of the Analysis Utilities package and saved to a new directory ''uid___A002_Xae00c5_X2a8d.ms.tsys.plots'' (Figure 2).
The System Temperature (Tsys) calibration gives a first-order correction for the atmospheric opacity as a function of time and frequency, and associates weights with each visibility that persist through imaging. The MS dataset contains Tsys measurements; the task {{gencal_6.5.4}} is used to generate a calibration table. Since the dataset obtained with TDM, the data in the channels near the both edges of the spectrum window (~10 channels) are flagged with {{flagdata_6.5.4}}. The plots for checking are created by the ''checkCalTable'' subroutine of the Analysis Utilities package and saved to a new directory ''uid___A002_Xae00c5_X2a8d.ms.tsys.plots'' (Figure 2).
Line 587: Line 324:
</source>
</source>


In the observations, we do not the atmospheric calibration of the phase calibrator between the scientific scans, because a long suspension of scientific observations has a bad influence on science. Therefore, we apply the Tsys calibration table, which is created from the data of the atmosphere calibration at the flux calibrator or the Sun, to the phase calibrator, as follows.
For normal ALMA observations, in between the science target scans, we would perform two adjacent scans of the phase calibrator: one with ATM intent, and one with PHASE intent. We would then apply the Tsys calibration table for the phase calibrator to itself.
 
<pre style="background-color: #E0FFFF;">
The Sun is very dynamic, and the time variations of solar structures are one of the main scientific topics of solar physics. To observe such structures, the suspension of science scans are problematic. Therefore, to return to the Sun as quickly as possible, we do NOT perform ATM scans of the phase calibrator in between science scans, instead only observing PHASE intent scans (9,12,15,18).
 
However, the antenna temperature (Tant) of the Sun is not negligible and must be calibrated. To do this, we perform ATM scans of the blank sky near the Sun, 2 degrees from the Sun's center. Field 0 is used for all calibration scans, and is the center of the mosaic, as we will see later (Fields 3~150 are used for science scans only). The first subscan of science scans (11,14,17) are listed as Field 0, but include the "OBSERVE_TARGET#OFF_SOURCE" intent, which slews to the blank sky.
 
In this Science Verification observation, one scan of the phase calibrator has ATM intent (scan 7) before any science scans begin, but more recent solar observations do not include this. Scan 7 can be ignored.
</pre>
 
'''For details, see [https://arxiv.org/abs/1704.03236v2 Shimojo+ 2017], in particular Figure 4.'''
 
Therefore, we calibrate the Tsys of the phase calibrator ''(field='2')'' by applying the calibrations for the Sun ''(gainfield='0')'', as follows:


<source lang="python">
<source lang="python">
#In CASA
#In CASA
applycal(vis=mso, field='2', spw='5,7,9,11', gaintable=mso+'.tsys', gainfield='2', interp='linear,linear', calwt=True, flagbackup=False)
applycal(vis=mso, field='2', spw='5,7,9,11', gaintable=mso+'.tsys', gainfield='0', interp='linear,linear', calwt=True, flagbackup=False)
</source>
</source>


Line 598: Line 347:
===Tsys+Tant calibration of the Sun===
===Tsys+Tant calibration of the Sun===


The standard method of Tsys calibration cannot be applied to the data of the Sun, because the antenna temperature (Tant) of the Sun cannot be neglected for estimating the system equivalent flux density (SEFD). To estimate correct SEFD at the Sun, the solar observing sequence includes some special measurements, like the measurement of zero-signal level. The subroutines for creating and applying Tsys+Tant calibration tables are prepared by the ALMA solar development team. They are part of the script we executed in the section [[#Import Tools and Scripts]]. For the Tsys+Tant calibration of the solar data, you will execute only the following command:
The standard method of Tsys calibration cannot be applied to the data of the Sun, because the antenna temperature (Tant) of the Sun cannot be neglected for estimating the system equivalent flux density (SEFD). To estimate correct SEFD at the Sun, the solar observing sequence includes some special measurements, like the measurement of zero-signal level (scan 1). The subroutines for creating and applying Tsys+Tant calibration tables are prepared by the ALMA solar development team. They are part of the script we executed in the section [[#Import Tools and Scripts]]. For the Tsys+Tant calibration of the solar data, you will execute only the following command:


<source lang="python">
<source lang="python">
Line 606: Line 355:


<font color="red">
<font color="red">
This process takes long time, about one night or day.
This process takes a long time. Depending on your machine, it may be between a couple hours and one day.
</font>
</font>


If you have already carried out the process before and there are the Tsys+Tant calibration tables (The directory name of the table includes “tsystant”), you can skip the generating of the tables using the following command.
If you have already carried out the process before and there are the Tsys+Tant calibration tables (The directory name of the table includes “tsystant”), you can skip the generating of the tables using the following command.
<font color="red">
Bug in Sun_reduction_util_6.5.4.py: underneath "if exisTbl == True:", ''scan = sciScan[i]'' should be ''scan = str(sciScan[i])''.
</font>


<source lang="python">
<source lang="python">
Line 617: Line 370:


===Flagging after Tsys and Tsys+Tant calibration===
===Flagging after Tsys and Tsys+Tant calibration===
Since we completed the Tsys and Tsys+Tant calibration, the data that are not used the image synthesis are flagged. At first, the data of auto-correlation and atmosphere calibration are flagged as follows:
Since we completed the Tsys and Tsys+Tant calibration, we can now flag the data that are not used for image synthesis. At first, the data of auto-correlation and atmosphere calibration are flagged with {{flagdata_6.5.4}} as follows:
<!-- unnecessary because we split out intents? -->


<source lang="python">
<source lang="python">
Line 625: Line 379:
</source>
</source>


Next is the spectrum windows that are not used in later:
Next are the spectral windows that are not used later:
<!-- unnecessary because we split out spws? -->


<source lang="python">
<source lang="python">
Line 633: Line 388:
</source>
</source>


The data in the channels near the both edges of the spectrum window (~10 channels) are flagged:
The channels near both edges of the science spectral windows (~10 channels) are flagged:
<!-- unnecessary because we already flagged these channels from the Tsys table before it was applied? -->


<source lang="python">
<source lang="python">
Line 640: Line 396:
</source>
</source>


The some sub-scans at the start and end of the scientific scans are used to measure the sky with the setting of optimized for the Sun. Since the data are used only for estimating the antenna temperatures at the Sun, we flag the data:
Some sub-scans at the start and end of the scientific scans are used to measure the sky with the settings optimized for the Sun. Since the data are used only for estimating the antenna temperatures at the Sun, we now flag the data:
<!-- unnecessary because we split out intents? -->


<source lang="python">
<source lang="python">
Line 653: Line 410:
       subNum = subInf[1]
       subNum = subInf[1]
       flagdata(vis=mso, timerange=subInf[3][1], mode='manual', flagbackup=False)       
       flagdata(vis=mso, timerange=subInf[3][1], mode='manual', flagbackup=False)       
       flagdata(vis=mso, timerange=subInf[3][2],mode='manual', flagbackup=False)       
       flagdata(vis=mso, timerange=subInf[3][2], mode='manual', flagbackup=False)       
       flagdata(vis=mso, timerange=subInf[3][subNum-1], mode='manual', flagbackup=False)       
       flagdata(vis=mso, timerange=subInf[3][subNum-1], mode='manual', flagbackup=False)       
       flagdata(vis=mso, timerange=subInf[3][subNum], mode='manual', flagbackup=False)
       flagdata(vis=mso, timerange=subInf[3][subNum], mode='manual', flagbackup=False)
</source>
</source>


As mentioned in the section [[#Tsys calibration of the visibilities of the calibrators]], the Tsys values of DA54 antennas are significantly large. Therefore, we flag the data of the antenna as follows:
As mentioned in the section [[#Tsys calibration of the calibrators]], the Tsys values of DA54 antennas are significantly large. Therefore, we flag the data of the antenna as follows:


<source lang="python">
<source lang="python">
Line 665: Line 422:
</source>
</source>


Then, we store the current flagging state for each dataset using the {{flagmanager_6.5.4}}:
Now we store the current flagging state using {{flagmanager_6.5.4}}:


<source lang="python">
<source lang="python">
Line 672: Line 429:
</source>
</source>


Now split out the CORRECTED data column, retaining spectral windows 5,7,9, and 11. This will get rid of the extraneous spectral windows.
Now {{split_6.5.4}} out the CORRECTED data column, retaining spectral windows 5,7,9,11. This will get rid of the extraneous spectral windows.


<source lang="python">
<source lang="python">
#In CASA
#In CASA
split(vis=mso, outputvis=mss, datacolumn='corrected', intent ='*BANDPASS*,*FLUX*,*PHASE*,*TARGET*', spw='5,7,9,11', keepflags=True)
split(vis=mso, outputvis=mss, datacolumn='corrected', intent ='*BANDPASS*,*FLUX*,*PHASE*,*TARGET*', spw='5,7,9,11')
</source>
</source>
'''CAUTION: When we run split, the spws that we choose to keep will be re-indexed to 0,1,2,3.'''


==Additional Data Inspection==
==Additional Data Inspection==


We will do some additional inspection.  
Now let's inspect the split MS with {{plotms_6.5.4}}.


The solar data frequently include zero values in some channels. To avoid their influence, we flag the data using the “clip” mode with “clipzeros=True” option.
The solar data frequently include zero values in some channels. This is best seen with all plot averaging turned off, though it may take {{plotms_6.5.4}} a moment to load. We can restrict the Y axis to a small number since we are looking for pure zero values.
 
[[Image:sunspot_zero_values.png|300px|thumb|right|'''Fig. 4.''' Amp vs Time. Pure zero values should be flagged.]]
 
<source lang="python">
#In CASA
plotms(vis=mss,
      xaxis='time',
      yaxis='amp',
      coloraxis='field',
      plotrange=[0,0,-0.1,0.1]
      )
</source>
 
To avoid their influence, we flag the data using the “clip” mode of {{flagdata_6.5.4}} with “clipzeros=True” option.


<source lang="python">
<source lang="python">
Line 690: Line 463:
</source>
</source>


The valleys appear in the plot of amplitude vs channel of the Spw #3. It shows the effect of the lines of the earth atmosphere. We need to flag the data.
Now let's plot Amp vs Channel for our two calibrators.
 
[[Image:sunspot_atm_lines.png|300px|thumb|right|'''Fig. 5.''' Amp vs Freq. Atmospheric lines should be flagged.]]
 
<source lang="python">
#In CASA
plotms(vis=mss,
      field='1,2',
      spw='3',
      xaxis='channel',
      yaxis='amp',
      avgantenna=True,
      iteraxis='field'
      )
</source>
 
Valleys appear in spw 3. It shows the effect of the lines of the Earth's atmosphere. We need to flag the data.


<source lang="python">
<source lang="python">
Line 702: Line 491:
<source lang="python">
<source lang="python">
#In CASA
#In CASA
flagmanager(vis=mss, mode='save', versionname='Initial')
flagmanager(vis=mss, mode='save', versionname='initial')
</source>
</source>


==Set up the Flux Calibration Model==
==Set up the Flux Calibration Model==
It is very rare that a useful planet for the flux calibration is located near the Sun. Therefore, we usually use a quasar near the Sun as a flux calibrator. In this observations, there is no quasar that flux is enough for the observations with the Mixer De-tuning mode, except nrao350. Therefore, the flux calibrator is the same as the phase calibrator in this observation.
It is very rare that a useful planet for the flux calibration is located near the Sun. Therefore, we usually use a quasar near the Sun as a flux calibrator. In these observations, there is no quasar with enough flux with the Mixer De-tuning mode except nrao530. Therefore, the flux calibrator is the same as the phase calibrator in these observation. We obtain the flux density and spectral index of nrao530 from the ALMA calibrator database.
To fill the model data column for nrao350 with a model for the flux density, we obtain the flux density and spectral index of the nro350 from the ALMA calibrator database.


<source lang="python">
<source lang="python">
#In CASA
#In CASA
intentSources=es.getIntentsAndSourceNames(mss)
intentSources = es.getIntentsAndSourceNames(mss)
ampCalId = intentSources['CALIBRATE_AMPLI']['id'] + intentSources['CALIBRATE_FLUX']['id']
ampCalId = intentSources['CALIBRATE_AMPLI']['id'] + intentSources['CALIBRATE_FLUX']['id']
calFieldNames = intentSources['CALIBRATE_AMPLI']['name'] + intentSources['CALIBRATE_FLUX']['name']
calFieldNames = intentSources['CALIBRATE_AMPLI']['name'] + intentSources['CALIBRATE_FLUX']['name']


amp_cal_name=calFieldNames[1]
amp_cal_name = calFieldNames[1]
spwInfo=es.getSpwInfo(mss)
spwInfo = es.getSpwInfo(mss)
obs_freq="%fGHz"%(spwInfo[0]['refFreq']/1e9)
obs_freq = "%fGHz"%(spwInfo[0]['refFreq']/1e9)


date=aU.getObservationStartDate(mss)
date = aU.getObservationStartDate(mss)
date_obs=date.split()[0]
date_obs = date.split()[0]
spw1_flux=aU.getALMAFlux(sourcename=amp_cal_name, date=date_obs,frequency=obs_freq)
spw1_flux = aU.getALMAFlux(sourcename=amp_cal_name, date=date_obs, frequency=obs_freq)
</source>
</source>


Then, fill the model data column for nrao350 with a model.
You will see the following output in terminal:
 
<pre style="background-color: #fffacd;">
Result using spectral index of -0.676873 for 231.000 GHz from 3.265 Jy at 97.500 GHz = 1.821049 +- 0.084762 Jy
</pre>
 
Then, fill the model data column for nrao530 by providing this model to {{setjy_6.5.4}}.


<source lang="python">
<source lang="python">
#In CASA
#In CASA
setjy(vis=mss, field='2', spw='0,1,2,3', standard='manual', fluxdensity=[spw1_flux['fluxDensity'], 0, 0, 0], spix= spw1_flux['spectralIndex'], reffreq=obs_freq)
setjy(vis=mss, field='2', spw='0,1,2,3', standard='manual', fluxdensity=[spw1_flux['fluxDensity'],0,0,0], spix=spw1_flux['spectralIndex'], reffreq=obs_freq, usescratch=True)
</source>
</source>
You will see the following output in the logger:
<pre style="background-color: #fffacd;">
Flux density as a function of frequency (channel 0 of each spw):
  Frequency (GHz)    Flux Density (Jy, Stokes I)
    230.992        1.82116
    232.992        1.81056
    245.008        1.74997
    247.008        1.74037
</pre>


==Creating the Bandpass Calibration Table==
==Creating the Bandpass Calibration Table==
''The bandpass and gain calibrations are the same as the calibrations of non-solar data basically. Please refer to the other ALMA tutorials of the CASA guide too. The calibration method presented below is the same as that used in the QA2 process of the ALMA observatory.''
The bandpass and gain calibrations are essentially the same as for non-solar data. First, we choose a reference antenna, then run {{gaincal_6.5.4}} on the bandpass calibrator to determine phase-only gain solutions. We will use ''solint='int''' for the solution interval, which means that one gain solution will be determined for every integration time.


At first, we determine the reference antenna, and run {{gaincal_6.5.4}} on the bandpass calibrator to determine phase-only gain solutions. We will use solint='int' for the solution interval, which means that one gain solution will be determined for every integration time.
[[Image:sunspot_pre-bandpass_gaincal.png|300px|thumb|right|'''Fig. 7.''' One example of the pre-bandpass gaincal solution plots. Note the phase wrapping between +/-180 deg.]]


<source lang="python">
<source lang="python">
#In CASA
#In CASA
ref_ant = 'DA41'
ref_ant = 'DA41'
gaincal(vis=mss, caltable=mss+'.ap_pre_bandpass', field='1', scan='5', solint='int', refant=ref_ant, calmode='p')
gaincal(vis=mss, caltable=mss+'.pre_bandpass', field='1', scan='5', solint='int', refant=ref_ant, calmode='p')
</source>
</source>


The plots are created for the checking.
Create plots to check the calibration table. These will be saved to a new directory ''uid___A002_Xae00c5_X2a8d.ms.split.pre_bandpass.plots''.


<source lang="python">
<source lang="python">
#In CASA
#In CASA
es.checkCalTable(mss+'.ap_pre_bandpass', msName=mss, interactive=False)
es.checkCalTable(mss+'.pre_bandpass', msName=mss, interactive=False)
</source>
</source>


The plots look good. Next, create the bandpass calibration table:  
The plots look good. Next, create the {{bandpass_6.5.4}} calibration table. Since we apply the per-integration phase solutions on-the-fly, we can use ''solint='inf''' to average the scan in time, resulting in higher signal-to-noise for our bandpass solution. You will see many warning reported about "insufficient unflagged antennas," which are for the edge channels we flagged so they are safe to ignore.
 
[[Image:sunspot_bandpass.png|300px|thumb|right|'''Fig. 7.''' One example of the bandpass solution plots.]]


<source lang="python">
<source lang="python">
#In CASA
#In CASA
bandpass(vis=mss, caltable=mss+'.bandpass', field='1', scan='5', solint='inf', refant=ref_ant, solnorm=True, bandtype='B', gaintable=mss+'.ap_pre_bandpass')
bandpass(vis=mss, caltable=mss+'.bandpass', field='1', scan='5', solint='inf', refant=ref_ant, solnorm=True, bandtype='B', gaintable=mss+'.pre_bandpass')
</source>
</source>


Check the bandpass calibration table using the plots created from the following command.
Check the bandpass calibration table using the plots created from the following command, which are saved into a new directory ''uid___A002_Xae00c5_X2a8d.ms.split.bandpass.plots''.


<source lang="python">
<source lang="python">
Line 763: Line 570:
</source>
</source>


==Creating the Gain Calibration Table==
<br clear="all">
At first, we determine phase-only gain solutions of the calibrators, using the bandpass calibration table and solint='int' option, and create the plots for checking.
 
==Creating the Gain Calibration Tables==
At first, we determine phase-only gain solutions of the calibrators with {{gaincal_6.5.4}}, applying the bandpass calibration table on-the-fly, using ''solint='int''' for per-integration solutions. Then create the plots for checking, saved to ''uid___A002_Xae00c5_X2a8d.ms.split.phase_int.plots''.


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
gaincal(vis=mss, caltable=mss+'.phase_int', field='1~2', solint='int', refant=ref_ant,gaintype='G',calmode='p',minsnr=3.0, gaintable=mss+'.bandpass')
gaincal(vis=mss, caltable=mss+'.phase_int', field='1~2', solint='int', refant=ref_ant, gaintype='G', calmode='p', minsnr=3.0, gaintable=mss+'.bandpass')
 
es.checkCalTable(mss+'.phase_int', msName=mss, interactive=False)
es.checkCalTable(mss+'.phase_int', msName=mss, interactive=False)
</source>
</source>


Using the bandpass and phase calibration tables, we obtain the amplitude-only gain solutions on the scan time, '''solint='inf''''.
Now apply both the bandpass and short phase solution tables on-the-fly to solve for amplitude-only gain solutions over the full scan time with ''solint='inf'''. Plots are saved to ''uid___A002_Xae00c5_X2a8d.ms.split.ampli_inf.plots''.


<source lang='python'>
<source lang='python'>
# In CASA
# In CASA
gaincal(vis=mss, caltable=mss+'.ampli_inf', field='1~2', solint='inf', refant=ref_ant, gaintype='G', calmode='a', minsnr=3.0, gaintable=[mss+'.bandpass', mss+'.phase_int'])
gaincal(vis=mss, caltable=mss+'.ampli_inf', field='1~2', solint='inf', refant=ref_ant, gaintype='G', calmode='a', minsnr=3.0, gaintable=[mss+'.bandpass', mss+'.phase_int'])
es.checkCalTable(mss+'.ampli_inf', msName=mss, interactive=False)  
es.checkCalTable(mss+'.ampli_inf', msName=mss, interactive=False)  
</source>
</source>


Usually, the gaintype is set to 'T' for nominal ALMA observations . Nevetheless, the gaintype is set to 'G' for solar data, because the diffrence between XX and YY images is used for for estimating the noise level (see Section 4.2, Shimojo et al. 2017).
Usually, the ''gaintype'' is set to 'T' for nominal ALMA observations. '''For solar data, the ''gaintype'' is set to 'G', because the difference between the XX and YY images is used for estimating the noise level''' (see Section 4.2, Shimojo et al. 2017).


Next we use the flux calibrator (whose flux density was set in {{setjy_6.5.4}} above) to derive the flux density of the other calibrators. Note that the flux table REPLACES the amp.gcal in terms of future application of the calibration to the data, i.e. the flux table contains both the amp.gcal and flux scaling. Unlike the gain calibration steps, this is not an incremental table.
Next we use the flux calibrator (whose flux density was set in {{setjy_6.5.4}} above) to derive the flux density of the other calibrators with {{fluxscale_6.5.4}}. Note that the flux table REPLACES '.ampli_inf' in terms of future application of the calibration to the data, i.e., the flux table contains both the amplitude calibration and flux scaling. Unlike the two {{gaincal_6.5.4}} steps above, this is not an incremental table. Plots are saved to ''uid___A002_Xae00c5_X2a8d.ms.split.flux_inf.plots''.
 
[[Image:sunspot_flux_inf.png|300px|thumb|right|'''Fig. 7.''' Example of the flux-scaled amplitude calibration table.]]


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
fluxscaleDict = fluxscale(vis=mss, caltable=mss+'.ampli_inf', fluxtable=mss+'.flux_inf', reference='2')  
fluxscaleDict = fluxscale(vis=mss, caltable=mss+'.ampli_inf', fluxtable=mss+'.flux_inf', reference='2')
es.fluxscale2(caltable=mss+'.ampli_inf', removeOutliers=True, msName=mss, writeToFile=True, preavg=10000)
 
es.checkCalTable(mss+'.flux_inf', msName=mss, interactive=False)
</source>
</source>


Finally, we create the gain calibration table of phase-only gain solutions on the scan time, '''solint='inf''''.
You should see the following output:
 
<pre style="background-color: #fffacd;">
Found reference field(s): nrao530
Found transfer field(s):  J1924-2914
Flux density for J1924-2914 in SpW=0 (freq=2.3e+11 Hz) is: 3.31068 +/- 0.0174892 (SNR = 189.299, N = 60)
Flux density for J1924-2914 in SpW=1 (freq=2.32e+11 Hz) is: 3.30137 +/- 0.016152 (SNR = 204.394, N = 60)
Flux density for J1924-2914 in SpW=2 (freq=2.46e+11 Hz) is: 3.16067 +/- 0.0180585 (SNR = 175.024, N = 60)
Flux density for J1924-2914 in SpW=3 (freq=2.48e+11 Hz) is: 3.14243 +/- 0.0205412 (SNR = 152.981, N = 60)
Fitted spectrum for J1924-2914 with fitorder=1: Flux density = 3.22814 +/- 0.002994 (freq=238.864 GHz) spidx: a_1 (spectral index) =-0.712504 +/- 0.0278242 covariance matrix for the fit:  covar(0,0)=7.85242e-06 covar(0,1)=9.80161e-05 covar(1,0)=9.80161e-05 covar(1,1)=0.0374698
</pre>
 
Finally, we create the gain calibration table of phase-only gain solutions on the scan time, ''solint='inf'''. Plots are saved to ''uid___A002_Xae00c5_X2a8d.ms.split.phase_inf.plots''.
 
[[Image:sunspot_phase_inf.png|300px|thumb|right|'''Fig. 7.''' Example of the long solution interval phase calibration table.]]


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
gaincal(vis=mss, caltable=mss+'.phase_inf', field='1~2', solint='inf', refant=ref_ant, gaintype='G', calmode='p', minsnr=3.0, gaintable=mss+'.bandpass')
gaincal(vis=mss, caltable=mss+'.phase_inf', field='1~2', solint='inf', refant=ref_ant, gaintype='G', calmode='p', minsnr=3.0, gaintable=mss+'.bandpass')
es.checkCalTable(mss+'.phase_inf', msName=mss, interactive=False)
es.checkCalTable(mss+'.phase_inf', msName=mss, interactive=False)
</source>
</source>
<br clear="all">


==Applying the Calibration Tables==
==Applying the Calibration Tables==
We apply the calibration solutions to each source individually, using the gainfield parameter to specify which calibrator's solutions should be applied from each of the gaintable calibration tables.
We apply the calibration solutions to each source individually with {{applycal_6.5.4}}, using the ''gainfield'' parameter to specify which calibrator's solutions should be applied from each of the calibration tables. We can leave ''gainfield'' blank when referencing the bandpass table, because this table has information derived only from the bandpass calibrator which will be applied to all sources.
Applying the tables to the bandpass, flux and phase calibrators:
 
Apply to the bandpass calibrator: the bandpass table; the short ''solint'' phase table, referencing the BP cal; and the flux table, referencing the BP cal.


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
for i in ['1', '2']:
applycal(vis=mss, field='1', gaintable=[mss+'.bandpass', mss+'.phase_int', mss+'.flux_inf'], gainfield=['', '1', '1'],
      applycal(vis=mss, field=str(i), gaintable=[mss+'.bandpass', mss+'.phase_int', mss+'.flux_inf'], gainfield=['', i, i], interp='linear,linear', calwt=True, flagbackup=False)
interp='linear,linear', calwt=True, flagbackup=False)
</source>
</source>


Next, applying them to the solar data:
Apply to the phase calibrator: the bandpass table; the short ''solint'' phase table, referencing the phase cal; and the flux table, referencing the phase cal.


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
applycal(vis=mss, field='0, 3~150', gaintable=[mss+'.bandpass', mss+'.phase_inf', mss+'.flux_inf'], gainfield=['', '2', ''], interp='linear,linear', calwt=True, flagbackup=False)
applycal(vis=mss, field='2', gaintable=[mss+'.bandpass', mss+'.phase_int', mss+'.flux_inf'], gainfield=['', '2', '2'],
interp='linear,linear', calwt=True, flagbackup=False)
</source>
 
Apply to the Sun: the bandpass table; the long ''solint'' phase table, referencing the phase cal; and the flux table, referencing the phase cal.
 
<source lang='python'>
#In CASA
applycal(vis=mss, field='0, 3~150', gaintable=[mss+'.bandpass', mss+'.phase_inf', mss+'.flux_inf'], gainfield=['', '2', '2'],
interp='linear,linear', calwt=True, flagbackup=False)
</source>
</source>


Line 819: Line 660:
<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
split(vis=mss, outputvis=msc, datacolumn='corrected', keepflags=True)
split(vis=mss, outputvis=msc, datacolumn='corrected')
</source>
</source>


Line 825: Line 666:


[[File:sunspot_mos1.png|thumb|right|'''Fig. 3.''' The pattern of mosaic BEFORE the re-calculation of the direction.]]
[[File:sunspot_mos1.png|thumb|right|'''Fig. 3.''' The pattern of mosaic BEFORE the re-calculation of the direction.]]
[[File:sunspot_mos2.png|thumb|right|'''Fig. 4.''' The pattern of mosaic AFTER the re-calculation of the direction.]]
[[File:sunspot_mosaic_after_CASA_6.5.4.png|thumb|right|'''Fig. 4.''' The pattern of mosaic AFTER the re-calculation of the direction.]]


During most solar observations, the antennas are tracking a structure on the Sun according to the solar differential rotation. The image frame is fixed on the solar frame, but the frame is moving on the RA/Dec coordinate frame. If we do not do any measures, the pattern of pointing in the MOSAIC observation is shown in Figure 3 , and the shape of the pattern is a rhombus though the correct shape is the square.
During most solar observations, the antennas are tracking a structure on the Sun according to the solar differential rotation. The image frame is fixed on the solar frame, but the frame is moving on the RA/Dec coordinate frame. If we do not correct for this, the pattern of pointing in the mosaic is a rhombus as shown in Figure 3, while the correct shape should be a square.
 
Plot the mosaic to see this:<br>
<font color="red">
The following command does not correctly create the plot as shown in Figure 3 at this time. Figure 3 is taken from a previous version of this guide.
</font>
 
<source lang='python'>
#In CASA
au.plotmosaic(vis=msc, sourceid='Sun', figfile=msc+'.pointings.sun.before.png')
</source>


To correct the MOSAIC pattern, we re-calculate the coordinate of each field. At first, we modify the coordinate of the field “0” from the reference time using {{fixplanets_6.5.4}} task. The reference time has to be the time when the antennas are directed to the field “0”.
To correct the MOSAIC pattern, we re-calculate the coordinate of each field. First, we modify the coordinate of ''field='0' ''from the reference time using {{fixplanets_6.5.4}} task. The reference time has to be the time when the antennas are directed to ''field='0'''.


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
reftime = '2015/12/18/19:49:00'
reftime = '2015/12/18/19:49:00'
fixplanets(msc, field='0', fixuvw=False, refant=ref_ant, reftime=reftime)
fixplanets(vis=msc, field='0', fixuvw=False, refant=ref_ant, reftime=reftime)
</source>
</source>


We define that the modified coordinate of field ’0’ is as the reference coordinate, and re-calculate the coordinate of each field, as follows.
We define that the modified coordinate of ''field='0''' is the reference coordinate, and re-calculate the coordinate of each field, as follows:


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
import math
pi = math.pi
tb.open(msc+'/FIELD', nomodify=True)
tb.open(msc+'/FIELD', nomodify=True)
phsCenOff = tb.getcol("PHASE_DIR")
phsCenOff = tb.getcol("PHASE_DIR")
tb.close()
tb.close()
refRaDec = aU.rad2radec(phsCenOff[0][0][0], phsCenOff[1][0][0], prec=1, hmsdms=True, delimiter=' ')   
 
refRaDec = aU.rad2radec(phsCenOff[0][0][0], phsCenOff[1][0][0], prec=1, hmsdms=True, delimiter=' ')   
for i in range(3, 151):
for i in range(3, 151):
     raOff = phsCenOff[0][0][i] * 180./pi * 60. *60.
     raOff = phsCenOff[0][0][i] * 180. / pi * 60. *60.
     deOff = phsCenOff[1][0][i] * 180./pi * 60. *60.
     deOff = phsCenOff[1][0][i] * 180. / pi * 60. *60.
     offRaDec = aU.radec2deg(aU.radecOffsetToRadec(refRaDec, raOff, deOff, prec=1))
     offRaDec = aU.radec2deg(aU.radecOffsetToRadec(refRaDec, raOff, deOff, prec=1))
     offRaDecF = 'J2000 ' + aU.deg2radec(offRaDec[0], offRaDec[1], prec=1, hmsdms=True, delimiter=' ')         
     offRaDecF = 'J2000 ' + aU.deg2radec(offRaDec[0], offRaDec[1], prec=1, hmsdms=True, delimiter=' ')         
     fixplanets(msc, field =str(i), fixuvw=False, direction=offRaDecF)
     fixplanets(vis=msc, field=str(i), fixuvw=False, direction=offRaDecF)


tb.open(msc+'/FIELD', nomodify=False)
tb.open(msc+'/FIELD', nomodify=False)
Line 861: Line 716:
</source>
</source>


Moreover, the direction in the pointing table has a bad influence to the coordinate system of the image synthesis. We erase the pointing table as follows:
Moreover, the direction in the pointing table has a bad influence on the coordinate system for image synthesis. We erase the pointing table as follows:


<source lang='python'>
<source lang='python'>
#In CASA
#In CASA
tb.open(msc+'/POINTING', nomodify = False)
tb.open(msc+'/POINTING', nomodify=False)
a = tb.rownumbers()
a = tb.rownumbers()
tb.removerows(a)
tb.removerows(a)
tb.close()
</source>
</source>


As a result, the pattern of the pointing in the MOSAIC observation becomes the map, as shown in Figure 4.
Now plot the corrected mosaic:
 
<source lang='python'>
#In CASA
au.plotmosaic(vis=msc, sourceid='Sun', figfile=msc+'.pointings.sun.after.png')
</source>

Latest revision as of 12:35, 22 April 2024

Last checked on CASA Version 6.5.4

Getting Started

Details of these ALMA observations are provided at Sunspot_Band6. This portion of the guide will cover the calibration of the raw visibility data.

WARNING: The command in the #Tsys+Tant calibration of the Sun section of this guide takes a long time to run. Depending on your machine, it may be between a couple hours and one day.

Download and Unpack the Data

To follow this guide you must download the file Sunspot_Band6_UncalibratedData.tgz from Sunspot_Band6#Obtaining the Data.

Unpack and cd to the directory:

#In bash
tar -xvzf Sunspot_Band6_UncalibratedData.tgz
cd Sunspot_Band6_UncalibratedData

Confirm your version of CASA

This guide has been written for CASA release 6.5. Start CASA and confirm your version before proceeding.

# In CASA
from casatools import version
vernum = str(version()[0])+'.'+str(version()[1])
print("You are using CASA ver. "+vernum)
if float(vernum) < 6.5:
 print("YOUR VERSION OF CASA IS TOO OLD FOR THIS GUIDE.")
 print("PLEASE UPDATE IT BEFORE PROCEEDING.")
else:
 print("Your version of CASA is appropriate for this guide.")

Import Tools and Scripts

The "Analysis Utilities" package must be used for the calibration of raw solar visibility data. Therefore, before starting the tutorial, you need to import the package to your data-analysis environment. The documents and software of the package can be obtained from the Analysis_Utilities guide.

We also need script to download the script Sun_reduction_util_6.5.4.py and execute it inside our CASA session.

#In CASA
import analysisUtils as aU
es = aU.stuffForScienceDataReduction() 
execfile('Sun_reduction_util_6.5.4.py')

Initial Inspection

Import the Data into CASA

We start by defining the directory name of the ASDM and some directory names of the Measurement Sets (MS) for the calibration.

#In CASA
asdm = 'uid___A002_Xae00c5_X2a8d' 
mso = asdm + '.ms'           #MS original
mss = asdm + '.ms.split'     #MS split (science spws)
msc = mss  + '.cal'          #MS calibrated

The raw data have been provided to you in the ASDM format. It is the native format of the data produced by the ALMA observatory.

Before we can proceed to the calibration, we will need to convert those data to the CASA MS format. This is done simply with the task importasdm.

#In CASA
importasdm(asdm=asdm, vis=mso, asis='Antenna Station Receiver Source CalAtmosphere CalWVR CorrelatorMode SBSummary CalDevice')

To check if a bug fix is needed in the SYSCAL table times, the following commands are executed.

#In CASA
from casarecipes.almahelpers import fixsyscaltimes
fixsyscaltimes(vis=mso)

Inspect the Data

The usual first step is then to get some basic information about the data. We do this using the task listobs, which will output a detailed summary of each dataset supplied.

#In CASA
listobs(mso, listfile=asdm+'.listobs.txt')

The output will be sent to the CASA logger, or saved in a text file.

The first section shows information about the project and observing time:

   Observer: shimojo     Project: uid://A002/Xac494e/X3  
Observation: ALMA
Data records: 20923884       Total elapsed time = 3184.8 seconds
   Observed from   18-Dec-2015/19:15:42.3   to   18-Dec-2015/20:08:47.1 (UTC)

Scan List

The second section describes details of each scan. This is often easier to read in a text viewer with line wrapping turned off. We see that our Fields have the following associated intents:

  • Field 0 = Sun: observed for the measurement of zero-signal level [scan 1] and the calibration of atmosphere [scans 10,13,16]. Scientific observations of the Sun (Target) [scans 11,14,17] are listed many times, which represent a mosaic.
  • Field 1 = J1924-2914: observed for the calibrations of pointing [scan 2], sideband ratio [scan 3], atmosphere [scan 4], and bandpass [scan 5]
  • Field 2 = nrao530: observed for the calibration of pointing [scan 6], atmosphere [scan 7], flux [scan 8], and phase [scans 9,12,15,18]

Subscans are NOT necessarily listed in chronological order, but rather ordered by Field ID. For example, for scan 14, Fields 3~9 are observed AFTER Fields 80~150, but listed first. See Sunspot_Band6_Imaging_for_CASA_6.5.4#Flag_the_surplus_subscans_and_baselines.

   ObservationID = 0         ArrayID = 0
  Date        Timerange (UTC)          Scan  FldId FieldName             nRows     SpwIds   Average Interval(s)    ScanIntent
  18-Dec-2015/19:15:42.3 - 19:16:47.5     1      0 Sun                     236530  [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.016, 0.016, 0.016, 0.016, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_ATMOSPHERE#AMBIENT,CALIBRATE_ATMOSPHERE#HOT,CALIBRATE_ATMOSPHERE#OFF_SOURCE,CALIBRATE_ATMOSPHERE#REFERENCE,CALIBRATE_WVR#AMBIENT,CALIBRATE_WVR#HOT,CALIBRATE_WVR#OFF_SOURCE,CALIBRATE_WVR#REFERENCE]
              19:17:31.8 - 19:19:26.4     2      1 J1924-2914             1081497  [4,25,26,27,28,29,30,31,32,33,34,35,36]  [1.15, 0.016, 0.016, 0.016, 0.016, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_POINTING#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              19:20:14.2 - 19:21:17.7     3      1 J1924-2914              895435  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_SIDEBAND_RATIO#IMAGE,CALIBRATE_SIDEBAND_RATIO#SIGNAL,CALIBRATE_WVR#IMAGE,CALIBRATE_WVR#SIGNAL]
              19:22:06.0 - 19:22:23.1     4      1 J1924-2914              177382  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_ATMOSPHERE#AMBIENT,CALIBRATE_ATMOSPHERE#HOT,CALIBRATE_ATMOSPHERE#OFF_SOURCE,CALIBRATE_WVR#AMBIENT,CALIBRATE_WVR#HOT,CALIBRATE_WVR#OFF_SOURCE]
              19:23:11.6 - 19:28:27.9     5      1 J1924-2914             3244522  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_BANDPASS#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              19:29:14.5 - 19:31:09.2     6      2 nrao530                1081528  [4,25,26,27,28,29,30,31,32,33,34,35,36]  [1.15, 0.016, 0.016, 0.016, 0.016, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_POINTING#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              19:31:54.6 - 19:32:11.4     7      2 nrao530                 177413  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_ATMOSPHERE#AMBIENT,CALIBRATE_ATMOSPHERE#HOT,CALIBRATE_ATMOSPHERE#OFF_SOURCE,CALIBRATE_WVR#AMBIENT,CALIBRATE_WVR#HOT,CALIBRATE_WVR#OFF_SOURCE]
              19:32:59.1 - 19:35:36.6     8      2 nrao530                1622261  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_FLUX#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              19:36:21.9 - 19:36:52.1     9      2 nrao530                 324446  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_PHASE#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              19:37:41.4 - 19:37:57.8    10      0 Sun                     177413  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_ATMOSPHERE#AMBIENT,CALIBRATE_ATMOSPHERE#HOT,CALIBRATE_ATMOSPHERE#OFF_SOURCE,CALIBRATE_WVR#AMBIENT,CALIBRATE_WVR#HOT,CALIBRATE_WVR#OFF_SOURCE]
              19:39:23.9 - 19:49:36.0    11      0 Sun                    3760610  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#OFF_SOURCE,CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#OFF_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11      3 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11      4 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11      5 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]

*** Scan 11, Fields 6~76 trimmed for brevity ***

              19:39:23.9 - 19:49:36.0    11     77 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11     78 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:39:23.9 - 19:49:36.0    11     79 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:50:21.7 - 19:50:52.0    12      2 nrao530                 324446  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_PHASE#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              19:51:41.2 - 19:51:58.0    13      0 Sun                     177413  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_ATMOSPHERE#AMBIENT,CALIBRATE_ATMOSPHERE#HOT,CALIBRATE_ATMOSPHERE#OFF_SOURCE,CALIBRATE_WVR#AMBIENT,CALIBRATE_WVR#HOT,CALIBRATE_WVR#OFF_SOURCE]
              19:53:26.0 - 20:03:37.0    14      0 Sun                    3742537  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#OFF_SOURCE,CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#OFF_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      3 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      4 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      5 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      6 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      7 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      8 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14      9 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14     80 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14     81 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14     82 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]

*** Scan 14, Fields 83~147 trimmed for brevity ***

              19:53:26.0 - 20:03:37.0    14    148 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    149 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              19:53:26.0 - 20:03:37.0    14    150 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              20:04:22.7 - 20:04:53.7    15      2 nrao530                 324446  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_PHASE#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
              20:05:43.3 - 20:06:00.2    16      0 Sun                     177413  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576, 0.576] [CALIBRATE_ATMOSPHERE#AMBIENT,CALIBRATE_ATMOSPHERE#HOT,CALIBRATE_ATMOSPHERE#OFF_SOURCE,CALIBRATE_WVR#AMBIENT,CALIBRATE_WVR#HOT,CALIBRATE_WVR#OFF_SOURCE]
              20:06:47.8 - 20:07:30.5    17      0 Sun                     227168  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#OFF_SOURCE,CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#OFF_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              20:06:47.8 - 20:07:30.5    17     10 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              20:06:47.8 - 20:07:30.5    17     11 Sun                      18042  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              20:06:47.8 - 20:07:30.5    17     12 Sun                      18011  [4,5,6,7,8,9,10,11,12]  [1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_WVR#ON_SOURCE,OBSERVE_TARGET#ON_SOURCE]
              20:08:16.5 - 20:08:47.1    18      2 nrao530                 324446  [0,1,2,3,4,5,6,7,8,9,10,11,12]  [0.016, 0.016, 0.016, 0.016, 1.15, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01, 2.02, 1.01] [CALIBRATE_PHASE#ON_SOURCE,CALIBRATE_WVR#ON_SOURCE]
           (nRows = Total number of rows per scan) 

Field List

The third section shows the information of the observing targets. This shows that three Sources with 151 Fields were observed: The Sun, J1924-2914, and nrao530 (J1733-1304). The Sun as one Source has many Fields - this is how a mosaic is represented:

Fields: 151
  ID   Code Name                RA               Decl           Epoch   SrcId      nRows
  0    none Sun                 17:44:04.359217 -23.19.29.92535 ICRS    0        8499084
  1    none J1924-2914          19:24:51.055957 -29.14.30.12103 ICRS    1        5398836
  2    none nrao530             17:33:02.705760 -13.04.49.54800 ICRS    2        4178986
  3    none Sun                 17:44:03.784899 -23.20.37.96258 ICRS    0          36022
  4    none Sun                 17:44:04.746002 -23.20.37.96880 ICRS    0          36053
  5    none Sun                 17:44:05.707104 -23.20.37.97466 ICRS    0          36053

*** Sun Fields 6~147 trimmed for brevity ***

  148  none Sun                 17:44:15.209549 -23.18.24.35893 ICRS    0          18011
  149  none Sun                 17:44:16.170347 -23.18.24.36194 ICRS    0          18011
  150  none Sun                 17:44:17.131145 -23.18.24.36460 ICRS    0          18011

Spectral Window List

The fourth section shows the information of the spws in the dataset. From the scan list, we see that science (Target) scans are done with the spw 0~12. Of these, the spws with 128 channels are 5,7,9,11, which are used for image synthesis. Therefore, we will apply calibrations to spws 5,7,9,11 only. The data of spws 0,1,2,3 are the data from the square-law detectors of the basebands, which will be used for creating Tsys+Tant tables, and are archived as the auto-correlation data in the dataset.

Spectral Windows:  (37 unique spectral windows and 2 unique polarization setups)
  SpwID  Name                                       #Chans   Frame   Ch0(MHz)  ChanWid(kHz)  TotBW(kHz) CtrFreq(MHz) BBC Num  Corrs  
  0      BB_1#SQLD                                       1   TOPO  230000.000   2000000.000   2000000.0 230000.0000        1  XX  YY
  1      BB_2#SQLD                                       1   TOPO  232000.000   2000000.000   2000000.0 232000.0000        2  XX  YY
  2      BB_3#SQLD                                       1   TOPO  246000.000   2000000.000   2000000.0 246000.0000        3  XX  YY
  3      BB_4#SQLD                                       1   TOPO  248000.000   2000000.000   2000000.0 248000.0000        4  XX  YY
  4      WVR#NOMINAL                                     4   TOPO  184550.000   1500000.000   7500000.0 187550.0000        0  XX
  5      X241538345#ALMA_RB_06#BB_1#SW-01#FULL_RES     128   TOPO  230992.188    -15625.000   2000000.0 230000.0000        1  XX  YY
  6      X241538345#ALMA_RB_06#BB_1#SW-01#CH_AVG         1   TOPO  229984.375   1796875.000   1796875.0 229984.3750        1  XX  YY
  7      X241538345#ALMA_RB_06#BB_2#SW-01#FULL_RES     128   TOPO  232992.188    -15625.000   2000000.0 232000.0000        2  XX  YY
  8      X241538345#ALMA_RB_06#BB_2#SW-01#CH_AVG         1   TOPO  231984.375   1796875.000   1796875.0 231984.3750        2  XX  YY
  9      X241538345#ALMA_RB_06#BB_3#SW-01#FULL_RES     128   TOPO  245007.813     15625.000   2000000.0 246000.0000        3  XX  YY
  10     X241538345#ALMA_RB_06#BB_3#SW-01#CH_AVG         1   TOPO  245984.375   1796875.000   1796875.0 245984.3750        3  XX  YY
  11     X241538345#ALMA_RB_06#BB_4#SW-01#FULL_RES     128   TOPO  247007.813     15625.000   2000000.0 248000.0000        4  XX  YY
  12     X241538345#ALMA_RB_06#BB_4#SW-01#CH_AVG         1   TOPO  247984.375   1796875.000   1796875.0 247984.3750        4  XX  YY
  13     BB_1#SQLD                                       1   TOPO  219559.000   2000000.000   2000000.0 219559.0000        1  XX  YY
  14     BB_2#SQLD                                       1   TOPO  219559.000   2000000.000   2000000.0 219559.0000        2  XX  YY
  15     BB_3#SQLD                                       1   TOPO  219559.000   2000000.000   2000000.0 219559.0000        3  XX  YY
  16     BB_4#SQLD                                       1   TOPO  219559.000   2000000.000   2000000.0 219559.0000        4  XX  YY
  17     X241538345#ALMA_RB_06#BB_1#SW-01#FULL_RES     128   TOPO  220551.188    -15625.000   2000000.0 219559.0000        1  XX  YY
  18     X241538345#ALMA_RB_06#BB_1#SW-01#CH_AVG         1   TOPO  219543.375   1796875.000   1796875.0 219543.3750        1  XX  YY
  19     X241538345#ALMA_RB_06#BB_2#SW-01#FULL_RES     128   TOPO  220551.188    -15625.000   2000000.0 219559.0000        2  XX  YY
  20     X241538345#ALMA_RB_06#BB_2#SW-01#CH_AVG         1   TOPO  219543.375   1796875.000   1796875.0 219543.3750        2  XX  YY
  21     X241538345#ALMA_RB_06#BB_3#SW-01#FULL_RES     128   TOPO  220551.188    -15625.000   2000000.0 219559.0000        3  XX  YY
  22     X241538345#ALMA_RB_06#BB_3#SW-01#CH_AVG         1   TOPO  219543.375   1796875.000   1796875.0 219543.3750        3  XX  YY
  23     X241538345#ALMA_RB_06#BB_4#SW-01#FULL_RES     128   TOPO  220551.188    -15625.000   2000000.0 219559.0000        4  XX  YY
  24     X241538345#ALMA_RB_06#BB_4#SW-01#CH_AVG         1   TOPO  219543.375   1796875.000   1796875.0 219543.3750        4  XX  YY
  25     BB_1#SQLD                                       1   TOPO  221538.000   2000000.000   2000000.0 221538.0000        1  XX  YY
  26     BB_2#SQLD                                       1   TOPO  223538.000   2000000.000   2000000.0 223538.0000        2  XX  YY
  27     BB_3#SQLD                                       1   TOPO  237538.000   2000000.000   2000000.0 237538.0000        3  XX  YY
  28     BB_4#SQLD                                       1   TOPO  239538.000   2000000.000   2000000.0 239538.0000        4  XX  YY
  29     X0000000000#ALMA_RB_06#BB_1#SW-01#FULL_RES    128   TOPO  222530.188    -15625.000   2000000.0 221538.0000        1  XX  YY
  30     X0000000000#ALMA_RB_06#BB_1#SW-01#CH_AVG        1   TOPO  221514.562   1781250.000   1781250.0 221514.5625        1  XX  YY
  31     X0000000000#ALMA_RB_06#BB_2#SW-01#FULL_RES    128   TOPO  224530.188    -15625.000   2000000.0 223538.0000        2  XX  YY
  32     X0000000000#ALMA_RB_06#BB_2#SW-01#CH_AVG        1   TOPO  223514.562   1781250.000   1781250.0 223514.5625        2  XX  YY
  33     X0000000000#ALMA_RB_06#BB_3#SW-01#FULL_RES    128   TOPO  236545.813     15625.000   2000000.0 237538.0000        3  XX  YY
  34     X0000000000#ALMA_RB_06#BB_3#SW-01#CH_AVG        1   TOPO  237514.563   1781250.000   1781250.0 237514.5625        3  XX  YY
  35     X0000000000#ALMA_RB_06#BB_4#SW-01#FULL_RES    128   TOPO  238545.813     15625.000   2000000.0 239538.0000        4  XX  YY
  36     X0000000000#ALMA_RB_06#BB_4#SW-01#CH_AVG        1   TOPO  239514.563   1781250.000   1781250.0 239514.5625        4  XX  YY

The fifth sections lists the spectral windows (spws) used for each Source. We exclude this for brevity.

Antenna List

The sixth section shows antenna information. Thirty-one antennas were used for the dataset. Note that numbering in python always begins with "0", so the antennas have IDs 0-30. To see what the antenna configuration looked like at the time of this observation, we use the task plotants (Figure 1):

Antennas: 31:
  ID   Name  Station   Diam.    Long.         Lat.                Offset from array center (m)                ITRF Geocentric coordinates (m)        
                                                                     East         North     Elevation               x               y               z
  0    CM01  N602      7.0  m   -067.45.17.4  -22.53.22.3          8.8042     -527.8587       22.2034  2225080.354846 -5440132.955920 -2481524.789784
  1    CM02  J502      7.0  m   -067.45.17.7  -22.53.23.0          2.1073     -549.4461       22.1460  2225070.957857 -5440127.670516 -2481544.655003
  2    CM03  J503      7.0  m   -067.45.17.4  -22.53.23.2          9.2482     -555.0637       22.1304  2225076.734430 -5440122.931505 -2481549.824201
  3    CM06  N606      7.0  m   -067.45.17.1  -22.53.23.6         19.1995     -566.5684       22.1011  2225084.240492 -5440114.997537 -2481560.411621
  4    CM07  N601      7.0  m   -067.45.17.0  -22.53.22.5         21.0633     -532.5817       22.1865  2225090.999805 -5440126.600430 -2481529.134327
  5    CM08  J505      7.0  m   -067.45.18.0  -22.53.22.8         -7.2123     -541.3466       22.1685  2225063.532652 -5440134.133528 -2481537.202006
  6    CM09  N603      7.0  m   -067.45.17.7  -22.53.22.3         -0.0497     -527.8657       22.1913  2225072.154896 -5440136.294753 -2481524.791487
  7    CM10  J501      7.0  m   -067.45.17.4  -22.53.22.9         10.0863     -545.4959       22.1606  2225078.929507 -5440126.084513 -2481541.021572
  8    CM11  N604      7.0  m   -067.45.17.8  -22.53.23.7         -0.2657     -571.8966       22.0829  2225065.433709 -5440120.432522 -2481565.313207
  9    DA41  A004      12.0 m   -067.45.15.9  -22.53.28.0         52.6609     -704.4171       21.7726  2225094.796703 -5440052.421785 -2481687.277348
  10   DA49  A002      12.0 m   -067.45.16.3  -22.53.27.6         40.6333     -690.2503       21.8023  2225085.761255 -5440062.100754 -2481674.237730
  11   DA50  A038      12.0 m   -067.45.18.5  -22.53.29.4        -22.4285     -745.7518       22.0606  2225019.312629 -5440066.210860 -2481725.468896
  12   DA52  A018      12.0 m   -067.45.17.2  -22.53.28.1         16.8264     -706.6065       21.7531  2225061.301483 -5440065.182346 -2481689.286759
  13   DA54  A005      12.0 m   -067.45.14.8  -22.53.28.7         83.3315     -725.0764       21.7237  2225120.123751 -5440033.331382 -2481706.290680
  14   DA55  A047      12.0 m   -067.45.16.4  -22.53.30.3         38.4542     -775.2187       21.5966  2225071.160250 -5440032.158745 -2481752.434689
  15   DA57  A025      12.0 m   -067.45.18.7  -22.53.27.4        -26.4296     -685.5228       22.2053  2225024.528985 -5440089.533369 -2481670.039352
  16   DA59  A001      12.0 m   -067.45.16.9  -22.53.27.7         24.1880     -693.3966       21.7925  2225070.073933 -5440067.185105 -2481677.132442
  17   DA61  A006      12.0 m   -067.45.15.0  -22.53.28.0         79.0341     -702.0939       21.7778  2225119.549746 -5440043.278703 -2481685.139167
  18   DA62  A035      12.0 m   -067.45.16.6  -22.53.28.1         32.0366     -706.8051       21.7640  2225075.353584 -5440059.362225 -2481689.473934
  19   DA63  A039      12.0 m   -067.45.18.0  -22.53.29.6         -6.1070     -751.7850       22.0696  2225033.533416 -5440057.867935 -2481731.030436
  20   DA65  A030      12.0 m   -067.45.18.1  -22.53.27.2        -10.3876     -679.1070       21.8267  2225040.189208 -5440085.447750 -2481663.981462
  21   DV02  A007      12.0 m   -067.45.15.1  -22.53.27.3         74.0140     -681.2928       21.3255  2225117.808983 -5440052.282474 -2481665.800190
  22   DV07  A027      12.0 m   -067.45.19.0  -22.53.28.7        -35.0460     -726.6032       21.5989  2225010.293430 -5440077.487712 -2481707.648731
  23   DV11  A049      12.0 m   -067.45.14.6  -22.53.29.6         88.4465     -754.5446       20.1464  2225119.968124 -5440019.440510 -2481732.824666
  24   DV12  A017      12.0 m   -067.45.15.9  -22.53.26.8         51.3634     -665.5893       21.3604  2225099.169786 -5440066.540560 -2481651.346952
  25   DV13  A064      12.0 m   -067.45.14.7  -22.53.31.4         85.6567     -808.0278       21.0176  2225109.813656 -5440001.983403 -2481782.434609
  26   DV17  A031      12.0 m   -067.45.19.1  -22.53.27.1        -37.8149     -675.5186       21.7325  2225015.299768 -5440097.041914 -2481660.638998
  27   DV18  A009      12.0 m   -067.45.16.1  -22.53.26.1         48.2542     -644.4621       21.0152  2225099.282791 -5440075.029589 -2481631.749120
  28   DV20  A033      12.0 m   -067.45.19.4  -22.53.29.0        -47.3621     -735.6360       21.8836  2224997.663712 -5440079.140555 -2481716.080921
  29   DV22  A014      12.0 m   -067.45.15.1  -22.53.26.4         74.5090     -654.2102       20.9880  2225122.137554 -5440061.557701 -2481640.719077
  30   DV23  A003      12.0 m   -067.45.16.5  -22.53.27.0         35.5295     -672.6352       21.3390  2225083.469962 -5440069.979758 -2481657.829615
Fig. 1. plotants output for the dataset uid___A002_Xae00c5_X2a8d. The plot has been scaled to approximately equal X and Y spacing. The compact cluster of 7m "CM" antennas is near the top, and the more spread out 12m antennas are labeled "DA" and "DV."

Now let's plot the antennas with plotants:

# In CASA
plotants(vis=mso, figfile='plotants.png')

Flagging Unnecessary Data

Some scans in the data were used by the online system for pointing and sideband ratio calibration. These scans are no longer needed, and we can flag them easily with flagdata by selecting on 'intent'.

#In CASA
flagdata(vis=mso, mode='manual', intent='*POINTING*,*SIDEBAND_RATIO*', flagbackup=False)

The averaged data of each spectrum window is not used, so we flagged the averaged data as a follow.

#In CASA
flagdata(vis=mso, mode='manual', spw='6,8,10,12', flagbackup=False)

We will then store the current flagging state for each dataset using the flagmanager:

#In CASA
flagmanager(vis=mso, mode='save', versionname='priori1')

A priori calibration

Tsys calibration of the calibrators

Fig. 2. One example of the Tsys plots.

The System Temperature (Tsys) calibration gives a first-order correction for the atmospheric opacity as a function of time and frequency, and associates weights with each visibility that persist through imaging. The MS dataset contains Tsys measurements; the task gencal is used to generate a calibration table. Since the dataset obtained with TDM, the data in the channels near the both edges of the spectrum window (~10 channels) are flagged with flagdata. The plots for checking are created by the checkCalTable subroutine of the Analysis Utilities package and saved to a new directory uid___A002_Xae00c5_X2a8d.ms.tsys.plots (Figure 2).

#In CASA
gencal(vis=mso, caltable=mso+'.tsys', caltype='tsys')

flagdata(vis=mso+'.tsys', mode='manual', spw='5:0~9;116~127,7:0~9;116~127,9:0~9;116~127,11:0~9;116~127', flagbackup=False)

es.checkCalTable(mso+'.tsys', msName=mso, interactive=False)

The Tsys of DA54 antennas are significant large, from the plot. In the later part, we will flag the data of the antenna. Please ignore the data of the scan #1, because the scan is not done for ATM cal. The scan is used for the measurement of zero-signal level .

We will apply the Tsys calibration table to the data of the calibrators with the task applycal, which reads the specified gain calibration tables, applies them to the (raw) data column, and writes the calibrated results into the corrected column. For non-solar observations, we also apply the WVR (Water Vapor Radiometer) calibration table to data. However, we must NOT apply the WVR table to the solar data, because the WVR receivers at the Sun are saturated. We apply the Tsys calibration table to the data of the bandpass calibrator:

#In CASA
applycal(vis=mso, field='1', spw='5,7,9,11', gaintable=mso+'.tsys', gainfield='1', interp='linear,linear', calwt=True, flagbackup=False)

For normal ALMA observations, in between the science target scans, we would perform two adjacent scans of the phase calibrator: one with ATM intent, and one with PHASE intent. We would then apply the Tsys calibration table for the phase calibrator to itself.

The Sun is very dynamic, and the time variations of solar structures are one of the main scientific topics of solar physics. To observe such structures, the suspension of science scans are problematic. Therefore, to return to the Sun as quickly as possible, we do NOT perform ATM scans of the phase calibrator in between science scans, instead only observing PHASE intent scans (9,12,15,18).

However, the antenna temperature (Tant) of the Sun is not negligible and must be calibrated. To do this, we perform ATM scans of the blank sky near the Sun, 2 degrees from the Sun's center. Field 0 is used for all calibration scans, and is the center of the mosaic, as we will see later (Fields 3~150 are used for science scans only). The first subscan of science scans (11,14,17) are listed as Field 0, but include the "OBSERVE_TARGET#OFF_SOURCE" intent, which slews to the blank sky.

In this Science Verification observation, one scan of the phase calibrator has ATM intent (scan 7) before any science scans begin, but more recent solar observations do not include this. Scan 7 can be ignored.

For details, see Shimojo+ 2017, in particular Figure 4.

Therefore, we calibrate the Tsys of the phase calibrator (field='2') by applying the calibrations for the Sun (gainfield='0'), as follows:

#In CASA
applycal(vis=mso, field='2', spw='5,7,9,11', gaintable=mso+'.tsys', gainfield='0', interp='linear,linear', calwt=True, flagbackup=False)

You can use plotms to plot channel-averaged amplitudes as a function of time, comparing the DATA and CORRECTED columns after applying the Tsys correction. This way you can check that calibration has done what was expected, which is put the data onto the Kelvin temperature scale.

Tsys+Tant calibration of the Sun

The standard method of Tsys calibration cannot be applied to the data of the Sun, because the antenna temperature (Tant) of the Sun cannot be neglected for estimating the system equivalent flux density (SEFD). To estimate correct SEFD at the Sun, the solar observing sequence includes some special measurements, like the measurement of zero-signal level (scan 1). The subroutines for creating and applying Tsys+Tant calibration tables are prepared by the ALMA solar development team. They are part of the script we executed in the section #Import Tools and Scripts. For the Tsys+Tant calibration of the solar data, you will execute only the following command:

#In CASA
sol_ampcal_2(mso, mso+'.tsys', exisTbl=False, outCSV=True)

This process takes a long time. Depending on your machine, it may be between a couple hours and one day.

If you have already carried out the process before and there are the Tsys+Tant calibration tables (The directory name of the table includes “tsystant”), you can skip the generating of the tables using the following command.

Bug in Sun_reduction_util_6.5.4.py: underneath "if exisTbl == True:", scan = sciScan[i] should be scan = str(sciScan[i]).

#In CASA
sol_ampcal_2(mso, mso+'.tsys', exisTbl=True, outCSV=False)

Flagging after Tsys and Tsys+Tant calibration

Since we completed the Tsys and Tsys+Tant calibration, we can now flag the data that are not used for image synthesis. At first, the data of auto-correlation and atmosphere calibration are flagged with flagdata as follows:

#In CASA
flagdata(vis=mso, mode='manual', autocorr=True, flagbackup=False)
flagdata(vis=mso, mode='manual', intent='*ATMOSPHERE*', flagbackup=False)

Next are the spectral windows that are not used later:

#In CASA
flagdata(vis=mso, mode='manual', spw='0~4', flagbackup=False)
flagdata(vis=mso, mode='manual', spw='13~36', flagbackup=False)

The channels near both edges of the science spectral windows (~10 channels) are flagged:

#In CASA
flagdata(vis=mso, mode='manual', flagbackup=False, spw='5:0~9;116~127,7:0~9;116~127,9:0~9;116~127,11:0~9;116~127')

Some sub-scans at the start and end of the scientific scans are used to measure the sky with the settings optimized for the Sun. Since the data are used only for estimating the antenna temperatures at the Sun, we now flag the data:

#In CASA
mymsmd = createCasaTool(msmdtool)
mymsmd.open(mso)
sciScan = mymsmd.scansforintent('*OBSERVE_TARGET*')
mymsmd.done()

for i in range(0, len(sciScan)):
      subInf=aU.computeDurationOfScan(sciScan[i], vis=mso, returnSubscanTimes=True)
      subNum = subInf[1]
      flagdata(vis=mso, timerange=subInf[3][1], mode='manual', flagbackup=False)       
      flagdata(vis=mso, timerange=subInf[3][2], mode='manual', flagbackup=False)       
      flagdata(vis=mso, timerange=subInf[3][subNum-1], mode='manual', flagbackup=False)       
      flagdata(vis=mso, timerange=subInf[3][subNum], mode='manual', flagbackup=False)

As mentioned in the section #Tsys calibration of the calibrators, the Tsys values of DA54 antennas are significantly large. Therefore, we flag the data of the antenna as follows:

#In CASA
flagdata(vis=mso, antenna='DA54', mode='manual', flagbackup=False)

Now we store the current flagging state using flagmanager:

#In CASA
flagmanager(vis=mso, mode='save', versionname='priori2')

Now split out the CORRECTED data column, retaining spectral windows 5,7,9,11. This will get rid of the extraneous spectral windows.

#In CASA
split(vis=mso, outputvis=mss, datacolumn='corrected', intent ='*BANDPASS*,*FLUX*,*PHASE*,*TARGET*', spw='5,7,9,11')

CAUTION: When we run split, the spws that we choose to keep will be re-indexed to 0,1,2,3.

Additional Data Inspection

Now let's inspect the split MS with plotms.

The solar data frequently include zero values in some channels. This is best seen with all plot averaging turned off, though it may take plotms a moment to load. We can restrict the Y axis to a small number since we are looking for pure zero values.

Fig. 4. Amp vs Time. Pure zero values should be flagged.
#In CASA
plotms(vis=mss,
       xaxis='time',
       yaxis='amp',
       coloraxis='field',
       plotrange=[0,0,-0.1,0.1]
       )

To avoid their influence, we flag the data using the “clip” mode of flagdata with “clipzeros=True” option.

#In CASA
flagdata(vis=mss, mode='clip', clipzeros=True, flagbackup=False)

Now let's plot Amp vs Channel for our two calibrators.

Fig. 5. Amp vs Freq. Atmospheric lines should be flagged.
#In CASA
plotms(vis=mss,
       field='1,2',
       spw='3',
       xaxis='channel',
       yaxis='amp',
       avgantenna=True,
       iteraxis='field'
       )

Valleys appear in spw 3. It shows the effect of the lines of the Earth's atmosphere. We need to flag the data.

#In CASA
flagdata(vis=mss, spw='3:68~85', mode='manual', flagbackup=False)
flagdata(vis=mss, spw='3:20~40', mode='manual', flagbackup=False)

Then, we store the current flagging state for each dataset using the flagmanager:

#In CASA
flagmanager(vis=mss, mode='save', versionname='initial')

Set up the Flux Calibration Model

It is very rare that a useful planet for the flux calibration is located near the Sun. Therefore, we usually use a quasar near the Sun as a flux calibrator. In these observations, there is no quasar with enough flux with the Mixer De-tuning mode except nrao530. Therefore, the flux calibrator is the same as the phase calibrator in these observation. We obtain the flux density and spectral index of nrao530 from the ALMA calibrator database.

#In CASA
intentSources = es.getIntentsAndSourceNames(mss)
ampCalId = intentSources['CALIBRATE_AMPLI']['id'] + intentSources['CALIBRATE_FLUX']['id']
calFieldNames = intentSources['CALIBRATE_AMPLI']['name'] + intentSources['CALIBRATE_FLUX']['name']

amp_cal_name = calFieldNames[1]
spwInfo = es.getSpwInfo(mss)
obs_freq = "%fGHz"%(spwInfo[0]['refFreq']/1e9)

date = aU.getObservationStartDate(mss)
date_obs = date.split()[0]
spw1_flux = aU.getALMAFlux(sourcename=amp_cal_name, date=date_obs, frequency=obs_freq)

You will see the following output in terminal:

Result using spectral index of -0.676873 for 231.000 GHz from 3.265 Jy at 97.500 GHz = 1.821049 +- 0.084762 Jy

Then, fill the model data column for nrao530 by providing this model to setjy.

#In CASA
setjy(vis=mss, field='2', spw='0,1,2,3', standard='manual', fluxdensity=[spw1_flux['fluxDensity'],0,0,0], spix=spw1_flux['spectralIndex'], reffreq=obs_freq, usescratch=True)

You will see the following output in the logger:

Flux density as a function of frequency (channel 0 of each spw):
  Frequency (GHz)    Flux Density (Jy, Stokes I)
     230.992         1.82116
     232.992         1.81056
     245.008         1.74997
     247.008         1.74037

Creating the Bandpass Calibration Table

The bandpass and gain calibrations are essentially the same as for non-solar data. First, we choose a reference antenna, then run gaincal on the bandpass calibrator to determine phase-only gain solutions. We will use solint='int' for the solution interval, which means that one gain solution will be determined for every integration time.

Fig. 7. One example of the pre-bandpass gaincal solution plots. Note the phase wrapping between +/-180 deg.
#In CASA
ref_ant = 'DA41'
gaincal(vis=mss, caltable=mss+'.pre_bandpass', field='1', scan='5', solint='int', refant=ref_ant, calmode='p')

Create plots to check the calibration table. These will be saved to a new directory uid___A002_Xae00c5_X2a8d.ms.split.pre_bandpass.plots.

#In CASA
es.checkCalTable(mss+'.pre_bandpass', msName=mss, interactive=False)

The plots look good. Next, create the bandpass calibration table. Since we apply the per-integration phase solutions on-the-fly, we can use solint='inf' to average the scan in time, resulting in higher signal-to-noise for our bandpass solution. You will see many warning reported about "insufficient unflagged antennas," which are for the edge channels we flagged so they are safe to ignore.

Fig. 7. One example of the bandpass solution plots.
#In CASA
bandpass(vis=mss, caltable=mss+'.bandpass', field='1', scan='5', solint='inf', refant=ref_ant, solnorm=True, bandtype='B', gaintable=mss+'.pre_bandpass')

Check the bandpass calibration table using the plots created from the following command, which are saved into a new directory uid___A002_Xae00c5_X2a8d.ms.split.bandpass.plots.

#In CASA
es.checkCalTable(mss+'.bandpass', msName=mss, interactive=False)


Creating the Gain Calibration Tables

At first, we determine phase-only gain solutions of the calibrators with gaincal, applying the bandpass calibration table on-the-fly, using solint='int' for per-integration solutions. Then create the plots for checking, saved to uid___A002_Xae00c5_X2a8d.ms.split.phase_int.plots.

#In CASA
gaincal(vis=mss, caltable=mss+'.phase_int', field='1~2', solint='int', refant=ref_ant, gaintype='G', calmode='p', minsnr=3.0, gaintable=mss+'.bandpass')

es.checkCalTable(mss+'.phase_int', msName=mss, interactive=False)

Now apply both the bandpass and short phase solution tables on-the-fly to solve for amplitude-only gain solutions over the full scan time with solint='inf'. Plots are saved to uid___A002_Xae00c5_X2a8d.ms.split.ampli_inf.plots.

# In CASA
gaincal(vis=mss, caltable=mss+'.ampli_inf', field='1~2', solint='inf', refant=ref_ant, gaintype='G', calmode='a', minsnr=3.0, gaintable=[mss+'.bandpass', mss+'.phase_int'])

es.checkCalTable(mss+'.ampli_inf', msName=mss, interactive=False)

Usually, the gaintype is set to 'T' for nominal ALMA observations. For solar data, the gaintype is set to 'G', because the difference between the XX and YY images is used for estimating the noise level (see Section 4.2, Shimojo et al. 2017).

Next we use the flux calibrator (whose flux density was set in setjy above) to derive the flux density of the other calibrators with fluxscale. Note that the flux table REPLACES '.ampli_inf' in terms of future application of the calibration to the data, i.e., the flux table contains both the amplitude calibration and flux scaling. Unlike the two gaincal steps above, this is not an incremental table. Plots are saved to uid___A002_Xae00c5_X2a8d.ms.split.flux_inf.plots.

Fig. 7. Example of the flux-scaled amplitude calibration table.
#In CASA
fluxscaleDict = fluxscale(vis=mss, caltable=mss+'.ampli_inf', fluxtable=mss+'.flux_inf', reference='2')

es.checkCalTable(mss+'.flux_inf', msName=mss, interactive=False)

You should see the following output:

Found reference field(s): nrao530
Found transfer field(s):  J1924-2914
Flux density for J1924-2914 in SpW=0 (freq=2.3e+11 Hz) is: 3.31068 +/- 0.0174892 (SNR = 189.299, N = 60)
Flux density for J1924-2914 in SpW=1 (freq=2.32e+11 Hz) is: 3.30137 +/- 0.016152 (SNR = 204.394, N = 60)
Flux density for J1924-2914 in SpW=2 (freq=2.46e+11 Hz) is: 3.16067 +/- 0.0180585 (SNR = 175.024, N = 60)
Flux density for J1924-2914 in SpW=3 (freq=2.48e+11 Hz) is: 3.14243 +/- 0.0205412 (SNR = 152.981, N = 60)
Fitted spectrum for J1924-2914 with fitorder=1: Flux density = 3.22814 +/- 0.002994 (freq=238.864 GHz) spidx: a_1 (spectral index) =-0.712504 +/- 0.0278242 covariance matrix for the fit:  covar(0,0)=7.85242e-06 covar(0,1)=9.80161e-05 covar(1,0)=9.80161e-05 covar(1,1)=0.0374698

Finally, we create the gain calibration table of phase-only gain solutions on the scan time, solint='inf'. Plots are saved to uid___A002_Xae00c5_X2a8d.ms.split.phase_inf.plots.

Fig. 7. Example of the long solution interval phase calibration table.
#In CASA
gaincal(vis=mss, caltable=mss+'.phase_inf', field='1~2', solint='inf', refant=ref_ant, gaintype='G', calmode='p', minsnr=3.0, gaintable=mss+'.bandpass')

es.checkCalTable(mss+'.phase_inf', msName=mss, interactive=False)


Applying the Calibration Tables

We apply the calibration solutions to each source individually with applycal, using the gainfield parameter to specify which calibrator's solutions should be applied from each of the calibration tables. We can leave gainfield blank when referencing the bandpass table, because this table has information derived only from the bandpass calibrator which will be applied to all sources.

Apply to the bandpass calibrator: the bandpass table; the short solint phase table, referencing the BP cal; and the flux table, referencing the BP cal.

#In CASA
applycal(vis=mss, field='1', gaintable=[mss+'.bandpass', mss+'.phase_int', mss+'.flux_inf'], gainfield=['', '1', '1'],
interp='linear,linear', calwt=True, flagbackup=False)

Apply to the phase calibrator: the bandpass table; the short solint phase table, referencing the phase cal; and the flux table, referencing the phase cal.

#In CASA
applycal(vis=mss, field='2', gaintable=[mss+'.bandpass', mss+'.phase_int', mss+'.flux_inf'], gainfield=['', '2', '2'],
interp='linear,linear', calwt=True, flagbackup=False)

Apply to the Sun: the bandpass table; the long solint phase table, referencing the phase cal; and the flux table, referencing the phase cal.

#In CASA
applycal(vis=mss, field='0, 3~150', gaintable=[mss+'.bandpass', mss+'.phase_inf', mss+'.flux_inf'], gainfield=['', '2', '2'],
interp='linear,linear', calwt=True, flagbackup=False)

Finally, split out the CORRECTED data column.

#In CASA
split(vis=mss, outputvis=msc, datacolumn='corrected')

Re-calculation of the direction

Fig. 3. The pattern of mosaic BEFORE the re-calculation of the direction.
Fig. 4. The pattern of mosaic AFTER the re-calculation of the direction.

During most solar observations, the antennas are tracking a structure on the Sun according to the solar differential rotation. The image frame is fixed on the solar frame, but the frame is moving on the RA/Dec coordinate frame. If we do not correct for this, the pattern of pointing in the mosaic is a rhombus as shown in Figure 3, while the correct shape should be a square.

Plot the mosaic to see this:
The following command does not correctly create the plot as shown in Figure 3 at this time. Figure 3 is taken from a previous version of this guide.

#In CASA
au.plotmosaic(vis=msc, sourceid='Sun', figfile=msc+'.pointings.sun.before.png')

To correct the MOSAIC pattern, we re-calculate the coordinate of each field. First, we modify the coordinate of field='0' from the reference time using fixplanets task. The reference time has to be the time when the antennas are directed to field='0'.

#In CASA
reftime = '2015/12/18/19:49:00'
fixplanets(vis=msc, field='0', fixuvw=False, refant=ref_ant, reftime=reftime)

We define that the modified coordinate of field='0' is the reference coordinate, and re-calculate the coordinate of each field, as follows:

#In CASA
import math
pi = math.pi

tb.open(msc+'/FIELD', nomodify=True)
phsCenOff = tb.getcol("PHASE_DIR")
tb.close()

refRaDec = aU.rad2radec(phsCenOff[0][0][0], phsCenOff[1][0][0], prec=1, hmsdms=True, delimiter=' ')   
for i in range(3, 151):
    raOff = phsCenOff[0][0][i] * 180. / pi * 60. *60.
    deOff = phsCenOff[1][0][i] * 180. / pi * 60. *60.
    offRaDec = aU.radec2deg(aU.radecOffsetToRadec(refRaDec, raOff, deOff, prec=1))
    offRaDecF = 'J2000 ' + aU.deg2radec(offRaDec[0], offRaDec[1], prec=1, hmsdms=True, delimiter=' ')         
    fixplanets(vis=msc, field=str(i), fixuvw=False, direction=offRaDecF)

tb.open(msc+'/FIELD', nomodify=False)
tgt_refdir = tb.getcol("RefDir_Ref")
for id in range(3, len(tgt_refdir)):
    tb.putcell("RefDir_Ref", id, 21)
    tb.putcell("DelayDir_Ref", id, 21)
    tb.putcell("PhaseDir_Ref", id, 21)
tb.close()

Moreover, the direction in the pointing table has a bad influence on the coordinate system for image synthesis. We erase the pointing table as follows:

#In CASA
tb.open(msc+'/POINTING', nomodify=False)
a = tb.rownumbers()
tb.removerows(a)
tb.close()

Now plot the corrected mosaic:

#In CASA
au.plotmosaic(vis=msc, sourceid='Sun', figfile=msc+'.pointings.sun.after.png')