Fortran Coder

查看: 13950|回复: 6
打印 上一主题 下一主题

[派生类型] 向派生数据类型中赋值出现奇怪的错误,谢谢各位高手

[复制链接]

4

帖子

1

主题

0

精华

入门

F 币
41 元
贡献
22 点

规矩勋章

跳转到指定楼层
楼主
发表于 2018-2-8 11:23:12 | 显示全部楼层 回帖奖励 |倒序浏览 |阅读模式
各位高手中午好!先向各位问好!

我定义了一个type,用来存储一个特定的化学反应类型的各个数据,如下
[Fortran] 纯文本查看 复制代码
    type :: L2R2
        !===============================================================================
        ! a1*L1^CL1 + a2*L2^CL2 + e*E^(-1) <---> b1*R1^CR1 + b2*R2^CR2 , p = XXX V
        !
        ! for example:
        !
        ! No = 1
        ! idn= 'L2R2_001'
        ! Rx = 'O3 + 2 H^(+) + 2 e^(-) <---> O2 + H2O'   Pot = +2.075 V
        !
        ! a1 = 1 ; L1 = 'O3'  ; CL1 = 0
        ! a2 = 2 ; L2 = 'H'   ; CL1 = 1
        ! e  = 2 ;
        ! b1 = 1 ; R1 = 'O2'  ; CR1 = 0
        ! b2 = 1 ; R1 = 'H2O' ; CL1 = 0
        ! p  = 2.075
        !==============================================================================
        integer        :: No
        character(10)  :: idn
        character(100) :: Rx
        integer        :: a1
        character(20)  :: L1
        integer        :: CL1
        integer        :: a2
        character(20)  :: L2
        integer        :: CL2
        integer        :: e
        integer        :: b1
        character(20)  :: R1
        integer        :: CR1
        integer        :: b2
        character(20)  :: R2
        integer        :: CR2
        real           :: p        
    end type L2R2

然后,两种赋值方式
第一种:
[Fortran] 纯文本查看 复制代码
    type (L2R2) :: L2R2_003
        !==============================================================================
        ! No = 3   Rx = 'ZnO2^(2-) + 2 H2O + 2 e^(-) <---> Zn + 4 OH^(-)'   Pot = -1.285 V
        !==============================================================================
        L2R2_003%No   =  3
        L2R2_003%idn  =  'L2R2_003'
        L2R2_003%Rx   =  'ZnO2^(2-) + 2 H2O + 2 e^(-) <---> Zn + 4 OH^(-)'
        L2R2_003%a1   =  1
        L2R2_003%L1   =  'ZnO2'
        L2R2_003%CL1  =  -2
        L2R2_003%a2   =  2
        L2R2_003%L2   =  'H2O'
        L2R2_003%CL2  =  0
        L2R2_003%e    =  2
        L2R2_003%b1   =  1
        L2R2_003%R1   =  'Zn'
        L2R2_003%CR1  =  0
        L2R2_003%b2   =  4
        L2R2_003%R2   =  'OH'
        L2R2_003%CR2  =  -1
        L2R2_003%p    =  -1.285

这个没有任何问题。

问题是第二种,出错了:
[Fortran] 纯文本查看 复制代码
        L2R2_003 = (/ 3 , "L2R2_003" , "ZnO2^(2-) + 2 H2O + 2 e^(-) <---> Zn + 4 OH^(-)" , &
                      1 , "ZnO2" , -2 , 2 , "H2O" , 0 , 2 , 1 , "Zn" , 0 , 4 , "OH" , -1 , -1.285 /)

编译时候的错误是
  1. Compiling ..\..\..\..\Desktop\ZZZ\hello.f90
  2. ..\..\..\..\Desktop\ZZZ\hello.f90:104.25:

  3.         L2R2_003 = (/ 3 , "L2R2_003" , "ZnO2^(2-) + 2 H2O + 2 e^(-) <---> Zn + 4 OH^(-)" , &
  4.                          1
  5. Error: Element in INTEGER(4) array constructor at (1) is CHARACTER(1)
  6. Error(E42): Last command making (build\hello.o) returned a bad status
  7. Error(E02): Make execution terminated

  8. * Failed *
复制代码

所用环境就是Simply Fortran 2.17 Build 1758
实在不能理解这个错误发生的原因,还请明示!十分感谢

我自己的理解是:( / ...... / )这种方式,实际上是先造一个数组或列表出来,然后依次对应到派生类型中的各个变量里。
而fortran不允许数组内各元素类型不同,他会去认第一个元素的类型,一看,是整形,3,于是整个数组开始被建立,是个整形的数组。但是读到第二个元素时,却变成了"L2R2_003",是字符串,好,报错!

像L2R2_003%No=3,L2R2_003%idn='L2R2_003'这样的赋值方式,太繁琐。那么,有什么好办法么?

啰嗦啰嗦写了这么多,谢谢您耐心看完,再次感谢所有提出意见的人!
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

4

帖子

1

主题

0

精华

入门

F 币
41 元
贡献
22 点

规矩勋章

沙发
 楼主| 发表于 2018-2-8 14:41:32 | 显示全部楼层
fcode 发表于 2018-2-8 12:29
(/  /) 只能构造一个数组。(实际上我建议用 [] 来构造 )

要构造 L2R2 的结构体,需要用  L2R2()

非常感谢您的回复!
另外,从您帮我修改的帖子格式来看,您太细心了。
顺便学会了mw_shl_code=fortran,这是您建站的时候自定义的吧?
再次感谢您的回答

4

帖子

1

主题

0

精华

入门

F 币
41 元
贡献
22 点

规矩勋章

板凳
 楼主| 发表于 2018-2-8 17:02:35 | 显示全部楼层
pasuka 发表于 2018-2-8 15:17
需要用到面向对象编程的内容,类似于C++的构造函数
[mw_shl_code=fortran,true]module example
    implici ...

谢谢!十分感谢提供例子!
我这儿只是做一个静态的数据库,为其他子程序提供一些基础数据(其实就是常用电极的基础数据),所以不用class啦
再次感谢!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-6 19:07

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表