爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4027|回复: 2

[脚本编辑] grads 编程提取数据

[复制链接]

新浪微博达人勋

发表于 2022-5-14 20:54:29 | 显示全部楼层 |阅读模式

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

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

x
大家好,想请教一下,我提取某时刻的数据存为grd,但画出的图与nc数据的出图存在数值差异,能帮我看一下gs和ctl文件哪里出错了吗?谢谢

提取某时刻数据的gs文件:
'reinit'
var.1='u10';var.2='v10';var.3='d2m';var.4='t2m';var.5='msl';var.6='tp'
'set gxout fwrite'
'set fwrite D:\MERRA-2\huanliu\dimian\GD\grd\1990.grd'

'sdfopen G:\ERA5_HOUR\onsingle\2007.nc'
'set lon -180 180'
'set lat 0 90'
'set t 465'



'set z 1'
j=1
while(j<=6)

'd 'var.j''

j=j+1
endwhile
'close 1'
'disable fwrite'
;

grd的ctl文件:
dset D:\MERRA-2\huanliu\dimian\GD\grd\1990.grd
title huanliu hh
undef -9.99e+33
xdef 1440 linear -180 0.25
ydef 361 linear 0 0.25
zdef 1 linear 0 1
tdef 1 linear 00Z01JAN1990 1yr
vars 6
u10=>u10  0  t,y,x  10 metre U wind component
v10=>v10  0  t,y,x  10 metre V wind component
d2m=>d2m  0  t,y,x  2 metre dewpoint temperature
t2m=>t2m  0  t,y,x  2 metre temperature
msl=>msl  0  t,y,x  Mean sea level pressure
tp=>tp  0  t,y,x  Total precipitation
endvars

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

新浪微博达人勋

发表于 2022-5-15 17:08:19 | 显示全部楼层
本帖最后由 lqouc 于 2022-5-15 17:11 编辑

set lon和set lat那里建议改为set x 1 1440和set y 361 721 (此处数值或根据实际情况修改)
如果还是不行请附上nc数据的ctlinfo
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 成长值: 0
发表于 2022-5-17 02:33:53 | 显示全部楼层
差异多大
有没有可能NC转grd时候没有考虑scale factor 和offset?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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