主要程序代码
cleara10b0hab500w1100w250w30001fori2501k11h11w3w111w3w1w11000005w2k12h115w3w2115w3w2w201w1k13h000001w1w3000002w2w3
k21h11w305k13w105k1111w305k13w105k11w105k111000005w205k12k22h115w305k13w205k12115w305k13w205k12w205k1201w105k11k23h000001w105k11w305k13000002w205k12w305k13
k31h11w305k23w105k2111w305k23w105k21w105k211000005w205k22k32h115w305k23w205k22115w305k23w205k22w205k2201w105k21k33h000001w105k21w305k23000002w205k22w305k23
k41h11w3k33w1k3111w3k33w1k31w1k311000005w2k32k42h115w3k33w2k32115w3k33w2k32w2k3201w1k31k43h000001w1k31w3k33000002w2k32w3k33
w1w1k112k212k31k416w2w2k122k222k32k426w3w3k132k232k33k436
p1iw1p2iw2p3iw3ppiie
dplotppp1plotppp2plotppp3plotppp3
1
f2Qi
put请输入漏源强度mgs：Qui
put请输入计算风速msudi
put请输人计算精度mdZoi
put请输人地面粗糙长度m：Zo：xymeshgrid0d1000100d100定义解空间和计算精度by0008x100001x12计算Y轴向的基本扩散参数bz0006x100015x12计算z轴向的扩散参数byby01038Zo对Y轴向的扩散参数按地面粗糙长度进行修正。fz253013logx0550042logx1Zo035003logx按地面粗糙长度计算z轴向的扩散参数修正系数bzbz0fz对轴向的扩散参数按地面粗糙长度进行修正。高斯扩散浓度模拟计算tempy1yybyby2tempy22718282tempy1cQpiubybz1tempy2Csi
put请输人所有求解浓度mgm3所输人数以””和””作为开头和结束。CsCs1000000单位转换co
tourxycCs以XYC变量分别作为的轴绘制扩散浓度分布图shadi
gi
terpgridxlabel轴向距离mylabelY轴向距离mtitle气体扩散浓度分布图
2
fr