Fortran Coder

查看: 58|回复: 1

[绘图界面库] fortran读取nc文件,读出的数据跟原始数据不能对应

[复制链接]

25

帖子

8

主题

0

精华

熟手

F 币
162 元
贡献
103 点
发表于 2017-12-3 22:46:03 | 显示全部楼层 |阅读模式
按照网上的方法写了一个读取nc文件的小程序nc文件为 QQ图片20171203224312.png
Band1(1145,702)原始数据为 QQ图片20171203224429.png
用该程序得出的结果为 QQ图片20171203224521.png ,请问是哪里出了错误
程序如下:

! This program is for reading area.nc
      program read
      implicit none
      integer, parameter :: longitude=3600, latitude=1800
      integer ierr,ncid,varid,varidx,varidy,len_file,err
      real(8)  area(longitude,latitude)
      real lon(longitude),lat(latitude)
      character*299 barl_out
      include 'netcdf.inc'
      barl_out='/home/iga/data/xycao/05area/barley/op/barl_harvest_2005_01.nc'

      len_file=len_trim(barl_out)
      ierr=nf_open(trim(barl_out),nf_write,ncid)    !open nc file,get ID number (ncid)
      ierr=nf_inq_varid (ncid, 'longitude',varidx)
      ierr=nf_get_var_real(ncid,varidx,lon)
      ierr=nf_inq_varid (ncid, 'latitude',varidy)
      ierr=nf_get_var_real(ncid,varidy,lat)
      ierr=nf_inq_varid (ncid, 'Band1', varid)      !open 'Band1',get ID number
      ierr=nf_get_var_real (ncid,varid,area)
      print*,'ierr'
      print*, area(1145,702),area(1145,704)                                !check
      err=nf_close(ncid)
      end program read



回复

使用道具 举报

25

帖子

8

主题

0

精华

熟手

F 币
162 元
贡献
103 点
 楼主| 发表于 2017-12-4 09:31:26 | 显示全部楼层
@fcode大神,请问Band1为什么读不出来
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|QQ群|Fcode

GMT+8, 2017-12-18 15:08

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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