数组实参传递到虚参出错

chiangtp

subroutine的声明部分,                               請列出各"虛參"的type/dimension宣告
CALL lineplu_v1(...). CALL doubplu_v1(...), 請列出各"實參"的type/dimension宣告

chiangtp

藍原なみだ 发表于 2018-6-27 19:19
subroutine lineplu_v1的宣告：

REAL*8 ALPHA,AREA,B,CO2,COMG,CN2,CNMG,DS,DZ,DENSEA,DENSEP,DENSEW,

1. REAL(8) 是 REAL(KIND=8),  REAL*8 是 8-byte REAL, 兩者不一定相同 (compiler dependent)
    如果有問題, compile-time 應該就知道, 所以就不是這個問題
    建議你通通改為 REAL*8, 或是
    INTEGER, PARAMETER :: r8=SELECTED_REAL_KIND(P=15)
    REAL(KIND=r8) :: ...

2. 請問你的 compiler ?

3. 請告知 完整的 SUBROUTINE title
   SUBROUTIEN doubplu_v1(...)
   SUBROUTINE lineplu_v1(...)

4. 請先自己檢查 "虛參"與"實參"的對應 是否一致

5. ksrc=?, km1=?
    REAL(8), DIMENSION (1:km1) :: qwd, bwd
    CALL lineplu_v1(..., qwd(1:ksrc), bwd(1:ksrc), ...)

chiangtp

chiangtp 发表于 2018-6-27 20:28
1. REAL(8) 是 REAL(KIND=8),  REAL*8 是 8-byte REAL, 兩者不一定相同 (compiler dependent)
    如果有 ...

1. 虛參CO2M有問題,   對應的實參是 DOamb(1:kms)
in SUBROUTINE DOUBPLU_v1,  REAL*8 :: CO2M
in SUBROUTINE LINEPLU_v1,     REAL*8 :: CO2M(1000)

2. 實參 salamb 與 patm 的 type?

3. runtime時, 請確認: km1=1000, krsc=500

chiangtp

4. 實參: CALL ABC(..., zamb (1:kms), Tamb (1:kms), DOamb(1:kms), ..., qwd(1:ksrc), bwd(1:ksrc), ...), array-setction的用法, 沒有必要
   --->  CALL ABC(..., zamb, Tamb, DOamb, ..., qwd, bwd, ...)

chiangtp

應該就是錯在這裡: in SUBROUTINE DOUBPLU_v1
REAL*8 :: CO2M   --->  REAL*8 :: CO2M(1000)

chiangtp

藍原なみだ 发表于 2018-6-28 16:10
非常感谢你，但是我修改为CO2M(1000)以后运行，赋值仍然和以前一样......

實參 salamb 與 patm 的 type?