请教假相当位温图问题

tszhaoweiwen

之前出图一直有问题，在论坛高人指点下能画了假相当位温，但是和我参考文章图不一样，是不是假相当位温计算公式不一样？资料都是同时次1X1资料

我画图GS
'reinit'
'open C:\Users\sta288\Desktop\fnl328\fnl08.ctl'
'set t 3'
'set grads off'
'set grid off'
'set csmooth on'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set grads off'
'set lon 112.5'
'set lat 0 40'
'set lev 1000 300'
'set zlog on'
'set ylevs 1000 925 900 850 800 750 700 650 600 550 500 450 400 350 300 '
'define prs=lev'
'define es=(6.112*exp(17.67*(tmpprs-273.15)/(tmpprs-29.65)))'
'define q=RHprs*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tmpprs)-log(e)-4.805)'
'define theta=tmpprs*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set gxout shaded'
'set cthick 6'
'set clskip 1'
'set clevs 320 330 340 350 360'
'd eqt'
'run cbarn.gs 1 0'
'set ylpos 0 r'
'set gxout contour'
'set cthick 11'
'set cint 4'
'd eqt'
'set gxout contour'
'set clevs 302 304 306 308 310 312'
'd eqt'
'writehz 1 8.1 假相当位温 1 2.5 1 1.3 0.6 0 7'
'printim D:\eqt2.png white'


下面是图对比
这是文献的图（图片中明显有个302~310低值区）
这是我画图
是不是文献图假相当位温计算公式更精细？希望有人教教

river

计算假相当位温有很多经验公式，不同的公式结果就会有不同。但是整体形势应该差不多的，你的这个图和文献里的整体形势是一致的。要想了解更多的假相当位温的计算的话，可以看一些文献，有很多文献就是研究假相当位温的计算的。

tszhaoweiwen

river 发表于 2013-6-25 18:26

                               
登录/注册后可看大图

计算假相当位温有很多经验公式，不同的公式结果就会有不同。但是整体形势应该差不多的，你的这个图和文献里 ...

您好，多谢您回答，我们家园论坛有无类似的文献或者研究啊？我搜索了下 没有发现呢

river

tszhaoweiwen 发表于 2013-6-25 19:23

                               
登录/注册后可看大图

您好，多谢您回答，我们家园论坛有无类似的文献或者研究啊？我搜索了下 没有发现呢

我说的文献其实就是有发表过的相似的文章，你应该去知网什么的地方搜索

liutaoz21

好长。。。帮顶一下~改天再来学习。

edwardli

你确定文献里的经度是112.5？  你的cint是4，文献的是2。

tszhaoweiwen

edwardli 发表于 2013-6-26 07:51

                               
登录/注册后可看大图

你确定文献里的经度是112.5？  你的cint是4，文献的是2。

您好，文献前一副图是112.5  是说广州的 所以应该没有错 se cint 2 后都没有出现低值

edwardli

tszhaoweiwen 发表于 2013-6-26 08:38

                               
登录/注册后可看大图

您好，文献前一副图是112.5  是说广州的 所以应该没有错 se cint 2 后都没有出现低值

那就你的数据有问题了，t=3？

mxp1

找了这么久还是不知道用哪个

wu3015869

学习，学习！！！