This script takes as inputs several horizons and a seg-y dataset. The objective is to manipulate values of the seg-y dataset to create a new attribute volume. Calculations can include seismic values, time, header information, and horizon attributes.

This is the result: