(custom-field-function/define '(((name q_crit) (display " - 0.5 * (dx-velocity-dx ^ 2 + dy-velocity-dy ^ 2 + dz-velocity-dz ^ 2) - dx-velocity-dy * dy-velocity-dx - dx-velocity-dz * dz-velocity-dx - dy-velocity-dz * dz-velocity-dy") (syntax-tree ("-" ("-" ("-" ("*" -0.5 ("+" ("+" ("**" "dx-velocity-dx" 2) ("**" "dy-velocity-dy" 2)) ("**" "dz-velocity-dz" 2))) ("*" "dx-velocity-dy" "dy-velocity-dx")) ("*" "dx-velocity-dz" "dz-velocity-dx")) ("*" "dy-velocity-dz" "dz-velocity-dy"))) (code (field-- (field-- (field-- (field-* -0.5 (field-+ (field-+ (field-** (field-load "dx-velocity-dx") 2) (field-** (field-load "dy-velocity-dy") 2)) (field-** (field-load "dz-velocity-dz") 2))) (field-* (field-load "dx-velocity-dy") (field-load "dy-velocity-dx"))) (field-* (field-load "dx-velocity-dz") (field-load "dz-velocity-dx"))) (field-* (field-load "dy-velocity-dz") (field-load "dz-velocity-dy"))))) ((name x-stretch) (display "x-vorticity * dx-velocity-dx + y-vorticity * dx-velocity-dy + z-vorticity * dx-velocity-dz") (syntax-tree ("+" ("+" ("*" "x-vorticity" "dx-velocity-dx") ("*" "y-vorticity" "dx-velocity-dy")) ("*" "z-vorticity" "dx-velocity-dz"))) (code (field-+ (field-+ (field-* (field-load "x-vorticity") (field-load "dx-velocity-dx")) (field-* (field-load "y-vorticity") (field-load "dx-velocity-dy"))) (field-* (field-load "z-vorticity") (field-load "dx-velocity-dz"))))) ((name y-stretch) (display "x-vorticity * dy-velocity-dx + y-vorticity * dy-velocity-dy + z-vorticity * dy-velocity-dz") (syntax-tree ("+" ("+" ("*" "x-vorticity" "dy-velocity-dx") ("*" "y-vorticity" "dy-velocity-dy")) ("*" "z-vorticity" "dy-velocity-dz"))) (code (field-+ (field-+ (field-* (field-load "x-vorticity") (field-load "dy-velocity-dx")) (field-* (field-load "y-vorticity") (field-load "dy-velocity-dy"))) (field-* (field-load "z-vorticity") (field-load "dy-velocity-dz"))))) ((name z-stretch) (display "x-vorticity * dz-velocity-dx + y-vorticity * dz-velocity-dy + z-vorticity * dz-velocity-dz") (syntax-tree ("+" ("+" ("*" "x-vorticity" "dz-velocity-dx") ("*" "y-vorticity" "dz-velocity-dy")) ("*" "z-vorticity" "dz-velocity-dz"))) (code (field-+ (field-+ (field-* (field-load "x-vorticity") (field-load "dz-velocity-dx")) (field-* (field-load "y-vorticity") (field-load "dz-velocity-dy"))) (field-* (field-load "z-vorticity") (field-load "dz-velocity-dz"))))) ))