I believe data conversion is my main problem. I just can't seem to find a way to force Kixtart to treat the gap variable as a double precision number throughout my routine. It seems that during function calls, kixtart does its own internal data format conversions & this conversion is not fixed but rather based upon the actual values in the arrays.
