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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 格式问题
291/3返回列表123下一页
查看: 948  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:21 by vZfe6xYu34yj
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:09 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:33 by vZfe6xYu34yj
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:12 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 QUjhNVAcOSff 2026-09-13 7/350 2026-09-14 21:21 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
关闭
请填处理意见
关闭