Add check if contfrac=0 to avoid division by zero. contfrac should not be 0 but still it might be the case in some special cases. Set aoutvar=-2 which is a diagnostic variable to evaluate the interpolation. -2 shows that something special happend on this grid cell.
(No files)
Note: See TracChangeset
for help on using the changeset viewer.