Code:
/*Define plot Function*/
void Plot(int nx, int ny, double x[350], double y[350], double double T_new[350][350], int iplot){
char fname[50];
char title[100];
char var[100];
char zone1[100];
char zone2[100];
char blank[5];
sprintf(title,"TITLE= Unsteady Advection Upwind Without Diffusion_%d", iplot);
sprintf(var, "VARIABLES=\"X\" \"Y\" \"C\"");
sprintf(zone1, "ZONE I=%d, ", nx);
sprintf(zone2, " J=%d, ", ny);
sprintf(blank, "\n");
FILE *fp;
sprintf(fname,"Unsteady Advection Upwind Without Diffusion_%d.PLT",iplot);
fp=fopen(fname, "wb");
fwrite(title, sizeof(title[0]),sizeof(title)/sizeof(title[0]),fp);
fwrite(blank,sizeof(blank[0]),sizeof(blank)/sizeof(blank[0]),fp);
fwrite(var, sizeof(var[0]),sizeof(var)/sizeof(var[0]),fp);
fwrite(blank,sizeof(blank[0]),sizeof(blank)/sizeof(blank[0]),fp);
fwrite(zone1, sizeof(zone1[0]),sizeof(zone1)/sizeof(zone1[0]),fp);
fwrite(blank,sizeof(blank[0]),sizeof(blank)/sizeof(blank[0]),fp);
fwrite(zone2, sizeof(zone2[0]),sizeof(zone2)/sizeof(zone2[0]),fp);
fwrite(blank,sizeof(blank[0]),sizeof(blank)/sizeof(blank[0]),fp);
fwrite(x,sizeof(x[0]),sizeof(x)/sizeof(x[0]),fp);
fwrite(blank,sizeof(blank[0]),sizeof(blank)/sizeof(blank[0]),fp);
fwrite(y,sizeof(y[0]),sizeof(y)/sizeof(y[0]),fp);
fwrite(blank,sizeof(blank[0]),sizeof(blank)/sizeof(blank[0]),fp);
fwrite(T_new,sizeof(T_new[0][0]),sizeof(T_new)/sizeof(T_new[0][0]),fp);
fclose(fp);
}