自定义类型

hanshikai

[

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	PROGRAM AA

02	IMPLICIT NONE

03

04	TYPE::ND !自定义类型:节点

05	INTEGER ::ID !节点编号

06	REAL ::X,Y !节点坐标

07

ENDTYPE

08

09	TYPE::EL !自定义类型:单元

10	INTEGER::ID !单元编号

11	! TYPE(ND%X)::N1

12	! TYPE(ND)::N1,N2

13	TYPE(ND)::N !组成单元的节点,为ND类型,将包含三个数据,读取时,只用到X,Y这两个

14	! TYPE(ND)%ID::N !不能这样定义

15

ENDTYPE

16

17

18	TYPE(ND)::ALLND(1,1)

19	TYPE(EL)::ALLEL(1,1)

20

21	OPEN(11,FILE="T.DAT")

22	READ(11,*) ALLND%ID, ALLND%X, ALLND%Y

23	!READ(11,*) ALLEL

24	!READ(11,*) ALLEL%NUM,ALLEL%ND%X,ALLEL%ND%Y

25	READ(11,*)ALLEL%ID, ALLEL%N%X, ALLEL%N%Y

26	!READ(11,*)ALLEL

27	PRINT*,ALLND

28	PRINT*,ALLEL

29

30

END

数据:
第一行:6 0.33 0.12
第二行:3 5 7

运行的结果是第二行出来四个数,我期望的是3个,在第一个后面多出一个0

[被我注释掉的都是验证过不能用的]

fcode

           6  0.3300000      0.1200000
           3           0   5.000000       7.000000
请按任意键继续. . .

你希望的是什么样子？

hanshikai

上面注释有错啊.

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1	TYPE(ND)::N !组成单元的节点,为ND类型,将包含三个数据,读取时,只用到X,Y这两个

应该是用到两个ND的ID,EL的定义有问题

不过这个应该不是最后输出错误的原因

hanshikai

fcode 发表于 2014-5-22 17:21
6  0.3300000      0.1200000
           3           0   5.000000       7.000000
请按任意键 ...

第二行 3 5 7

fcode

PRINT*,ALLEL
改为
PRINT*,ALLEL%ID, ALLEL%N%X, ALLEL%N%Y

hanshikai

fcode 发表于 2014-5-22 17:28
PRINT*,ALLEL
改为
PRINT*,ALLEL%ID, ALLEL%N%X, ALLEL%N%Y

ALLE是包含四个元素的,因此直接输出会有四个
got it !