任意两点间斜剖面做法总结【利用GRADS自带函数】

river

liuyw 发表于 2017-9-10 18:21
嗯嗯，我采用的方法是先算出来假相当位温，生产一个dat文件，然后再调用它，利用您的那个单一变量的脚本 ...

没问题，这样是可以的。我觉得也可以尝试在一个脚本里完成。

liuyw

river 发表于 2017-9-10 19:30
没问题，这样是可以的。我觉得也可以尝试在一个脚本里完成。

一开始也有这个打算的，只是没想好怎么放在一块

river

liuyw 发表于 2017-9-10 22:05
一开始也有这个打算的，只是没想好怎么放在一块

先解决燃眉之急，后面回头再试好了

micolian

谢谢楼主分享！有个问题，就是在GS中有个-U，这个之前没有定义的啊

river

micolian 发表于 2017-9-11 14:15
谢谢楼主分享！有个问题，就是在GS中有个-U，这个之前没有定义的啊

这个是函数里的参数，不明白的话去看coll2gr 的用法

micolian

river 发表于 2017-9-11 16:08
这个是函数里的参数，不明白的话去看coll2gr 的用法

你好，我通过你的GS，能画出图形，但是有一些类似x axis labels overridden by set xaxis提示。但是，画水汽通量的图，特别是加了地形处理的，结果画不出来图形

micolian

river 发表于 2017-9-11 16:08
这个是函数里的参数，不明白的话去看coll2gr 的用法

可以加个QQ或微信吗？我的QQ42396561，方便的话留下QQ呀

river

micolian 发表于 2017-9-11 18:54
可以加个QQ或微信吗？我的QQ42396561，方便的话留下QQ呀

coll2gr()

This function creates a grid from a collection of station data. The syntax is:

coll2gr(cnum <,num>)

cnum     collection number
num      number of vertical levels in the output grid (default is 10)
If num is given as "-u", then the vertical levels will be a union of all the levels in the collection of station data.

Usage Notes

Before using coll2gr, you must create the collection of stations using the GrADS command collect.
coll2gr does not yet support time slices; currently, it will only work when the collection of stations is a collection of vertical profiles.
coll2gr produces an output grid that varies in X and Z; the dimension environment used when coll2gr is invoked must also be X and Z varying. The X axis of the output grid will contain the equally spaced station profiles and will span the range of the current X dimension environment. The Z axis of the output grid will span the range of the current Z dimension environment and will have either the specified number of levels or a union of the levels. Data points outside of the range of levels will be used for interpolating to within the range if appropriate.
The "-u" option will only work for data in pressure vertical coordinates. If your vertical coordinate is height, num should be the actual number of vertical levels in the collection of profiles. Note: This restriction has been removed as of version 2.1.0.
See the section of the User's Guide on Arbitrary Cross Sections for more information.

micolian

river 发表于 2017-9-11 23:07
coll2gr()

This function creates a grid from a collection of station data. The syntax is:

关于命令语句，我查了学习手册，已经知道用法。问题是这样的，我用的是NCEP1*1的资料，在用你第一个GS文件可以画出垂直速度，但第二个GS，特别是加了地形后，好多错误提示，所以想请教下，我的问题出在了哪里？

river

micolian 发表于 2017-9-12 09:58
关于命令语句，我查了学习手册，已经知道用法。问题是这样的，我用的是NCEP1*1的资料，在用你第一个GS文 ...

那个地形数据的层次和FNL资料的层次不一样，所以会报错。但是如果是正常的经纬度的剖面是可以正常使用的。我建议你把地形数据插值到和FNL资料一样的层次上再用吧。但是目前我没有发现很好地插值方法