24小时热门版块排行榜    

查看: 1590  |  回复: 10
本帖产生 1 个 程序强帖 ,点击这里进行查看

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★
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

11楼2011-06-14 11:42:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿黛拉 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268求调剂 +7 简单点0 2026-03-17 7/350 2026-03-21 08:32 by JourneyLucky
[考研] 313求调剂 +3 肆叁贰壹22 2026-03-19 3/150 2026-03-21 08:01 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 8/400 2026-03-20 15:58 by babero
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
信息提示
请填处理意见