Numa3DFunctions.c File Reference

#include <math.h>
#include <mathfunctions.h>
#include <physicalmodels/numa3d.h>

Go to the source code of this file.

Functions
int	Numa3DCalculateStandardAtmosphere_1 (void *p, double z, double *ExnerP, double *P, double *rho, double *T)
int	Numa3DCalculateStandardAtmosphere_2 (void *p, double z, double *ExnerP, double *P, double *rho, double *T)

Function Documentation

int Numa3DCalculateStandardAtmosphere_1	(	void *	p,
		double	z,
		double *	ExnerP,
		double *	P,
		double *	rho,
		double *	T
	)

Definition at line 5 of file Numa3DFunctions.c.

{
  Numa3D *physics = (Numa3D*) p;

  double R      = physics->R;
  double gamma  = physics->gamma;
  double g      = physics->g;

  /* reference quantities at zero altitude */
  double P0, T0; 
  P0   = physics->Pref;
  T0   = physics->Tref;

  double inv_gamma_m1 = 1.0/(gamma-1.0);
  double Cp = gamma * inv_gamma_m1 * R;

  double theta  = T0;
  *ExnerP = 1.0 - (g/(Cp*theta))*z;
  *P      = P0 * raiseto((*ExnerP),gamma*inv_gamma_m1);
  *rho    = (P0/(R*theta)) * raiseto((*ExnerP),inv_gamma_m1);
  *T      = (*rho) * theta;

  return(0);
}

int Numa3DCalculateStandardAtmosphere_2	(	void *	p,
		double	z,
		double *	ExnerP,
		double *	P,
		double *	rho,
		double *	T
	)

Definition at line 30 of file Numa3DFunctions.c.

{
  Numa3D *physics = (Numa3D*) p;

  double R      = physics->R;
  double gamma  = physics->gamma;
  double g      = physics->g;

  /* reference quantities at zero altitude */
  double P0, T0; 
  P0 = physics->Pref;
  T0 = physics->Tref;

  double BV = 0.01; /* Brunt-Vaisala frequency */
  double inv_gamma_m1 = 1.0/(gamma-1.0);
  double Cp = gamma * inv_gamma_m1 * R;

  double term   = BV*BV*z/g;
  double theta  = T0 * exp(term);
  *ExnerP = 1.0 + (g*g/(Cp*T0*BV*BV)) * (exp(-term) - 1.0);
  *P      = P0 * raiseto((*ExnerP),gamma*inv_gamma_m1);
  *rho    = (P0/(R*theta)) * raiseto((*ExnerP),inv_gamma_m1);
  *T      = (*rho) * theta;

  return(0);
}