Tools/ShowerUniversality/UnivParamNS/Atmosphere.h

Go to the documentation of this file.

#ifndef __Atmosphere_h_
#define __Atmosphere_h_


#include <stdio.h>
#include <stdlib.h>
#include <string>

namespace AtmosphereNS {

  const float c_cmperns = 29.9792458;
  const float hGroundRef = 1452.e2;

  class AtmModel {
  public:
    float ai[5], bi[5], ci[5];
    float hi[5]; // cm
    int atmflag; // 1...12 Montly models  10,...,50 Seasonal models
  };

  class Atmosphere {
  private:

    AtmModel amonth[12];
    AtmModel aseason[5];
    AtmModel acurrent;
    bool ok; // Initialized to false, only set to true when acurrent is set
    std::string atmType; // currently implemented: "monthly", "seasonal"

  public:
    Atmosphere(const bool UseGDAS = false, const std::string atmtype = "monthly");
    ~Atmosphere();
    AtmModel* GetCurrentAtmosphere()
    {
      return &acurrent;
    };
    int GetCurrentAtmosphereMonth()
    {
      return acurrent.atmflag;
    };
    void SetCurrentAtmosphere(AtmModel theAtm);
    void SetCurrentAtmosphere(int imonth);  // imonth=1,...,12
    void SetCurrentAtmosphere(float* ai, float* bi, float* ci, float* hi);
    void printCurrentAtmosphere(FILE* fp);
    void ClearCurrentAtmosphere();

    //Using current atmosphere
    float VertDepthToHeight(float vdepth);
    float SlantDepthToHeight(float sdepth, float theta);
    float HeightToVertDepth(float h);
    float HeightToSlantDepth(float h, float theta);
    float GetDensity_i(float h);
    float GetDensity(float h);
    float GetDensityAboveGround(float hground, float theta, float nRad);
    float GetDensityAboveGround_i(float hground, float theta, float nRad);
    float GetSlant0()
    {
      return acurrent.ai[4] * acurrent.ci[4] / acurrent.bi[4];
    };
    float GetC0()
    {
      return acurrent.ci[0];
    };

    float Get_DX_DL(float r, float psi, float SlantDepth, float theta, float hground, bool UseDL, bool IsDiffusive);
    float Get_DX_DL_i(float r, float psi, float SlantDepth, float theta, float hground, bool UseDL, bool IsDiffusive);

    float GetTimePF_delay(float x, float y, float z, float theta, float azi);

    float GetTimePF(float r, float psi, float SlantDepth, float theta, float hground);
    float GetTimePF_i(float r, float psi, float SlantDepth, float theta, float hground);
    float GetTimePF_h(float r, float psi, float hStart, float theta, float hground);

    float GetTimeCF_h(float r, float psi, float h1, float theta, float hground);
    float GetTimeCF(float r, float psi, float SlantDepth, float theta, float hground);
    float GetTimeCF_i(float r, float psi, float SlantDepth, float theta, float hground);

  };
}
#endif // __Atmosphere_h_