output.c

这是一个用SIMPLE法模拟的矩形流场的程序

#include "global.h"
void output()
{
	fprintf(filehandle1,"ZONE T=t%f\n I=%d, J=%d,F=POINT\n",time*dt,l,m);
	 for(j=1;j<=m;j++)
		{
			for(i=1;i<=l;i++)
			{
			if(i==1)fprintf(filehandle1,"%f,%f,%f,%f,%f,%f\n",xp[i],yp[j],f[i+1][j][1],(f[i][j][2]+f[i][j+1][2])/2,p[i][j],f[i][j][4]);
			else if(i==l)fprintf(filehandle1,"%f,%f,%f,%f,%f,%f\n",xp[i],yp[j],f[i][j][1],(f[i][j][2]+f[i][j+1][2])/2,p[i][j],f[i][j][4]);
					else if(j==1)fprintf(filehandle1,"%f,%f,%f,%f,%f,%f\n",xp[i],yp[j],(f[i][j][1]+f[i+1][j][1])/2,f[i][j+1][2],p[i][j],f[i][j][4]);
						else if(j==m)fprintf(filehandle1,"%f,%f,%f,%f,%f,%f\n",xp[i],yp[j],(f[i][j][1]+f[i+1][j][1])/2,f[i][j][2],p[i][j],f[i][j][4]);
							else fprintf(filehandle1,"%f,%f,%f,%f,%f,%f\n",xp[i],yp[j],(f[i][j][1]+f[i+1][j][1])/2,(f[i][j][2]+f[i][j+1][2])/2,p[i][j],f[i][j][4]);
							 
			}
		}

}