对流不稳定差分算法（Qse/p）

278803532

前阵子用中央差分算了风切变，最近算了稳定度，就随便用前差、后差、中央差三种格式算了下对流稳定度，画出来的图如下，这个程序算出来是每一层的，画高度-时间剖面图的程序看我垂直风切变的那个帖子。不过建议，对流不稳定还是用上层减下层好用些，但是如果想要的画，可以借鉴一下，看程序吧
*============================================================
*对流性不稳定，Qse/P。分别用三种差分格式。
*上下边界用与最近那层次差分来代替
*------------------------------------------------------------
'reinit'
'open d:\zdz\qse.ctl'
'set gxout fwrite'
'set fwrite E:\fyb\qsep\qse-p.grd'

id = 1
while( id <=22)
  'set t 'id''

  'set lon 0 180'
  'set lat 0 80'
*==============================================================
*  First method----forward difference
*==============================================================
  'set lev 1000'
  'define vv1 = qse(lev=1000)'
  'define vv2 = qse(lev=975)'
  'display (vv2-vv1)/(975-1000)'

  'set lev 975'
  'define vv1 = vv2'
  'define vv2 = qse(lev=950)'
  'display (vv2-vv1)/(950-975)'

  'set lev 950'
  'define vv1 = vv2'
  'define vv2 = qse(lev=925)'
  'display (vv2-vv1)/(925-950)'

  'set lev 925'
  'define vv1 = vv2'
  'define vv2 = qse(lev=900)'
  'display (vv2-vv1)/(900-925)'

xlxh1023

占个沙发，谢谢楼主分享

qxtlyf

谢谢楼主分享

nunu18

哇~~~~楼主自己写的程序？太强啦，谢谢分享

zxc

非常感谢楼主！先收藏了！以后用的着。

longlong

谢谢，无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）

涛声依旧否

学习一下！

shayingying

这个很有用 不过想问一下 这样一层一层的算完以后 能不能把不同层次的差分值赋给一个变量呀？ 也就是说一个变量里面包含不同的层次

萌萌_pYGjF

维稳么呢也是醉了

依恋。

虽然用不上。。也是一种读数据的思路。棒