dolandscape
[ub,x,y]=testhisread('test.ubot.dat');       
pslice(x,y,ub,[0 0.35],'Max. orbital vel (m/s)');shading('flat');...
title('mb10.ubot.dat:  North wind, 20 m/s')