版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

本帖产生 1 个程序强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

阿黛拉

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 173.9
散金: 515
红花: 6
帖子: 133
在线: 72.8小时
虫号: 1182523
注册: 2011-01-03
专业: 理论和计算化学

[求助] 哪位帮忙写份Fortran程序？

实现的目标如下
有一份文件包含（注省略号代表省去的内容）
。。。。
。。。。。
1f=。。。。
2f=。。。。。
。。。。。
117 f  = 0.554979 THz    3.487036 2PiTHz 18.512109 cm-1    2.295210 meV
         X       Y       Z                            dx       dy       dz
   10.946913 16.075862  8.898705 -0.004611 -0.000783 0.001817
   14.007941 14.566718  9.506160    0.000936 -0.001987 -0.016250
。。。。

118 f/i= 0.135429 THz    0.850926 2PiTHz 4.517429 cm-1    0.560090 meV
         X       Y       Z          dx       dy       dz
   10.946913 16.075862  8.898705    0.005213 -0.000239 -0.017779
   14.007941 14.566718  9.506160    0.002094 -0.001910 -0.017946
   。。。。。。

119 f/i= 0.159039 THz    0.999273 2PiTHz 5.304980 cm-1    0.657734 meV
         X       Y       Z                               dx       dy       dz
   10.946913 16.075862  8.898705    0.022673 0.022642 -0.016684
   14.007941 14.566718  9.506160    0.007804 0.011458 0.000320
8.454822 15.389561  9.528711 -0.152499 0.049474 -0.037482
   13.830276 10.899067 11.066884 -0.153476 0.086312 0.031470
   12.338448 15.061294  9.339663 -0.169789 0.066115 0.026027
   9.593558  9.592271  5.477447 -0.074791 0.052529 -0.022671
   13.693851 11.017968  7.269960 -0.101499 0.083451 0.019408
   9.904298  9.337579 13.105627 -0.167523 0.071767 -0.030243
   6.232093 11.731397 11.419762 -0.195444 0.039903 -0.103436
   10.015420 13.425170  5.570057 -0.090249 0.056042 -0.030760
   6.064181 11.852228  7.590656 -0.105687 0.040804 -0.073813
   10.337222 13.175130 13.204681 -0.238876 0.033485 -0.011847
   7.633679  7.764385  9.311316 -0.111119 0.040341 -0.055429
   9.462183  7.612316  7.318255 -0.077203 0.052896 -0.022021
   6.351238 13.692965  9.553844 -0.167824 0.053177 -0.090210
   9.607166  7.489229 11.137700 -0.128770 0.065609 -0.030289
   11.433621  7.344679  9.147421 -0.099001 0.063935 -0.000888
   12.007007 11.053752 13.065257 -0.170872 0.075687 0.008814
   5.930483  9.881390  9.455522 -0.137995 0.032910 -0.103045
   11.702241 11.302622  5.446557 -0.083672 0.057803 -0.002193
   7.916422 11.708652  5.607436 -0.081477 0.048928 -0.041310
   14.026017 12.894285  9.217672 -0.147413 0.078933 0.039205
   8.229457 11.458811 13.243556 -0.231964 0.030838 -0.050940
   10.306047 15.272435  7.514793 -0.123688 0.051121 -0.020875
   10.448729 15.124387 11.347733 -0.213204 0.057002 0.013285
   13.551900  9.050957  9.123624 -0.118530 0.089605 0.021981
   7.836243  9.541040 11.362086 -0.166481 0.046336 -0.064171
   11.792534  9.107562 11.181748 -0.137181 0.069631 0.006275
   12.243843 13.090462 11.283023 -0.186831 0.071512 0.018850
   12.093316 13.222266  7.317261 -0.116961 0.064358 0.005111
   7.667685  9.667356  7.378393 -0.091330 0.042996 -0.047419
   8.268422 13.511407 11.475978 -0.205406 0.043669 -0.043708
   11.643939  9.237999  7.210266 -0.090099 0.067539 0.000851
   8.101048 13.650800  7.486080 -0.112671 0.051988 -0.042587
   11.858891 11.176573  9.250229 -0.132036 0.071164 0.004425
   9.739673  9.454709  9.279567 -0.118119 0.058373 -0.021995
   9.877664 11.444483  7.432702 -0.107815 0.061457 -0.019708
   8.037331 11.580078  9.417518 -0.140308 0.050828 -0.046809
   10.170866 13.296060  9.386659 -0.146412 0.057393 -0.005684
   10.027708 11.313744 11.241268 -0.169636 0.056018 -0.015943

