met_sample.c File Reference

Sample meteorological data at given geolocations. More...

#include "mptrac.h"

Go to the source code of this file.

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

Detailed Description

Sample meteorological data at given geolocations.

Definition in file met_sample.c.

Function Documentation

main()

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

Definition at line 31 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;
 
  /* Check arguments... */
  if (argc < 3)
    ERRMSG("Give parameters: <ctl> <sample.tab> <atm_in>");
 
  /* 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: