READ_MESH tet1.msh DIM 3
FUNCTION f(x,y,z) MESH tet1.msh

vec_f[i] = 1*vec_f_x[i] + 2*vec_f_y[i] + 3*vec_f_z[i]
# PRINT_FUNCTION f
# PRINT_VECTOR vec_f_x vec_f_y vec_f_z vec_f

PRINT f(0.1,0.1,0.1)