grid anduct.B.cgd
 solution anduct.B.cfl
 units inches
 zone 1
 subset j all k all i last
 variables p0 p0inf; M
 integrate area output cf.B.lis
 quit