perturbation.c File Reference

Extract gravity wave perturbations. More...

#include "libgps.h"

Go to the source code of this file.

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

Detailed Description

Extract gravity wave perturbations.

Definition in file perturbation.c.

Function Documentation

main()

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

Definition at line 32 of file perturbation.c.

                {
 
  gps_t *gps;
 
  char metbase[LEN];
 
  /* Allocate... */
  ALLOC(gps, gps_t, 1);
 
  /* Check arguments... */
  if (argc < 4)
    ERRMSG("Give parameters: <ctl> <out.nc> <gps1.nc> [<gps2.nc> ...]");
 
  /* Get control parameters... */
  const double dt_met = scan_ctl(argc, argv, "DT_MET", -1, "21600", NULL);
  const double gauss_dx = scan_ctl(argc, argv, "GAUSS_DX", -1, "-999", NULL);
  const double gauss_dy = scan_ctl(argc, argv, "GAUSS_DY", -1, "-999", NULL);
  const double grid_zmin = scan_ctl(argc, argv, "GRID_ZMIN", -1, "0", NULL);
  const double grid_zmax = scan_ctl(argc, argv, "GRID_ZMAX", -1, "40", NULL);
  const int grid_nz = (int) scan_ctl(argc, argv, "GRID_NZ", -1, "-1", NULL);
  const double ham_dz = scan_ctl(argc, argv, "HAM_DZ", -1, "-999", NULL);
  const double ham_dz2 = scan_ctl(argc, argv, "HAM_DZ2", -1, "-999", NULL);
  scan_ctl(argc, argv, "METBASE", -1, "", metbase);
  const double prof_zmax_min =
    scan_ctl(argc, argv, "PROF_ZMAX_MIN", -1, "35", NULL);
  const double prof_zmin_max =
    scan_ctl(argc, argv, "PROF_ZMIN_MAX", -1, "5", NULL);
  const int poly_dim = (int) scan_ctl(argc, argv, "POLY_DIM", -1, "5", NULL);
  const double poly_zmin = scan_ctl(argc, argv, "POLY_ZMIN", -1, "0", NULL);
  const double poly_zmax = scan_ctl(argc, argv, "POLY_ZMAX", -1, "40", NULL);
 
  /* Read individual GPS-RO data files... */
  for (int iarg = 3; iarg < argc; iarg++) {
    FILE *in;
    if (!(in = fopen(argv[iarg], "r")))
      continue;
    else {
      fclose(in);
      read_gps_prof(argv[iarg], gps, prof_zmin_max, prof_zmax_min);
    }
  }
 
  /* Check number of profiles... */
  if (gps->nds <= 0)
    ERRMSG("No profiles found!");
 
  /* Grid profile... */
  if (grid_nz > 0)
    grid_gps(gps, grid_zmin, grid_zmax, grid_nz);
 
  /* Get tropopause... */
  tropopause(gps);
 
  /* Get perturbations from horizontal Gaussian mean... */
  if (gauss_dx > 0 && gauss_dy > 0)
    gauss(gps, gauss_dx, gauss_dy);
 
  /* Get perturbations from vertical Hamming filter... */
  if (ham_dz > 0)
    hamming_low_pass(gps, ham_dz);
 
  /* Use vertical Hamming filter to reduce noise... */
  if (ham_dz2 > 0)
    hamming_high_pass(gps, ham_dz2);
 
  /* Use meteo data for detrending... */
  if (metbase[0] != '-')
    detrend_met(gps, metbase, dt_met);
 
  /* Remove polynomial fit from perturbation profile... */
  if (poly_dim > 0)
    poly(gps, poly_dim, poly_zmin, poly_zmax);
 
  /* Write GPS-RO data file... */
  write_gps(argv[2], gps);
 
  /* Free... */
  free(gps);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: