爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3483|回复: 0

[求助] fortran读取nc文件,经纬度与原nc文件不对应[已经解决]

[复制链接]

新浪微博达人勋

发表于 2017-12-2 17:47:16 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 cxyykl 于 2017-12-2 19:50 编辑

对不起,不知道 怎么删贴,给大家添麻烦了,已经解决了,谢谢nc文件中纬度范围 QQ图片20171204161138.png
用程序输出的纬度:-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文件中的经纬度?拜托给个解答,谢谢大家!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

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