radiation.c

Optically Thin Radiation Model

We test the calculation of the radiation using the optically-thin model at different temperautres.

We define variables that are necessary for the radiation model to compile.

#define NGS 2
const char* gas_species[2] = {"CO2", "H2O"};

#include "common.h"
#include "utils.h"
#include "radiation.h"

int main (void) {

The properties for the radiation model are gathered in the following struct.

  OpticallyThinProperties otp;
  otp.T = 1400.;
  otp.P = 101325.;
  otp.xH2O = 0.21;
  otp.xCO2 = 0.1;

We set the radiation model, which is set to no_radiation by default.

  divq_rad = optically_thin;

We loop over different temperatures and compute the \nabla\cdot\mathbf{q}_{rad} contribution.

  for (otp.T = 300.; otp.T <= 3000.; otp.T += 50.) {
    fprintf (stderr, "divq_rad = %f\n", divq_rad(&otp));
  }
}