clear;

load data/katrina_condensate;

x=squeeze(longitude(113,:,1));
y=squeeze(latitude(:,135,1))';
z=squeeze(height(113,135,:))';
%[X,Y,Z]=meshgrid(x,y,z);

[X,Y,Z]=meshgrid([1:270],[1:225],[1:59]);
hiso=patch(isosurface(X,Y,Z,condensate,2));

isonormals(x,y,z,condensate,hiso);
set(hiso,'EdgeColor','none')
%hcap = patch(isocaps(x,y,z,condensate,1),'FaceColor','interp','EdgeColor','none');

%axis tight
box on
%camproj perspective
%camzoom(1.25)
view(0,90)

%colormap(jet)
 
box on

camproj perspective;
 
%camva(34);
 
%campos([165 -20 65]);
 
%camtarget([100 40 -5])
 
%camlight left;
 
%lighting gouraud


-can


view(-37.5,30);
---can

%contourf(x,y,condensate(:,:,24),'edgecolor','none');
%h=colorbar;
%set(get(h,'title'),'string','g/kg');
%---can

v=shiftdim(condensate,1);

%p=patch(isosurface(X,Y,Z,rr,0.1));
%isonormals(X,Y,Z, p)
%set(p, 'FaceColor', [0.8 0.8 0.8], 'EdgeColor', 'none');


