% figure with c_V and c_P
h1=subplot(2,1,1);
plot(mhd5(:,2),mhd5(:,13),'b');hold on;plot(mhd5(:,2),opal(:,13),'r');plot(mhd5(:,2),isg(:,13),'g');
h2=subplot(2,1,2);
plot(mhd5(:,2),mhd5(:,14),'b');hold on;plot(mhd5(:,2),opal(:,14),'r');plot(mhd5(:,2),isg(:,14),'g');
