elephant-sparselizard-axi/diff.was

FUNCTION u_fino(y)  FILE_PATH axi-fino-100.dat

MESH FILE_PATH u.vtk  DIMENSIONS 2 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)

# 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