1
grid run.cgd
solution run.cfl
units fss
subset i all j all k all
zone 1
zone 2
zone 3
zone 4
cut at z 0.2  
variables M
genplot surface output mach.m overwrite
contour increment 0.05
plot color contours edge data mach.m
quit