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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 格式问题
291/3返回列表123下一页
查看: 826  |  回复: 28
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

贺仪

铁杆木虫 (著名写手)

[求助] 格式问题

write(19,100) ikpt,(kpt(i),i=1,3),npw                     
100    FORMAT('k-point #',I3,':  (',3f7.4,')    npw=',I6)
这个句子哪里有问题？
回复此楼

» 猜你喜欢
嘟啦啦嘟啦啦嘟
1楼 2011-04-25 22:22:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

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

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

赞一下(1人) 回复此楼
2楼2011-04-26 08:33:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★
微尘、梦想(金币+1): 谢谢参与！ 2011-04-26 17:41:23
gfortran 4.3.3测试，没有问题
赞一下(1人) 回复此楼
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的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
291/3返回列表123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