shear = 1:.5:40;
conc = 0.05:.005:.65;
shear = shear';
conc = conc';

con_size_sub = size(conc);
shear_size_sub =size(shear);
con_size = con_size_sub(1);
shear_size = shear_size_sub(1);


vis = zeros(shear_size,con_size);
con_mat = zeros(shear_size,con_size);
shear_mat = zeros(shear_size,con_size);

for x = 1:shear_size
    for y = 1:con_size
        vis(x,y) = (400.7 *conc(y) + 1.0) * shear(x)^-0.5;
        con_mat(x,y) = conc(y);
        shear_mat(x,y) = shear(x);
    end
end
%aviobj = avifile('test5.avi');
frame_num = 1;
for ang = 115:2:475

    mesh(con_mat,shear_mat,vis, 'FaceColor', [0.804 0.878 0.969], 'LineStyle', 'none','Marker','none');
    hold on ;
    title('A Comparison of Modeled and Measured Effective Viscosity for Ash-poor Slurries');
    xlabel('Concentration (Sed:Water)');
    ylabel('Shear Rate (1/sec)');
    zlabel('Effective Viscosity (cP)');
    xlim([0 0.60]);
    set(GCA, 'CameraPosition',[5.33 192.7 620]);
    set(GCA, 'CameraPositionMode', 'manual')
    set(GCA, 'CameraTarget' , [0.382 19.5 136.1]);
    set(gca, 'CameraTargetMode', 'manual');
    set(gca, 'CameraViewAngle', 10);
    set(gca, 'CameraViewAngleMode', 'manual');
    view([ang,25]);
    plot3(ash_sand(:,2),ash_sand(:,1),ash_sand(:,3),'o', 'MarkerEdgeColor',[0.478 0.063 0.894], 'MarkerFaceColor',[0.478 0.063 0.894]);
    plot3(ash_silt(:,2),ash_silt(:,1),ash_silt(:,3),'o','MarkerEdgeColor',[0.0 0.749 0.749], 'MarkerFaceColor',[0.0 0.749 0.749]);
    plot3(ash_silt_sand(:,2),ash_silt_sand(:,1),ash_silt_sand(:,3),'o','MarkerEdgeColor',[0.871 0.49 0.0], 'MarkerFaceColor',[0.871 0.49 0.0]);
    plot3(sand(:,2),sand(:,1),sand(:,3),'o','MarkerEdgeColor',[0.847 0.161 0.0], 'MarkerFaceColor',[0.847 0.161 0.0]);
    plot3(silt(:,2),silt(:,1),silt(:,3),'o','MarkerEdgeColor',[0.23 0.44 0.33], 'MarkerFaceColor',[0.23 0.44 0.33]);
    plot3(silt_sand(:,2),silt_sand(:,1),silt_sand(:,3),'o','MarkerEdgeColor',[1.0 0.6 0.78], 'MarkerFaceColor',[1.0 0.6 0.78]);
    hold off
    %frame =getframe;
    %aviobj = addframe(aviobj,frame);
    M(frame_num) =getframe;
    frame_num = frame_num +1;
end