CRadtran/brightness_temp.cpp at main · rtlines/CRadtran · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
/*****************************************************************************
 * Source Code:
 *    brightness_temp.cpp
 * Author:
 *    Alexander Marvin
 * Summary:
 *    Calculates blackbody radiation, using Planck's equation
 ****************************************************************************/

#include <iostream>
#include <cstdlib>
#include <cmath>
#include "brightness_temp.h"

using namespace std;

long double planckv(double freq, double temp)
{
  long double B;
  long double c = SPEED_OF_LIGHT;
  long double h = PLANCK_CONST*pow(10,-34);
  long double kB = BOLTZMAN_CONST*pow(10,-23);
  B = 2*h*pow(freq, 3)/pow(c, 2)*1/(exp(h*freq/(kB*temp)) - 1);
  return B;
};

long double planckw(double wave, double temp)
{
  long double B, nu;
  long double c = SPEED_OF_LIGHT;
  long double h = PLANCK_CONST*pow(10,-34);
  long double kB = BOLTZMAN_CONST*pow(10,-23);
  nu = c/wave;
  //B = 2*h*pow(c,2)/pow(wave,5)*1/(exp(h*c/(wave*kB*temp)) - 1);
  B = 2*h*pow(nu,3)/( pow(c,2)*(exp(h*nu/(kB*temp))) - 1);
  return B;
};