>> 
>> 
>> clear
>> load GOMtempNov98
>> whos
  Name        Size                   Bytes  Class

  C        8563x1                    68504  double array
  Err      8563x1                    68504  double array
  I        1213x1                     9704  double array
  J        2194x1                    17552  double array
  T        8563x1                    68504  double array
  Val      8563x1                    68504  double array
  amap        1x64                     512  double array
  ans         1x1                        8  double array
  blue      256x3                     6144  double array
  ele     16749x4                   535968  double array
  h           1x1                        8  double array
  map       513x3                    12312  double array
  nod      8563x3                   205512  double array
  nodll    8563x3                   205512  double array
  p           1x4                       32  double array

Grand total is 158410 elements using 1267280 bytes

>> 
>> 
>> h=polygonsurf(ele(:,2:4),nodll(:,2:3),Val);
>> shading interp
>> C=(Val-20)/(100-20)+1;I=find(Val<20);C(I)=nan;
>> clf
>> h=polygonsurf(ele(:,2:4),nodll(:,2:3),C);
>> colorbar
>> shading interp
>> colormap(blue)
>> colorbar
>> clf
>> g=GOMplot3D_2;
>> hold
Current plot held
>> hold on
>> h=polygonsurf(ele(:,2:4),nodll(:,2:3),C);
>> colormap([map;blue]);caxis([-1 2]);
>> shading interp
>> set(h,'backfacelighting','lit')
>> set(h,'zdata',get(h,'xdata')*0-100);
>> get(h)
	AlphaDataMapping = scaled
	CData = [ (3 by 16749) double array]
	CDataMapping = scaled
	FaceVertexAlphaData = []
	FaceVertexCData = [ (50247 by 1) double array]
	EdgeAlpha = [1]
	EdgeColor = none
	EraseMode = normal
	FaceAlpha = [1]
	FaceColor = interp
	Faces = [ (16749 by 3) double array]
	LineStyle = -
	LineWidth = [0.5]
	Marker = none
	MarkerEdgeColor = auto
	MarkerFaceColor = none
	MarkerSize = [6]
	Vertices = [ (50247 by 3) double array]
	XData = [ (3 by 16749) double array]
	YData = [ (3 by 16749) double array]
	ZData = [ (3 by 16749) double array]
	FaceLighting = flat
	EdgeLighting = none
	BackFaceLighting = lit
	AmbientStrength = [0.3]
	DiffuseStrength = [0.6]
	SpecularStrength = [0.9]
	SpecularExponent = [10]
	SpecularColorReflectance = [1]
	VertexNormals = [ (50247 by 3) double array]
	NormalMode = auto

	BeingDeleted = off
	ButtonDownFcn = 
	Children = []
	Clipping = on
	CreateFcn = 
	DeleteFcn = 
	BusyAction = queue
	HandleVisibility = on
	HitTest = on
	Interruptible = on
	Parent = [3.00146]
	Selected = off
	SelectionHighlight = on
	Tag = 
	Type = patch
	UIContextMenu = []
	UserData = []
	Visible = on

>> A=Err(ele(:,2:4));
>> whos A
  Name      Size                   Bytes  Class

  A     16749x3                   401976  double array

Grand total is 50247 elements using 401976 bytes

>> A=A';
>> A=A(:);
>> whos A
  Name      Size                   Bytes  Class

  A     50247x1                   401976  double array

Grand total is 50247 elements using 401976 bytes

>> set(h,'FaceVertexAlphaData',A,'FaceAlpha','interp');
>> set(gcf,'AlphaMap',amap);
>> set(gca,'alim',[0.2 1.4]);
>> exit