|
- 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
|