24小时热门版块排行榜    

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

贺仪

铁杆木虫 (著名写手)

[求助] 格式问题

write(19,100) ikpt,(kpt(i),i=1,3),npw                     
100    FORMAT('k-point #',I3,':  (',3f7.4,')    npw=',I6)
这个句子哪里有问题?
回复此楼
嘟啦啦嘟啦啦嘟
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖


微尘、梦想(金币+1): 谢谢参与应助! 2011-04-26 17:41:07
这句话,语法上没有问题啊,至少这里可以编译通过……

你是不是不想输出结果中有那对括号??
CODE:
100    FORMAT('k-point #',I3,':  ',(3f7.4),'    npw=',I6)

2楼2011-04-26 08:33:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)


微尘、梦想(金币+1): 谢谢参与! 2011-04-26 17:41:23
gfortran 4.3.3测试,没有问题
3楼2011-04-26 08:38:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by snoopyzhao at 2011-04-26 08:33:16:
这句话,语法上没有问题啊,至少这里可以编译通过……

你是不是不想输出结果中有那对括号??
CODE:
100    FORMAT('k-point #',I3,':  ',(3f7.4),'    npw=',I6)


是的
但是我这里用ifort总是提示出错
嘟啦啦嘟啦啦嘟
4楼2011-04-26 12:29:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

这是程序:
      program WAVECAR
      integer   nkpt,nband,emax,ikpt,npw,iband
      real*8  A kpt
      dimension   A(3,3),kpt(3)           

      open(unit=12,file="WAVECAR",status="old",form="unformatted"
      open(unit=19,file="wave",status="new"
      read(12,*) nkpt,nband,emax,((A(i,j),i=1,3),j=1,3)
      write(19,*)
      write(19,*) 'nkpt  =',nkpt
      write(19,*) 'nband =',nband
      write(19,*) 'emax  =',emax
      write(19,*) 'A='
      write(19,'(3X,3(1X,f8.3))') (A(i,1),i=1,3)
      write(19,'(3X,3(1X,f8.3))') (A(i,2),i=1,3)
      write(19,'(3X,3(1X,f8.3))') (A(i,3),i=1,3)
      read(12,*) code   
                                                
      do ikpt=1, nkpt                                                                                                      
          read(12,*) npw, (kpt(i),i=1,3)
          write(19,*)                                             
          write(19,100) ikpt,(kpt(i),i=1,3),npw
          100  FORMAT('k-point #',I3,':  ',(3f7.4),'    npw=',I6)                        
          write(19,*) 'kpt =',kpt                                 
          write(19,*) 'npw =',npw                                 
          write(19,*) "  band       energy        weight"         
c         
          ALLOCATE(coef(npw))                                                
          ALLOCATE(ac(npw,200))                                            
          ic = 1                                                            
             DO iband = 1, nband                                                                             
               read(12,*) eval, fweight, (coef(i),i=1,npw)                        
                                                                        
c               
               write(19,'(5X,I3,5X,f8.4,5x,f8.4)') iband, dreal(eval), fweight
c               
                                                                        
                IF(iband >= nbandmin.and.iband.le.nbandmax) then               
                     DO i = 1,npw                                                     
                       ac(i,ic) = coef(i)                                            
                     ENDDO                                                            
                       ic = ic+1;                                                      
                     end if                                                           
                                                                                 
            ENDDO                                                              
           DEALLOCATE(coef)                                                   
       ENDDO                                                               
close(unit=12)  
close(unit=19)
c                                                                                   
      stop                                             
      end                                               
c----------------------- end ---------------------------
嘟啦啦嘟啦啦嘟
5楼2011-04-26 12:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by 贺仪 at 2011-04-26 12:31:45:
这是程序:
      program WAVECAR
      integer   nkpt,nband,emax,ikpt,npw,iband
      real*8  A kpt
      dimension   A(3,3),kpt(3)           

      open(unit=12,file="WAVECAR",sta ...

那个笑脸是一个  )
嘟啦啦嘟啦啦嘟
6楼2011-04-26 12:32:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by holmescn at 2011-04-26 08:38:43:
gfortran 4.3.3测试,没有问题

我班程序贴出来了 麻烦你在看看,谢谢咯
嘟啦啦嘟啦啦嘟
7楼2011-04-26 12:33:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

引用回帖:
Originally posted by 贺仪 at 2011-04-26 12:33:48:
我班程序贴出来了 麻烦你在看看,谢谢咯

数据文件?
8楼2011-04-26 14:40:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

引用回帖:
Originally posted by 贺仪 at 2011-04-26 12:31:45:
这是程序:
      program WAVECAR
      integer   nkpt,nband,emax,ikpt,npw,iband
      real*8  A kpt
      dimension   A(3,3),kpt(3)           

      open(unit=12,file="WAVECAR",sta ...

你这程序编译不过去啊,coef 是啥没有定义啊……
9楼2011-04-26 14:44:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by snoopyzhao at 2011-04-26 14:44:50:
你这程序编译不过去啊,coef 是啥没有定义啊……

这个另外有定义,但是你编译时,前面没有错误吗?
嘟啦啦嘟啦啦嘟
10楼2011-04-26 19:28:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:27 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我,科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 02:52 by vs90ilomwc
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 02:45 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 3/150 2026-02-08 02:32 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:27 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:25 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:12 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:07 by vs90ilomwc
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见