climatology.c File Reference

Prepare atmospheric data file from climatological data. More...

#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

Prepare atmospheric data file from climatological data.

Definition in file climatology.c.

Function Documentation

main()

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

Definition at line 27 of file climatology.c.

                {
 
  static atm_t atm;
  static ctl_t ctl;
 
  double clk[NCL], sfeps[NSF];
 
  /* Check arguments... */
  if (argc < 3)
    ERRMSG("Give parameters: <ctl> <atm>");
 
  /* Read control parameters... */
  read_ctl(argc, argv, &ctl);
  double t0 = scan_ctl(argc, argv, "T0", -1, "0", NULL);
  double t1 = scan_ctl(argc, argv, "T1", -1, "0", NULL);
  double dt = scan_ctl(argc, argv, "DT", -1, "1", NULL);
  double z0 = scan_ctl(argc, argv, "Z0", -1, "0", NULL);
  double z1 = scan_ctl(argc, argv, "Z1", -1, "90", NULL);
  double dz = scan_ctl(argc, argv, "DZ", -1, "1", NULL);
  double clz = scan_ctl(argc, argv, "CLZ", -1, "0", NULL);
  double cldz = scan_ctl(argc, argv, "CLDZ", -1, "0", NULL);
  for (int icl = 0; icl < ctl.ncl; icl++)
    clk[icl] = scan_ctl(argc, argv, "CLK", icl, "0", NULL);
  double sfz = scan_ctl(argc, argv, "SFZ", -1, "0", NULL);
  double sfp = scan_ctl(argc, argv, "SFP", -1, "0", NULL);
  double sft = scan_ctl(argc, argv, "SFT", -1, "0", NULL);
  for (int isf = 0; isf < ctl.nsf; isf++)
    sfeps[isf] = scan_ctl(argc, argv, "SFEPS", isf, "1", NULL);
 
  /* Set atmospheric grid... */
  for (double t = t0; t <= t1; t += dt)
    for (double z = z0; z <= z1; z += dz) {
      atm.time[atm.np] = t;
      atm.z[atm.np] = z;
      if ((++atm.np) >= NP)
        ERRMSG("Too many atmospheric grid points!");
    }
 
  /* Interpolate climatological data... */
  climatology(&ctl, &atm);
 
  /* Set cloud layer... */
  atm.clz = clz;
  atm.cldz = cldz;
  for (int icl = 0; icl < ctl.ncl; icl++)
    atm.clk[icl] = clk[icl];
 
  /* Set surface layer... */
  atm.sfz = sfz;
  atm.sfp = sfp;
  atm.sft = sft;
  for (int isf = 0; isf < ctl.nsf; isf++)
    atm.sfeps[isf] = sfeps[isf];
 
  /* Write data to disk... */
  write_atm(NULL, argv[2], &ctl, &atm);
 
  return EXIT_SUCCESS;
}

Here is the call graph for this function: