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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10110|回复: 8

[脚本编辑] 新手求助!在画水汽通量散度场时等值线只显示一部分,这是为何?该如何处理?

[复制链接]

新浪微博达人勋

发表于 2021-9-27 20:23:21 | 显示全部楼层 |阅读模式

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

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

x
一、给出现有成果图。(450hpa水汽通量散度场)
二、给出我的gs程序。
'reinit'
'open E:\Master_thesis\GrADS\YFY_practice\llj\ctl\sandu\450hpa.ctl'
'set z 1'
'set t 1 12'
'set x 1 241'
'set y 1 161'
'define hdivg450=hdivg(q450*u450/9.8,q450*v450/9.8)*1e5'
'set t 1'
'set parea 1.5 10.0 1.5 8.0'
'set map 1 1 3'
'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 1'
'set ylint 1'
'set lon 81 98'
'set lat 27 32'
'set grads off'
'set grid off'
'set xlpos -10'
'set ylpos -10'
'set font 5'
'set gxout contour'
'set ccolor 1'
'set clopts 1 5 0.1'
'd hdivg450'
'run axis.gs -type b  -position o -label on -start 81 -end 98 -suffix `3. -interval 2  -size 0.2 -lsize 0.15 -lthick 0.3 '
'run axis.gs -type L  -position o -start 27 -end 32 -suffix `3. -interval 1 -lsize 0.15 -lthick 0.3'
'set line 2'
'draw shp E:\Master_thesis\GrADS\thesis\arcgis\yarlungfenqv'
'printim E:\Master_thesis\GrADS\thesis\png\sandu\450_1mo.png white'


三、设置的范围都对,水汽通量图也能正常出图,排查了几遍,但是还是不知道这个是什么问题引起的,请各位大神指点!感恩!
450_1mo.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-9-27 22:23:49 | 显示全部楼层
是不是原始资料的范围没这么大啊
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-9-28 08:44:49 | 显示全部楼层
计算的范围一定要比显示的范围大才行,你看下你的数据经纬度范围是多大
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-9-28 09:31:58 | 显示全部楼层
river 发表于 2021-9-27 22:23
是不是原始资料的范围没这么大啊

您好,我把我的过程和您说一下:
一、我先利用U、V和shum三个数据生成可识别的dat数据,gs程序如下:
'reinit'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\U_1979_2018.nc'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\V_1979_2018.nc'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\shum_1979_2018.nc'
'set gxout fwrite'
'set fwrite E:\Master_thesis\GrADS\thesis\dat\sandu\450hpa.dat'
'set lon 60 120'
'set lat 5 45'
'set z 1 27'
'set t 1 480'
'set lev 450'
'define u450=ave(u.1,t+0,t=480,12)'
'define v450=ave(v.2,t+0,t=480,12)'
'define q450=ave(q.3*1000,t+0,t=480,12)'
'set t 1 12'
'd u450'
'd v450'
'd q450'
'disable fwrite'
return

(注:我的原始资料范围确认为经度60~120,纬度5~45,如图所示)
二、为dat数据配备ctl
DSET  E:\Master_thesis\GrADS\YFY_practice\llj\dat\sandu\450hpa.dat
UNDEF  -9999.0
XDEF   241  LINEAR  60   0.25
YDEF   161  LINEAR  5   0.25
ZDEF   1    LINEAR  1    450
TDEF   12  linear 00z0jan1979 1mo
VARS  3
u450 0  99 u450 water all
v450 0  99 v450 water all
q450 0  99 q450 water all
ENDVARS
三、运行我提问时附上的gs程序
麻烦老师再帮我看看,是不是我哪里出错了?感谢!
U.png
shum.png
V.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-9-28 09:37:04 | 显示全部楼层
灵心流萤 发表于 2021-9-28 08:44
计算的范围一定要比显示的范围大才行,你看下你的数据经纬度范围是多大

我之前考虑过会不会是这个问题,但是我查了一下,我的原始资料范围确认为经度60~120,纬度5~45,如图所示,所以不知道是哪里出错了。
U.png
shum.png
V.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-9-28 15:03:24 | 显示全部楼层
我又重新梳理了一遍,不知道为何,图又能都显示出来了,但还是谢谢各位前辈。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-9-28 16:37:47 | 显示全部楼层
周一的07 发表于 2021-9-28 15:03
我又重新梳理了一遍,不知道为何,图又能都显示出来了,但还是谢谢各位前辈。

知道问题了,在写ctl的时候,我把数据的路径输错了。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-10-4 09:39:12 | 显示全部楼层
感谢分享!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2022-3-18 10:13:30 | 显示全部楼层
{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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