extract.c File Reference

Extract radiance data for retrievals. More...

#include "libairs.h"

Go to the source code of this file.

Functions
double	gph2z (double gph)
int	main (int argc, char *argv[])

Variables
int	airs_chan [L1_NCHAN]

Detailed Description

Extract radiance data for retrievals.

Definition in file extract.c.

Function Documentation

gph2z()

double gph2z

(

double

gph

)

Definition at line 163 of file extract.c.

              {
 
  double a = 3.086e-3;
 
  return G0 / a - sqrt(gsl_pow_2(G0 / a) - 2 * G0 * gph / a);
}

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 51 of file extract.c.

                {
 
  static airs_rad_gran_t airs_rad_gran;
  static airs_ret_gran_t airs_ret_gran;
 
  static airs_l1_t l1;
  static airs_l2_t l2;
 
  /* Check arguments... */
  if (argc != 4)
    ERRMSG("Give parameters: <airs_l1_file> <airs_l2_file> <out.nc>");
 
  /* Check Level-1 filename... */
  if (argv[1][0] != '-') {
 
    /* Read data... */
    printf("Read AIRS Level-1 file: %s\n", argv[1]);
    airs_rad_rdr(argv[1], &airs_rad_gran);
 
    /* Flag bad data... */
    for (int track = 0; track < AIRS_RAD_GEOTRACK; track++)
      for (int xtrack = 0; xtrack < AIRS_RAD_GEOXTRACK; xtrack++)
        for (int ichan = 0; ichan < L1_NCHAN; ichan++)
          if ((airs_rad_gran.state[track][xtrack] != 0)
              || (airs_rad_gran.ExcludedChans[airs_chan[ichan]] > 2)
              || (airs_rad_gran.CalChanSummary[airs_chan[ichan]] & 8)
              || (airs_rad_gran.CalChanSummary[airs_chan[ichan]] & (32 + 64))
              || (airs_rad_gran.CalFlag[track][airs_chan[ichan]] & 16))
            airs_rad_gran.radiances[track][xtrack][airs_chan[ichan]]
              = GSL_NAN;
 
    /* Copy data to struct... */
    for (int track = 0; track < AIRS_RAD_GEOTRACK; track++)
      for (int xtrack = 0; xtrack < AIRS_RAD_GEOXTRACK; xtrack++) {
        l1.time[track][xtrack]
          = airs_rad_gran.Time[track][xtrack] - 220838400.;
        l1.lon[track][xtrack]
          = airs_rad_gran.Longitude[track][xtrack];
        l1.lat[track][xtrack]
          = airs_rad_gran.Latitude[track][xtrack];
        l1.sat_z[track]
          = airs_rad_gran.satheight[track];
        l1.sat_lon[track]
          = airs_rad_gran.sat_lon[track];
        l1.sat_lat[track]
          = airs_rad_gran.sat_lat[track];
        for (int ichan = 0; ichan < L1_NCHAN; ichan++) {
          l1.nu[ichan]
            = airs_rad_gran.nominal_freq[airs_chan[ichan]];
          l1.rad[track][xtrack][ichan]
            = airs_rad_gran.radiances[track][xtrack][airs_chan[ichan]] *
            0.001f;
        }
      }
 
    /* Write netCDF file... */
    write_l1(argv[3], &l1);
  }
 
  /* Check Level-2 filename... */
  if (argv[2][0] != '-') {
 
    /* Read data... */
    printf("Read AIRS Level-2 file: %s\n", argv[2]);
    airs_ret_rdr(argv[2], &airs_ret_gran);
 
    /* Flag bad data... */
    for (int track = 0; track < AIRS_RET_GEOTRACK; track++)
      for (int xtrack = 0; xtrack < AIRS_RET_GEOXTRACK; xtrack++)
        for (int lay = 1; lay < AIRS_RET_STDPRESSURELAY; lay++)
          if (airs_ret_gran.GP_Height[track][xtrack][lay] <= -9000.
              || airs_ret_gran.TAirStd[track][xtrack][lay] <= -9000.) {
            airs_ret_gran.GP_Height[track][xtrack][lay] = GSL_NAN;
            airs_ret_gran.TAirStd[track][xtrack][lay] = GSL_NAN;
          }
 
    /* Save data in struct... */
    for (int track = 0; track < AIRS_RET_GEOTRACK; track++)
      for (int xtrack = 0; xtrack < AIRS_RET_GEOXTRACK; xtrack++)
        for (int lay = 1; lay < AIRS_RET_STDPRESSURELAY; lay++) {
          l2.time[track][xtrack]
            = airs_ret_gran.Time[track][xtrack] - 220838400.;
          l2.z[track][xtrack][lay - 1]
            = airs_ret_gran.GP_Height[track][xtrack][lay] / 1000.;
          l2.lon[track][xtrack]
            = airs_ret_gran.Longitude[track][xtrack];
          l2.lat[track][xtrack]
            = airs_ret_gran.Latitude[track][xtrack];
          l2.p[lay - 1]
            = airs_ret_gran.pressStd[lay];
          l2.t[track][xtrack][lay - 1]
            = airs_ret_gran.TAirStd[track][xtrack][lay];
        }
 
    /* Convert geopotential heights to geometric heights... */
    for (int track = 0; track < L2_NTRACK; track++)
      for (int xtrack = 0; xtrack < L2_NXTRACK; xtrack++)
        for (int lay = 0; lay < L2_NLAY; lay++)
          l2.z[track][xtrack][lay]
            = gph2z(l2.z[track][xtrack][lay]);
 
    /* Write netCDF file... */
    write_l2(argv[3], &l2);
  }
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function:

Variable Documentation

airs_chan

int airs_chan[L1_NCHAN]

Initial value:

= { 54, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82,
  2035, 2036, 2040, 2041, 2052, 2053, 2054, 2055,
  2067, 2075, 2076, 2077, 2078, 2079, 2080, 2081,
  2082, 2086, 2088, 2089, 2091, 2092, 2093
}

Definition at line 33 of file extract.c.