版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

zhaohuxian

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3090.5
散金: 2
帖子: 314
在线: 94.6小时
虫号: 474887
注册: 2007-12-06
性别: MM
专业: 凝聚态物性 II ：电子结构

[求助] 请大侠帮忙看一下下面的程序

program angle

!Define the parameters

   implicit real*8(a-h,o-z)
   parameter (x=10000)
   parameter (y=10000)
   parameter (radtodeg=57.29577951d0)

integer i, j, atomkind, lower, upper
integer iatom(x)

   dimension rv(3,3),temp(6)
   dimension dx(x,y), dy(x,y), dz(x,y)

   character(len=132) line, system
   character*4 lable(x), latom(x,y)

   double precision lattice

   logical eof, con, out

eof=.false.
con=.false.
out=.false.

! Check the required files

inquire(file='OUTCAR',exist=out)
inquire(file='CONTCAR',exist=con)

if (.not. out) then
print *, "OUTCAR: file not exist."
stop
end if

if (.not. con) then
print *, "CONTCAR: file not exist."
stop
end if

! Read the species in OUTCAR

open(unit=out,file='OUTCAR',status='OLD')
atomkind=0
do while(.not.eof)
   read(out,*,end=10) line
   if(index(line,'VRHFIN').ne.0) then
      atomkind=atomkind+1
      upper=index(line,':')-1
      lower=index(line,'=')+1
      lable(atomkind)=line(lower:upper)
   end if
end do
close(out)
10 continue

   ! Read POSCAR to process lattice information and amount of each kind

open(unit=con,file='CONTCAR',status='OLD')
read(con,*,end=20) system
read(con,*,end=20) lattice
read(con,*,end=20) rv(1,1), rv(2,1), rv(3,1)
read(con,*,end=20) rv(1,2), rv(2,2), rv(3,2)
read(con,*,end=20) rv(1,3), rv(2,3), rv(3,3)
read(con,*,end=20), (iatom(i), i=1,atomkind)

! Process the label list
do i=1,atomkind
   do j=1,iatom(i)
      latom(j,i)=lable(i)
   end do
end do

do i=1,3
   temp(i)=0.0
   do j=1,3
      temp(i)=temp(i)+rv(j,i)**2
   end do
   temp(i)=sqrt(temp(i))
end do
a=temp(1)
b=temp(2)
c=temp(3)

do i=1,3
   temp(3+i)=0.0
end do

do j=1,3
   temp(4)=temp(4)+rv(j,2)*rv(j,3)
   temp(5)=temp(5)+rv(j,1)*rv(j,3)
   temp(6)=temp(6)+rv(j,1)*rv(j,2)
end do
temp(4)=temp(4)/(temp(2)*tem(3))
temp(5)=temp(5)/(temp(1)*tem(3))
temp(6)=temp(6)/(temp(1)*tem(2))

alpha=radtodeg*acos(temp(4))
beta=radtodeg*acos(temp(5))
gamma=radtodeg*acos(temp(6))

20 continue

   ! Open file for record, and write

open (11,file=opsition)
write(11,300) system
write(11,300) "lattice lengths: "
write(11,300) "a=", a
write(11,300) "b=", b
write(11,300) "c=", c
write(11,300) "lattice angles: "
write(11,300) "alpha=", alpha
write(11,300) "beta=", beta
write(11,300) "gamma=", gamma

close(con)

! Read OUTCAR to process the coordinates

open(unit=out,file='OUTCAR',status='OLD')
do while(.not.eof)
   read(out,*,end=30) line
   if (index(line,'POSITION').ne.0) then
   read(out,*,end=30) line

      do i=1,atomkind
      do j=1,iatom(i)
         read(out,*,end=30) dx(j,i), dy(j,i), dz(j,i)
         write(11,300) latom(j,i), dx(j,i), dy(j,i), dz(j,i)
      end do
      end do

   end if
end do
close(out)
30 continue

300 format(a4,1x,3f15.9,1x,a4,1x,i4,2(1x,a2),1x,f8.4,1x,i4)

stop
   end

编译通过了，但是云心的时候报错：

--------------------Configuration: adf - Win32 Debug--------------------
Linking...
main.obj : error LNK2001: unresolved external symbol _TEM@4
Debug/adf.exe : fatal error LNK1120: 1 unresolved externals
Error executing link.exe.

类似的另外一个程序运行没有问题

回复此楼

» 猜你喜欢

346求调剂[0856] 已经有7人回复
一志愿山大07化学 332分四六级已过本科山东双非求调剂！已经有3人回复
310求调剂已经有3人回复
机械专硕299求调剂至材料已经有4人回复
070300化学319求调剂已经有7人回复
08工科 320总分求调剂已经有6人回复
一志愿天津大学化学工艺专业（081702）315分求调剂已经有12人回复
307求调剂已经有3人回复
265求调剂已经有3人回复
301求调剂已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-03-07 21:10:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

