Changeset 496

Timestamp:

02/22/05 11:21:13 (20 years ago)

Author:

kil064

Message:

write header and data into seaprate files
write abcissa to data file

File:

• trunk/src/SDAsciiWriter.cc (modified) (5 diffs)

trunk/src/SDAsciiWriter.cc

@@ -37 +37 @@
 #include <casa/Quanta/Quantum.h>
 #include <casa/Quanta/MVAngle.h>
-
+#include <casa/Utilities/Assert.h>
 
 #include <casa/iostream.h>
@@ -94 +94 @@
    const Unit RAD(String("rad"));
 
-// Open file
-
-   String fName(fileName);
-   if (fileName.length()==0) fName = String("ascii.txt");
+// Open and write header file
+
+
+   String rootName(fileName);
+   if (rootName.length()==0) rootName = String("ascii");
+   {
+      String fName = String(rootName) + String("_header.txt");
+      cout << "Writing header to " << fName << endl;
+      ofstream of(fName.chars(), ios::trunc);
+      std::string summary = sdTable.summary(true);
+      of << summary;
+      of.close();
+   }
+
+// Open data file
+
+   String fName = rootName + String(".txt");
    ofstream of(fName.chars(), ios::trunc);
 
@@ -116 +129 @@
       const Array<Bool>& mask = dataIn.getMask();
 
+// Get abcissa
+
+      std::vector<double> abcissa = sdTable.getAbcissa(Int(iRow));
+      const uInt n = abcissa.size();
+
 // Iterate through data in this row by spectra
 
@@ -122 +140 @@
       while (!itData.pastEnd()) {
          const IPosition& pos = itData.pos();
+         AlwaysAssert(itData.vector().nelements()==n,AipsError);
 
 // FreqID
@@ -144 +163 @@
                pos(asap::PolAxis) << " " <<
                src(iRow) <<  " " << formatDirection(lonLat) << " " << 
-               sdTable.getTime(iRow,True) << " " << itData.vector().nelements() << " ";
+               sdTable.getTime(iRow,True) << " " << n << " ";
+
+// Write abcissa
+
+         of.setf(ios::fixed, ios::floatfield);
+         of.precision(4);
+         for (uInt i=0; i<n; i++) {
+            of << abcissa[i] << " ";
+         }
+
 // Write data
 
          const Vector<Float>& data = itData.vector();
          const Vector<Bool>& mask = itMask.vector();
-         const uInt n = data.nelements();
          for (uInt i=0; i<n; i++) {
             of << data[i] << " ";