Changeset 1940 for branches/asap4casa3.1.0/external-alma/atnf/PKSIO

Timestamp:

09/30/10 14:28:10 (14 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: s=sd.scantable('sdfile')

s.save('output.ms',format='MS2')

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

Filling empty arrays by dummy values for POLARIZATION and SPECTRAL_WINDOW
subtables.

File:

• branches/asap4casa3.1.0/external-alma/atnf/PKSIO/PKSMS2writer.cc (modified) (6 diffs)

branches/asap4casa3.1.0/external-alma/atnf/PKSIO/PKSMS2writer.cc

@@ -732 +732 @@
   // Extend the DATA_DESCRIPTION subtable.
   while (Int(cDataDescription.nrow()) < IFno) {
-    cDataDescription.addRow();
+    cDataDescription.addRow(1,True);
   }
   Int n = IFno - 1;
@@ -966 +966 @@
 {
   // Extend the POLARIZATION subtable.
+  Int nr = cPolarization.nrow() ;
   while (Int(cPolarization.nrow()) < IFno) {
-    cPolarization.addRow();
+    cPolarization.addRow(1,True);
   }
   Int n = IFno - 1;
@@ -985 +986 @@
   corrType(1) = Stokes::YY;
   }
+  // initialization
+  for ( Int i = nr ; i < n ; i++ )
+    cPolarizationCols->corrType().put( i, corrType ) ;
   cPolarizationCols->corrType().put(n, corrType);
 
@@ -1092 +1096 @@
 {
   // Extend the SPECTRAL_WINDOW subtable.
+  Int nr = cSpectralWindow.nrow() ;
   while (Int(cSpectralWindow.nrow()) < IFno) {
-    cSpectralWindow.addRow();
+    cSpectralWindow.addRow(1,True);
   }
   Int n = IFno - 1;
@@ -1111 +1116 @@
     freqs(i) = refFreq + (i - refChan)*freqInc;
   }
+  // initialization
+  for ( Int i = nr ; i < n ; i++ ) 
+    cSpWindowCols->chanFreq().put( i, freqs ) ;
   cSpWindowCols->chanFreq().put(n, freqs);
 
   Vector<Double> chanWidths(nChan, freqInc);
+  // initialization
+  for ( Int i = nr ; i < n ; i++ ) {
+    cSpWindowCols->chanWidth().put( i, chanWidths ) ;
+    cSpWindowCols->effectiveBW().put( i, chanWidths ) ;
+  }
   cSpWindowCols->chanWidth().put(n, chanWidths);
 
@@ -1120 +1133 @@
 
   Vector<Double> resolution(nChan, fabs(freqInc));
+  // initialization
+  for ( Int i = nr ; i < n ; i++ ) 
+    cSpWindowCols->resolution().put( i, resolution ) ;
   cSpWindowCols->resolution().put(n, resolution);