Volume:

• rass.f - the same problem as for surface case
• rags.f - if the value of rags4 is equal to 0 the inverted value causes division by zero, because output value rarags is similar to cbrags (k), in such case we set rarags equal to cbrags (k)