三层蒸发模型编程示例.doc

《三层蒸发模型编程示例.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、#include#includevoidmain(){floatP[18]={0,0,0.8,0,0,0.2,0,0,0,0,0,0.3,0,0,0,7.4,0.1,0};floatEp[18]={7.9,7.4,5.9,6.1,6.2,5.8,5.0,5.2,5.4,6.9,6.7,4.1,5.8,4.0,4.3,5.9,4.2,6.3};floatEu[18]={0},El[18]={0},Ed[18]={0},E[18],WU[18]={0},WL[18]={14.9},WD[18]={40.0},C=1.0/6;intWUM=2

2、0,WLM=60,WDM=40;for(inti=0;i<=17;i++){if(WU[i]+P[i]>=Ep[i]){Eu[i]=Ep[i];El[i]=0;Ed[i]=0;}elseEu[i]=WU[i]+P[i];if(WU[i]+P[i]=C*WLM){El[i]=(Ep[i]-Eu[i])*WL[i]/WLM;Ed[i]=0;}if(WU[i]+P[i]=C*(Ep[i]-Eu[i])){El[i]=C*(Ep[i]-Eu[i]);Ed[i]=0;}if(WU[i]+P[i]

3、i]&&WL[i]=0)WU[i+1]=P[i]+WU[i]-Eu[i];elseWU[i+1]=0;if((WL[i]-El[i])>=0)WL[i+1]=WL[i]-El[i];elseWL[i+1]=0;if((WD[i]-Ed[i])>=0)WD[i+1]=WD[i]-Ed[i];elseWD[i+1]=0;E[i]=Eu[i]+El[i]+Ed[i];if((WL[i]-El[i])>=0)WL[i+1]

4、=WL[i]-El[i];elseWL[i+1]=0;if((WD[i]-Ed[i])>=0)WD[i+1]=WD[i]-Ed[i];elseWD[i+1]=0;E[i]=Eu[i]+El[i]+Ed[i];printf("1970.8.%d：Eu=%.3f,El=%.3f,Ed=%.3f,E=%.3f，WU=%.3f,WL=%.3f,WD=%.3f",8+i,Eu[i],El[i],Ed[i],E[i],WU[i],WL[i],WD[i]);}system("pause");}