pure function surfaceTensionAux input Real p; input Real T; input Real x[:]; output SurfaceTension sigma "surface tension"; end surfaceTensionAux;