nadir.c File Reference

Create observation geometry for a nadir sounder. More...

#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

Create observation geometry for a nadir sounder.

Definition in file nadir.c.

Function Documentation

main()

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

Definition at line 39 of file nadir.c.

                {
 
  static ctl_t ctl;
  static obs_t obs;
 
  /* Print usage information... */
  USAGE;
 
  /* Check arguments... */
  if (argc < 3)
    ERRMSG("Missing or invalid command-line arguments.\n\n"
           "Usage: nadir <ctl> <obs> [KEY VALUE ...]\n\n"
           "Use -h for full help.");
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
  const double t0 = scan_ctl(argc, argv, "T0", -1, "0", NULL);
  const double t1 = scan_ctl(argc, argv, "T1", -1, "0", NULL);
  const double dt = scan_ctl(argc, argv, "DT", -1, "1", NULL);
  const double obsz = scan_ctl(argc, argv, "OBSZ", -1, "700", NULL);
  const double lat0 = scan_ctl(argc, argv, "LAT0", -1, "-8.01", NULL);
  const double lat1 = scan_ctl(argc, argv, "LAT1", -1, "8.01", NULL);
  const double dlat = scan_ctl(argc, argv, "DLAT", -1, "0.18", NULL);
 
  /* Create measurement geometry... */
  for (double t = t0; t <= t1; t += dt)
    for (double lat = lat0; lat <= lat1; lat += dlat) {
      obs.time[obs.nr] = t;
      obs.obsz[obs.nr] = obsz;
      obs.vplat[obs.nr] = lat;
      if ((++obs.nr) >= NR)
        ERRMSG("Too many rays!");
    }
 
  /* Write observation data... */
  write_obs(NULL, argv[2], &ctl, &obs, 0);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: