| 查看: 1860 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
sci_papers银虫 (正式写手)
|
[交流]
【求助】帮忙改写展宽程序 已有3人参与
|
||
|
大家好,我想用高斯展宽程序,现在得到了一个,但是感觉不是很满意,再次贴出,麻烦大家帮我看看怎么修改,而且我也不太懂编程,所以大家修改后麻烦也全部贴上.谢谢,我现在的程序如下: program CONVOLUTE implicit none integer noofmodes,arbitrary parameter(arbitrary=3500) double precision wavenr(arbitrary),IRint(arbitrary) integer weigenv ! wavenumber start and finish double precision wbegin,wend character*80 fname character*30 tc ! Filename: inputfilename read(5,*) tc,fname ! wbegin: 0.0 read(5,*) tc,wbegin ! wend: 3500.0 read(5,*) tc,wend ! print the data print *,'Filename: ',fname print *,'wbegin: ',wbegin print *,'wend: ',wend print *,'calling subroutine...' call readfreq(fname,wavenr,IRint,weigenv,wbegin) call convolution(wbegin,wend,wavenr,IRint,weigenv) stop ! end program end subroutine readfreq(fname,wavenr,IRint,weigenv,wbegin) ! print *,'beginning of subroutine...' ! read in eigenfrequencies and IR intensities implicit none integer noofmodes,arbitrary parameter(arbitrary=3500) double precision wavenr(arbitrary),IRint(arbitrary) double precision energy(2,arbitrary),wbegin character*80 fname ! print *,'opening the file' integer i,reason,weigenv double precision wavetmp, inttmp open (unit=3,file=fname,status='OLD') i=0 ! print *,'entering do statement' do read(3,*,IOSTAT=reason) wavetmp,inttmp if (reason>0) then print *,'reason>0, i.e. somethings wrong' else if (reason<0) then print *,'reason<0, reached the end of file' close(3) return else if(wavetmp.gt.wbegin) then i=i+1 ! print *,'eigenfreq: ',wavetmp,'IR intensity',inttmp wavenr(i)=wavetmp IRint(i)=inttmp print *,'eigenfreq: ',wavenr(i),' IR intensity: ',IRint(i),' i = ',i weigenv=i endif end if end do return ! print *,'end of subroutine' end subroutine subroutine convolution(wbegin,wend,wavenr,IRint,weigenv) implicit none integer bignumber,arbitrary parameter(bignumber=10000000,arbitrary=3500) double precision wbegin,wend,wavenr(arbitrary),IRint(arbitrary) integer noofscanpoints,n,m,k,weigenv parameter(noofscanpoints=100000) double precision deltaw,w1(bignumber),conv(bignumber) double precision mygamma,wave2ev parameter(wave2ev=8065.46) ! mygamma=4.0d-3*wave2ev mygamma=10.0 ! mygamma=4.0d-3 in eV unit, convert to wavenumber by multiple wave2ev deltaw=(wend-wbegin)/dfloat(noofscanpoints) do n=1,noofscanpoints w1(n)=wbegin+n*deltaw conv(n)=0.0 enddo do m=1,weigenv do k=1,noofscanpoints ! print *,'m=',m,'k=',k conv(k)=conv(k)+IRint(m)/((w1(k)-wavenr(m))**2+mygamma**2) ! print *,'w1=',w1(k),' conv(k)=',conv(k) enddo enddo open(unit=1,file='IETSconv.dat',status='unknown') do k=1,noofscanpoints write(1,*) w1(k),conv(k) enddo close(1) return end subroutine 现在麻烦大家: 1. 麻烦检查下这个程序的高斯展宽是否正确. 自己没用过,真不知道.谢谢 2. 这个程序需要提供jobinput.txt文件,然后这个文件中含有相关数值: inputfile anth.txt wbegin 1000.0 wend 1800.0 这些数值也包括了需要展宽的波段数,比如:1000.0-1800.0。我觉得这个也挺好,因为有时候需要分析的波段不是全部波段, 然后这个展宽程序自动从anth.txt文件中自动读取这个范围的数据展宽.然后会输出一个展宽后的数据文件IETSconv.dat.当然我说这些对于编程的人来说就太啰嗦了,因为大家一看就知道什么意思.不好意思 3. 我现在想要的就是,自己输入需要读取的文件,也输出保存的文件名.这样比较方便.不想这样读取inputfile,因为如果处理数据多的话,每次都要修改inputfile文件,很麻烦.当然在这儿,如果我只输入一个文件,最后数据自动保存为同名的dat类型文件则更好. 比如./gauvib.exe < anth.txt > anth.dat 或者./gauvib.exe < anth.txt 然后得到anth.dat数据文件更好 4. 既然想不用inputfile文件,那这个文件中的波段要在程序中实现了.在程序中加入可以设置波段的语句. 比如: wbegin=1000.0 wend=1800.0 5. 如果谁有这样的或者更好的高斯或者洛伦茨展宽程序可以帮忙给我,因为要用,自己又不会. 再次先感谢大家.深深的希望大家帮个忙.看怎么修改程序. |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高手帮忙写个简单程序
已经有6人回复
- 帮忙编一段程序
已经有3人回复
- 求助一段代换程序
已经有6人回复
- 求助一段程序……
已经有3人回复
- 大家帮忙看一个matlab的小程序，看看哪里出错了
已经有5人回复
- 求助两个完成方程式的题目，麻烦高手指点下，谢了！
已经有7人回复
- 大家帮我看看这个方法怎么操作吧，谢谢~
已经有5人回复
- 求助一道题。。希望大家帮帮我，谢谢大家！！
已经有6人回复
- 大家来帮个忙，看看这是什么东西
已经有6人回复
- 【求助】大家帮忙修改一下程序
已经有17人回复
- 求能修改语言的软件
已经有5人回复
- 【求助】帮忙看一个程序
已经有16人回复
- 【求助】求高手帮忙改写数据实时动态显示程序~
已经有8人回复
- 【求助】求助MOLvib程序
已经有3人回复
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
6楼2010-05-15 21:32:38
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
★ ★
resonant(金币+2):新政策摸索中，感谢回帖交流。 2010-05-14 20:43:10
sci_papers(金币+2):谢谢,swizard,我也试过了,最后也是生成一个相同的文件,如果处理很多数据的话,也比较麻烦.所以就想得到,输入一个然后就输出一个相同文件名的文件.这样感觉在Linux下更好一点.所以还是想麻烦各位帮看看我的这个程序. 2010-05-14 22:46:06
sci_papers(金币+1):我不会编程哦,哎,惭愧.所以麻烦大家帮忙看看,要不看怎么在这个现有的程序上改. 2010-05-14 22:51:15
resonant(金币+2):新政策摸索中，感谢回帖交流。 2010-05-14 20:43:10
sci_papers(金币+2):谢谢,swizard,我也试过了,最后也是生成一个相同的文件,如果处理很多数据的话,也比较麻烦.所以就想得到,输入一个然后就输出一个相同文件名的文件.这样感觉在Linux下更好一点.所以还是想麻烦各位帮看看我的这个程序. 2010-05-14 22:46:06
sci_papers(金币+1):我不会编程哦,哎,惭愧.所以麻烦大家帮忙看看,要不看怎么在这个现有的程序上改. 2010-05-14 22:51:15
|
为什么不直接用 swizard （http://www.sg-chem.net/swizard/）？ 那个有更多的选项啥的…… 如果一定要自己写，你可以参考它的说明书中的相应的方法，呵呵…… |
2楼2010-05-14 20:38:32
sci_papers
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 532.7
- 散金: 40
- 帖子: 328
- 在线: 210.7小时
- 虫号: 811122
- 注册: 2009-07-18
- 性别: GG
- 专业: 能源化工
3楼2010-05-15 18:58:18
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
4楼2010-05-15 19:12:25
200