新人求助，混编中参数输入出现系列问题

ptxyzs2 · 发表于 2021-11-29 16:54:36

本帖最后由 ptxyzs2 于 2021-11-29 17:05 编辑

VS2019，x86+win32的fortran与c的混编下，出于某些原因没有使用iso_c_binding来混编，而是通过在fortran代码dailymoist.f中编写interface来调用c语言编写的trace_gas_model函数。
代码是对上一个版本代码的更新，interface的编写是根据上一个版本编写的，调用的c语言代码只是增加了一个float *efscltef参数，因此对原interface的修改也仅仅是增加了一个对efscltef参数的声明。上一个版本的代码运行是没有问题的，但是新代码不知道为什么在参数这里出现了多个相关的不同错误：
(a) error #6633: The type of the actual argument differs from the type of the dummy argument.;
(b) error #6784: The number of actual arguments cannot be greater than the number of dummy arguments.;
(c) error #8284: If the actual argument is scalar, the dummy argument shall be scalar unless the actual argument is of type character or is an element of an array that is not assumed shape, pointer, or polymorphic. [DN2LYR, DN2OLYR]。

个人认为问题首先出在error #6784实参、虚参的数目不一致上，但是我多次对过，数目是一样的43个，十分不理解为什么会出现这个问题。而其他的类型不一样、数组与标量的问题，我一一对过，也是一致的，所以实在不知道问题出在哪里，不知道接下来从何入手，还请大家指教。提供的代码比较长，实在是不好意思了

出问题的dailymoist.f文件中，trace_gas_model函数的interface编写如下，共有43个参数， dN2lyr, dN2Olyr为数组：

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

c ... fortran代码，trace_gas_model的interface编写
      INTERFACE
        SUBROUTINE trace_gas_model(newminrl, ammonium, nitrate, texture,
     &                             sand, silt, clay, afiel, bulkd, maxt,
     &                             ppt, snow, avgwfps,stormf, basef,
     &                             frlechd, stream, inorglch, critflow, 
     &                             wfluxout, newCO2, co2_conc, efscltef, time, 
     &                             NOflux, Nn2oflux, Dn2oflux, Dn2flux, CH4,
     &                             isdecid, isagri, aglivc, rleavc, btolai, 
     &                             crpstore, forstore, nit_amt, nreduce, curday, 
     &                             pHscale, dN2lyr, dN2Olyr, month)
          !MS$ATTRIBUTES ALIAS:'_trace_gas_model' :: trace_gas_model
          DOUBLE PRECISION newminrl
          DOUBLE PRECISION ammonium
          DOUBLE PRECISION nitrate(*)
          INTEGER          texture
          REAL             sand
          REAL             silt
          REAL             clay
          REAL             afiel
          REAL             bulkd
          REAL             maxt
          REAL             ppt
          REAL             snow
          REAL             avgwfps
          REAL             stormf
          REAL             basef
          REAL             frlechd(*)
          REAL             stream(*)
          DOUBLE PRECISION inorglch
          REAL             critflow
          REAL             wfluxout(*)
          REAL             newCO2
          REAL             co2_conc(*)
          REAL             efscltef
          REAL             time
          DOUBLE PRECISION NOflux
          DOUBLE PRECISION Nn2oflux
          DOUBLE PRECISION Dn2oflux
          DOUBLE PRECISION Dn2flux
          DOUBLE PRECISION CH4
          INTEGER          isdecid
          INTEGER          isagri
          REAL             aglivc
          REAL             rleavc
          REAL             btolai
          REAL             crpstore
          REAL             forstore
          DOUBLE PRECISION nit_amt
          REAL             nreduce
          INTEGER          curday
          REAL             pHscale
          DOUBLE PRECISION dN2lyr(*)
          DOUBLE PRECISION dN2Olyr(*)
          INTEGER          month
        END SUBROUTINE trace_gas_model

dailymoist.f中调用该trace_gas_model的代码如下，共有43个参数：

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

