爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3322|回复: 7

[程序设计] 求解MATLAB编程

[复制链接]

新浪微博达人勋

发表于 2016-9-24 15:14:53 | 显示全部楼层 |阅读模式

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

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

x
addpath(genpath('netcdf_tools'))  
nc=netcdf('sst.mon.mean.nc','nowrite');
lat = nc{'lat'}(:);
lon = nc{'lon'}(:);
time = nc{'time'}(:);
tp=1273:1284;
tt=time(tp);
T = nc{'sst'}(31,-10:10,50:70);

我这个样子运行是可以的,但是当time=tt(1)时,就不会得到T值了。为什么?
因为tt(1)太大了?tt(1)=38716;(31是time,-10:10是经度范围;50:70是纬度范围)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-9-24 16:34:20 | 显示全部楼层
时间是在第三维吧?我不太清楚,建议楼主看一下nc的数据信息,时间在哪一维
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2016-9-24 16:37:42 | 显示全部楼层
为什么不用ncread呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-9-24 16:45:01 | 显示全部楼层
oucdaqi 发表于 2016-9-24 16:34
时间是在第三维吧?我不太清楚,建议楼主看一下nc的数据信息,时间在哪一维

这是我用ncdump看到的:nc{'sst'} = ncfloat('time', 'lat', 'lon'); %% 97718400 elements.
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-9-24 16:45:29 | 显示全部楼层

不会用,老师教的就是用这个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2016-9-24 18:30:42 | 显示全部楼层
liushujie 发表于 2016-9-24 16:45
不会用,老师教的就是用这个

那你学啊…………………………
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-9-24 21:57:46 | 显示全部楼层
二爷名声在外 发表于 2016-9-24 18:30
那你学啊…………………………

老师没有让学啊~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-9-29 14:11:29 | 显示全部楼层
ncdisp 和 ncread 比较方便吧
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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