gsiadd function

Internal functions: Parallel operations of single and multiple datasets