Changeset 1493

Timestamp:

01/27/09 13:57:36 (16 years ago)

Author:

Takeshi Nakazato

Message:

New Development: No

JIRA Issue: Yes CAS-1043

Ready to Release: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: No

Module(s): Module Names change impacts.

Description: Fixed frequency and velocity conversion from/to

channel:

1. frequency conversion is done in NROReader, doesn't require conversion using casacore.
2. velocity obtained by

scantable.set_unit('km/s')
scantable._getabcissa()

doesn't contain source velocity (SRCVELOCITY)
contribution, so that the velocity information
should be interpreted as an offset from the
velocity center.

File:

• branches/alma/src/STFiller.cpp (modified) (3 diffs)

branches/alma/src/STFiller.cpp

@@ -492 +492 @@
   else if ( strncmp( eq, "J2000", 5 ) == 0 ) 
     header_->equinox = 2000.0 ;
-  char *vref = nheader->getVREF() ;
+//   char *vref = nheader->getVREF() ;
 //   if ( strncmp( vref, "LSR", 3 ) == 0 ) {
 //     strcat( vref, "K" ) ;
 //   }
-  header_->freqref = vref ;
+//   header_->freqref = vref ;
   nreader_->getData( 0 ) ;
   header_->reffreq = nreader_->getData()->FREQ0 ;
@@ -556 +556 @@
   table_->setHeader( *header_ ) ;
 
-  // set velocity reference 
-//   table_->frequencies().setFrame( vref ) ;
-
   // DEBUG
   //cout << "STFiller::openNRO() Velocity Definition = " << nheader->getVDEF() << endl ;
@@ -643 +640 @@
       RecordFieldPtr<Double> svelCol(rec, "SRCVELOCITY") ;
       *svelCol = h->getURVEL() ;   // [m/s]
-      //*svelCol = h->getURVEL() + d->VRAD ;  // [m/s]
-      //*svelCol = d->VRAD ;  // [m/s]
-      //*svelCol = 0.0 ; 
       RecordFieldPtr<Int> fitCol(rec, "FIT_ID") ;
       *fitCol = -1 ;