user, 'drytol', drytol_default) q = current_data.q h=ma.masked_where(q[:,:,0]<drytol,q[:,:,0]) m=solid_frac(current_data) rho = density(current_data) var = ma.masked_where(h<drytol,gmod(current_data)*rho*h) return var...