【已解决】我编了一个K指数的程序，可是运行出来的图好奇怪，程序和图如下

莎朗嘿哟YY

生成的ctl:
dset E:/fnl/fnl_20120617_00_00_c
index E:/fnl/fnl_20120617_00_00_c.idx
undef 9.999E+20
title E:/fnl/fnl_20120617_00_00_c
*  produced by grib2ctl v0.9.12.5p33k
dtype grib 3
options yrev
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 1 linear 00Z17jun2012 1mo
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 121



配套的gs
'reinit'

'open e:/fnl/1.ctl'

'enable print e:\fnl\61700k.gmf'
'set lon 50 160'

'set lat 10 70'

'set lev 850'

'define T850=tmpprs-273.16'

'define RH=rhprs'

'define Td850=T850-((14.55+0.114*T850)*(1-0.01*rh) + pow((2.5+0.007*T850)*(1-0.01*rh),3) + (15.9+0.117*T850)*pow((1-0.01*rh),14))'

'set lev 700'

'define T700=tmpprs-273.16'

'define RH=rhprs'

'define Td700=T700-((14.55+0.114*T700)*(1-0.01*rh) + pow((2.5+0.007*T700)*(1-0.01*rh),3) + (15.9+0.117*T700)*pow((1-0.01*rh),14))'

'set lev 500'

'define T500=tmpprs-273.16'

'define K=T850-T500+Td850-(T700-Td700)'

'set cint 2.5'

'set cthick 6'

'd k'
'draw title K on 17,JUN 00'
'print'
'disable print'
;

出来的图如下：

请大神帮我看一下！谢谢！！！

kongfeng0824

请问ctl是你自己写的还是？

river

看着挺好的呀

传说中的谁

没记错的话这段gs应该是@传说中的谁 编的，楼主只是改了两三个参数

mofangbao

传说中的谁 发表于 2013-4-24 08:59

                               
登录/注册后可看大图

没记错的话这段gs应该是@传说中的谁 编的，楼主只是改了两三个参数

哈哈 这都能看出来

topmad

mofangbao 发表于 2013-4-24 09:08

                               
登录/注册后可看大图

哈哈 这都能看出来

经典的笨式风格...

莎朗嘿哟YY

kongfeng0824 发表于 2013-4-24 05:07

                               
登录/注册后可看大图

请问ctl是你自己写的还是？

ctl是 用grads生成的

莎朗嘿哟YY

river 发表于 2013-4-24 07:24

                               
登录/注册后可看大图

看着挺好的呀

前辈，这图没问题吗？好密集 啊

莎朗嘿哟YY

传说中的谁 发表于 2013-4-24 08:59

                               
登录/注册后可看大图

没记错的话这段gs应该是@传说中的谁 编的，楼主只是改了两三个参数

嗯嗯，是看到你的程序然后借用的~可是我的这个图出来好密集，你出来的图就好清爽~不知道能否告诉我我的程序里面是不是有什么地方要改，谢谢！！！！

莎朗嘿哟YY

topmad 发表于 2013-4-24 09:10

                               
登录/注册后可看大图

经典的笨式风格...

是挺笨的。。。。以前一直没接触过K指数。。。