120 f/i= 0.255221 THz    1.603598 2PiTHz 8.513242 cm-1    1.055508 meV
         X       Y       Z                                  dx       dy       dz
   10.946913 16.075862  8.898705 -0.017367 0.005145 0.000110
   14.007941 14.566718  9.506160 -0.016450 0.007712 0.003796
   8.454822 15.389561  9.528711 -0.152499 0.049474 -0.037482
   13.830276 10.899067 11.066884 -0.153476 0.086312 0.031470
   12.338448 15.061294  9.339663 -0.169789 0.066115 0.026027
   9.593558  9.592271  5.477447 -0.074791 0.052529 -0.022671
   13.693851 11.017968  7.269960 -0.101499 0.083451 0.019408
   9.904298  9.337579 13.105627 -0.167523 0.071767 -0.030243
   6.232093 11.731397 11.419762 -0.195444 0.039903 -0.103436
   10.015420 13.425170  5.570057 -0.090249 0.056042 -0.030760
   6.064181 11.852228  7.590656 -0.105687 0.040804 -0.073813
   10.337222 13.175130 13.204681 -0.238876 0.033485 -0.011847
   7.633679  7.764385  9.311316 -0.111119 0.040341 -0.055429
   9.462183  7.612316  7.318255 -0.077203 0.052896 -0.022021
   6.351238 13.692965  9.553844 -0.167824 0.053177 -0.090210
   9.607166  7.489229 11.137700 -0.128770 0.065609 -0.030289
   11.433621  7.344679  9.147421 -0.099001 0.063935 -0.000888
   12.007007 11.053752 13.065257 -0.170872 0.075687 0.008814
   5.930483  9.881390  9.455522 -0.137995 0.032910 -0.103045
   11.702241 11.302622  5.446557 -0.083672 0.057803 -0.002193
   7.916422 11.708652  5.607436 -0.081477 0.048928 -0.041310
   14.026017 12.894285  9.217672 -0.147413 0.078933 0.039205
   8.229457 11.458811 13.243556 -0.231964 0.030838 -0.050940
   10.306047 15.272435  7.514793 -0.123688 0.051121 -0.020875
   10.448729 15.124387 11.347733 -0.213204 0.057002 0.013285
   13.551900  9.050957  9.123624 -0.118530 0.089605 0.021981
   7.836243  9.541040 11.362086 -0.166481 0.046336 -0.064171
   11.792534  9.107562 11.181748 -0.137181 0.069631 0.006275
   12.243843 13.090462 11.283023 -0.186831 0.071512 0.018850
   12.093316 13.222266  7.317261 -0.116961 0.064358 0.005111
   7.667685  9.667356  7.378393 -0.091330 0.042996 -0.047419
   8.268422 13.511407 11.475978 -0.205406 0.043669 -0.043708
   11.643939  9.237999  7.210266 -0.090099 0.067539 0.000851
   8.101048 13.650800  7.486080 -0.112671 0.051988 -0.042587
   11.858891 11.176573  9.250229 -0.132036 0.071164 0.004425
   9.739673  9.454709  9.279567 -0.118119 0.058373 -0.021995
   9.877664 11.444483  7.432702 -0.107815 0.061457 -0.019708
   8.037331 11.580078  9.417518 -0.140308 0.050828 -0.046809
   10.170866 13.296060  9.386659 -0.146412 0.057393 -0.005684
   10.027708 11.313744 11.241268 -0.169636 0.056018 -0.015943
实现的想法是将带有f/i  的这些行下的 dx的坐标乘以一个小数N加到X坐标上，以此类推dy*N+Y,dz*N+Z，然后输出加和后的 X Y Z 坐标即可。
特点：1  此处的X  Y  Z坐标行数是40行，不同的文件行数可能不同。2  XYZ坐标，同一文件都是相同的。只是对应的f/i不同dx  dy dz 不同。

[ Last edited by 阿黛拉 on 2011-6-9 at 09:13 ]

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有6人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
心脉受损已经有3人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复

进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。

1楼2011-06-09 08:35:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

【答案】应助回帖

★ ★ ★
jjdg(金币+3): 辛苦了 2011-06-14 12:40:47
余泽成(程序强帖+1): 2011-06-15 21:04:28

引用回帖:

Originally posted by 阿黛拉 at 2011-06-14 09:04:31:
你好
输出三组
这三组是dx*N+X=X#       dy*N+Y=Y#    dz*N+Z=Z# X# Y# Z# 是所想要的新数据。
比如          X       Y       Z                      dx       dy       dz
   10 ...

好吧，大概的程序如下：

CODE:

program outcar

implicit none

real, dimension(:,:), allocatable :: coor, diff

character(len=128) :: line

integer :: num, ios, i

real :: frac

open(unit=12, file='OUTCARf', status='old')

open(unit=13, file='Result', status='new')

do

   read(12,'(a)', iostat=ios) line

   if (ios /= 0) exit

   if (index(line, 'f/i') /= 0) then

       read (12, '(a)') line

       exit

   end if

end do

num = 0

do

   read (12, '(a)') line

!  write(*,*) trim(line)

   if (len_trim(line) == 0) exit

   num = num + 1

end do

write (*,*) num

rewind 12

allocate(coor(num,3), diff(num,3))

write(*,*) 'please input the frac: '

read(*,*) frac

do

   read(12,'(a)', iostat=ios) line

   if (ios /= 0) exit

   if (index(line, 'f/i') /= 0) then

       read (12, '(a)') line

       do i=1,num

          read (12, *) coor(i,1:3), diff(i,1:3)

          write(13,'(3F11.6)') coor(i,1:3) + diff(i,1:3) * frac

       end do

       write(13,*)

   end if

end do

end program outcar

赞一下(3人)

回复此楼

11楼2011-06-14 11:42:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

【答案】应助回帖

★
阿黛拉(金币+2): 2011-06-10 18:07:56
余泽成(金币+1): 谢谢参与应助！ 2011-06-15 21:03:38

引用回帖:

Originally posted by 阿黛拉 at 2011-06-09 08:35:59:
实现的目标如下
有一份文件包含（注省略号代表省去的内容）
。。。。
。。。。。
1f=。。。。
2f=。。。。。
。。。。。
117 f = 0.554979 THz 3.487036 2PiTHz 18.512109 cm-1 2.295210 me ...

给个完整的输入数据文件

赞一下(1人)

回复此楼

3楼2011-06-10 10:46:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

阿黛拉

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 173.9
散金: 515
红花: 6
帖子: 133
在线: 72.8小时
虫号: 1182523
注册: 2011-01-03
专业: 理论和计算化学

引用回帖:

Originally posted by snoopyzhao at 2011-06-10 10:46:07:
给个完整的输入数据文件

这是完整的文件。谢谢你。

赞一下

回复此楼

进退有度，才不至进退维谷；宠辱皆忘，方可以宠辱不惊。

4楼2011-06-10 18:12:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yongma2008

木虫 (著名写手)

论文？苦！主。

程序强帖: 1
应助: 34 (小学生)
金币: 2089.6
散金: 1670
红花: 15
帖子: 1610
在线: 398.5小时
虫号: 1091013
注册: 2010-09-06
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★
jjdg(金币+1): 感谢参与 2011-06-11 00:30:41
阿黛拉(金币+5): 2011-06-11 10:30:50

你传的文件我打不开，所以没调试，程序主体在这里，你根据情况再改改参数，调试一下吧。不保证没问题啊，调试了之后才能知道。

program test
implicit none
integer       ::i,j
integer       ::NF,Nx,M
real(kind=8) ::N,x(1000,100),y(1000,100),z(1000,100),dx(1000,100),dy(1000,100),dz(1000,100)
character    ::inputfile*80,outputfile*80

NF=120  !number of the f values
Nx=40 !number of x values
M=5    !number of the useless lines on the top of the input file
N=0.1 !the parameter in formula x=x+N*dx

write(*,*) "Name of input file:"
read(*,*) inputfile
open(unit=7,file=inputfile)
do i=1,M
read(7,*)
enddo

do i=1,NF,1
read(7,*)
read(7,*)
do j=1,Nx,1
read(7,*) x(i,j),y(i,j),z(i,j),dx(i,j),dy(i,j),dz(i,j)
x(i,j)=x(i,j)+N*dx(i,j)
y(i,j)=y(i,j)+N*dy(i,j)
z(i,j)=z(i,j)+N*dz(i,j)
enddo
read(7,*)
enddo

write(*,*) "Name of the output file:"
read(*,*) outputfile
open(unit=17,file=outputfile)
do i=1,NF,1
write(17,*) 'f=',i
do j=1,Nx,1
write(17,999) x(i,j),y(i,j),z(i,j)
enddo
enddo
999 format(3(f10.6,1x))
end program

赞一下(1人)

回复此楼

6楼2011-06-10 20:49:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答