planck.c File Reference

Convert brightness temperature to radiance. More...

#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

Convert brightness temperature to radiance.

Definition in file planck.c.

Function Documentation

main()

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

Definition at line 39 of file planck.c.

                {
 
  /* Print usage information... */
  USAGE;
 
  /* Check arguments... */
  if (argc != 3 && argc != 7)
    ERRMSG("Missing or invalid command-line arguments.\n\n"
           "Usage: planck <t> <nu>\n"
           "   or: planck <t0> <t1> <dt> <nu0> <nu1> <dnu>\n\n"
           "Use -h for full help.");
 
  /* Calculate single value... */
  if (argc == 3) {
 
    /* Read arguments... */
    const double t = atof(argv[1]);
    const double nu = atof(argv[2]);
 
    /* Compute Planck function... */
    printf("%.10g\n", PLANCK(t, nu));
  }
 
  /* Calculate table... */
  else if (argc == 7) {
 
    /* Read arguments... */
    const double t0 = atof(argv[1]);
    const double t1 = atof(argv[2]);
    const double dt = atof(argv[3]);
    const double nu0 = atof(argv[4]);
    const double nu1 = atof(argv[5]);
    const double dnu = atof(argv[6]);
 
    /* Write header... */
    printf("# $1 = brightness temperature [K]\n"
           "# $2 = wavenumber [cm^-1]\n"
           "# $3 = radiance [W/(m^2 sr cm^-1)]\n");
 
    /* Compute Planck function... */
    for (double t = t0; t <= t1; t += dt) {
      printf("\n");
      for (double nu = nu0; nu <= nu1; nu += dnu)
        printf("%.10g %.4f %.10g\n", t, nu, PLANCK(t, nu));
    }
  }
 
  return EXIT_SUCCESS;
}