casa5

//# ATAtmosphere.cc: Model of atmospheric opacity

//# Copyright (C) 2004

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

//#

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

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

//# the Free Software Foundation; either version 2 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 Library General Public

//# License for more details.

//#

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

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

//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.

//#

//# Correspondence concerning AIPS++ should be addressed as follows:

//#        Internet email: aips2-request@nrao.edu.

//#        Postal address: AIPS++ Project Office

//#                        National Radio Astronomy Observatory

//#                        520 Edgemont Road

//#                        Charlottesville, VA 22903-2475 USA

//#

//#---------------------------------------------------------------------------

​

// own includes

#include "ATAtmosphere.h"

​

// casa includes

#include <casa/Utilities/Assert.h>

#include <casa/Quanta.h>

​

// std includes

//#include <cmath>

​

using namespace casacore;

namespace casa{

​

/**

 * Default Constructor (apart from optional parameters).

 * The class set up this way will assume International Standard Atmosphere (ISA) conditions,

 * except for humidity. The latter is assumed to be 50%, which seems more realistic for 

 * Australian telescopes than 0%. 

 * @param[in] wvScale water vapour scale height (m), default is 1540m to match MIRIAD's model

 * @param[in] maxAlt maximum altitude of the model atmosphere (m), plane parallel layers are spread linearly up to

 *            this height, default is 10000m to match MIRIAD.

 * @param[in] nLayers number of plane parallel layers in the model (essentially for a numberical integration),

 *            default is 50 to match MIRIAD.

 **/

ATAtmosphere::ATAtmosphere(Double wvScale, Double maxAlt, Int nLayers) :

   itsHeights(nLayers), itsTemperatures(nLayers), 

   itsDryPressures(nLayers), itsVapourPressures(nLayers),

   itsGndTemperature(288.), itsPressure(101325.), itsGndHumidity(0.5), 

   itsLapseRate(0.0065), itsWVScale(wvScale), itsMaxAlt(maxAlt), itsObsHeight(200.)

{

  recomputeAtmosphereModel();

}

​

/**

 * Constructor with explicitly given parameters of the atmosphere 

 * @param[in] temperature air temperature at the observatory (K)