c ... trace_gas_model调用
      call trace_gas_model(newminrl, ammonium, nitrate, texture, 
     &                     sand, silt, clay, afiel(1), bulkd, maxt,
     &                     ppt(curday), snow, avgwfps, stormf, basef,
     &                     frlechd, stream, inorglch, critflow,
     &                     wfluxout, newCO2, co2_conc, efscltef, time, 
     &                     NOflux, Nn2oflux, Dn2oflux, Dn2flux, CH4, 
     &                     isdecid, isagri, aglivc, rleavc, btolai, 
     &                     crpstg(N),forstg(N),nit_amt,nreduce,curday,
     &                     pHscalar(month), dN2lyr, dN2Olyr, month)

dailymoist.f中与trace_gas_model相关的local variables声明如下，dN2lyr, dN2Olyr为数组：

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

c ... 相关的local variable声明
c ... LOCAL VARIABLES
c ... SWMAXLYR为parameter :: 43
      double precision newminrl, inorglch
      real             avgwfps
      real             frlechd(MAXIEL)
      real             critflow
      real             wfluxout(SWMAXLYR)
      real             newCO2
      real             co2_conc(SWMAXLYR)
      real             dstr, efscltef
      double precision Nn2oflux, Dn2oflux, Dn2flux, NOflux
      double precision CH4, nit_amt
      integer          isdecid, isagri
      double precision dN2lyr(SWMAXLYR), dN2Olyr(SWMAXLYR)

被调用c代码声明如下：

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

/* c source code, tgmodel.c */
    void trace_gas_model(double *newminrl, double *ammonium, double nitrate[], int *texture,
                         float *sand, float *silt, float *clay, float *afiel, float *bulkd, float *maxt,
                         float *ppt, float *snow, float *avgwfps, float *stormf, float *basef,  
                         float frlechd[], float stream[], double *inorglch, float *critflow,
                         float wfluxout[], float *newCO2, float co2PPM[], float *efscltef, float *time,
                         double *NOflux, double *Nn2oflux, double *Dn2oflux, double *Dn2flux, double *CH4,
                         int *isdecid, int *isagri, float *aglivc, float *rleavc, float *btolai,  
                         float *crpstore, float *forstore, double *nit_amt, float *nreduce, int *jday,
                         float *pHscale, double dN2lyr[], double dN2Olyr[], int *month)

如果需要的话，这是上一个版本interface的写法，可以见到我仅仅是多加了一个efscltef的说明：

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

        SUBROUTINE trace_gas_model(newminrl, ammonium, nitrate,
     &                             texture, sand, silt, clay, afiel,
     &                             bulkd, maxt, ppt, snow, avgwfps,
     &                             stormf, basef, frlechd, stream,
     &                             inorglch, critflow, wfluxout,
     &                             newCO2, co2_conc, time, NOflux,
     &                             Nn2oflux, Dn2oflux, Dn2flux, CH4,
     &                             isdecid, isagri, aglivc, rleavc,
     &                             btolai, crpstore, forstore, nit_amt,
     &                             nreduce, curday, pHscale, dN2lyr,
     &                             dN2Olyr, month)
          !MS$ATTRIBUTES ALIAS:'_trace_gas_model' :: trace_gas_model
          DOUBLE PRECISION newminrl
          DOUBLE PRECISION ammonium
          DOUBLE PRECISION nitrate(*)
          INTEGER          texture
          REAL             sand
          REAL             silt
          REAL             clay
          REAL             afiel
          REAL             bulkd
          REAL             maxt
          REAL             ppt
          REAL             snow
          REAL             avgwfps
          REAL             stormf
          REAL             basef
          REAL             frlechd(*)
          REAL             stream(*)
          DOUBLE PRECISION inorglch
          REAL             critflow
          REAL             wfluxout(*)
          REAL             newCO2
          REAL             co2_conc(*)
          REAL             time
          DOUBLE PRECISION NOflux
          DOUBLE PRECISION Nn2oflux
          DOUBLE PRECISION Dn2oflux
          DOUBLE PRECISION Dn2flux
          DOUBLE PRECISION CH4
          INTEGER          isdecid
          INTEGER          isagri
          REAL             aglivc
          REAL             rleavc
          REAL             btolai
          REAL             crpstore
          REAL             forstore
          DOUBLE PRECISION nit_amt
          REAL             nreduce
          INTEGER          curday
          REAL             pHscale
          DOUBLE PRECISION dN2lyr(*)
          DOUBLE PRECISION dN2Olyr(*)
          INTEGER          month
        END SUBROUTINE trace_gas_model