本文主要是介绍(三)NCL读取经度:西经10度 (东经350度)-东经20度,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
题目
nc文件中lon的范围是0-360
我只会两种方法
方法1:
{lon|(/350,352.5,355,357.5,0,2.5,5,7.5,10,12.5,15,17.5,20/)}
begin
f = addfile("D:/cygwin/home/wunan/hgt.day.mean.200001-200004.nc","r")hgt = f->hgt
hgtF = short2flt(f->hgt)
time = f->timetime_1 = cd_calendar(time,2)
timeInd1 = ind(time_1.eq.20000101)
timeInd2 = ind(time_1.eq.20000107)
hgtF_sub = hgtF(time|timeInd1:timeInd2,{level|(/250,500,850/)},{lat|0:90},{lon|(/350,352.5,355,357.5,0,2.5,5,7.5,10,12.5,15,17.5,20/)})printVarSummary(hgtF_sub)end
方法2:直接350:20!!!!
begin
f = addfile("D:/cygwin/home/wunan/hgt.day.mean.200001-200004.nc","r")hgt = f->hgt
hgtF = short2flt(f->hgt)
time = f->timetime_1 = cd_calendar(time,2)
timeInd1 = ind(time_1.eq.20000101)
timeInd2 = ind(time_1.eq.20000107)
hgtF_sub = hgtF(time|timeInd1:timeInd2,{level|(/250,500,850/)},{lat|0:90},{lon|350:20})printVarSummary(hgtF_sub)end
结果
这篇关于(三)NCL读取经度:西经10度 (东经350度)-东经20度的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!