mc_atmprof.h

Go to the documentation of this file.

/* ============================================================================

   Copyright (C) 2019  Konrad Bernloehr

   This file is part of the eventio/hessio library and other packages.

   The eventio/hessio 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, see <http://www.gnu.org/licenses/>.

============================================================================ */

/* ========================================================= */

#ifndef HAVE_MC_ATMPROF

# define HAVE_MC_ATMPROF 1

#ifdef __cplusplus
extern "C" {
#endif

struct atmospheric_profile
{
   int atmprof_id;     
   char *atmprof_fname;
   double obslev;      
   unsigned n_alt;     
   double *alt_km;     
   double *rho;        
   double *thick;      
   double *refidx_m1;  
   int have_lay5_param;
   double hlay[6];     
   double aatm[5];     
   double batm[5];     
   double catm[5];     
   double datm[5];     
   double thickl[6];   
   double htoa;        
};
typedef struct atmospheric_profile AtmProf;

/* Make the common profile available */
AtmProf *get_common_atmprof(void);
/* Set the common profile from a separate copy. */
void set_common_atmprof(AtmProf *atmprof);
void show_atmprof (AtmProf *atmprof);

void atmegs_(int *nlay, double *hlay,  double *aatm, 
   double *batm, double *catm, double *datm, double *htoa);
void atmegs_default(void);

double rhofc (double *height);
double thickc (double *height);
double refidc (double *height);
double refim1c (double *height);
double heighc (double *thick);

#ifdef __cplusplus
}
#endif

#endif