casa6

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

 * ALMA - Atacama Large Millimeter Array

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

 * (in the framework of the ALMA collaboration).

 * All rights reserved.

 * 

 * 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

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

​

#include <string>

#include <vector>

#include <iostream>

#include <fstream>

#include <sstream>

#include <math.h>

​

using namespace std;

​

#include <atmosphere/ATM/ATMPercent.h>

#include <atmosphere/ATM/ATMPressure.h>

#include <atmosphere/ATM/ATMNumberDensity.h>

#include <atmosphere/ATM/ATMMassDensity.h>

#include <atmosphere/ATM/ATMTemperature.h>

#include <atmosphere/ATM/ATMLength.h>

#include <atmosphere/ATM/ATMInverseLength.h>

#include <atmosphere/ATM/ATMOpacity.h>

#include <atmosphere/ATM/ATMAngle.h>

#include <atmosphere/ATM/ATMHumidity.h>

#include <atmosphere/ATM/ATMFrequency.h>

#include <atmosphere/ATM/ATMWaterVaporRadiometer.h>

#include <atmosphere/ATM/ATMWVRMeasurement.h>

#include <atmosphere/ATM/ATMProfile.h>

#include <atmosphere/ATM/ATMSpectralGrid.h>

#include <atmosphere/ATM/ATMRefractiveIndex.h>

#include <atmosphere/ATM/ATMRefractiveIndexProfile.h>

#include <atmosphere/ATM/ATMSkyStatus.h>

​

int main()

{  

  // Initialize the Atmospheric model

​

  unsigned int      atmType = 1;       // 1=tropical (ALMA site), 2=midlatSummer, 3=midlatWinter 

  atm::Temperature  T(273. ,"K" );     // Ground temperature

  atm::Pressure     P(55000. ,"Pa");   // Ground Pressure

  atm::Humidity     H(8.,"%" );       // Ground Relative Humidity (indication)

  atm::Length       Alt(5000,"m" );    // Altitude of the site

  atm::Length       WVL(   3.0,"km");  // Water vapor scale height

  double            TLR=  -6.5      ;  // Tropospheric lapse rate (must be in K/km)

  atm::Length       topAtm(  48.0,"km");  // Upper atm. boundary for calculations

  atm::Pressure     Pstep(  10.0,"mb"); // Primary pressure step (10.0 mb)

  double   PstepFact=         1.2;     // Pressure step ratio between two consecutive layers