爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4125|回复: 3

[分享资料] 【求助】用Grads将nc格式转dat,绘图出错~~

[复制链接]

新浪微博达人勋

发表于 2015-1-10 21:10:35 | 显示全部楼层 |阅读模式

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

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

x
一堆soda资料,cdf格式,命名有规律:
SODA_2.2.4_187101.cdf
SODA_2.2.4_187102.cdf
SODA_2.2.4_187103.cdf
……
SODA_2.2.4_200812.cdf

用GrADS将1950年1月~2008年12月(共708个月),区域是120°E~80°W,30°S~30°N的temp,taux两个变量提取出来放到一个.dat文件中
提取数据所用.gs文件如下:
'reinit'
'set gxout fwrite'
'set fwrite I:\SODA2.2.4\pic\SODA_2.2.4.dat'
yy=1950
while(yy<=2008)
  mm=1
  while(mm<=12)
  if(mm<=9)
  'sdfopen I:\SODA2.2.4\SODA_2.2.4_'yy'0'mm'.cdf'
  else
  'sdfopen I:\SODA2.2.4\SODA_2.2.4_'yy''mm'.cdf'
  endif
  'set lon 120 280'
  'set lat -30 30'
  'set z 40'
  'set t 1'
  'd temp'
  'd taux'
  'd tauy'
  'close 1'
  mm=mm+1
  endwhile
yy=yy+1
endwhile
'disable fwrite'
读出来后,写.ctl文件是:
dset     I:\SODA2.2.4\SODA_2.2.4.dat
undef    -9.99e+08
xdef     321 linear 120 0.5
ydef     121 linear -30 0.5
zdef     1   levels 0 1
tdef     708 linear jan1950 1mo
vars     3
temp 1  99  sst
taux 1  99  taux wind
endvars
画出来的图却是这样的:
在论坛上也搜索过很多大神前辈们的指导文章,受益不少,但实在是解决不动了,这是数据提取的时候出现问题了么?求赐教~
pivvv.jpg
ctl.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-10 21:11:16 | 显示全部楼层
选择“脚本编辑”项无法提交帖子。。。只好选了“分享资料”,请见谅~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-10 21:59:27 | 显示全部楼层
把你的经纬度设置改成用x y 的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-13 23:20:26 | 显示全部楼层
变量不是2个么 vars怎么是3
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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