利用NCEP资料画剖面图的一点儿经验

承德气象

river 发表于 2012-12-7 19:42

                               
登录/注册后可看大图

你不是都说出来了么。道理都是一样的，按照我那个gs，就是在出图的时候设置改一下就行啊。先固定经度
se ...

非常感谢！

river

承德气象 发表于 2012-12-8 21:23

                               
登录/注册后可看大图

非常感谢！

不客气，互相交流

随缘

*--------计算假相当位温----------
'define br=(0.1158*log(rhum.3)-1.5332787)*air.6-789.92816'                              
'define cr=789.92816*air.6'                                                         
'define tc2=(-br-sqrt(br*br-4*cr))/2.0'                                               
'define pc2=lev*pow(tc2/air.6,3.5005574)'                                                        
'define etc2=6.1078*pow(273.16/tc2,5.1347779)*exp(3148.0973*(tc2-273.16)/(126.06334*tc2))'   
'define wc2=0.622*etc2/(pc2-etc2)'                                                        
'define ed2=tc2*pow(1000.0/(pc2-etc2),0.28765924)'                                          
'define lc2=2500.79-2.3697*(tc2-273.16)'                                                  
'define eqt=ed2*exp(wc2*lc2/1.0048/tc2)'
楼主，这个thetase的计算公式是根据哪个来的？方便将计算公式附上么？计算结果楼主可否验证过？

river

随缘 发表于 2012-12-9 21:43

                               
登录/注册后可看大图

*--------计算假相当位温----------
'define br=(0.1158*log(rhum.3)-1.5332787)*air.6-789.92816'        ...

这个公式我记得是物理量诊断里有的，而且这是最常用的计算方法吧，和坛子里其他人分享的是一样的。这个计算完单位是开尔文，如果你理解成摄氏度的话肯定就觉得不对了

baggio

l楼主，太谢谢了！急需的东西啊！

蓝冰草

很有用，剖面图总是要画的，谢谢楼主分享

残阳飞血

好东西    多谢楼主   学习了

koala54

谢谢楼主分享，学到很多，赞一个！

gynxz

非常感谢！

流水2012

楼主，你那个多时次的固定纬度的剖面图，我要输出水汽通量散度
*---------输出图形------------
'set lat 15 35'
'set lon 113'
'set zlog on'
'set gxout shaded'
'set cint 0.3'
'set cmax 0'
'd qdiv'
'run cbarn.gs'
'set ylpos 0 r
'set gxout contour'
'set cint 0.3'
'set cthick 11'
'set clskip 1'
'd qdiv'
'writehz 1 8.1 沿东经113度剖面图 1 2.5 1 1.0 0.6 0 3'
'printim E:/201206/113E'i'.gif white'
*'print'
'c'
i=i+1
endwhile
*'disable print'
;
这样写怎么还是错误啊？