fortran读取nc文件，读出的数据跟原始数据不能对应

桂圆莲子汤 · 发表于 2017-12-3 22:46:03

按照网上的方法写了一个读取nc文件的小程序nc文件为

Band1（1145,702）原始数据为

用该程序得出的结果为

，请问是哪里出了错误
程序如下：

! 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

桂圆莲子汤 · 发表于 2017-12-4 09:31:26

@fcode大神，请问Band1为什么读不出来

yanghw · 发表于 2018-3-8 00:36:36

ierr=nf_get_var_real (ncid,varid,area)这是读单精度变量的，你用读双精度的试试，我一般直接用nf_get_var就可以了。

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