snoopyzhao

至尊木虫 (职业作家)

程序强帖: 16
应助: 157 (高中生)
贵宾: 0.02
金币: 18844.7
红花: 29
帖子: 3803
在线: 1422.4小时
虫号: 183750
注册: 2006-02-13
专业: 污染生态化学

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
jjdg(金币+2): 感谢应助 2012-03-08 01:56:40

引用回帖:

3楼: Originally posted by zhaohuxian at 2012-03-07 21:29:07:
不好意思，我也是第一次写程序，我是用CVF 6.5编译通过的
不知道哪里有问题，谢谢

你的 cvf 好强啊……

你看一下下面的编译错误提示：

CODE:

C:\Documents and Settings\Administrator\Desktop>gfortran -Wall -Wextra zf.f90

zf.f90:15.19:

      dimension dx(x,y), dy(x,y), dz(x,y)

                   1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:15.21:

      dimension dx(x,y), dy(x,y), dz(x,y)

                     1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:15.18:

      dimension dx(x,y), dy(x,y), dz(x,y)

                  1

Error: The module or main program array 'dx' at (1) must have constant shape

zf.f90:15.28:

      dimension dx(x,y), dy(x,y), dz(x,y)

                            1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:15.30:

      dimension dx(x,y), dy(x,y), dz(x,y)

                              1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:15.27:

      dimension dx(x,y), dy(x,y), dz(x,y)

                           1

Error: The module or main program array 'dy' at (1) must have constant shape

zf.f90:15.37:

      dimension dx(x,y), dy(x,y), dz(x,y)

                                     1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:15.39:

      dimension dx(x,y), dy(x,y), dz(x,y)

                                       1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:15.36:

      dimension dx(x,y), dy(x,y), dz(x,y)

                                    1

Error: The module or main program array 'dz' at (1) must have constant shape

zf.f90:11.22:

        integer iatom(x)

                      1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:11.24:

        integer iatom(x)

                        1

Error: The module or main program array 'iatom' at (1) must have constant shape

zf.f90:18.24:

      character*4 lable(x), latom(x,y)

                        1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:18.26:

      character*4 lable(x), latom(x,y)

                          1

Error: The module or main program array 'lable' at (1) must have constant shape

zf.f90:18.34:

      character*4 lable(x), latom(x,y)

                                  1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:18.36:

      character*4 lable(x), latom(x,y)

                                    1

Error: Expression at (1) must be of INTEGER type, found REAL

zf.f90:18.38:

      character*4 lable(x), latom(x,y)

                                      1

Error: The module or main program array 'latom' at (1) must have constant shape

zf.f90:45.18:

        open(unit=out,file='OUTCAR',status='OLD')

                  1

Error: UNIT tag at (1) must be of type INTEGER

zf.f90:48.15:

          read(out,*,end=10) line

               1

Error: UNIT specification at (1) must be an INTEGER expression or a CHARACTER va

riable

zf.f90:56.14:

        close(out)

              1

Error: UNIT tag at (1) must be of type INTEGER

zf.f90:61.18:

        open(unit=con,file='CONTCAR',status='OLD')

                  1

Error: UNIT tag at (1) must be of type INTEGER

zf.f90:62.13:

        read(con,*,end=20) system

             1

Error: UNIT specification at (1) must be an INTEGER expression or a CHARACTER va

riable

zf.f90:63.13:

        read(con,*,end=20) lattice

             1

Error: UNIT specification at (1) must be an INTEGER expression or a CHARACTER va

riable

zf.f90:64.13:

        read(con,*,end=20) rv(1,1), rv(2,1), rv(3,1)

             1

Error: UNIT specification at (1) must be an INTEGER expression or a CHARACTER va

riable

zf.f90:65.13:

        read(con,*,end=20) rv(1,2), rv(2,2), rv(3,2)

             1

Error: UNIT specification at (1) must be an INTEGER expression or a CHARACTER va

riable

zf.f90:66.13:

        read(con,*,end=20) rv(1,3), rv(2,3), rv(3,3)

             1

Error: UNIT specification at (1) must be an INTEGER expression or a CHARACTER va

riable

Fatal Error: Error count reached limit of 25.

赞一下(1人)

回复此楼

» 本帖已获得的红花（最新10朵）

zhaohuxian

5楼2012-03-07 23:17:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

普通回帖

lurencyj

木虫 (著名写手)

应助: 159 (高中生)
金币: 2869.2
散金: 520
红花: 8
沙发: 10
帖子: 1244
在线: 148.3小时
虫号: 888093
注册: 2009-10-29
性别: GG
专业: 凝聚态物性I:结构、力学和

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jjdg(金币+1): 感谢应助 2012-03-08 01:56:17
zhaohuxian(金币+30): ★有帮助 2 2012-03-08 09:01:51

不知道楼主是怎么写的程序，gfortran编译，一大堆Error！

赞一下

回复此楼

很女子很弓虽大

2楼2012-03-07 21:23:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaohuxian

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3090.5
散金: 2
帖子: 314
在线: 94.6小时
虫号: 474887
注册: 2007-12-06
性别: MM
专业: 凝聚态物性 II ：电子结构

引用回帖:

楼: Originally posted by lurencyj at 2012-03-07 21:23:48:
不知道楼主是怎么写的程序，gfortran编译，一大堆Error！

不好意思，我也是第一次写程序，我是用CVF 6.5编译通过的
不知道哪里有问题，谢谢

赞一下

回复此楼

3楼2012-03-07 21:29:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

uboat

银虫 (初入文坛)

应助: 3 (幼儿园)
金币: 397
帖子: 44
在线: 25.4小时
虫号: 627033
注册: 2008-10-15
性别: GG
专业: 生物电子学与生物信息处理

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg(金币+1): 感谢应助 2012-03-08 01:56:28

temp(4)=temp(4)/(temp(2)*tem(3))
temp(5)=temp(5)/(temp(1)*tem(3))
temp(6)=temp(6)/(temp(1)*tem(2))

看最后的错误提示应该是这里的错误
里面没有定义tem的变量
楼主好好的看看

赞一下

回复此楼

4楼2012-03-07 22:25:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lwk205

至尊木虫 (著名写手)

应助: 7 (幼儿园)
金币: 11155.9
散金: 198
红花: 16
帖子: 1655
在线: 1223.7小时
虫号: 648089
注册: 2008-11-07
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg(金币+1): 感谢应助 2012-03-08 01:56:44

没有定义tem为数组

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

zhaohuxian

6楼2012-03-08 01:11:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lurencyj

木虫 (著名写手)

应助: 159 (高中生)
金币: 2869.2
散金: 520
红花: 8
沙发: 10
帖子: 1244
在线: 148.3小时
虫号: 888093
注册: 2009-10-29
性别: GG
专业: 凝聚态物性I:结构、力学和

引用回帖:

6楼: Originally posted by lwk205 at 2012-03-08 01:11:45:
没有定义tem为数组

记得去找房东

