met_sample.c File Reference

Sample meteorological data at given geolocations. More...

#include "mptrac.h"

Go to the source code of this file.

Functions
void	usage (void)
	Print command-line help. More...
int	main (int argc, char *argv[])

Detailed Description

Sample meteorological data at given geolocations.

Definition in file met_sample.c.

Function Documentation

usage()

void usage

(

void

)

Print command-line help.

Definition at line 185 of file met_sample.c.

        {
 
  printf("\nMPTRAC met_sample tool.\n\n");
  printf("Sample meteorological data at atmospheric particle locations.\n");
  printf("\n");
  printf("Usage:\n");
  printf("  met_sample <ctl> <sample.tab> <atm_in> [KEY VALUE ...]\n");
  printf("\n");
  printf("Arguments:\n");
  printf("  <ctl>         Control file.\n");
  printf("  <sample.tab>  Output table.\n");
  printf("  <atm_in>      Atmospheric input file with sample locations.\n");
  printf("  [KEY VALUE]  Optional control parameters.\n");
  printf("\nFurther information:\n");
  printf("  Manual: https://slcs-jsc.github.io/mptrac/\n");
}

main()

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

Definition at line 39 of file met_sample.c.

                {
 
  ctl_t ctl;
 
  clim_t *clim;
 
  atm_t *atm;
 
  met_t *met0, *met1;
 
  dd_t *dd;
 
  FILE *out;
 
  double h2o, h2ot, o3, lwc, rwc, iwc, swc, cc, p0, p1, ps, ts, zs, us, vs,
    ess, nss, shf, lsm, sst, pbl, pt, pct, pcb, cl, plcl, plfc, pel,
    cape, cin, o3c, pv, t, tt, u, v, w, z, zm, zref, zt, time_old = -999,
    p_old = -999, lon_old = -999, lat_old = -999;
 
  /* Print usage information... */
  USAGE;
 
  /* Check arguments... */
  if (argc < 3)
    ERRMSG("Missing or invalid command-line arguments.\n\n"
           "Usage: met_sample <ctl> <sample.tab> <atm_in> [KEY VALUE ...]\n\n"
           "Use -h for full help.");
 
  /* Allocate... */
  ALLOC(clim, clim_t, 1);
  ALLOC(atm, atm_t, 1);
  ALLOC(met0, met_t, 1);
  ALLOC(met1, met_t, 1);
  ALLOC(dd, dd_t, 1);
 
  /* Read control parameters... */
  mptrac_read_ctl(argv[1], argc, argv, &ctl);
  const int geopot =
    (int) scan_ctl(argv[1], argc, argv, "SAMPLE_GEOPOT", -1, "0", NULL);
  const int grid_time =
    (int) scan_ctl(argv[1], argc, argv, "SAMPLE_GRID_TIME", -1, "0", NULL);
  const int grid_z =
    (int) scan_ctl(argv[1], argc, argv, "SAMPLE_GRID_Z", -1, "0", NULL);
  const int grid_lon =
    (int) scan_ctl(argv[1], argc, argv, "SAMPLE_GRID_LON", -1, "0", NULL);
  const int grid_lat =
    (int) scan_ctl(argv[1], argc, argv, "SAMPLE_GRID_LAT", -1, "0", NULL);
 
  /* Read climatological data... */
  mptrac_read_clim(&ctl, clim);
 
  /* Read atmospheric data... */
  if (!mptrac_read_atm(argv[3], &ctl, atm))
    ERRMSG("Cannot open file!");
 
  /* Create output file... */
  LOG(1, "Write meteorological data file: %s", argv[2]);
  if (!(out = fopen(argv[2], "w")))
    ERRMSG("Cannot create file!");
 
  /* Write header... */
  MET_HEADER;
 
  /* Loop over air parcels... */
  for (int ip = 0; ip < atm->np; ip++) {
 
    /* Get meteorological data... */
    mptrac_get_met(&ctl, clim, atm->time[ip], &met0, &met1, dd);
 
    /* Set reference pressure for interpolation... */
    INTPOL_INIT;
    double pref = atm->p[ip];
    if (geopot) {
      zref = Z(pref);
      p0 = met0->p[0];
      p1 = met0->p[met0->np - 1];
      for (int it = 0; it < 24; it++) {
        pref = 0.5 * (p0 + p1);
        intpol_met_time_3d(met0, met0->z, met1, met1->z, atm->time[ip], pref,
                           atm->lon[ip], atm->lat[ip], &zm, ci, cw, 1);
        if (zref > zm || !isfinite(zm))
          p0 = pref;
        else
          p1 = pref;
      }
      pref = 0.5 * (p0 + p1);
    }
 
    /* Interpolate meteo data... */
    INTPOL_TIME_ALL(atm->time[ip], pref, atm->lon[ip], atm->lat[ip]);
 
    /* Make blank lines... */
    if (ip == 0 || (grid_time && atm->time[ip] != time_old)
        || (grid_z && atm->p[ip] != p_old)
        || (grid_lon && atm->lon[ip] != lon_old)
        || (grid_lat && atm->lat[ip] != lat_old))
      fprintf(out, "\n");
    time_old = atm->time[ip];
    p_old = atm->p[ip];
    lon_old = atm->lon[ip];
    lat_old = atm->lat[ip];
 
    /* Write data... */
    fprintf(out,
            "%.2f %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g"
            " %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g"
            " %g %g %g %g %g %g %g %g %g %g %g %g %g %g 1 1 1\n",
            atm->time[ip], Z(atm->p[ip]), atm->lon[ip], atm->lat[ip],
            atm->p[ip], t, u, v, w, h2o, o3, z, pv, ps, ts, zs, us, vs,
            ess, nss, shf, lsm,
            sst, pt, zt, tt, h2ot, lwc, rwc, iwc, swc, cc, cl, pct, pcb, plcl,
            plfc, pel, cape, cin, RH(atm->p[ip], t, h2o), RHICE(atm->p[ip], t,
                                                                h2o),
            TDEW(atm->p[ip], h2o), TICE(atm->p[ip], h2o),
            nat_temperature(atm->p[ip], h2o,
                            clim_zm(&clim->hno3, atm->time[ip], atm->lat[ip],
                                    atm->p[ip])), clim_zm(&clim->hno3,
                                                          atm->time[ip],
                                                          atm->lat[ip],
                                                          atm->p[ip]),
            clim_oh(&ctl, clim, atm->time[ip], atm->lon[ip], atm->lat[ip],
                    atm->p[ip]), clim_zm(&clim->h2o2, atm->time[ip],
                                         atm->lat[ip], atm->p[ip]),
            clim_zm(&clim->ho2, atm->time[ip], atm->lat[ip], atm->p[ip]),
            clim_zm(&clim->o1d, atm->time[ip], atm->lat[ip], atm->p[ip]), pbl,
            o3c);
  }
 
  /* Close file... */
  fclose(out);
 
  /* Free... */
  free(clim);
  free(atm);
  free(met0);
  free(met1);
  free(dd);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: