火焰通过平板传播自编程序VC

1、等温-火花点火程序Setting.c#include global_var.hvoid default_setting()缺省设定一些全局变量。包括计算域的尺寸以及网格密度、需要解的变量方程组、采用的坐标系、迭代参次数、每次迭代扫描次数、松弛因子、时间步长等等。计 算中的一些策略设定，如是否块修正计算域也包括在内。int ij;dt=O.OOOOOl; 时间步长 0.1/12 妙 nt=5800;xl=0.1; 求解域x方向宽度=6.0米 yl=0.002; 求解域y方向半宽度=0.06米11=5001;ml=101;z=3;for(i=l;i = nt;i+)printi=O;/printl

2、=l;/print50=l;/print100=l;/print150=l;/print200=l;/print250=l;/print300=l;/print350=l;/print400=l;/print450=l;/print500=l;/print550=l;/print600=l;/print650=l;/print700=l;/print750=l;/print800=l;/print850=l;/print900=l;/print950=l;/print1000=l;/for(i=l;ill;i+)/1fclose(filehandleul600);fclose(filehand

3、levl600);fclose(filehandlepl600);fclose(filehandleT1600);fclose(filehandlerhol600);fclose(filehandleyfl600);fclose(filehandleomegattl600);fclose(filehandleWml600);fclose(filehandleEccl600);if(tt=1650)filehandleul650=fopen(e:/simpele/result/ul650.txt,1,w);filehandlevl650=fopen(,e:/simpele/result/vl65

4、0.txt,w)；filehandlepl650=fopen(e:/simpele/result/pl650.txt,w);filehandleT1650=fopen(e:/simpele/result/T1650.txt,w);filehandlerhol650=fopen(e:/simpele/result/rhol650.txt,w);filehandleyfl650=fopen(,e:/simpele/result/yfl650.txtJw);filehandleomegattl650 = fopen(”e:/simpele/result/omegattl650.txt”Jw”);fi

5、lehandleWml650=fopen(e:/simpele/result/Wml650.txt,w);filehandleEccl650=fopen(e:/simpele/result/Eccl650.txt,w);for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleul650,%e ,fijl);if(j = = ml) fprintf(filehandleul650,n);fprintf(filehandlevl650,%e ;fi皿 2);if(j=ml) fprintf(filehandlevl650,n);fprintf(file

6、handlepl650,%e ,fij5);if(j = = ml) fprintfilehandleplGSO/Xn);fprintf(filehandleT1650；%e ,fij4);if(j = = ml) fprintf(filehandleT1650,n);fprintf(filehandlerhol650,1,%e ,fij6);if(j = = ml) fprintf(filehandlerhol650,n);fprintf(filehandleyfl650,%e ,fij3);if(j=ml) fprintf(filehandleyfl650,n);fprintf(fileh

7、andleomegattl650,%e ,omegaij)；if(j = = ml) fprintf(filehandleomegattl650,n);fprintf(filehandleWml650,%e ,Wmij);if(j = = ml) fprintf(filehandleWml650,n,)；fprintf(filehandleEccl650,%e ,Eccij);if(j=ml) fprintf(filehandleEccl650,n);)/for(i=l;ill;i+)/)fclose(filehandleul650);fclose(filehandlevl650);fclos

8、e(filehandlepl650);fclose(filehandleT1650);fclose(filehandlerhol650);fclose(filehandleyfl650);fclose(filehandleomegattl650);fclose(filehandleWml650);fclose(filehandleEccl650);)if(tt=1700)filehandleul700 = fopen(”e:/simpele/result/ul700.txt”Jw”)；filehandlevl700=fopen(e:/simpele/result/vl700.txt,w);fi

9、lehandlepl700=fopen(e:/simpele/result/pl700.txtJw);filehandleT1700=fopen(e:/simpele/result/T1700.txt,w);filehandlerhol700=fopen(e:/simpele/result/rhol700.txt,1iw);filehandleyfl700=fopen(e:/simpele/result/yfl700.txt,w)；filehandleomegattl700=fopen(e:/simpele/result/omegattl700.txt,w); filehandleWml700

10、=fopen(e:/simpele/result/Wml700.txt,，w); filehandleEccl700=fopen(e:/simpele/result/Eccl700.txt,w); for(i=l;i = ll;i+)for(j=l;j=ml;j + 4-)fprintf(filehandleul700,%e ,fijl);if(j = = ml) fprintf(filehandleul700,n);fprintf(filehandlevl700,%eif(j=ml) fprintf(filehandlevl700,n);fprintf(filehandlepl700,%e皿

11、 5);if(j = = ml) fprintf(filehandlepl700,n);fprintf(filehandleT1700,H%e ,fij4);if(j=ml) fprintfilehandleTlYOO/Xn);fprintf(filehandlerhol700,%e ,fij6);if(j = = ml) fprintf(filehandlerhol700,n);fprintf(filehandleyfl700,%e皿 3);if(j = = ml) fprintf(filehandleyfl700,n);fprintf(filehandleomegattl700,%e ,o

12、megai|j)； if(j = = ml) fprintfilehandleomegattiyOO/Xn)；fprintf(filehandleWml700：,%e ,Wmij);if(j = = ml) fprintf(filehandleWml700,n);fprintf(filehandleEccl700,%e ,Eccij);if(j= = ml) fprintf(filehandleEccl700,n);)/for(i=l;ill;i+)/)fclose(filehandleul700);fclose(filehandlevl700);fclose(filehandlepl700)

13、;fclose(filehandleT1700);fclose(filehandlerhol700);fclose(filehandleyfl700);fclose(filehandleomegattl700);fclose(filehandleWml700);fclose(filehandleEccl700);) if(tt=1750)(filehandleul750=fopen(e:/simpele/result/ul750.txtJ,w);filehandlevl750=fopen(,e:/simpele/result/vl750.txt,w);filehandlepl750=fopen

14、(e:/simpele/result/pl750.txtJw);filehandleT1750=fopen(e:/simpele/result/T1750.txt,w);filehandlerhol750=fopen(,e:/simpele/result/rhol750.txt,w);filehandleyfl750=fopen(e:/simpele/result/yfl750.txt,w)；filehandleomegattl750=fopen(e:/simpele/result/omegattl750.txt,w);filehandleWml750=fopen(e:/simpele/res

15、ult/Wml750.txt,w);filehandleEccl750=fopen(e:/simpele/result/Eccl750.txt,w);for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleul750,%e ,fijl);if(j = = ml) fprintfilehandleulYSO/An);fprintf(filehandlevl750,%e ,fij2);if(j=ml) fprintf(filehandlevl750,n);fprintf(filehandlepl750,%e ,fij5);if(j = = ml) fp

16、rintf(filehandlepl750,n)；fprintf(filehandleT1750,%e H,fij4);if(j=ml) fprintf(filehandleT1750,n);fprintf(filehandlerhol750,%e ,fij6); if(j=ml) fprintf(filehandlerhol750,n);fprintf(filehandleyfl750,%e ,fij3);if(j = = ml) fprintf(filehandleyfl750,n);fprintf(filehandleomegattl750,%e ,omegaij)； if(j = =

17、ml) fprintf(filehandleomegattl750,n);fprintf(filehandleWml750,%e ,Wmij);if(j = = ml) fprintf(filehandleWml750,n)；fprintf(filehandleEccl750,%e ,Eccij);if(j=ml) fprintf(filehandleEccl7501n);/for(i=l;ill;i+)/)fclose(filehandleul750);fclose(filehandlevl750);fclose(filehandlepl750);fclose(filehandleT1750

18、);fclose(filehandlerhol750);fclose(filehandleyfl750);fclose(filehandleomegattl750);fclose(filehandleWml750);fclose(filehandleEccl750);)if(tt=1800)(filehandleul800=fopen(e:/simpele/result/ul800.txt,w);filehandlevl800=fopen(,e:/simpele/result/vl800.txt,w);filehandlepl800=fopen(e:/simpele/result/pl800.

19、txt,w);filehandleT1800 = fopen(”e:/simpele/result/T1800.txt,“w”)；filehandlerhol800=fopen(e:/simpele/result/rhol800.txt,1,w);filehandleyfl800=fopen(e:/simpele/result/yfl800.txt,w);filehandleomegattl800=fopen(e:/simpele/result/omegattl800.txt,w);filehandleWml800=fopen(e:/simpele/result/Wml800.txt,w);f

20、ilehandleEccl800=fopen(e:/simpele/result/Eccl800.txt,w);for(i=l;i=ll;i+)for(j=l;j=ml;j +) fprintf(filehandleul800,%e ,fijl)；if(j = = ml) fprintf(filehandleul800,n);fprintf(filehandlevl800,%e ,fij2);if(j=ml) fprintf(filehandlevl800,n);fprintf(filehandlepl800,%e ,fij5);if(j=ml) fprintf(filehandlepl800

21、,n);fprintf(filehandleT1800,%e 1i皿 4);if(j = = ml) fprintf(filehandleT1800,n);fprintf(filehandlerhol800,%e ,fij6);if(j=ml) fprintf(filehandlerhol800,n);fprintf(filehandleyfl800,%e ,fij3);if(j=ml) fprintf(filehandleyfl800,n);fprintf(filehandleomegattl800,%e ,omegaij)； if(j = = ml) fprintfilehandleome

22、gattlSOO/Xn);fprintf(filehandleWml800,%e ,Wmij)；if(j=ml) fprintf(filehandleWml800/n);fprintf(filehandleEccl800,%e ,Eccij)； if(j= = ml) fprintf(filehandleEccl800,n);)/for(i=l;ill;i+)/)fclose(filehandleul800);fclose(filehandlevl800);fclose(filehandlepl800);fclose(filehandleT1800);fclose(filehandlerhol

23、800);fclose(filehandleyfl800);fclose(filehandleomegattl800);fclose(filehandleWml800);fclose(filehandleEccl800);)if(tt= = 1850)filehandleul850=fopen(e:/simpele/result/ul850.txt,w);filehandlevl850=fopen(e:/simpele/result/vl850.txt,w)； filehandlepl850 = fopen(e:/simpele/result/pl850.txt”Jw) filehandleT

24、1850=fopen(e:/simpele/result/T1850.txt,w); filehandlerhol850=fopen(e:/simpele/result/rhol850.txt,w);filehandleyfl850=fopen(e:/simpele/result/yfl850.txt,w)；filehandleomegattl850=fopen(e:/siinpele/result/omegattl850.txt1,w)； filehandleWml850=fopen(,e:/simpele/result/Wml850.txt,w);filehandleEccl850=fop

25、en(,e:/simpele/result/Eccl850.txt,w);for(i=l;i=ll;i+)for(j=l;j=ml;j +) fprintf(filehandleul850,%e ,fijl)； if(j=ml) fprintf(filehandleul850,n);fprintf(filehandlevl850,%e ;fi皿 2);if(j=ml) fprintf(filehandlevl850,n);fprintf(filehandlepl850,%e ,fij5);if(j=ml) fprintf(filehandlepl850,n);fprintf(filehandl

26、eT1850,%e ,fij4);if(j = = ml) fprintf(filehandleT1850,n);fprintf(filehandlerhol850,%e ,fij6);if(j = = ml) fprintf(filehandlerhol850,n);fprintf(filehandleyfl850,%e ,fij3);if(j=ml) fprintf(filehandleyfl850,n);fprintf(filehandleomegattl850,%e ,omegaij)； if(j = = ml) fprintfilehandleomegattlSSO/Xn);fpri

27、ntf(filehandleWml850,%e ,Wmij);if(j=ml) fprintfilehandleWmlSSO/An);fprintf(filehandleEccl850,%e ,Eccij);if(j=ml) fprintf(filehandleEccl850,n);/for(i=l;ill;i+) /)fclose(filehandleul850);fclose(filehandlevl850);fclose(filehandlepl850);fclose(filehandleT1850);fclose(filehandlerhol850);fclose(filehandle

28、yfl850);fclose(filehandleomegattl850);fclose(filehandleWml850);fclose(filehandleEccl850);if(tt=1900)filehandleul900=fopen(,e:/simpele/result/ul900.txt,1,w);filehandlevl900=fopen(,e:/simpele/result/vl900.txt,w)；filehandlepl900=fopen(e:/simpele/result/pl900.txt,w);filehandleT1900=fopen(e:/simpele/resu

29、lt/T1900.txt,w);filehandlerhol900=fopen(e:/simpele/result/rhol900.txt,w);filehandleyfl900=fopen(e:/simpele/result/yfl900.txt,w)；filehandleomegattl900=fopen(,e:/simpele/result/omegattl900.txt,w);filehandleWml900=fopen(e:/simpele/result/Wml900.txt,w);filehandleEccl900=fopen(e:/simpele/result/Eccl900.t

30、xt1A/);for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleul900,%e ,fijl);if(j=ml) fprintf(filehandleul900,n);fprintf(filehandlevl900,%e ,fij2);if(j = = ml) fprintf(filehandlevl900,n);fprintf(filehandlepl900,%e n,fij5);if(j = = ml) fprintfilehandleplGOO/Xn);fprintf(filehandleT1900,%e ,fij4);if(j = =

31、 ml) fprintf(filehandleT1900,n);fprintf(filehandlerhol900,%e ,fij6);if(j = = ml) fprintf(filehandlerhol900,n);fprintf(filehandleyfl900,%e ,fij3);if(j=ml) fprintf(filehandleyfl900,n);fprintf(filehandleomegattl900,%e ,omegaij)；if(j = = ml) fprintf(filehandleomegattl900,n);fprintf(filehandleWml900,%e ,

32、Wmij);if(j = = ml) fprintf(filehandleWml900,n,)；fprintf(filehandleEccl900,%e ,Eccij);if(j=ml) fprintf(filehandleEccl900Jn);)/for(i=l;ill;i + +)/)fclose(filehandleul900);fclose(filehandlevl900);fclose(filehandlepl900);fclose(filehandleT1900);fclose(filehandlerhol900);fclose(filehandleyfl900);fclose(f

33、ilehandleomegattl900);fclose(filehandleWml900);fclose(filehandleEccl900);if(tt=1950)filehandleul950 = fopen(”e:/simpele/result/ul950.txt“Jw”)；filehandlevl950=fopen(e:/simpele/result/vl950.txt1w);filehandlepl950=fopen(e:/simpele/result/pl950.txtJw);filehandleT1950=fopen(e:/simpele/result/T1950.txt,w)

34、;filehandlerhol950=fopen(e:/simpele/result/rhol950.txt,w);filehandleyfl950=fopen(,e:/simpele/result/yfl950.txt,w)；filehandleomegattl950=fopen(e:/simpele/result/omegattl950.txt,w);filehandleWml950=fopen(e:/simpele/result/Wml950.txt,w);filehandleEccl950=fopen(e:/simpele/result/Eccl950.txt,w);for(i=l;i

35、=ll;i+)for(j=l;j=ml;j+)fprintf(filehandleul950,%e ,fijl);if(j = = ml) fprintf(filehandleul950,n);fprintf(filehandlevl950,%e ,fij2);if(j=ml) fprintf(filehandlevl950,n);fprintf(filehandlepl950,%e ,fij5);if(j = = ml) fprintf(filehandlepl950,n)；fprintf(filehandleT1950,H%e ,fij4);if(j=ml) fprintf(filehan

36、dleT1950,n)；fprintf(filehandlerhol950,%e ,fij6);if(j = = ml) fprintf(filehandlerhol950,n);fprintf(filehandleyfl950,%e皿 3);if(j=ml) fprintf(filehandleyfl950,n)；fprintf(filehandleomegattl950,%e ,omegai|j)；if(j = = ml) fprintf(filehandleomegattl950,n)；fprintf(filehandleWml950J%e ,Wmij)；if(j = = ml) fpr

37、intf(filehandleWml950,n);fprintf(filehandleEccl950,%e ,Eccij);if(j=ml) fprintf(filehandleEccl950,n);)/for(i=l;ill;i+)/11fclose(filehandleul950);fclose(filehandlevl950);fclose(filehandlepl950);fclose(filehandleT1950);fclose(filehandlerhol950);fclose(filehandleyfl950);fclose(filehandleomegattl950);fcl

38、ose(filehandleWml950);fclose(filehandleEccl950);addl.c#include global_var.hvoid addl()int i;double tempc,tempw;for(i=l;i=(int)(nt/z);i +)(tempc=alllocationi*z+l-alllocation(i-l)*z+l;tempw=blllocationi*z+l-blllocation(i-l)*z+l;Sfci=tempc/(z*dt);Sfwi=tempw/(z*dt);)if(tt=2000)(filehandleu2000=fopen(e:/

39、simpele/result/u2000.txt,w);filehandlev2000=fopen(e:/simpele/result/v2000.txtJw);filehandlep2000=fopen(e:/simpele/result/p2000.txt,w);filehandleT2000=fopen(,e:/simpele/result/T2000.txt,w);filehandlerho2000=fopen(e:/simpele/result/rho2000.txt,1,w);filehandleyf2000=fopen(e:/siinpele/result/yf2000.txt,

40、w);filehandleomegatt2000= fopen(e:/simpele/result/omegatt2000.txt”Jw”);filehandleWm2000=fopen(e:/simpele/result/Wm2000.txt,w);filehandleEcc2000= fopen(”e:/simpele/result/Ecc2000.txt“Jw”)；for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleu2000,|1%e 1i皿 1);if(j = = ml) fprintf(filehandleu2000,n);fpri

41、ntf(filehandlev2000,%e ,fij2);if(j = = ml) fprintf(filehandlev2000,n);fprintf(filehandlep2000,%e皿 5);if(j = = ml) fprintf(filehandlep2000,n)；fprintf(filehandleT2000J%e fi皿 4);if(j = = ml) fprintf(filehandleT2000,n)；fprintf(filehandlerho2000；,%e n,fij6);if(j = = ml) fprintf(filehandlerho2000,n);fprin

42、tf(filehandleyf2000,%e ,fij3);if(j = = ml) fprintf(filehandleyf2000,n);fprintf(filehandleomegatt2000,%e ,omegaij);if(j = = ml) fprintf(filehandleomegatt2000,(,n);fprintf(filehandleWm2000,%e ,Wmij);if(j=ml) fprintf(filehandleWm2000,1,nn);fprintf(filehandleEcc2000,%e ,Eccij);if(j=ml) fprintf(filehandl

43、eEcc2000,n)；/for(i=l;ill;i+)/)fclose(filehandleu2000);fclose(filehandlev2000);fclose(filehandlep2000);fclose(filehandleT2000);fclose(filehandlerho2000);fclose(filehandleyf2000);fclose(filehandleomegatt2000);fclose(filehandleWm2000);fclose(filehandleEcc2000);if(tt=2200)filehandleu2200=fopen(e:/simpel

44、e/result/u2200.txt,w)；filehandlev2200=fopen(,e:/simpele/result/v2200.txt1,w);filehandlep2200=fopen(,e:/simpele/result/p2200.txtJw);filehandleT2200=fopen(e:/simpele/result/T2200.txt,w);filehandlerho2200=fopen(e:/simpele/result/rho2200.txt,w);filehandleyf2200=fopen(e:/simpele/result/yf2200.txt,w);file

45、handleomegatt2200=fopen(,e:/simpele/result/omegatt2200.txt,Jw)； filehandleWm2200=fopen(e:/simpele/result/Wm2200.txt,w);filehandleEcc2200=fopen(e:/simpele/result/Ecc2200.txt1w);for(i=l;i=ll;i+)for(j=l;j=ml;j+)fprintf(filehandleu2200,%e ,fijl); if(j=ml) fprintf(filehandleu2200,n);fprintf(filehandlev22

46、00,%e ,fij2); if(j = = ml) fprintf(filehandlev2200,n)；fprintf(filehandlep2200,n%e ,fij5); if(j=ml) fprintf(filehandlep2200Jn);fprintf(filehandleT2200,%e ,fij4); if(j=ml) fprintf(filehandleT22001n);fprintf(filehandlerho2200,%e ,fij6);if(j = = ml) fprintf(filehandlerho2200,n);fprintf(filehandleyf2200,

47、%e ,fij3);if(j=ml) fprintf(filehandleyf2200ln);fprintf(filehandleomegatt2200,%e ,omegaij)； if(j = = ml) fprintf(filehandleomegatt2200/,n);fprintf(filehandleWm2200,%e ,Wmij);if(j=ml) fprintf(filehandleWm2200J,n);fprintf(filehandleEcc2200,%e ,Eccij);if(j=ml) fprintf(filehandleEcc2200,n);/for(i=l;ill;i

48、+)/)fclose(filehandleu2200);fclose(filehandlev2200);fclose(filehandlep2200);fclose(filehandleT2200);fclose(filehandlerho2200);fclose(filehandleyf2200);fclose(filehandleomegatt2200);fclose(filehandleWm2200);fclose(filehandleEcc2200);if(tt=2400)filehandleu2400=fopen(,e:/simpele/result/u2400.txtJ,w);fi

49、lehandlev2400=fopen(,e:/simpele/result/v2400.txt,w)；filehandlep2400=fopen(e:/simpele/result/p2400.txt,w);filehandleT2400=fopen(e:/simpele/result/T2400.txt,w);filehandlerho2400=fopen(e:/simpele/result/rho2400.txt,w);filehandleyf2400=fopen(e:/simpele/result/yf2400.txt,w)；filehandleomegatt2400=fopen(e:

50、/simpele/result/omegatt2400.txt1w)； filehandleWm2400=fopen(e:/simpele/result/Wm2400.txt,w);filehandleEcc2400=fopen(e:/simpele/result/Ecc2400.txt,1,w)； for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleu2400,%e ,fijl)；if(j=ml) fprintf(filehandleu2400,n);fprintf(filehandlev2400,%e ,fij2);if(j=ml) fpr

51、intf(filehandlev2400,n);fprintf(filehandlep2400J%e ,fij5);if(j = = ml) fprintf(filehandlep2400,n);fprintf(filehandleT2400J,%e ,fij4);if(j=ml) fprintf(filehandleT2400,n);fprintf(filehandlerho2400,%e ,fij6);if(j = = ml) fprintf(filehandlerho2400,n);fprintf(filehandleyf2400,n%e ,fi|j3);if(j = = ml) fpr

52、intf(filehandleyf2400,n);fprintf(filehandleomegatt2400,%e ,omegaij);if(j = = ml) fprintf(filehandleomegatt2400Jn”)；fprintf(filehandleWm2400J%e ,Wmij);if(j = = ml) fprintf(filehandleWm24001n);fprintf(filehandleEcc2400,%e ,Eccij);if(j=ml) fprintf(filehandleEcc24001n);)/for(i=l;ill;i+)/1fclose(filehand

53、leu2400);fclose(filehandlev2400);fclose(filehandlep2400);fclose(filehandleT2400);fclose(filehandlerho2400);fclose(filehandleyf2400);fclose(filehandleomegatt2400);fclose(filehandleWm2400);fclose(filehandleEcc2400);)if(tt=2600)filehandleu2600 = fopen(”e:/simpele/result/u2600.txt”Jw”)；filehandlev2600=f

54、open(e:/simpele/result/v2600.txt,w);filehandlep2600=fopen(e:/simpele/result/p2600.txtJw);filehandleT2600=fopen(e:/simpele/result/T2600.txt,w);filehandlerho2600=fopen(e:/simpele/result/rho2600.txt,w);filehandleyf2600=fopen(e:/simpele/result/yf2600.txt,w);filehandleomegatt2600=fopen(e:/siinpele/result

55、/omegatt2600.txt,w);filehandleWm2600=fopen(,e:/simpele/result/Wm2600.txt,w);filehandleEcc2600=fopen(1,e:/simpele/result/Ecc2600.txt,w);for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleu2600,%eif(j = = ml) fprintf(filehandleu2600Jn);fprintf(filehandlev2600,%e ,fij2);if(j=ml) fprintf(filehandlev2600

56、,n);fprintf(filehandlep2600,%e ,fij5);if(j = = ml) fprintf(filehandlep2600,n);fprintf(filehandleT2600,%e 1i皿 4);if(j = = ml) fprintf(filehandleT2600,n);fprintf(filehandlerho2600,%e ,fij6);if(j=ml) fprintf(filehandlerho2600,n);fprintf(filehandleyf2600,%e ”而皿 3);if(j = = ml) fprintf(filehandleyf2600,n

57、);fprintf(filehandleomegatt2600,%e ,omegaij)； if(j二二ml) fprintf(filehandleomegatt26001n)；fprintf(filehandleWm2600,%e ,Wmij);if(j=ml) fprintf(filehandleWm2600,n)；fprintf(filehandleEcc2600,%e ,Eccij);if(j= = ml) fprintf(filehandleEcc2600,n);/for(i=l;ill;i+)/(/)fclose(filehandleu2600);fclose(filehandle

58、v2600);fclose(filehandlep2600);fclose(filehandleT2600);fclose(filehandlerho2600);fclose(filehandleyf2600);fclose(filehandleomegatt2600);fclose(filehandleWm2600);fclose(filehandleEcc2600);)if(tt=2800)(filehandleu2800=fopen(e:/simpele/result/u2800.txt,w);filehandlev2800=fopen(e:/simpele/result/v2800.t

59、xt,w);filehandlep2800=fopen(,e:/simpele/result/p2800.txt1w);filehandleT2800=fopen(e:/simpele/result/T2800.txt1w);filehandlerho2800=fopen(e:/simpele/result/rho2800.txt,w);filehandleyf2800=fopen(e:/simpele/result/yf2800.txt,w);filehandleomegatt2800=fopen(e:/simpele/result/omegatt2800.txt,w); filehandl

60、eWm2800=fopen(e:/simpele/result/Wm2800.txt,w);filehandleEcc2800= fopen(”e:/simpele/result/Ecc2800.txt”Jw”)；for(i=l;i=ll;i+)for(j=l;j=ml;j +)fprintf(filehandleu2800/%e ,fijl);if(j=ml) fprintf(filehandleu2800,n);fprintf(filehandlev2800,%e皿 2);if(j = = ml) fprintf(filehandlev2800,n);fprintf(filehandlep