风垂直切变（V/Z）中央差分格式算法

278803532 · 发表于 2014-5-19 22:18:54

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

前两天本来做风垂直切变的，搜了一下论坛，发现大家都想着用CDIFF函数来做，因为要用到位势高度，本人就自己动手一层层算了，
说明一下，*============================================================
*风垂直切变程序，v/z，由于位势高度和几何高度差别很小，用位势高度
*代替几何高度。分别用三种差分格式。
*上下边界用与最近那层次差分来代替
*==============================================================
'set lev 1000'
  'define vv1 = UGRDprs(lev=1000)'
  'define vv2 = UGRDprs(lev=975)'
  'define hh1 = HGTprs(lev=1000)'
  'define hh2 = HGTprs(lev=975)'
  'display (vv2-vv1)/(hh2-hh1)'

  'set lev 975'
  'define vv3 = UGRDprs(lev=950)'
  'define hh3 = HGTprs(lev=950)'
  'display (vv3-vv1)/(hh3-hh1)'

  'set lev 950'
  'define vv1 = vv2'
  'define vv2 = vv3'
  'define hh1 = hh2'
  'define hh2 = hh3'
  'define vv3 = UGRDprs(lev=925)'
  'define hh3 = HGTprs(lev=925)'
  'display (vv3-vv1)/(hh3-hh1)'

发三个层次上来看看吧，因为我不知道怎么发那个可以滚动的让你们粘贴，所以还是看附件，附件里面lon7.gs，7.gs,7.ct这个三个程序是用来高度-时间垂直剖面图的，也算一点福利吧。大家看图了

Mon.Chichi · 发表于 2014-5-19 22:24:01

学习学习，留着以后有用！

278803532 · 发表于 2014-5-19 22:33:56

忘记说明了，本来想用三种差分格式的，后面我就只用了中央差分的哈，前差和后差有人写了的话，可以讨论讨论

season100 · 发表于 2014-5-20 00:03:02

本帖最后由 season100 于 2014-5-20 00:08 编辑

Fnl数据，不用Z循环么 set z' zz

tszhaoweiwen · 发表于 2014-5-20 08:25:54

很好的学习资料

wangjb · 发表于 2014-5-20 08:54:07

下载学习了,谢谢楼主!

278803532 · 发表于 2014-5-20 10:13:12

season100 发表于 2014-5-20 00:03
Fnl数据，不用Z循环么 set z' zz

我没有用z循环，有兴趣可以试试的，因为用循环处理边界有些麻烦，但是感觉的也行的，

xinhu123456 · 发表于 2014-5-20 16:08:41

很好的东西,正好用上.

kumaei · 发表于 2014-5-20 16:14:58

捧场

green · 发表于 2014-5-22 20:05:43

收藏下来，以后学习，谢谢楼主！！！

		自动登录	找回密码
密码			立即注册

[脚本编辑] 风垂直切变（V/Z）中央差分格式算法

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

评分