def compute_eqmc_bar(self, daylengths) :
bufs = self.get_buffer_group()
derived = self.get_derived_buffers()
temps = bufs.get('Tair')
rh = derived.get('rh')
t_min = temps.min()
t_max = temps.max()
rh_min = rh.min()
rh_max = rh.max()
return fm.eqmc_bar(daylengths, t_max, t_min, rh_max, rh_min)
def compute_eqmc_bar(self, daylengths):
bufs = self.get_buffer_group()
derived = self.get_derived_buffers()
temps = bufs.get('Tair')
rh = derived.get('rh')
t_min = temps.min()
t_max = temps.max()
rh_min = rh.min()
rh_max = rh.max()
return fm.eqmc_bar(daylengths, t_max, t_min, rh_max, rh_min)
def test_eqmc_bar_units(self) :
"""checks that eqmc_bar returns expected values"""
rh_min = 15*u.pct
rh_max = 60*u.pct
t_min = 60*iu.deg_F
t_max = 85*iu.deg_F
daylength = 15*u.hour
testval = fm.eqmc_bar(daylength, t_max, t_min, rh_max, rh_min)
emax = fm.eqmc(t_min,rh_max)
emin = fm.eqmc(t_max,rh_min)
refval = (daylength*emin + (1*u.day-daylength)*emax)/(1*u.day)
self.assertEqual(testval,refval)
self.assertEqual(testval.unit, u.pct)
def test_eqmc_bar_units(self):
"""checks that eqmc_bar returns expected values"""
rh_min = 15 * u.pct
rh_max = 60 * u.pct
t_min = 60 * iu.deg_F
t_max = 85 * iu.deg_F
daylength = 15 * u.hour
testval = fm.eqmc_bar(daylength, t_max, t_min, rh_max, rh_min)
emax = fm.eqmc(t_min, rh_max)
emin = fm.eqmc(t_max, rh_min)
refval = (daylength * emin +
(1 * u.day - daylength) * emax) / (1 * u.day)
self.assertEqual(testval, refval)
self.assertEqual(testval.unit, u.pct)
