def test_thermo_domain(): '''Can we call qsat, etc on a multi-dim climlab state temperature object?''' state = climlab.column_state(num_lev = 30, num_lat=3) T = state.Tatm p = T.domain.lev.points thermo.clausius_clapeyron(T) thermo.qsat(T, p) thermo.pseudoadiabat(T, p) thermo.blackbody_emission(T)
def test_thermo(): '''Basic single value tests for the thermodynamic routines.''' assert np.isclose(thermo.potential_temperature(250., 500.), 304.783) assert np.isclose(thermo.theta(250., 500.), 304.783) assert np.isclose(thermo.temperature_from_potential(300., 500.), 246.076) assert np.isclose(thermo.T(300., 500.), 246.076) assert np.isclose(thermo.clausius_clapeyron(300.), 35.345) assert np.isclose(thermo.qsat(300., 1000.), 0.02227839) assert np.isclose(thermo.estimated_inversion_strength(300., 290.), 5.3605345) assert np.isclose(thermo.EIS(300., 290.), 5.3605345) assert np.isclose(thermo.blackbody_emission(300.), 459.3)
def test_thermo(): '''Basic single value tests for the thermodynamic routines.''' assert np.isclose(thermo.potential_temperature(250., 500.), 304.783) assert np.isclose(thermo.theta(250., 500.), 304.783) assert np.isclose(thermo.temperature_from_potential(300., 500.), 246.076) assert np.isclose(thermo.T(300., 500.), 246.076) assert np.isclose(thermo.clausius_clapeyron(300.), 35.345) assert np.isclose(thermo.qsat(300., 1000.), 0.02227839) assert np.isclose(thermo.estimated_inversion_strength(300., 290.), 2.58025) assert np.isclose(thermo.EIS(300., 290.), 2.58025) assert np.isclose(thermo.blackbody_emission(300.), 459.3)
def compute_emission(self): return self.emissivity * blackbody_emission(self.Tatm)
def _compute_emission(self): return self.emissivity * blackbody_emission(self.Tatm)
def _compute_emission_sfc(self): return self.emissivity_sfc * blackbody_emission(self.Ts)