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