[ 发自手机版 http://muchong.com/3g ]

回复此楼

很女子很弓虽大

7楼2012-03-08 03:10:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaohuxian

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3090.5
散金: 2
帖子: 314
在线: 94.6小时
虫号: 474887
注册: 2007-12-06
性别: MM
专业: 凝聚态物性 II ：电子结构

送鲜花一朵

引用回帖:

楼: Originally posted by lwk205 at 2012-03-08 01:11:45:
没有定义tem为数组

谢谢，但是我不知道为啥不能给金币了，sorry

赞一下

回复此楼

8楼2012-03-08 09:15:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaohuxian

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3090.5
散金: 2
帖子: 314
在线: 94.6小时
虫号: 474887
注册: 2007-12-06
性别: MM
专业: 凝聚态物性 II ：电子结构

送鲜花一朵

引用回帖:

楼: Originally posted by snoopyzhao at 2012-03-07 23:17:32:
你的 cvf 好强啊……

你看一下下面的编译错误提示：

C:\Documents and Settings\Administrator\Desktop>gfortran -Wall -Wextra zf.f90
zf.f90:15.19:

dimension dx(x,y), dy(x,y), dz(x,y)
...

谢谢您的热心帮助，根据帖子里的错误信息，我改了程序：
program angle
!Define the parameters

   parameter (x=1000)
   parameter (y=1000)
   parameter (radtodeg=57.29577951d0)

integer i, j, atomkind, lower, upper
integer iatom(x)

   dimension rv(3,3),temp(6)
   dimension dx(x,y), dy(x,y), dz(x,y)

   character(len=132) line, system
   character*4 lable(x), latom(x,y)

   double precision lattice

   logical eof, con, out

eof=.false.
con=.false.
out=.false.

! Check the required files

inquire(file='OUTCAR',exist=out)
inquire(file='CONTCAR',exist=con)

if (.not. out) then
print *, "OUTCAR: file not exist."
stop
end if

if (.not. con) then
print *, "CONTCAR: file not exist."
stop
end if

! Read the species in OUTCAR

open(unit=12,file='OUTCAR',status='OLD')
atomkind=0
do while(.not.eof)
   read(12,*,end=10) line
   if(index(line,'VRHFIN').ne.0) then
      atomkind=atomkind+1
      upper=index(line,':')-1
      lower=index(line,'=')+1
      lable(atomkind)=line(lower:upper)
   end if
end do
close(12)
10 continue

   ! Read POSCAR to process lattice information and amount of each kind

open(unit=12,file='CONTCAR',status='OLD')
read(12,*,end=20) system
read(12,*,end=20) lattice
read(12,*,end=20) rv(1,1), rv(2,1), rv(3,1)
read(12,*,end=20) rv(1,2), rv(2,2), rv(3,2)
read(12,*,end=20) rv(1,3), rv(2,3), rv(3,3)
read(12,*,end=20), (iatom(i), i=1,atomkind)

! Process the label list
do i=1,atomkind
   do j=1,iatom(i)
      latom(j,i)=lable(i)
   end do
end do

do i=1,3
   temp(i)=0.0
   do j=1,3
      temp(i)=temp(i)+rv(j,i)**2
   end do
   temp(i)=sqrt(temp(i))
end do
a=temp(1)
b=temp(2)
c=temp(3)

do i=1,3
   temp(3+i)=0.0
end do

do j=1,3
   temp(4)=temp(4)+rv(j,2)*rv(j,3)
   temp(5)=temp(5)+rv(j,1)*rv(j,3)
   temp(6)=temp(6)+rv(j,1)*rv(j,2)
end do
temp(4)=temp(4)/(temp(2)*temp(3))
temp(5)=temp(5)/(temp(1)*temp(3))
temp(6)=temp(6)/(temp(1)*temp(2))

alpha=radtodeg*acos(temp(4))
beta=radtodeg*acos(temp(5))
gamma=radtodeg*acos(temp(6))

20 continue

   ! Open file for record, and write

open (11,file=opsition)
write(11,300) system
write(11,300) "lattice lengths: "
write(11,300) "a=", a
write(11,300) "b=", b
write(11,300) "c=", c
write(11,300) "lattice angles: "
write(11,300) "alpha=", alpha
write(11,300) "beta=", beta
write(11,300) "gamma=", gamma

close(con)

! Read OUTCAR to process the coordinates

open(unit=16,file='OUTCAR',status='OLD')
do while(.not.eof)
   read(16,*,end=30) line
   if (index(line,'POSITION').ne.0) then
   read(16,*,end=30) line

      do i=1,atomkind
      do j=1,iatom(i)
         read(1,*,end=30) dx(j,i), dy(j,i), dz(j,i)
         write(11,300) latom(j,i), dx(j,i), dy(j,i), dz(j,i)
      end do
      end do

   end if
end do
close(16)
30 continue

300 format(a4,1x,3f15.9,1x,a4,1x,i4,2(1x,a2),1x,f8.4,1x,i4)

   stop
end

赞一下

回复此楼

9楼2012-03-08 10:54:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zhaohuxian

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3090.5
散金: 2
帖子: 314
在线: 94.6小时
虫号: 474887
注册: 2007-12-06
性别: MM
专业: 凝聚态物性 II ：电子结构

送鲜花一朵

引用回帖:

但是运行的时候debug出现了这样的信息，想问一下为什么，是第一次程序，谢谢热心帮助
Loaded 'ntdll.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\kernel32.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\imagehlp.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\msvcrt.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\dbghelp.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\version.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\advapi32.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\rpcrt4.dll', no matching symbolic information found.
Loaded 'C:\WINDOWS\system32\secur32.dll', no matching symbolic information found.
The program 'E:\SSC_date\temp\ADF\Debug\main.exe' has exited with code 0 (0x0).

赞一下

回复此楼

10楼2012-03-08 10:55:39

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zhaohuxian 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 08工科 320总分求调剂 +6	梨花珞晚风 2026-03-17	6/300	2026-03-21 03:40 by JourneyLucky
[考研] 南京大学化学376求调剂 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 294求调剂材料与化工专硕 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 288求调剂 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 287求调剂 +7	晨昏线与星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 一志愿西南交通专硕材料355 本科双非求调剂 +5	西南交通专材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 一志愿南京航空航天大学大学，080500材料科学与工程学硕 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 材料考研调剂 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐阳w
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 无际的草原
[考研] 材料工程专硕调剂 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 334求调剂 +3	志存高远意在机� 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考博] 26博士申请 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3	橘颂. 2026-03-17	4/200	2026-03-17 21:43 by 有只狸奴
[考研] 308求调剂 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿南京大学，080500材料科学与工程，调剂 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3	听风就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 热情沙漠
[考研] 304求调剂 +5	素年祭语 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 本科南京大学一志愿川大药学327 +3	麦田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明

24小时热门版块排行榜

[求助] 请大侠帮忙看一下下面的 程序

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

【答案】应助回帖

» 本帖已获得的红花（最新10朵）

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

» 本帖已获得的红花（最新10朵）

[求助] 请大侠帮忙看一下下面的程序