ncl中调用fortran程序生成so文件后无法调用

松下风

用fortran写了一个读取micaps第四类格式的程序（fortran中可以使用），准备在ncl中调用，成功生成了so文件，最后一步ncl脚本中调用时报错（已排除平台问题，ncl官网中的实例能够运行成功）报错如下：
warning:Could not find Init() in external file /home/jxlc/ncl/ncl-fortran/readm490.so, file not loaded
warning:error at line 2 in file test.ncl
ncl脚本第二行：external EX01 "./readm490.so"附上fortran程序，stub文件和ncl脚本
fortran程序：

subroutine readm4(filename,x,y,&
r1h,diamond,num,title,yy,&
mm,dd,starthr,ybhr,level,xspace,yspace,&
startlon,endlon,startlat,endlat,xnum,ynum,&
cint,cstart,cend,smooth,jiacu)
    character(16) filename
    integer x,y,num,yy,mm,dd,starthr,ybhr,level,xnum,ynum
    integer i,j,k,m
    character(8) diamond
    character(40) title
    real r1h(x,y),xspace,yspace,startlon,endlon,startlat,endlat,cint,cstart,cend,smooth,jiacu
    open(10,file=filename,status='old',ERR=111)
    read(10,*) diamond,num,title
    read(10,*) yy,mm,dd,starthr,ybhr,level,xspace,yspace,startlon,endlon,startlat,endlat,xnum,ynum
    read(10,*) cint,cstart,cend,smooth,jiacu
    do i=1,y
      read(10,*) (r1h(j,i),j=1,x)
    enddo
    close(10)
    111 continue
    return
end
相应的readm490.stub文件内容：
C NCLFORTSTART
subroutine readm4(filename,x,y,&
r1h,diamond,num,title,yy,&
mm,dd,starthr,ybhr,level,xspace,yspace,&
startlon,endlon,startlat,endlat,xnum,ynum,&
cint,cstart,cend,smooth,jiacu)
    character(16) filename
    integer x,y,num,yy,mm,dd,starthr,ybhr,level,xnum,ynum
    integer i,j,k,m
    character(8) diamond
    character(40) title
    real r1h(x,y),xspace,yspace,startlon,endlon,startlat,endlat,cint,cstart,cend,smooth,jiacu
C NCLEND
ncl脚本内容如下：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
external READ "/home/jxlc/ncl/ncl-fortran/readm490.so"
begin
;
; Calculate three values of a quadraticequation
;
  filename="16070108.001"
  x=351
  y=301
  r1h=new((/x,y/),float)
  diamond="20160808"
  num1=1
  yy=1
  mm=1
  dd=1
  starthr=1
  ybhr=1
  level=1
  xnum=1
  ynum=1
  title="20160808"
  xspace=0.1
  yspace=0.1
  startlon=0.1
  endlon=0.1
  startlat=0.1
  endlat=0.1
  cint=0.1
  cstart=0.1

  cend=0.1
  smooth=0.1
  jiacu=0.1
  READ::readm4(filename,x,y,r1h,diamond,num1,title,yy,mm,dd,starthr,ybhr,level,xspace,yspace,startlon,endlon,startlat,endlat,xnum,ynum,cint,cstart,cend,smooth,jiacu) ; Call the NCL version of

end

东西有点乱，请大神指导

robindct

刚才的回复有点问题，重新更正了一下

                               
登录/注册后可看大图

我是神奇分割线

                               
登录/注册后可看大图

有两个问题你需要解决：
1、stub文件必须使用f77格式声明；
2、stub声明部分， 对于数组，ncl向fortran传递时会发生转置，所以要注意(x,y) --> (y,x) 的变换。对于字符串，ncl有string和character（字符型数组）两种，所以要注意ncl和fortran数据类型对应关系。建议使用ncl传递string，在stub声明character*(*)。如果ncl一定要传递character，那么一定要注意数组的size必须一致。

我手头没有Micaps第四类数据，建议楼主发一个样例文件，以便测试

guoyunqian

非常感谢！

robindct

guoyunqian 发表于 2016-8-11 10:22
非常感谢！

更正了一点小错误，敬请参考

松下风

添加个micaps第四类文件供大神测试

松下风

robindct 发表于 2016-8-11 09:57
刚才的回复有点问题，重新更正了一下

我是神奇分割线

感谢回复，收获良多，能不能对f77格式不太熟悉，能不能贴一个stub样本，非常感谢

robindct

根据你的stub修改，仅供参考～～

C NCLFORTSTART
      subroutine readm4(filename,x,y,
     $ r1h,diamond,num,title,yy,
     $ mm,dd,starthr,ybhr,level,xspace,yspace,
     $ startlon,endlon,startlat,endlat,xnum,ynum,
     $ cint,cstart,cend,smooth,jiacu)
        character*16 filename
        integer x,y,num,yy,mm,dd,starthr,ybhr,level,xnum,ynum
        integer i,j,k,m
        character*8 diamond
        character*40 title
        real r1h(x,y)
        real xspace,yspace,startlon,endlon,startlat,
     $   endlat,cint,cstart,cend,smooth,jiacu
C NCLEND

松下风

robindct 发表于 2016-8-11 11:52
根据你的stub修改，仅供参考～～

C NCLFORTSTART

非常感谢，根据你的stub，问题解决了！