爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 13115|回复: 14

[分享资料] grads画剖面图skip的使用疑问

[复制链接]

新浪微博达人勋

发表于 2013-6-26 09:23:54 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 寄居蟹 于 2013-6-26 14:57 编辑

用GRADS画剖面图,画垂直速度时,使用语句d u;w时,风速比较密。但是使用skip(u,4);w时,水平的风速没那么密了,但是垂直层高度上也变稀疏了。如下图
我希望既能在水平方向上不那么密,垂直方向上也不稀疏。怎样能只改变水平方向的疏密,同时垂直方向的不改变

求教各位这种应该怎么改?谢谢!!感激不尽!!

d skip(u,4);w*10

d skip(u,4);w*10

d u;w*10

d u;w*10
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-6-26 09:42:14 | 显示全部楼层
你用fortran重新写一下不就好了。
或者skip(u,4,1)w
总之去看看手册好了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-6-26 12:26:42 | 显示全部楼层
把你ctl贴出来看一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-6-26 14:57:44 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-6-27 16:55:47 | 显示全部楼层
基础知识嘛,skip函数分别对u 和w控制不就行了么。
skip()

skip (expr, skipi, skipj)

Sets alternating values of expr to the missing data value.

expr     A valid grid expression that may have 1 or 2 varying dimensions
skipi    Skip factor in the I dimension of expr
skipj    Skip factor in the J dimension of expr
Usage Notes

This function is often used while displaying wind arrows or barbs to thin the number of arrows or barbs. It is not necessary to use the skip function on both the U and V wind components; it is sufficient to populate only one component with missing data values to suppress the plotting of the wind arrow or barb.
Examples

Suppose you have a time series of 3-hourly data, but you want to display values at 6-hourly time steps:

set x 1
set y 1
set t 1 last
d skip(var,2)

To display every other grid point in both the X and Y direction:

d skip(u,2,2);v

To display every grid point in the Y direction, but every 5th grid point in the X direction:

d skip(u,5,1);v

This example script "d_uv.gs" written by Wesley Ebisuzaki automatically sets the skip factor based on the plot dimensions.
*
* This function does a d skip(ugrd,n);v
* where n is automatically set to an appropriate value
*
* usage: d_uv ugrd vgrd
*
* v1.1 w. ebisuzaki
* v1.2 4/6/98 revised empirical formula for skip
*
function duv(arg)
u = subwrd(arg,1)
v = subwrd(arg,2)

* get lat/lon info
'query dims'
lons = sublin(result,2)
lats = sublin(result,3)
dx = subwrd(lons,13) - subwrd(lons,11)
dy = subwrd(lats,13) - subwrd(lats,11)

* Determine skip factor
dn = dx
if (dy > dx) ; dn = dy ; endif
skip = dn / 50 + 0.5
if (skip < 1) ; skip=1 ; endif

* Display the plot
'd skip('u','skip');'v
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-3-29 21:27:04 | 显示全部楼层
感谢分享,一样的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-3-29 21:27:11 | 显示全部楼层
感谢分享,一样的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-3-30 19:20:27 | 显示全部楼层
楼主风场的剖面图怎么画
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-3-30 19:20:52 | 显示全部楼层
是不是要先做纬向平均啊?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-27 21:58:56 | 显示全部楼层
你的图为什么往右偏
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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