爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2836|回复: 4

[分享资料] grads提取nc的dat文件画不出图来

[复制链接]

新浪微博达人勋

发表于 2014-4-1 17:20:42 | 显示全部楼层 |阅读模式

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

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

x
资料是ncep的0.5°x0.5°的,每天每小时(hourly)的地面气温,一共有一个月的,现在我要把云南地区的提取出来,用grads的fwrite,可是画了图之后不知道是哪里错了,图一眼看上去就是错的。
这是gs文件
'reinit'
'set gxout fwrite'
'set fwrite e:\lunwen\data\07t2m.dat'
'sdfopen e:\lunwen\data\NRA_T2M_1979-07.nc'
'set lon  97 107'
'set lat 20 30'
t=1
while(t<=744)
'set t 't
'd t2m'
t=t+1
endwhile
'disable fwrite'
'reinit'


这是dat的ctl
dset e:\lunwen\data\07t2m.dat
title yunnan temperature
undef -9.99e+33
xdef 21 linear 97 0.5
ydef 21 linear 20 0.5
zdef 1 levels 1000
tdef 744 linear 01jun1979 1dy
vars 1
t2m 0 -999 2meter air temperation
endvars

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

新浪微博达人勋

发表于 2014-4-1 18:28:29 | 显示全部楼层
tt=1试试。。。。。
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2014-4-1 22:24:39 | 显示全部楼层
改成这个试试
mt=1
while(mt<=744)
'set t 'mt''
'd t2m'
mt=mt+1
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-4-2 21:29:01 | 显示全部楼层
胖胖 发表于 2014-4-1 22:24
改成这个试试
mt=1
while(mt

好的。解决了,谢谢你啊~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-28 23:20:20 | 显示全部楼层
过来学习中
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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