湿位涡剖面图 画图出现问题 求教 用FNL资料

heyufei

GS脚本如下 画图边界层有问题
’reinit'
'open E:\aa\fnl\fnl_20130706_06_00.ctl'
'enable print E:\aa\0606.jpg'
'set z 2 20'
'set zlog on'
'set lon 100 130'
'set lat 20 45'
'define tc=tmpprs-273.15'
'define rh=rhprs'
'define prs=lev'  
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=rh*qs/100'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tc+273.16)-log(e)-4.805)'
'define theta=(tc+273.16)*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

'define u=UGRDprs'
'define v=VGRDprs'
'define vo=hcurl(u,v)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=u(z-1)-u(z+1)'
'define dv=u(z-1)-u(z+1)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp'  
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'

'set lon 114'
'set lat 20 45'
'set grid off'
'set grads off'
'set gxout contour'
'set csmooth on'
'set cthick 7'
'set cint 0.2'
*'set clopts -1 -1 0.10'
'd pv1*10e5'
*'d pv2*10e6'
*'d pv*10e5'
'print'
'printim E:\aa\0606.jpg white'
'disable print'

river

大姐，先把你的经纬度范围设置的大点儿

nunu18

这个我之前也遇到过，主要是你的经纬度设置小啦，你先把范围扩大运行'd pv1*10e5'，执行c命令，再缩小到你要的范围，再'd pv1*10e5'，就可以啦

写论文的兔子

受教了~

mengnannan1314

我也出现这种问题，也按下面的方法把经纬度设置大了，可是没变化呀？是什么情况呢

天边的彩虹

请问，顶部出现空白，怎么设置啊

clare

天边的彩虹 发表于 2015-6-23 21:34
请问，顶部出现空白，怎么设置啊

同理，计算的时候所有维度均设置大一些；画图时则缩小到所需范围。

fujianbeige

受教了，谢谢

oydl1985

1. 'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
   
2. 'define dy=2.0*6370949.0*3.14159/180.0'
   

复制代码

请问下楼主这里为什么要乘2呀？有的帖子没有乘2，不知道相信哪个

js250768097

谢谢这个对我帮助很大