Changeset 338 for trunk/src/SDReader.cc

Timestamp:

01/31/05 20:54:05 (19 years ago)

Author:

kil064

Message:

add brightness unit to constructor

File:

• trunk/src/SDReader.cc (modified) (4 diffs)

trunk/src/SDReader.cc

@@ -32 +32 @@
 #include <casa/OS/Path.h>
 #include <casa/OS/File.h>
+#include <casa/Quanta/Unit.h>
 #include <atnf/PKSIO/PKSreader.h>
 
@@ -45 +46 @@
   cursor_ = 0;
 }
-SDReader::SDReader(const std::string& filename, int whichIF, int whichBeam) :
+SDReader::SDReader(const std::string& filename, const std::string& unit,
+                   int whichIF, int whichBeam) :
   reader_(0),
   table_(new SDMemTable()) {
   cursor_ = 0;
-  open(filename, whichIF, whichBeam);
+  open(filename, unit, whichIF, whichBeam);
 }
 
@@ -86 +88 @@
 
 
-void SDReader::open(const std::string& filename, int whichIF, int whichBeam) {
+void SDReader::open(const std::string& filename, const std::string& unit, 
+                    int whichIF, int whichBeam) {
   if (reader_) delete reader_; reader_ = 0;
   Bool   haveBase, haveSpectra, haveXPol;
@@ -146 +149 @@
 // Determine Telescope and set units...
 
-  Bool throwIt = False;
-  Instrument inst = SDMemTable::convertInstrument (header_.antennaname, throwIt);
-  header_.fluxunit = "Jy";
-  if (inst==ATMOPRA || inst==TIDBINBILLA) {
-     header_.fluxunit = "K";
-  }
-  cerr << "Assuming brightness units are " << header_.fluxunit << endl;
+  String sUnit(unit);
+  if (sUnit.empty()) {
+     Bool throwIt = False;
+     Instrument inst = SDMemTable::convertInstrument (header_.antennaname, throwIt);
+     header_.fluxunit = "Jy";
+     if (inst==ATMOPRA || inst==TIDBINBILLA) {
+        header_.fluxunit = "K";
+     }
+     cerr << "Assuming brightness unit is " << header_.fluxunit << endl;
+   } else {
+      Unit u(sUnit);
+      if (u==Unit("Jy") || u==Unit("K")) {
+         header_.fluxunit = unit;
+         cerr << "Setting brightness unit to " << header_.fluxunit << endl;
+      } else {
+         throw(AipsError("Specified brightness unit is illegal - must be consistent with Jy or K"));
+      }
+   }
 //
   header_.nif = nIF_;