Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/external-alma/atnf/PKSIO/NROFITSDataset.cc

    r1757 r1869  
    179179    + sizeof( double ) * ARYNM * 10 * 3    // FQCAL, CHCAL, CWCAL
    180180    + sizeof( char ) * 180 ;               // CDMY1
     181
     182  refFreq_.resize( ARYNM, 0.0 ) ;
    181183}
    182184
     
    11651167  }
    11661168  else {
    1167     sprintf( record_->LAVST, "%d%d%d%d%d%d.000", itmp[0], itmp[1], itmp[2], itmp[3], itmp[4], itmp[5] ) ;
     1169    sprintf( record_->LAVST, "%4d%02d%02d%02d%02d%02d.000", itmp[0], itmp[1], itmp[2], itmp[3], itmp[4], itmp[5] ) ;
    11681170  }
    11691171  // DEBUG
     
    14601462//   //cout << endl ;
    14611463  //
     1464
     1465
     1466  // Update IPTIM since it depends on the row for NROFITS
     1467  int integ ;
     1468  status = readTable( integ, "INTEG", same_, i ) ;
     1469  if ( !status ) {
     1470    IPTIM = (double)integ ;
     1471  }
     1472
    14621473  return status ;
    14631474}
     
    15131524    // debug
    15141525    //cout << "i=" << i << endl ;
    1515     if ( i == 16 ) {
    1516       ofstream ofs0( "spgrid0.dat" ) ;
    1517       for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
    1518         ofs0 << xout[ii] << "," ;
    1519       ofs0 << endl ;
    1520       for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
    1521         ofs0 << setprecision(16) << record->FREQ0+yout[ii] << "," ;
    1522       ofs0 << endl ;
    1523       ofs0.close() ;
    1524     }
     1526//     if ( i == 16 ) {
     1527//       ofstream ofs0( "spgrid0.dat" ) ;
     1528//       for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
     1529//         ofs0 << xout[ii] << "," ;
     1530//       ofs0 << endl ;
     1531//       for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
     1532//         ofs0 << setprecision(16) << record->FREQ0+yout[ii] << "," ;
     1533//       ofs0 << endl ;
     1534//       ofs0.close() ;
     1535//     }
    15251536    //
    15261537    Vector<Double> z( nchan ) ;
     
    27152726  return ib ;
    27162727}
     2728
     2729double NROFITSDataset::getStartIntTime( int i )
     2730{
     2731  double v ;
     2732  readTable( v, "MJDST", same_, i ) ;
     2733  return v/86400.0 ;
     2734}
     2735
     2736// double NROFITSDataset::toLSR( double v, double t, double x, double y )
     2737// {
     2738//   return v ;
     2739// }
Note: See TracChangeset for help on using the changeset viewer.