#pragma ident "$Id: DCBDataReader.hpp 3143 2012-06-19 16:19:50Z snelsen $"

/**
 * @file DCBReader.hpp
 * Class to read DCB data from CODE.
 */

#ifndef GPSTK_DCBDATAREADER_HPP
#define GPSTK_DCBDATAREADER_HPP

//============================================================================
//
//  This file is part of GPSTk, the GPS Toolkit.
//
//  The GPSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 2.1 of the License, or
//  any later version.
//
//  The GPSTk 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 Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GPSTk; if not, write to the Free Software Foundation,
//  Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110, USA
//
//  Wei Yan - Chinese Academy of Sciences  2009, 2010
//
//============================================================================


#include <string>
#include <map>

#include "Exception.hpp"
#include "FFTextStream.hpp"
#include "StringUtils.hpp"
#include "SatID.hpp"


namespace gpstk
{

      /** @addtogroup formattedfile */
      //@{

      /** This is a class to read and DCB(Differences of Code Biases) data file
       *  from CODE.
       *
       * You can find DCB data at:
       *
       *    ftp.unibe.ch/aiub/BSWUSER50/ORB   - daily P1-P2
       *    ftp.unibe.ch/aiub/CODE            - monthly P1-P2 and P1-C1
       *
       *
       *  You should use different objects to load different DCB files. A typical 
       *  way to use these classes follows:
       *
       * @code
       *      // Declare some Antenna objects
       *   DCBDataReader dcbP1P2("P1P21002_ALL.DCB");
       *   DCBDataReader dcbP1C1("P1C11002.DCB");
       *   
       *   double p1p2Sat1 = dcbP1P2.getDCB(1, SatID::systemGPS);
       *   double p1c1Sat1 = dcbP1C1.getDCB(1, SatID::systemGPS);
       *   
       *   double p1p2ALGO = dcbP1P2.getDCB("ALGO");
       *
       * @endcode
       *
       * @sa DCBDataReader.hpp
       */
   class DCBDataReader : public FFTextStream
   {
   public:
         /// Default constructor
      DCBDataReader()
      {};
         
         /** Common constructor. It will always open file for read and will
          *  load DCB data in one pass.
          *
          * @param fn   DCB data file to read
          *
          */
      DCBDataReader(const char* fn)
         : FFTextStream(fn, std::ios::in)
      { loadData(); };


         /** Common constructor. It will always open file for read and will
          *  load DCB data in one pass.
          *
          * @param fn   DCB data file to read
          *
          */
      DCBDataReader(const std::string& fn)
         : FFTextStream(fn.c_str(), std::ios::in)
      { loadData(); };


         /// Method to open AND load DCB data file.
      virtual void open(const char* fn);


         /// Method to open AND load DCB data file. It doesn't
         /// clear data previously loaded.
      virtual void open(const std::string& fn);


         /// Get DCB data of a satellite
         /// @param    sat   the satellite you desired
         /// @return         P1-P2 or P1-C1 depend what you have loaded
      double getDCB( const SatID& sat);

         /// Get DCB data of a satellite
         /// @param    prn    the satellite id you desired
         /// @param    system the satellite system you desired
         /// @return          P1-P2 or P1-C1 depend what you have loaded
      double getDCB(const int& prn,
         const SatID::SatelliteSystem& system = SatID::systemGPS);


         /// Get DCB data of a receiver
         /// @param    station    the receiver name you desired
         /// @param    system     the satellite system you desired
         /// @return              P1-P2
      double getDCB(const std::string& station,
         const SatID::SatelliteSystem& system = SatID::systemGPS);


         /// Destructor
      virtual ~DCBDataReader() {};


   private:

         // Map holding satellite DCB data
      typedef std::map< SatID, double > SatDCBData;
      
         // Map holding receiver DCB data
      typedef std::map< std::string, double > ReceiverDCBData;
         
         /// A structure used to store daily DCB data
      struct DailyDCBData
      {
         SatDCBData        satDCB;

         ReceiverDCBData   gpsDCB;

         ReceiverDCBData   glonassDCB; 
      };

         /// Object holding all of the DCB data
      DailyDCBData allDCB;


         /// Method to store ocean tide harmonics data in this class' data map
      virtual void loadData()
         throw( FFStreamError, gpstk::StringUtils::StringException );


   };  // End of class 'DCBDataReader'

       //@}

}  // End of namespace gpstk

#endif  // GPSTK_DCBDATAREADER_HPP




