- 积分
- 39
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-3-22
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 cxyykl 于 2017-12-2 19:50 编辑
对不起,不知道 怎么删贴,给大家添麻烦了,已经解决了,谢谢nc文件中纬度范围
:
用程序输出的纬度:-179.9~0
程序如下:
! This program is for reading area.nc
program read
implicit none
integer, parameter :: longitude=3600, latitude=1800
integer ierr,ncid,varid,len_file
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_nowrite,ncid) !open nc file,get ID number (ncid)
ierr=nf_inq_varid (ncid, 'longitude',varid) !open longitude
ierr=nf_get_var_real(ncid,varid,lon)
ierr=nf_inq_varid (ncid, 'latitude',varid) !open latitude
ierr=nf_get_var_real(ncid,varid,lat)
ierr=nf_inq_varid (ncid, 'Band1', varid) !open 'Band1',get ID number
ierr=nf_get_var_real (ncid,varid,area)
print*,'ierr'
print*, lon(1),lat(1)
print*, area(1145,702),area(1145,704) !check
end program read
求助:为什么程序不读取原nc文件中的经纬度?拜托给个解答,谢谢大家!
|
|