casa6

//# Importmiriad: miriad dataset to MeasurementSet conversion

//# Copyright (C) 1997,2000,2001,2002,2013,2015

//# Associated Universities, Inc. Washington DC, USA.

//#

//# This program is free software; you can redistribute it and/or modify

//# it under the terms of the GNU General Public License as published by

//# the Free Software Foundation; either version 2 of the License, or

//# (at your option) any later version.

//#

//# This program is distributed in the hope that it will be useful,

//# but WITHOUT ANY WARRANTY; without even the implied warranty of

//# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

//# GNU General Public License for more details.

//# 

//# You should have received a copy of the GNU General Public License

//# along with this program; if not, write to the Free Software

//# Foundation, Inc., 675 Masve, Cambridge, MA 02139, USA.

//#

​

//#Includes

#include <miriad/Filling/Importmiriad.h>

​

#include <casacore/casa/Inputs/Input.h>

#include <casacore/casa/OS/File.h>

#include <casacore/casa/Utilities/GenSort.h>

​

#include <casacore/casa/Arrays/Cube.h>

#include <casacore/casa/Arrays/Matrix.h>

#include <casacore/casa/Arrays/Vector.h>

#include <casacore/casa/Arrays/ArrayMath.h>

#include <casacore/casa/Arrays/ArrayUtil.h>

#include <casacore/casa/Arrays/ArrayLogical.h>

#include <casacore/casa/Arrays/MatrixMath.h>

​

​

#include <casacore/measures/Measures.h>

#include <casacore/measures/Measures/MPosition.h>

#include <casacore/measures/Measures/MeasData.h>

#include <casacore/measures/Measures/Stokes.h>

​

#include <casacore/tables/Tables.h>

#include <casacore/tables/Tables/TableInfo.h>

​

#include <casacore/ms/MeasurementSets.h> 

​

#include <casacore/mirlib/maxdimc.h>

#include <casacore/mirlib/miriad.h>

#include <sstream>

​

using namespace casa;

​

// helper functions

​

String show_version_info()

{

  return "============================================================\n"

   "Importmiriad - last few updates:\n"

   " Mar 2013 - make it process ATCA/CABB data \n"

   " Jul 2015 - deal with multiple zoom setups \n"

   "           ...\n"

   "============================================================\n";

}

​

​

// Convert fits date string of form dd/mm/yy to mjd seconds

​

Double date2mjd(const String& date)

{

  Int day,month,year;

  if (date[2] == '/') {     // old FITS style (dd/mm/yy)

    sscanf(date.chars(),"%2d/%2d/%2d",&day,&month,&year);

    year+=1900;

    if (year<1950) year+=100;   // yuck !

  } else {      // YEAR-2000 (ISO) convention (ccyy-mm-ddThh:mm:ss.sss)

    // cerr<< "Parsing new Year-2000 notation" << endl;

    sscanf(date.chars(),"%4d-%2d-%2d",&year,&month,&day);

    //sscanf(date,"%4d-%2d-%2dT%2d:%2d:%f",&year,&month,&day,&hour,&min,&sec);

  }

  MVTime mjd_date(year,month,(Double)day);

  return mjd_date.second();

}

​

// apply CARMA line calibration, the 'linecal' method

// check MIRIADs 'uvcal options=linecal' for the other approach

​

void linecal(int ndata, float *data, float phi1, float phi2)

{

  float x,y,c,s;

  if (ndata <= 0) return;

​