pro visc
openr,1,'visc_he.dat'
x=dblarr(19)
y=dblarr(9)
z=dblarr(9,19)
readf,1,x
readf,1,y
y=alog10(y)
readf,1,z
close,1

surface,z,y,x,ytitle='Temperature/K',xtitle='Log(P/bar)',$
	ztitle='Viscosity / 10!E-4!N Pa.s',charsize=1.8,zrange=[0,0.30]
end