def test_atmosphere_yearly_cycle(self):
"Model 'yearly_cycle'"
model = PISM.AtmosphereCosineYearlyCycle(self.grid)
model.init(self.geometry)
one_year = 365 * 86400.0
model.update(self.geometry, 0, one_year)
summer_peak = config.get_number(
"atmosphere.fausto_air_temp.summer_peak_day") / 365.0
ts = np.linspace(0, one_year, 13)
cycle = np.cos(2.0 * np.pi * (ts / one_year - summer_peak))
T = (self.T_summer - self.T_mean) * cycle + self.T_mean
P = np.zeros_like(T) + self.P
check_model(model, T=self.T_mean, P=self.P, ts=ts, Ts=T, Ps=P)
