FUNCTION u_fino(y)  FILE_PATH 3d-fino-20.dat

MESH FILE_PATH u.vtk  DIMENSIONS 3 READ_SCALAR u1 AS u_lizard_vtk
# extract the profile at 95% of the radius
u_lizard(y) := u_lizard_vtk(0.95*0.5, y, 0)

# write the two profiles and their differences
PRINT_FUNCTION FORMAT %e FILE_PATH diff-fino-lizard.dat \
  u_fino u_lizard u_lizard(y)-u_fino(y) (u_lizard(y)-u_fino(y))/u_lizard(y) \
  MIN 2e-2 MAX 2 NSTEPS 100
