请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4704|回复: 1

[混合编程] 求助:nc文件批量读取

[复制链接]

新浪微博达人勋

发表于 2019-10-4 19:22:33 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 突破的渴望 于 2019-10-4 19:23 编辑

请问各位大神 我求的是ncep数据的nc文件中特定经纬度(lon>=120 & lon<=290;lat>=-30 & lat<=30)下 高度为第一层 (level=1)的多年平均,请问这个代码错在哪里? 平均的应该是怎么写呢?写在哪里呢?谢谢大家~~~
  1. clc;  %清屏
  2. clear; %清空
  3. datadir='C:\Users\xuchu\Desktop\air\'; %指定批量数据所在的文件夹
  4. filelist=dir([datadir,'*.nc']); %指定批量数据的类型
  5. a=filelist(1).name; %查看你要读取的文件的编号。filelist(1).name在window下为第一个标号数据
  6. b=filelist(2).name; %查看你要读取的文件的编号。filelist(2).name在window下为第二个标号数据
  7. k=length(filelist);
  8. for s=1:k
  9.    filename=[datadir,filelist(s).name];
  10.    lon=ncread(filename,'lon');
  11.    lat=ncread(filename,'lat');
  12.    c=find(lon>=120 & lon<=290);
  13.    d=find(lat>=-30 & lat<=30);
  14.    lon_num=length(c);
  15.    lat_num=length(d);
  16.    level_lim=ncread(filename,'level',[1], [1]);
  17.    level_num=length(level_lim);
  18.    ai  = ncread(filename,'air',[level_lim c d],[level_num lon_num lat_num]); %读入变量air
  19. end
复制代码

air.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2019-10-5 15:46:58 | 显示全部楼层
问题已解决~~~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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