cal_dlnT_dlnP.cpp

Go to the documentation of this file.

// canoe
#include <air_parcel.hpp>
 
// snap
#include "thermodynamics.hpp"
 
Real Thermodynamics::calDlnTDlnP(AirParcel const& qfrac, Real latent[]) const {
  // calculate gammad
  Real gammad = GetGammad(qfrac);
 
  Real q_gas = 1.;
  for (int n = 0; n < NCLOUD; ++n) q_gas -= qfrac.c[n];
 
  Real f_sig = 1.;
  // vapor
  for (int n = 1; n <= NVAPOR; ++n)
    f_sig += qfrac.w[n] * (cp_ratio_mole_[n] - 1.);
  // cloud
  for (int n = 0; n < NCLOUD; ++n)
    f_sig += qfrac.c[n] * (cp_ratio_mole_[1 + NVAPOR + n] - 1.);
  Real cphat_ov_r = gammad / (gammad - 1.) * f_sig / q_gas;
 
  // vapor
  Real xd = q_gas;
  for (int n = 1; n <= NVAPOR; ++n) xd -= qfrac.w[n];
 
  Real c1 = 0., c2 = 0., c3 = 0.;
  for (int iv = 1; iv <= NVAPOR; ++iv) {
    c1 += qfrac.w[iv] / xd * latent[iv];
    c2 += qfrac.w[iv] / xd * latent[iv] * latent[iv];
    c3 += qfrac.w[iv] / xd;
  }
 
  return (1. + c1) / (cphat_ov_r + (c2 + c1 * c1) / (1. + c3));
}