kernel.c File Reference

Calculate kernel functions. More...

#include "jurassic.h"

Go to the source code of this file.

Functions
void	call_kernel (ctl_t *ctl, const char *wrkdir, const char *obsfile, const char *atmfile, const char *kernelfile)
	Perform kernel calculations in a single directory. More...
int	main (int argc, char *argv[])

Detailed Description

Calculate kernel functions.

Definition in file kernel.c.

Function Documentation

call_kernel()

void call_kernel	(	ctl_t *	ctl,
		const char *	wrkdir,
		const char *	obsfile,
		const char *	atmfile,
		const char *	kernelfile
	)

Perform kernel calculations in a single directory.

Definition at line 96 of file kernel.c.

                          {
 
  static atm_t atm;
  static obs_t obs;
 
  /* Read observation geometry... */
  read_obs(wrkdir, obsfile, ctl, &obs);
 
  /* Read atmospheric data... */
  read_atm(wrkdir, atmfile, ctl, &atm);
 
  /* Get sizes... */
  const size_t n = atm2x(ctl, &atm, NULL, NULL, NULL);
  const size_t m = obs2y(ctl, &obs, NULL, NULL, NULL);
 
  /* Check sizes... */
  if (n == 0)
    ERRMSG("No state vector elements!");
  if (m == 0)
    ERRMSG("No measurement vector elements!");
 
  /* Allocate... */
  gsl_matrix *k = gsl_matrix_alloc(m, n);
 
  /* Compute kernel matrix... */
  kernel(ctl, &atm, &obs, k);
 
  /* Write matrix to file... */
  write_matrix(wrkdir, kernelfile, ctl, k, &atm, &obs, "y", "x", "r");
 
  /* Free... */
  gsl_matrix_free(k);
}

Here is the call graph for this function:

main()

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

Definition at line 43 of file kernel.c.

                {
 
  static ctl_t ctl;
 
  char dirlist[LEN];
 
  /* Check arguments... */
  if (argc < 5)
    ERRMSG("Give parameters: <ctl> <obs> <atm> <kernel>");
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
 
  /* Get dirlist... */
  scan_ctl(argc, argv, "DIRLIST", -1, "-", dirlist);
 
  /* Set flags... */
  ctl.write_matrix = 1;
 
  /* Single kernel calculation... */
  if (dirlist[0] == '-')
    call_kernel(&ctl, NULL, argv[2], argv[3], argv[4]);
 
  /* Work on directory list... */
  else {
 
    /* Open directory list... */
    FILE *in;
    if (!(in = fopen(dirlist, "r")))
      ERRMSG("Cannot open directory list!");
 
    /* Loop over directories... */
    char wrkdir[LEN];
    while (fscanf(in, "%4999s", wrkdir) != EOF) {
 
      /* Write info... */
      LOG(1, "\nWorking directory: %s", wrkdir);
 
      /* Call forward model... */
      call_kernel(&ctl, wrkdir, argv[2], argv[3], argv[4]);
    }
 
    /* Close dirlist... */
    fclose(in);
  }
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: