noise_pert.c File Reference

Estimate noise based on perturbation data. More...

#include "libairs.h"

Go to the source code of this file.

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

Detailed Description

Estimate noise based on perturbation data.

Definition in file noise_pert.c.

Function Documentation

main()

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

Definition at line 28 of file noise_pert.c.

                {
 
  static pert_t *pert;
  static wave_t wave;
 
  FILE *out;
 
  char pertname[LEN];
 
  double mu, nedt = -1e99;
 
  /* Check arguments... */
  if (argc < 4)
    ERRMSG("Give parameters: <ctl> <pert.nc> <noise.tab>");
 
  /* Read control parameters... */
  scan_ctl(argc, argv, "PERTNAME", -1, "4mu", pertname);
  int bsize = (int) scan_ctl(argc, argv, "BSIZE", -1, "-999", NULL);
  const double maxvar = scan_ctl(argc, argv, "MAXVAR", -1, "-999", NULL);
 
  /* Allocate... */
  ALLOC(pert, pert_t, 1);
 
  /* Read perturbation data... */
  read_pert(argv[2], pertname, pert);
 
  /* Set block size... */
  if (bsize < 0)
    bsize = pert->nxtrack;
 
  /* 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 = longitude [deg]\n"
          "# $2 = latitude [deg]\n"
          "# $3 = mean brightness temperature [K]\n"
          "# $4 = noise estimate [K]\n\n");
 
  /* Loop over granules... */
  for (int itrack = 0; itrack < pert->ntrack; itrack += bsize) {
 
    /* Convert retrieval data to wave struct... */
    pert2wave(pert, &wave, itrack, itrack + bsize,
              pert->nxtrack / 2 - bsize / 2, pert->nxtrack / 2 + bsize / 2);
 
    /* Estimate noise... */
    double nedt_old = nedt;
    noise(&wave, &mu, &nedt);
 
    /* Write output... */
    if (maxvar <= 0
        || fabs(200 * (nedt - nedt_old) / (nedt + nedt_old)) < maxvar)
      fprintf(out, "%g %g %g %g\n", wave.lon[wave.nx / 2][wave.ny / 2],
              wave.lat[wave.nx / 2][wave.ny / 2], mu, nedt);
  }
 
  /* Close file... */
  fclose(out);
 
  /* Free... */
  free(pert);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: