Hello,
Normally this set-up provides a black surface but for this input-file the Facecolorsentences seems inactive since I get the same plot if I remove the two sentences. I like the way of the plot but the coloors are not right. A drayscale would be best
load thermoeqall
x3=thermoeqall(:,1)
y3=thermoeqall(:,2)
z3=thermoeqall(:,3)
xv = linspace(min(x3), max(x3), 1000);
yv = linspace(min(y3), max(y3), 1000);
[X,Y] = meshgrid(xv, yv);
Z = griddata(x3, y3, z3, X, Y);
A = mesh (X, Y, Z);
set(A(1),'FaceColor','k')
meshc(X,Y,Z)
grid on
xlabel('Temperature (K)','Fontsize',8)
ylabel('Pressure (Bars)','Fontsize',8)
zlabel('Hydrate Gibbs free energy (kJ/mole)','Fontsize',10)
set(gca,'fontsize',10)