Read the input parameters Raise an error if the number of input time series is less than 2 From the RATE extension header of the first time series Get the number of rows of the RATE extension (nRates) Get some mandatory keywords (e.g. TSTART, TSTOP, TIMEDEL, TIMEZERO) Do for input time series 2 to NTimeSeries Check keyword values against those of the first input time series (nRates, TSTART, TSTOP, TIMEDEL, TIMEZERO) If not the same then raise an error and exit End Do Create the output dataset Propagate the keywords from the first time series to the output time series for both the primary and RATE extension headers. Do for each input time series (dataset) Read the RATE extension from file i Read the values of E_MIN, E_MAX from the header Store these values as E_MINi, E_MAXi in the header of the RATE extension of the output file Get columns RATE, ERROR, BACKV, BACKE, FRACEXP Store corresponding arrays in matrix(NtimeSeries,nRates) End Do Write matrix in the RATE extension of the output dataset Release output dataset