retrieval.c File Reference

JURASSIC retrieval processor. More...

#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

JURASSIC retrieval processor.

Definition in file retrieval.c.

Function Documentation

main()

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

Definition at line 35 of file retrieval.c.

                {
 
  static atm_t atm_i, atm_apr;
  static ctl_t ctl;
  static obs_t obs_i, obs_meas;
  static ret_t ret;
 
  FILE *dirlist;
 
  /* MPI task distribution (optional)... */
  int ntask = -1;
  int rank = 0;
  int size = 1;
 
#ifdef MPI
  MPI_Init(&argc, &argv);
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
  MPI_Comm_size(MPI_COMM_WORLD, &size);
#endif
 
  /* Check arguments... */
  if (argc < 3)
    ERRMSG("Give parameters: <ctl> <dirlist>");
 
  /* Measure CPU-time... */
  TIMER("total", 1);
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
  read_ret(argc, argv, &ctl, &ret);
 
  /* Initialize look-up tables... */
  tbl_t *tbl = read_tbl(&ctl);
 
  /* Open directory list... */
  if (!(dirlist = fopen(argv[2], "r")))
    ERRMSG("Cannot open directory list!");
 
  /* Loop over directories... */
  while (fscanf(dirlist, "%4999s", ret.dir) != EOF) {
 
    /* Distribute directories with MPI (optional)... */
    if ((++ntask) % size != rank)
      continue;
 
    /* Write info... */
    if (size > 1) {
      LOG(1, "\nRetrieve in directory %s on rank %d of %d...\n",
          ret.dir, rank + 1, size);
    } else
      LOG(1, "\nRetrieve in directory %s...\n", ret.dir);
 
    /* Read atmospheric data... */
    read_atm(ret.dir, "atm_apr.tab", &ctl, &atm_apr);
 
    /* Read observation data... */
    read_obs(ret.dir, "obs_meas.tab", &ctl, &obs_meas);
 
    /* Run retrieval... */
    double chisq;
    optimal_estimation(&ret, &ctl, tbl, &obs_meas, &obs_i, &atm_apr, &atm_i,
                       &chisq);
 
    /* Measure CPU-time... */
    TIMER("total", 2);
  }
 
  /* Write info... */
  LOG(1, "\nRetrieval done...");
 
  /* Measure CPU-time... */
  TIMER("total", 3);
 
  /* Free... */
  tbl_free(&ctl, tbl);
 
#ifdef MPI
  MPI_Finalize();
#endif
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: