brightness.c File Reference

Convert radiance to brightness temperature. More...

#include "jurassic.h"

Go to the source code of this file.

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

Detailed Description

Convert radiance to brightness temperature.

Definition in file brightness.c.

Function Documentation

main()

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

Definition at line 39 of file brightness.c.

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