求助！急！高斯公式的fortran源代码，新手跪求大神们给予帮助

Kelvina_mu

各位大神们，小女初来乍到，还不太懂用，还请大家谅解！只是上面这个公式的fortran码我急求，但是自己写不出来，才疏学浅，望莫笑话。。上面主要是输入各项数据最后能得到q的值就行。Q、y、z、σ、μ都是要根据实际情况来输入的。谢谢啦~~~！！！！！
祝大家新年快乐，马年行大运，大吉大利！！！
貌似我刚开的账号还不会怎么奖励，希望大家帮帮忙

又是那隻貓

是个妹子@版主@管理员

言深深

额，应该输入的是xyz吧？否则的话计算器都可以求出来了？
要是不会fortran干嘛非要fortran，个人感觉，excel都能求的···
ps：fortran程序的话不是很难，建议楼主思考一下，再想想，你的公式明白否？到底输入什么，得到什么？？

Kelvina_mu

言深深 发表于 2014-2-3 11:22

                               
登录/注册后可看大图

额，应该输入的是xyz吧？否则的话计算器都可以求出来了？
要是不会fortran干嘛非要fortran，个人 ...

因为只是考虑x风向（u风）所以浓度分布公式中只剩下了yz两个如果是xyz的话就是无风点源的了。其实这个程序最后要用matlab作动态图表现出来，可是我现在还什么都不会，比较没头没脑，所以想先求一段fortran程序看一看。。。是不是有点儿无厘头了？？？

Kelvina_mu

又是那隻貓 发表于 2014-2-3 11:11

                               
登录/注册后可看大图

是个妹子@版主@管理员

谢谢你啊

言深深

Kelvina_mu 发表于 2014-2-3 11:53

                               
登录/注册后可看大图

因为只是考虑x风向（u风）所以浓度分布公式中只剩下了yz两个如果是xyz的话就是无风点源的了。其实这个程序 ...

这样做的意义不大，建议你直接用matlab，或者找会matlab的朋友帮忙指导，否则fortran编出来也没有用···

金色牛神

PROGRAM Gauss
       IMPLICIT NONE
       REAL Q,SIGMAz,SIGMAy,y,z,miu,result
       PRINT *,请依次输入Q,y,z,SIGMAy,SIGMAz,miu
       result=Q/(6.283*SIGMAy*SIGMAz*miu)*EXP(-y*y/2/SIGMAy**2-z*z/SIGMAz**2)
       print *,result
END

把上述代码拷进去，大致调试一下就可以了，可其实这种计算用MATLAB最方便了，有空可以接触一下……

金色牛神

呀，刚才那段代码在回帖时出现了差错，第4，6，8，9那四行代码不该换行……

Kelvina_mu

言深深 发表于 2014-2-3 14:06

                               
登录/注册后可看大图

这样做的意义不大，建议你直接用matlab，或者找会matlab的朋友帮忙指导，否则fortran编出来也没有用··· ...

谢谢~

Kelvina_mu

金色牛神 发表于 2014-2-3 15:25

                               
登录/注册后可看大图

PROGRAM Gauss
       IMPLICIT NONE
       REAL Q,SIGMAz,SIGMAy,y,z,miu,result

十分感谢！啊！