爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9570|回复: 6

我想把nc数据中所需要的范围内的数据提取出来,代码有点问题

[复制链接]

新浪微博达人勋

发表于 2014-9-14 18:15:26 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 Rolle0524 于 2014-9-14 21:48 编辑

下载了gcm的nc格式的数据,不过这里面是所有经纬度范围内的数据,我只想提取出自己想要的那部分数据,并输出为txt,GrADS里面可以实现吗?我的文件路径是:E:\DATA\wvap-morton_Amon_CNRM-CM5_historical_r1i1p1_native.nc
我要的范围是经度在90.27°~91.05°,纬度在30.5°~30.92°范围内
所有时间段的数据都要。
(该文件一共有1872个time,月数据)
经度范围在nc数据中对应的是89.30~92.11
纬度范围在nc数据中对应的是 29.42~32.22                                    看了别人的帖子之后,我改了下别人的代码:sdfopen E:\DATA\wvap-morton_Amon_CNRM-CM5_historical_r1i1p1_native.nct=1
while(t<=1872)
set t t
set lon 89.30 92.11
set lat 29.42 32.22
set gxout fwrite
set fwrite E:\DATA\CNRM-CM5.dat
d uwndt=t+1
endwhile
disable fwrite
输入到while(t<=1872)这句的时候,提示unknown  command,为什么会这样?本人新手,不过急着用这个结果,跪求大神指导下代码。


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

新浪微博达人勋

发表于 2014-9-14 19:08:49 | 显示全部楼层
《=改成 .le.  ?同时新手,
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-9-14 20:13:33 | 显示全部楼层
循环变量用 i 吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-9-14 21:12:26 | 显示全部楼层
先看一下ctl文件,看一下变量名是不是一致的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-9-14 22:26:35 | 显示全部楼层
{:eb303:}{:eb303:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2014-9-14 22:36:27 | 显示全部楼层
可以使用兰溪的grads2ascii来实现
http://bbs.06climate.com/forum.p ... &extra=page%3D4
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-12-9 20:03:26 | 显示全部楼层
{:5_217:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

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

本版积分规则

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

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

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