casa6

/*******************************************************************************

 * ALMA - Atacama Large Millimiter Array

 * (c) Instituto de Estructura de la Materia, 2009

 *

 * This library 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 (at your option) any later version.

 *

 * This library 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 this library; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307  USA

 *

 * "@(#) $Id: ATMRefractiveIndexProfile.cpp Exp $"

 *

 * who       when      what

 * --------  --------  ----------------------------------------------

 * pardo     24/03/09  created

 */

​

#include "ATMRefractiveIndexProfile.h"

​

#include <iostream>

#include <math.h>

#include <string>

#include <vector>

​

​

​

ATM_NAMESPACE_BEGIN

​

// Constructors

​

RefractiveIndexProfile::RefractiveIndexProfile(const Frequency &freq,

                                               const AtmProfile &atmProfile) :

  AtmProfile(atmProfile), SpectralGrid(freq)

{

  mkRefractiveIndexProfile();

}

​

RefractiveIndexProfile::RefractiveIndexProfile(const SpectralGrid &spectralGrid,

                                               const AtmProfile &atmProfile) :

  AtmProfile(atmProfile), SpectralGrid(spectralGrid)

{

  mkRefractiveIndexProfile();

}

​

RefractiveIndexProfile::RefractiveIndexProfile(const RefractiveIndexProfile & a) : AtmProfile(a), SpectralGrid(a)

{

  //   std::cout<<"Enter RefractiveIndexProfile copy constructor version Fri May 20 00:59:47 CEST 2005"<<endl;

​

  // level AtmProfile

​

  // type_ = a.type_;

  // prLimit_ = a.prLimit_;

  /*

   v_hx_.reserve(a.v_hx_.size());

   v_px_.reserve(a.v_px_.size());