The calculation in routine time is straightforward. First, the boundary conditions are checked. Then the integration routine is called and the result is passed back to the IDL function calc_t.