高空颠簸dutton指数的gs文件

saharast

'reinit'
'sdfopen d:/monthly/uwnd.mon.mean.nc'
'sdfopen d:/monthly/vwnd.mon.mean.nc'
'sdfopen d:/monthly/hgt.mon.mean.nc'
'set mproj scaled'
'set vpage 0 11 0 8.5'
'set parea 2 8 1 7'
'set mproj scaled'
'draw rec 0 0 10.9 8.49'
'draw rec 0 0 10.9 8.49'
'draw rec 2 1 8 7 '
'set lev 200'
'set lon 100 150'
'set lat 15 60'


i=433
m=1984

while (i<=793)
'set t 'i
m=1984+((i-433)/12)

'define a=uwnd(lev-50)-uwnd(lev+50)'
'define b=vwnd.2(lev-50)-vwnd.2(lev+50)'
'define c=mag(a,b)'
'define d=hgt.3(lev-50)-hgt.3(lev+50)'
'define e=(c/(0.001*d))'

'define dux=cdiff(uwnd,x)'
'define dvy=cdiff(vwnd.2,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'define p=dux/dx/637800'
'define q=dvy/dy/637800'
'define r=mag(p,q)*100000'
'define Dutton=(1.25*r+0.25*e*e+10.5)'

'set gxout shaded'
'set cint 5'
'set cmin 25'
'd Dutton'

'draw title 200.Dutton.'i'.'m''
'enable print d:/gradsst/200.Dutton.'i'.'m''
'print'
'disable print'
'c'
i=i+12
endwhile
;

wheelar

不错，先收下了，以前都是通过Fortran计算的

saharast

除dutton指数以外，其它颠簸预报指数如布朗指数、湍流指数TI1、湍流指数TI2等，你晓得它们判定颠簸的阈值吗？查了好多资料，没找到

syl136

你好，请问你也是研究航空气象方向的么？能不能交流一下

cms001

这是用的月平均资料吗，2.5*2.5是不是有点粗啊

awan

很强大，先收藏了，到时用到就来学习