libiasi.h

Go to the documentation of this file.

/*
  This file is part of the IASI Code Collection.
  
  the IASI Code Collections is free software: you can redistribute it
  and/or modify it under the terms of the GNU General Public License
  as published by the Free Software Foundation, either version 3 of
  the License, or (at your option) any later version.
  
  The IASI Code Collection 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
  General Public License for more details.
  
  You should have received a copy of the GNU General Public License
  along with the IASI Code Collection. If not, see
  <http://www.gnu.org/licenses/>.
  
  Copyright (C) 2019-2026 Forschungszentrum Juelich GmbH
*/
 
#include <netcdf.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_fft_complex.h>
#include <gsl/gsl_multifit.h>
#include <gsl/gsl_poly.h>
#include <gsl/gsl_sort.h>
#include <gsl/gsl_spline.h>
#include "coda.h"
#include "jurassic.h"
 
/* ------------------------------------------------------------
   Dimensions...
   ------------------------------------------------------------ */
 
#define L1_NCHAN 33
 
#define L1_NTRACK 1800
 
#define L1_NXTRACK 60
 
#define L2_NLAY 27
 
#define L2_NTRACK 1800
 
#define L2_NXTRACK 60
 
#define IASI_L1_NCHAN 8700
 
#define IASI_NXTRACK 30
 
#define IASI_PM 4
 
#define IASI_IDefNsfirst1b 2581
 
#define IASI_IDefNslast1b 11041
 
#define IASI_IDefSpectDWn1b 25
 
#define PERT_NTRACK 132000
 
#define PERT_NXTRACK 360
 
#define WX 300
 
#define WY 33000
 
/* ------------------------------------------------------------
   Macros...
   ------------------------------------------------------------ */
 
#define CODA(cmd) {                                     \
    int coda_result=(cmd);                              \
    if(coda_result!=0)                                  \
      ERRMSG("%s", coda_errno_to_string(coda_errno));   \
  }
 
#define NC(cmd) {                                    \
  int nc_result=(cmd);                               \
  if(nc_result!=NC_NOERR)                            \
    ERRMSG("%s", nc_strerror(nc_result));            \
}
 
/* ------------------------------------------------------------
   Structs...
   ------------------------------------------------------------ */
 
typedef struct {
 
  size_t ntrack;
 
  double time[L1_NTRACK][L1_NXTRACK];
 
  double lon[L1_NTRACK][L1_NXTRACK];
 
  double lat[L1_NTRACK][L1_NXTRACK];
 
  double sat_z[L1_NTRACK];
 
  double sat_lon[L1_NTRACK];
 
  double sat_lat[L1_NTRACK];
 
  double nu[L1_NCHAN];
 
  float rad[L1_NTRACK][L1_NXTRACK][L1_NCHAN];
 
} iasi_l1_t;
 
typedef struct {
 
  size_t ntrack;
 
  double time[L2_NTRACK][L2_NXTRACK];
 
  double z[L2_NTRACK][L2_NXTRACK][L2_NLAY];
 
  double lon[L2_NTRACK][L2_NXTRACK];
 
  double lat[L2_NTRACK][L2_NXTRACK];
 
  double p[L2_NLAY];
 
  double t[L2_NTRACK][L2_NXTRACK][L2_NLAY];
 
} iasi_l2_t;
 
typedef struct {
 
  int ntrack;
 
  int nxtrack;
 
  double time[PERT_NTRACK][PERT_NXTRACK];
 
  double lon[PERT_NTRACK][PERT_NXTRACK];
 
  double lat[PERT_NTRACK][PERT_NXTRACK];
 
  double dc[PERT_NTRACK][PERT_NXTRACK];
 
  double bt[PERT_NTRACK][PERT_NXTRACK];
 
  double pt[PERT_NTRACK][PERT_NXTRACK];
 
  double var[PERT_NTRACK][PERT_NXTRACK];
 
} pert_t;
 
typedef struct {
 
  long ntrack;
 
  float IDefSpectDWn1b[L1_NTRACK];
 
  int32_t IDefNsfirst1b[L1_NTRACK];
 
  int32_t IDefNslast1b[L1_NTRACK];
 
  double Time[L1_NTRACK][IASI_NXTRACK];
 
  double Loc[L1_NTRACK][IASI_NXTRACK][IASI_PM][2];
 
  float Wavenumber[IASI_L1_NCHAN];
 
  short int Radiation[L1_NTRACK][IASI_NXTRACK][IASI_PM][IASI_L1_NCHAN];
 
  unsigned int Sat_z[L1_NTRACK];
 
} iasi_raw_t;
 
typedef struct {
 
  int ntrack;
 
  double freq[IASI_L1_NCHAN];
 
  double Time[L1_NTRACK][L1_NXTRACK];
 
  double Longitude[L1_NTRACK][L1_NXTRACK];
 
  double Latitude[L1_NTRACK][L1_NXTRACK];
 
  float Rad[L1_NTRACK][L1_NXTRACK][IASI_L1_NCHAN];
 
  double Sat_z[L1_NTRACK];
 
  double Sat_lon[L1_NTRACK];
 
  double Sat_lat[L1_NTRACK];
 
} iasi_rad_t;
 
typedef struct {
 
  int nx;
 
  int ny;
 
  double time;
 
  double z;
 
  double lon[WX][WY];
 
  double lat[WX][WY];
 
  double x[WX];
 
  double y[WY];
 
  double temp[WX][WY];
 
  double bg[WX][WY];
 
  double pt[WX][WY];
 
  double var[WX][WY];
 
} wave_t;
 
/* ------------------------------------------------------------
   Functions...
   ------------------------------------------------------------ */
 
void add_var(
  int ncid,
  const char *varname,
  const char *unit,
  const char *longname,
  int type,
  int dimid[],
  int *varid,
  int ndims);
 
void background_poly(
  wave_t * wave,
  int dim_x,
  int dim_y);
 
void background_poly_help(
  double *xx,
  double *yy,
  int n,
  int dim);
 
void background_smooth(
  wave_t * wave,
  int npts_x,
  int npts_y);
 
void gauss(
  wave_t * wave,
  double fwhm);
 
void hamming(
  wave_t * wave,
  int nit);
 
void iasi_read(
  int format,
  char *filename,
  iasi_rad_t * iasi_rad);
 
void iasi_read_native(
  char *filename,
  iasi_rad_t * iasi_rad);
 
void iasi_read_netcdf(
  char *filename,
  iasi_rad_t * iasi_rad);
 
void median(
  wave_t * wave,
  int dx);
 
void noise(
  wave_t * wave,
  double *mu,
  double *sig);
 
void pert2wave(
  pert_t * pert,
  wave_t * wave,
  int track0,
  int track1,
  int xtrack0,
  int xtrack1);
 
void read_pert(
  char *filename,
  char *pertname,
  pert_t * pert);
 
void variance(
  wave_t * wave,
  double dh);
 
double wgs84(
  double lat);
 
void write_l1(
  char *filename,
  iasi_l1_t * l1);
 
void write_l2(
  char *filename,
  iasi_l2_t * l2);