CoordTransform.hxx

// Copyright (C) 2003-2007, ENPC - INRIA - EDF R&D
// Author(s): Vivien Mallet
//
// This file is part of AtmoData library, a tool for data processing in
// atmospheric sciences.
//
// AtmoData is developed in the INRIA - ENPC joint project-team CLIME and in
// the ENPC - EDF R&D joint laboratory CEREA.
//
// AtmoData 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.
//
// AtmoData 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.
//
// For more information, visit the AtmoData home page:
//      http://cerea.enpc.fr/polyphemus/atmodata.html


#ifndef ATMODATA_FILE_COORDTRANSFORM_HXX

#include <cmath>

namespace AtmoData
{
  
  template<class T>
  class LaeaToLonlat
  {
  protected:
    const T Earth_radius_;
    const T pi_;
    const T limit_;

    T lat_origin_;
    T lon_origin_;

  public:
    LaeaToLonlat(T lon_origin, T lat_origin)  throw();
    void operator() (const T x, const T y,
                     T& lon, T& lat);
  };


  template<class T>
  class MM5LccIndToLonlat
  {
  protected:
    const int jmx_;
    const int imx_;
    const double jx_;
    const double ix_;
    const double phic_;
    const double lambdac_;
    const double phi1_;
    const double phi2_;
    const double ds_;
    const int ratio_;
    const double Earth_radius_;
    const double pi_;

  public:
    MM5LccIndToLonlat(int jmx, int imx, double jx, double ix,
                      double phic, double lambdac,
                      double phi1, double phi2, double ds,
                      int ratio)  throw();
    void operator() (const T j, const T i,
                     T& lon, T& lat);
  };
  

  template<class T>
  class LonlatToMM5LccInd
  {
  protected:
    const int jmx_;
    const int imx_;
    const double jx_;
    const double ix_;
    const double phic_;
    const double lambdac_;
    const double phi1_;
    const double phi2_;
    const double ds0_;
    const int ratio_;
    const double Earth_radius_;
    const double pi_;

  public:
    LonlatToMM5LccInd(int jmx, int imx, double jx, double ix,
                      double phic, double lambdac,
                      double phi1, double phi2, double ds0,
                      int ratio)  throw();
    void operator() (const T lon, const T lat,
                     T& j, T& i);
  };


  template<class T>
  class MM5MercIndToLonlat
  {
  protected:
    const int jmx_;
    const int imx_;
    const double jx_;
    const double ix_;
    const double phic_;
    const double lambdac_;
    const double phi1_;
    const double ds_;
    const int ratio_;
    const double Earth_radius_;
    const double pi_;

  public:
    MM5MercIndToLonlat(int jmx, int imx, double jx, double ix,
                       double phic, double lambdac,
                       double phi1, double ds, int ratio)  throw();
    void operator() (const T j, const T i,
                     T& lon, T& lat);
  };


  template<class T>
  class LonlatToMM5MercInd
  {
  protected:
    const int jmx_;
    const int imx_;
    const double jx_;
    const double ix_;
    const double phic_;
    const double lambdac_;
    const double phi1_;
    const double ds0_;
    const int ratio_;
    const double Earth_radius_;
    const double pi_;

  public:
    LonlatToMM5MercInd(int jmx, int imx, double jx, double ix,
                       double phic, double lambdac,
                       double phi1, double ds0, int ratio)  throw();
    void operator() (const T lon, const T lat,
                     T& j, T& i);
  };


  template<class T>
  class MM5StereIndToLonlat
  {
  protected:
    const int jmx_;
    const int imx_;
    const double jx_;
    const double ix_;
    const double phic_;
    const double lambdac_;
    const double phi1_;
    const double ds_;
    const int ratio_;
    const double Earth_radius_;
    const double pi_;

  public:
    MM5StereIndToLonlat(int jmx, int imx, double jx, double ix,
                        double phic, double lambdac,
                        double phi1, double ds, int ratio)  throw();
    void operator() (const T j, const T i,
                     T& lon, T& lat);
  };


  template<class T>
  class LonlatToMM5StereInd
  {
  protected:
    const int jmx_;
    const int imx_;
    const double jx_;
    const double ix_;
    const double phic_;
    const double lambdac_;
    const double phi1_;
    const double ds0_;
    const int ratio_;
    const double Earth_radius_;
    const double pi_;

  public:
    LonlatToMM5StereInd(int jmx, int imx, double jx, double ix,
                        double phic, double lambdac,
                        double phi1, double ds0, int ratio)  throw();
    void operator() (const T lon, const T lat,
                     T& j, T& i);
  };


  template<class T>
  class LonlatToWRFLccInd
  {
  protected:
    const int imx_;
    const int jmx_;
    const double lambdar_;
    const double lambda0_;
    const double lambda1_;
    const double phir_;
    const double max_phi0_;
    const double max_phi1_;
    const double phi1_;
    const double phi2_;
    const double dsi0_;
    const double dsj0_;
    const double Earth_radius_;
    const double pi_;

  public:
    LonlatToWRFLccInd(int imx, int jmx,
                      double lambdar, double phir,
                      double lambda0, double lambda1,
                      double max_phi0, double max_phi1,
                      double phi1, double phi2,
                      double dsi0, double dsj0)  throw();
    void operator() (const T lon, const T lat,
                     T& i, T& j);
  };

  
  template<class T>
  class LonlatToWRFMercInd
  {
  protected:
    const int imx_;
    const int jmx_;
    const double lambdac_;
    const double phic_;
    const double phi1_;
    const double dsi0_;
    const double dsj0_;
    const double Earth_radius_;
    const double pi_;

  public:
    LonlatToWRFMercInd(int imx, int jmx,
                       double lambdac, double phic,
                       double phi1,
                       double dsi0, double dsj0)  throw();
    void operator() (const T lon, const T lat,
                     T& i, T& j);
  };


  template<class T>
  class LonlatToWRFStereInd
  {
  protected:
    const int imx_;
    const int jmx_;
    const double lambdac_;
    const double phic_;
    const double phi1_;
    const double dsi0_;
    const double dsj0_;
    const double Earth_radius_;
    const double pi_;

  public:
    LonlatToWRFStereInd(int imx, int jmx,
                        double lambdac, double phic,
                        double phi1,
                        double dsi0, double dsj0)  throw();
    void operator() (const T lon, const T lat,
                     T& j, T& i);
  };


}  // namespace AtmoData.

#define ATMODATA_FILE_COORDTRANSFORM_HXX
#endif

---