| 查看: 849 | 回复: 28 | ||
| 本帖产生 1 个 程序强帖 ,点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | ||
[求助]
格式问题
|
||
|
write(19,100) ikpt,(kpt(i),i=1,3),npw 100 FORMAT('k-point #',I3,': (',3f7.4,') npw=',I6) 这个句子哪里有问题? |
» 猜你喜欢
酰胺脱乙酰基
已经有13人回复
有院领导为了换新车,用横向课题经费买了俩车
已经有8人回复
博士延得我,科研能力直往上蹿
已经有8人回复
同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗
已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人)
已经有5人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优
已经有6人回复
面上项目申报
已经有3人回复
遇见不省心的家人很难过
已经有22人回复


4楼2011-04-26 12:29:01
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2011-04-26 08:33:16
holmescn
金虫 (正式写手)
- 程序强帖: 37
- 应助: 1 (幼儿园)
- 金币: 1918.8
- 散金: 275
- 红花: 1
- 帖子: 699
- 在线: 102.6小时
- 虫号: 913482
- 注册: 2009-11-26
- 性别: GG
- 专业: 凝聚态物性 II :电子结构
3楼2011-04-26 08:38:43
|
这是程序: 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













回复此楼
