interpolate.c File Reference

Interpolate atmospheric data to another spatial grid. More...

#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

Interpolate atmospheric data to another spatial grid.

Definition in file interpolate.c.

Function Documentation

main()

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

Definition at line 39 of file interpolate.c.

                {
 
  static atm_t atm_in, atm_pts;
  static ctl_t ctl;
 
  double k[NW], q[NG];
 
  /* Print usage information... */
  USAGE;
 
  /* Interpolate atmospheric data... */
 
  /* Check arguments... */
  if (argc < 5)
    ERRMSG("Missing or invalid command-line arguments.\n\n"
           "Usage: interpolate <ctl> <atm_in> <atm_pts> <atm_out> [KEY VALUE ...]\n\n"
           "Use -h for full help.");
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
 
  /* Read atmospheric data... */
  read_atm(NULL, argv[2], &ctl, &atm_in, 0);
  read_atm(NULL, argv[3], &ctl, &atm_pts, 0);
 
  /* Interpolate atmospheric data... */
  for (int ip = 0; ip < atm_pts.np; ip++) {
    intpol_atm(&ctl, &atm_in, atm_pts.z[ip],
               &atm_pts.p[ip], &atm_pts.t[ip], q, k);
    for (int ig = 0; ig < ctl.ng; ig++)
      atm_pts.q[ig][ip] = q[ig];
    for (int iw = 0; iw < ctl.nw; iw++)
      atm_pts.k[iw][ip] = k[iw];
  }
 
  /* Save interpolated data... */
  write_atm(NULL, argv[4], &ctl, &atm_pts, 0);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: