求助大气可降水量变量不显示

houzhaocan

用的欧洲中心era5的资料，积分比湿求大气可降水量，显示不出来，求助
我的GS文件：
'reinit'
'sdfopen d:\1921\E5\pw.nc'
'set lon 90 150'
'set lat 0 60'
'set t 1'
'set mpdset cnworld'
'set map 1 1 4'
'set parea 1 10.2 1.5 8'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set ylint 5'
'set xlint 5'
'set grads off'
'set grid off'
'define w=q/9.8'
'define pw=vint(lev(z=1),w,200)/1000'
'set gxout shaded'
'set lon 100 135'
'set lat 15 45'
'set cint 5'
'set ccolor rainbow'
'd pw'
'run C:\OpenGrADS-2.2\Contents\Resources\Scripts\cbarm_wql.gs 1.0 1'
'set gxout vector'
'set lon 100 135'
'set lat 15 45'
'set z 7'
'set arrscl 0.1 10 '
'set arrowhead 0.06 '
'set ccolor 1'
'd skip(u*2.5,5,5);skip(v*2.5,5,5)'
'run C:\OpenGrADS-2.2\Contents\Resources\Scripts\drawfigstr.gs  TopLeft (a) 3'
'printim d:\1921\08时大气可降水量及850风场.png white'
;

river

houzhaocan 发表于 2022-4-12 23:27
抱歉，不好意思，之前运行都没有报错，但是出的图不对，所以不知道怎么回事。改过之后出现这样的报错，警 ...

你的地面气压就一层数据吧？再加一个限定条件，把sp.2 改成 sp.2(z=1) 试试

river

我觉得是你地面气压的问题，lev(z=1) 这不就是固定值1000吗？而函数要求的是单位为mb的地面气压值，地面格点上的气压不可能是恒定的。你看下官网说明http://cola.gmu.edu/grads/gadoc/gadocindex.html。
还有一点是vint算法里应该已经 /9.8 了，所以我觉得是不是不需要再除了

houzhaocan

river 发表于 2022-4-10 22:32
我觉得是你地面气压的问题，lev(z=1) 这不就是固定值1000吗？而函数要求的是单位为mb的地面气压值，地面格 ...

修改过后，显示的数据还是错误的
'reinit'
'sdfopen d:\1921\E5\pw.nc'
'sdfopen d:\1921\E5\pres.nc'
'set lon 100 135'
'set lat 15 45'
'set t 1'
'set mpdset cnworld'
'set map 1 1 4'
'set parea 1 10.2 1 8.4'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set ylint 5'
'set xlint 5'
'set grads off'
'set grid off'
'set gxout shaded'
'define pw=vint(sp/100,q,200)/1000'
'set rgb 16 255  255  255'
'set rgb 17 247  252  240'
'set rgb 18 224  243  219'
'set rgb 19 204  235  197'
'set rgb 20 168  221  181'
'set rgb 21 123  204  196'
'set rgb 22  78  179  211'
'set rgb 23  43  140  190'
'set rgb 24   8  104  172'
'set rgb 25   8   64  129'
'set rgb 26   0   32   62'
'set rbcols 16 17 18 19 20 21 22 23 24 25 26 '
*'set clevs   0 5 10 15 20 25 30 35'
*'set cint 5'
'd pw'
'run C:\OpenGrADS-2.2\Contents\Resources\Scripts\cbarm_wql.gs 1.0 1'
'set gxout vector'
'set z 7'
'set arrscl 0.1 10 '
'set arrowhead 0.06 '
'set ccolor 1'
'd skip(u*2.5,5,5);skip(v*2.5,5,5)'
'run C:\OpenGrADS-2.2\Contents\Resources\Scripts\drawfigstr.gs  TopLeft (b) 3'
'printim d:\1921\08时大气可降水量及850风场.png white'
;

river

houzhaocan 发表于 2022-4-11 08:56
修改过后，显示的数据还是错误的
'reinit'
'sdfopen d:\1921\E5\pw.nc'

你这个脚本竟然没有报错？？？

sp这个变量是哪个资料的？pw.nc里面有这个变量？

q是什么单位，你计算出来的可降水量单位是mm还是kg?

houzhaocan

river 发表于 2022-4-11 13:15
你这个脚本竟然没有报错？？？

sp这个变量是哪个资料的？pw.nc里面有这个变量？

sp是我下载的地面表面压力（Surface pressure），是pres.nc文件的，q是比湿，kg/kg，我想计算mm

houzhaocan

river 发表于 2022-4-11 13:15
你这个脚本竟然没有报错？？？

sp这个变量是哪个资料的？pw.nc里面有这个变量？

函数不是要求单位为mb的地面气压值吗？我就下载了地面的表面压力数据，但是感觉不太会

river

houzhaocan 发表于 2022-4-12 19:24
函数不是要求单位为mb的地面气压值吗？我就下载了地面的表面压力数据，但是感觉不太会

肯定是不太会有问题才来提问的对吧。但是你提问要抓住要害，不是说随便弄个脚本往这儿一放等着别人给你挑错啊。
既然是个脚本程序，运行起来的错误提示呢？从错误提示一步一步来解决才是正确思路啊，一句不会就能都推给别人吗······
你的地面气压资料是第二个打开的文件，所以你要用角标告诉GrADS，我sp这个数据在第二个文件里，所以你脚本里的sp都要写成 sp.2

houzhaocan

river 发表于 2022-4-12 22:33
肯定是不太会有问题才来提问的对吧。但是你提问要抓住要害，不是说随便弄个脚本往这儿一放等着别人给你挑 ...

抱歉，不好意思，之前运行都没有报错，但是出的图不对，所以不知道怎么回事。改过之后出现这样的报错，警告超出文件限制

栤唥de嶶笑

参考一下我这个帖子，里面给出了vint函数的详细解析。
http://bbs.06climate.com/forum.p ... &extra=page%3D1