obs2spec.c File Reference

Converter for spectra. More...

#include <omp.h>
#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

Converter for spectra.

Definition in file obs2spec.c.

Function Documentation

main()

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

Definition at line 32 of file obs2spec.c.

                {
 
  static ctl_t ctl;
 
  obs_t *obs;
 
  FILE *out;
 
  /* Check arguments... */
  if (argc < 4)
    ERRMSG("Give parameters: <ctl> <obs> <spec.tab>");
 
  /* Allocate... */
  ALLOC(obs, obs_t, 1);
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
 
  /* Read observation geometry... */
  read_obs(".", argv[2], &ctl, obs);
 
  /* Write info... */
  LOG(1, "Write spectra: %s", argv[3]);
 
  /* Create file... */
  if (!(out = fopen(argv[3], "w")))
    ERRMSG("Cannot create file!");
 
  /* Write header... */
  fprintf(out,
          "# $1 = time (seconds since 2000-01-01T00:00Z)\n"
          "# $2 = observer altitude [km]\n"
          "# $3 = observer longitude [deg]\n"
          "# $4 = observer latitude [deg]\n"
          "# $5 = view point altitude [km]\n"
          "# $6 = view point longitude [deg]\n"
          "# $7 = view point latitude [deg]\n"
          "# $8 = tangent point altitude [km]\n"
          "# $9 = tangent point longitude [deg]\n"
          "# $10 = tangent point latitude [deg]\n"
          "# $11 = channel frequency [cm^-1]\n"
          "# $12 = channel radiance [W/(m^2 sr cm^-1)]\n"
          "# $13 = channel transmittance [1]\n");
 
  /* Write data... */
  for (int ir = 0; ir < obs->nr; ir++) {
    fprintf(out, "\n");
    for (int id = 0; id < ctl.nd; id++)
      fprintf(out, "%.2f %g %g %g %g %g %g %g %g %g %.4f %g %g\n",
              obs->time[ir], obs->obsz[ir], obs->obslon[ir], obs->obslat[ir],
              obs->vpz[ir], obs->vplon[ir], obs->vplat[ir], obs->tpz[ir],
              obs->tplon[ir], obs->tplat[ir], ctl.nu[id], obs->rad[id][ir],
              obs->tau[id][ir]);
  }
 
  /* Close file... */
  fclose(out);
 
  /* Free... */
  free(obs);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: