SPH无网格程序

SPH无网格程序，可以调试使用的工程文件，C++和Fortran语音混合编程！

/**************************************************************************************************/
void export_vart（int mint ntdouble t）
{
int inreglength;
char filename[10]filename0[10];
FILE *fp2;
double c0u0p0d0e0;

itoa（mfilename010）;//把整数转换为字符串。
length=strlen（filename0）;//求字符串长度。
if（length==1） strcpy（filename“0000“）;//将一个字符串复制到另一字符串中。
if（length==2） strcpy（filename“000“）;
if（length==3） strcpy（filename“00“）;
if（length==4） strcpy（filename“0“）;
if（length==5） strcpy（filename““）;
strcat（filenamefilename0）;//将一个字符串附加到另一个字符串上。
strcat（filename“.plt“）;
fp2=fopen（filename“w“）;//打开流!
fprintf（fp2“%s\n““VARIABLES = X U P D E H S A U0 P0 D0 E0“）;
fprintf（fp2“%s\n““TEXTCS=frameC=REDS=LOCALX=60.0Y=2.0“）;
fprintf（fp2“%s\n““HU=POINTLS=1AN=LEFTBX=FILLEDBXM=20LT=0.4“）;
fprintf（fp2“%s\n““BXO=WHITEBXF=WHITEF=TIMES-BOLDH=20A=0MFC=\“\““）;
fprintf（fp2“%s%12.5le%s\n““T=\“Time = “t/1.0e-6“ ms\““）;

double u1=0.0;//Input Velocity1!
double u2=0.0;//Input Velocity2!
double c1=1290.0;
double c2=6519.97863934307;
double p1=1.0e6;//Input Pressure1!
double p2=3.5e10;//Input Pressure2!
double p3=2.570826e10;//2.461595124e10;//Input Pressure3!
double d1=15750.0;//15.75;//Input Density1!
double d2=2470.0;//Input Density2!
double e1=15.8730158730159;
double e2=7.08502024291498e6;
double k1=5.0k2=3.0;
double a2=2.32261477582487*1.0e9/pow（1000.0k2）;//Input a2!

double d3=d1*（p1*（k1-1.0）+p3*（k1+1.0）+2.0*c1*c1*d1）/（p3*（k1-1.0）+p1*（k1+1.0）+2.0*c1*c1*d1）;
double d4=pow（p3/a21.0/k2）;
double u3=-sqrt（（p3-p1）*（1.0/d1-1.0/d3））;
double e3=（p3-c1*c1*（d3-d1））/（k1-1.0）/d3;
double e4=p3/d4/（k2-1.0）;
double D=d3*u3/（d3-d1）;
double XJ=u3*t;
double XD=D*t;
double XRH=c2*t;
double XRT=（c2+（k2+1.0）/2.0*u3-（k2-1.0）/2.0*u2）*t;

for（nreg=1;nreg<=mreg;nreg++）
{
if（openreg[nreg]==0） continue;
fprintf（fp2“%s%s%s%3d%s\n““ZONE T=“namereg[nreg]“ I=“iip[nreg]“ F=POINT“）;
for（i=1;i<=iip[nreg];i++）
{
if（rp[nm[i][nreg]][nt].x<=XD && nreg==1）
{
	u0=u1;
	p0=p1;
	d0=d1;
	e0=e1;
}
if（rp[nm[i][nreg]][nt].x>XD && nreg==1）
{
	u0=u3;
	p0=p3;
	d0=d3;
	e0=e3;
}
if（rp[nm[i][nreg]][nt].x<=XRT && nreg==2）
{
	u0=u3;
	p0=p3;
	d0=d4;
	e0=e4;
}
if（rp[nm[i][nreg]][nt].x>XRT && rp[nm[i][nreg]][nt].x{
	c0=（k2-1.0）/（k2+1.0）*（rp[nm[i][nreg]][nt].x/t-u2+2.0/（k2-1.0）*c2）;
	u0=（k2-1.0）/（k2+1.0）*（2.0/（k2-1.0）*rp[nm[i][nreg]][nt].x/t+u2-2.0/（k2-1.0）*c2）;
	p0=p2*pow（c0/c22.0*k2/（k2-1.0））;
	d0=d2*pow（p0/p21.0/k2）;
	e0=p0/d0/（k2-1.0）;
}
if（rp[nm[i][nreg]][nt].x>=XRH && nreg==2）
{
	u0=u2;
	p0=p2;
	d0=d2;
	e0=e2;
}

fprintf（fp2“%15.5le%15.5le%15.5le%15.5le%15.5le%15.5le%5d%5d%15.5le%15.5le%15.5le%15.5le\n“
	    rp[nm[i][nreg]][nt].xvp[nm[i][nreg]][nt].x
        pp[nm[i][nreg]][nt]  dp[nm[i][nreg]][nt]  
		ep[nm[i][nreg]][nt]  hp[nm[i][nreg]][nt]  
		sp[nm[i][nreg]]      ap[nm[i][nreg]]      u0p0d0e0）;
}
}
fclose（fp2）;//关闭流!
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2454  2008-12-03 15:31  SPH_YJW\artificial.for

     文件        592  2008-12-03 15:31  SPH_YJW\density.for

     文件        341  2008-12-03 15:31  SPH_YJW\eos.for

     文件       1239  2008-12-03 14:36  SPH_YJW\exact.for

     文件       1184  2009-03-05 19:21  SPH_YJW\hupdate.for

     文件       5890  2008-12-07 10:31  SPH_YJW\intergral.for

     文件       1514  2008-06-15 09:43  SPH_YJW\internal.for

     文件       1264  2008-02-28 17:29  SPH_YJW\kernel.for

     文件       3051  2009-10-14 15:44  SPH_YJW\M4_exact.cpp

     文件       3011  2009-10-14 15:44  SPH_YJW\M5_exact.cpp

     文件       2745  2008-12-07 08:59  SPH_YJW\model.for

     文件       1192  2009-10-15 14:10  SPH_YJW\neighbor.for

     文件       2141  2008-12-03 15:59  SPH_YJW\output.for

     文件        777  2009-10-14 08:46  SPH_YJW\param.inc

     文件       2454  2008-12-03 15:31  SPH_YJW\R4\artificial.for

     文件       3432  2009-10-15 13:11  SPH_YJW\R4\Debug\BuildLog.htm

     文件       2048  2009-10-15 13:11  SPH_YJW\R4\Debug\SPH.exe.embed.manifest

     文件        592  2008-12-03 15:31  SPH_YJW\R4\density.for

     文件        489  2009-10-14 19:27  SPH_YJW\R4\eos.for

     文件       2184  2009-10-15 14:03  SPH_YJW\R4\exact.for

     文件       1184  2009-03-05 19:21  SPH_YJW\R4\hupdate.for

     文件       4315  2009-10-15 14:12  SPH_YJW\R4\intergral.for

     文件       1514  2008-06-15 09:43  SPH_YJW\R4\internal.for

     文件       1264  2008-02-28 17:29  SPH_YJW\R4\kernel.for

     文件       3051  2009-10-14 15:44  SPH_YJW\R4\M4_exact.cpp

     文件       3011  2009-10-14 15:44  SPH_YJW\R4\M5_exact.cpp

     文件       2511  2009-10-15 13:09  SPH_YJW\R4\model.for

     文件       1155  2009-10-15 14:07  SPH_YJW\R4\neighbor.for

     文件       5040  2009-10-14 22:31  SPH_YJW\R4\output.for

     文件        777  2009-10-14 08:46  SPH_YJW\R4\param.inc

............此处省略100个文件信息