Changeset 363 for trunk/src

Timestamp:

02/03/05 16:45:28 (20 years ago)

Author:

kil064

Message:

gain elevation poly moved to SDAttr

File:

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

trunk/src/SDMath.cc

@@ -1160 +1160 @@
   ROScalarColumn<Float> elev(tab, "ELEVATION");
   Vector<Float> x = elev.getColumn();
-  x *= Float(180 / C::pi);
+  x *= Float(180 / C::pi);                        // Degrees
 //
   const uInt nC = coeffs.nelements();
@@ -1186 +1186 @@
         msg = String("user");
      } else {
-        if (inst==ATPKSMB) {
-        } else if (inst==ATPKSHOH) {
-        } else if (inst==TIDBINBILLA) {
-           pPoly = new Polynomial<Float>(3);
-           coeff.resize(3);
-           coeff(0) = 3.58788e-1;
-           coeff(1) = 2.87243e-2;
-           coeff(2) = -3.219093e-4;
-        } else if (inst==ATMOPRA) {
-        } else {
-        }
+        SDAttr sdAttr;
+        coeff = sdAttr.gainElevationPoly(inst);
+        pPoly = new Polynomial<Float>(3);
         msg = String("built in");
      }
@@ -1203 +1195 @@
         pPoly->setCoefficients(coeff);
      } else {
-        throw(AipsError("There is no known gain-el polynomial known for this instrument"));
+        throw(AipsError("There is no known gain-elevation polynomial known for this instrument"));
      }
 //
@@ -1303 +1295 @@
 //
    Vector<Float> JyPerK = sdAtt.JyPerK(inst, dateObs, freqs);
+   cerr << "Jy/K = " << JyPerK << endl;
    Vector<Float> factors = cFac * JyPerK;
    if (toKelvin) factors = Float(1.0) / factors;