请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5033|回复: 1

[图形美化] 地形图叠加出现错误

[复制链接]

新浪微博达人勋

发表于 2020-5-5 16:29:33 | 显示全部楼层 |阅读模式

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

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

x
以下为中国地形的ctl:
dset D:\ncep\dx.nc
title
undef -9.99e+33
dtype netcdf
xdef 2011 levels 70 70.0333 70.0667 70.1 70.1333等

ydef 1201 levels 15 15.0333 15.0667 15.1等
zdef 1 linear 0 1
tdef 1 linear 00Z01JAN0001 1mn
vars 1
z=>z  0  y,x  z
endvars


地形的gs:
'reinit'
'reset'
'sdfopen d:\ncep\dx.nc'
'set grads off '
'set grid off'
'set lon 73 137'
'set lat 15 55'
define_colors
'set gxout shaded'
'set clevs 3000'
'set ccols 0 16 '
'd z'
'draw title (a)'
'printim d:\ncep\dx.png x1200 y1000 white'
'disable print'

能出图

热源的ctl:
dset d:\ncep\q1pj.dat
undef  -9.99e+33
xdef 144 linear   0 2.5
ydef 69 linear  -85 2.5
zdef 1 linear 1000   1000
tdef 12 linear 0Z1mar1979 1mo
vars 1
q1n  0  99 *Monthly Mean of Q1
endvars


热源的gs:
'reinit'
'reset'
'open d:\ncep\q1pj.ctl'
'set grads off '
'set grid off'
'set lon 73 137'
'set lat 15 55'
'set z 1'
'set gxout contour'
'set clevs  -60  -40  -20  0  20  40  60  80  100  120'
'set ccolor 1 '
'set cthick 4'
'd ave(q1n,t=1,t=3)'
'draw title (a)'
'printim d:\ncep\q1.png x1200 y1000 white'
'disable print'

能出图

两个数据合到一起后的gs:
'reinit'
'reset'
'set grads off '
'set grid off'
'sdfopen d:\ncep\dx.nc'
'open d:\ncep\q1pj.ctl'
'set lon 73 137'
'set lat 15 55'
'set z 1'
define_colors
'set gxout shaded'
'set clevs 3000'
'set ccols 0 16 '
'd z.1'
'set gxout contour'
'set clevs  -60  -40  -20  0  20  40  60  80  100  120'
'set ccolor 1 '
'set cthick 4'
'd ave(q1n.2,t=1,t=3)'
'draw title (a)'
'printim d:\ncep\z.png x1200 y1000 white'
'disable print'

grads报错:
Data Request Error: Invalid grid coordinates
  World coordinates convert to non-integer  grid coordinates
    Variable = q1n.2  Dimension = 3
Operation Error:  Error from ave function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = ave(q1n.2,t=1,t=3)

请问这个gs应该怎么改?因为地形本来就好像只有两维,没有高度层和时间层,而热源有一层高度和12个时间层?有没有什么办法将等值线直接叠加到阴影图上?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-5-5 17:13:49 | 显示全部楼层
是是这样可以不'd z.1(t=1)'
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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