prof.c File Reference

Extract vertical profiles. More...

#include "libgps.h"

Go to the source code of this file.

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

Detailed Description

Extract vertical profiles.

Definition in file prof.c.

Function Documentation

main()

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

Definition at line 28 of file prof.c.

                {
 
  gps_t *gps;
 
  FILE *out;
 
  int ids, iz;
 
  /* Allocate... */
  ALLOC(gps, gps_t, 1);
 
  /* Check arguments... */
  if (argc < 4)
    ERRMSG("Give parameters: <ctl> <gps.nc> <prof.tab>");
 
  /* Get control parameters... */
  ids = (int) scan_ctl(argc, argv, "IDS", -1, "0", NULL);
 
  /* Read gps data... */
  read_gps(argv[2], gps);
 
  /* Check profile number... */
  if (ids < 0 || ids > gps->nds)
    ERRMSG("Profile not available!");
 
  /* Create output file... */
  printf("Write profile data: %s\n", argv[3]);
  if (!(out = fopen(argv[3], "w")))
    ERRMSG("Cannot create file!");
 
  /* Write header... */
  fprintf(out,
          "# $1 = time [sec]\n"
          "# $2 = altitude [km]\n"
          "# $3 = longitude [deg]\n"
          "# $4 = latitude [deg]\n"
          "# $5 = pressure [hPa]\n"
          "# $6 = temperature [K]\n"
          "# $7 = water vapor vmr [ppm]\n"
          "# $8 = temperature perturbation [K]\n"
          "# $9 = tropopause height [km]\n\n");
 
  /* Write data... */
  for (iz = 0; iz < gps->nz[ids]; iz++)
    if (gps->z[ids][iz] > 0)
      fprintf(out, "%.2f %g %g %g %g %g %g %g %g\n",
              gps->time[ids], gps->z[ids][iz], gps->lon[ids][iz],
              gps->lat[ids][iz], gps->p[ids][iz], gps->t[ids][iz],
              gps->wv[ids][iz], gps->pt[ids][iz], gps->th[ids]);
 
  /* Close file... */
  fclose(out);
 
  /* Free... */
  free(gps);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: