求助：怎样计算垂直变化？

彩虹依旧

也在学习中，同问？

dcx

karenlk 发表于 2015-6-5 10:51
cdiff函数用在z方向时不可以define，必须直接出图

'reinit'
'open E:/FNL/fnl_20140723_00_00.ctl'

'set parea 1 10 1 8'
'set grads off'
'set grid off'
'set parea 1 10 1 8'
'set csmooth on'
'set xlint 5'
'set ylint 5'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set lat 0 70'
'set lon 80 150'
'set z 1 21'
'set time 06Z24JUL2014'
'set map 1 1 2'
'set mpdset cnworld'
'define prs=lev'
'define ek=17.67*(tmpprs-273.16)/(tmpprs-29.65)'
'define es=6.112*exp(ek)'  
*比湿（克/千克）
'define qs=0.622*es/(prs-0.378*es)'
'define rh=rhprs*0.01*qs'
'define e=prs*rh/(0.62197+rh)+1e-10'
'define tk=55.0+2840.0/(3.5*log(tmpprs)-log(e)-4.805)'
'define pot=tmpprs*pow((1000/prs),(0.2854*(1.0-0.28*rh)))'
*位温
'set z 19'

'set lat 15 55'
'set lon 90 140'


'd cdiff(pot,z)'
'cbarn 1 0'
'printim fnl_weiwen.png x1024 y768 white'
请问你说的是这样么？提示错误：
error from sdiff:specified dimension non varying
operation error:error from cdiff function
应该怎么改啊大神

karenlk

dcx 发表于 2015-9-10 15:01
'reinit'
'open E:/FNL/fnl_20140723_00_00.ctl'

你是想说出什么样的图？是xy还是垂直剖面图

dcx

karenlk 发表于 2015-9-10 15:27
你是想说出什么样的图？是xy还是垂直剖面图

xy图，垂直稳定度的，d位温/dz

karenlk

dcx 发表于 2015-9-10 20:26
xy图，垂直稳定度的，d位温/dz

因为你要输出的xy的图，必须要z固定，而cdiff（pot，z）又需要z取一定的范围，这样就产生矛盾，不能使用这个函数。实际当中计算的时候可以采用上下层相减的方法，例如'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'然后再相除，得到结果

dcx

karenlk 发表于 2015-9-11 10:26
因为你要输出的xy的图，必须要z固定，而cdiff（pot，z）又需要z取一定的范围，这样就产生矛盾，不能使用 ...

大神，如果要画垂直剖面图怎么画啊，是画高度纬度剖面（固定经度）好还是画高度时间剖面比较好啊。按您的建议，我画了高度纬度图，但是不知道代表什么意思。高度时间图应该怎么画哪

                               
登录/注册后可看大图

karenlk

dcx 发表于 2015-9-11 12:15
大神，如果要画垂直剖面图怎么画啊，是画高度纬度剖面（固定经度）好还是画高度时间剖面比较好啊。按您的 ...

画垂直剖面图的时候你可以在计算的时候把z的范围扩大，然后出图的时候再截取你想要的范围，这样就不会出现你到100hPa的时候出现空白的情况了。
至于你想画什么样的剖面图，你可以根据你所研究的多看看文献的例子。资料都是xyzt四个变量，xy取区域平均，得到的就是zt的剖面。但是你算的时候就要把资料好好的处理。
我会的也不是很多，大家有问题互相切磋

dcx

karenlk 发表于 2015-9-11 15:05
画垂直剖面图的时候你可以在计算的时候把z的范围扩大，然后出图的时候再截取你想要的范围，这样就不会出 ...

好的，我再看一看，谢谢你

lnsqxt

也在学习中，同问？

那XX

Malik 发表于 2013-9-24 16:13
哦哦哦，因为直接弹出来了我不知道什么错所以没说清楚，不好意思

可能是中央差函数我不太会用，我改成 ...

能不能贴一下gs,谢谢