三层蒸发模型编程示例

#includestdio.h#includestdlib.hvoidmain(){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=20,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]Ep[i]&&WL[i]=C*WLM){El[i]=(Ep[i]-Eu[i])*WL[i]/WLM;Ed[i]=0;}if(WU[i]+P[i]Ep[i]&&WL[i]C*WLM&&WL[i]=C*(Ep[i]-Eu[i])){El[i]=C*(Ep[i]-Eu[i]);Ed[i]=0;}if(WU[i]+P[i]Ep[i]&&WL[i]C*(Ep[i]-Eu[i])){El[i]=WL[i];Ed[i]=C*(Ep[i]-Eu[i])-El[i];}if((P[i]+WU[i]-Eu[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]=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\n,8+i,Eu[i],El[i],Ed[i],E[i],WU[i],WL[i],WD[i]);}system(pause);}