noise.c File Reference

Get noise estimates from radiance data. More...

#include "libiasi.h"

Go to the source code of this file.

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

Detailed Description

Get noise estimates from radiance data.

Definition in file noise.c.

Function Documentation

main()

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

Definition at line 28 of file noise.c.

                {
 
  static iasi_rad_t *iasi_rad;
 
  static wave_t wave;
 
  static FILE *out;
 
  static double mu, nesr, sigma;
 
  static int ichan, itrack, ix, iy, format;
 
  /* Check arguments... */
  if (argc < 4)
    ERRMSG("Give parameters: <ctl> <iasi_l1_file> <noise.tab>");
 
  /* Read control parameters... */
  format = (int) scan_ctl(argc, argv, "FORMAT", -1, "1", NULL);
 
  /* Allocate... */
  ALLOC(iasi_rad, iasi_rad_t, 1);
 
  /* Read IASI data... */
  printf("Read IASI data: %s\n", argv[2]);
  iasi_read(format, argv[2], iasi_rad);
 
  /* Create file... */
  printf("Write noise data: %s\n", argv[3]);
  if (!(out = fopen(argv[3], "w")))
    ERRMSG("Cannot create file!");
 
  /* Write header... */
  fprintf(out,
          "# $1 = track index\n"
          "# $2 = channel index\n"
          "# $3 = wavenumber [1/cm]\n"
          "# $4 = mean BT [K]\n"
          "# $5 = NEDT [K]\n" "# $6 = NESR [W/(m^2 sr cm^-1)]\n");
 
  /* Analyze blocks of data... */
  for (itrack = 0; itrack < iasi_rad->ntrack; itrack += 60) {
 
    /* Write empty line... */
    fprintf(out, "\n");
 
    /* Loop over channels... */
    for (ichan = 0; ichan < IASI_L1_NCHAN; ichan++) {
 
      /* Set wave struct... */
      wave.nx = L1_NXTRACK;
      wave.ny = 0;
      for (iy = itrack; iy < GSL_MIN(itrack + 60, iasi_rad->ntrack); iy++) {
        for (ix = 0; ix < wave.nx; ix++)
          wave.temp[ix][wave.ny] = BRIGHT(iasi_rad->Rad[iy][ix][ichan],
                                          iasi_rad->freq[ichan]);
        wave.ny++;
      }
 
      /* Check number of data points... */
      if (wave.ny >= 55) {
 
        /* Get noise... */
        noise(&wave, &mu, &sigma);
 
        /* Get NESR... */
        nesr = PLANCK(mu + sigma, iasi_rad->freq[ichan])
          - PLANCK(mu, iasi_rad->freq[ichan]);
 
        /* Write output... */
        if (gsl_finite(sigma))
          fprintf(out, "%d %d %.4f %g %g %g\n", itrack, ichan,
                  iasi_rad->freq[ichan], mu, sigma, nesr);
      }
    }
  }
 
  /* Close file... */
  fclose(out);
 
  /* Free... */
  free(iasi_rad);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: