public void Run(Clock clock, IClimateCH4CycleState state, IDimensions dimensions) { // create shortcuts for commonly accessed data var s = state; var t = clock.Current; if (clock.IsFirstTimestep) { s.ch4decay = 1.0 / s.lifech4; s.acch4[t] = 1222.0; } else { // Calculate CH4 concentrations s.acch4[t] = s.acch4[t - 1] + 0.3597 * s.globch4[t] - s.ch4decay * (s.acch4[t - 1] - s.ch4pre); if (s.acch4[t] < 0) throw new ApplicationException("ch4 atmospheric concentration out of range"); } }
public void Run(Clock clock, IClimateCH4CycleState state, IDimensions dimensions) { // create shortcuts for commonly accessed data var s = state; var t = clock.Current; if (clock.IsFirstTimestep) { s.ch4decay = 1.0 / s.lifech4; s.acch4[t] = 1222.0; } else { // Calculate CH4 concentrations s.acch4[t] = s.acch4[t - 1] + 0.3597 * s.globch4[t] - s.ch4decay * (s.acch4[t - 1] - s.ch4pre); if (s.acch4[t] < 0) { throw new ApplicationException("ch4 atmospheric concentration out of range"); } } }