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

登录注册

应《网络安全法》要求，自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.

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

回复此楼

» 猜你喜欢

心脉受损已经有5人回复
博士读完未来一定会好吗已经有15人回复
Springer期刊投稿求助已经有4人回复
读博已经有3人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复

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

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 (幼儿园)
金币: 11154.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 的主题更新

返回列表

24小时热门版块排行榜

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

» 猜你喜欢

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

【答案】应助回帖

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

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

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