等熵位涡 求指导错误

平流层的萝卜

Joucly 发表于 2013-10-22 21:32

                               
登录/注册后可看大图

set lev 的问题···  grads里面不能经度  维度  高度  都是变化的

哦？是这个问题吗？我写的也是三个都变化，但是也没什么问题啊，只要最后d的时候设置一个z就好。。。。{:eb512:}

Joucly

平流层的萝卜 发表于 2013-10-22 21:37

                               
登录/注册后可看大图

哦？是这个问题吗？我写的也是三个都变化，但是也没什么问题啊，只要最后d的时候设置一个z就好。 ...

我是这样改的呢···

我把层次分开一个一个计算···  图就出来了

平流层的萝卜

Joucly 发表于 2013-10-22 21:32

                               
登录/注册后可看大图

set lev 的问题···  grads里面不能经度  维度  高度  都是变化的

应该是set lev 1000 150时，z=1的时候是1000，这个时候dp就没法算了因为z-1=0不存在，是这个问题吧？楼主可以把最底层设成975试试

ayzqs

呼呼，我去试一试

Joucly

Joucly 发表于 2013-10-22 21:39

                               
登录/注册后可看大图

我是这样改的呢···

我把层次分开一个一个计算···  图就出来了

恩· ··我把每一层单独分开时·· 没有求1000  是从975开始

平流层的萝卜

同学，为什么我的isen.gs运行不了？

Joucly

平流层的萝卜 发表于 2013-10-30 23:11

                               
登录/注册后可看大图

同学，为什么我的isen.gs运行不了？

什么意思？
你像我一样复制在后面 试试呢

平流层的萝卜

Joucly 发表于 2013-10-31 09:58

                               
登录/注册后可看大图

什么意思？
你像我一样复制在后面 试试呢

试了几次，还是不行。

脚本也是照着isen.gs下边的用法来的，没错吧：
'open h:\paper\fnl.ctl'
'enable print H:\paper\mpv\pv_function_test.gmf'
i=1
while(i<=13)
     'set t 'i''
    'set lat 15 55'
    'set lon 70 150'
    "set lev 950 150"
   'define pp=lev'
    "define coriol=2*7.29e-5*sin(lat*3.1415/180)"
    "define dudy=cdiff(ugrdprs,y)/(111177*cdiff(lat,y))"
    "define dvdx=cdiff(vgrdprs,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))"
    "define dt=tmpprs(z-1)*pow(1000/PP(z-1),0.286)-tmpprs(z+1)*pow(1000/PP(z+1),0.286)"
    "define dp=100*(PP(z-1)-PP(z+1))"
    "define dtdp=dt/dp"
    'define part1=isen(dvdx,tmpprs,PP,320)'
    'define part2=isen(dudy,tmpprs,PP,320)'
    'define part3=isen(dtdp,tmpprs,PP,320)'
    "define pv320=-9.8*(coriol+part1-part2)*part3"
    "set z 1"
    "d pv320"
    'print'
    'c'
i=i+1
endwhile
'disable print'
;script文件夹里明显有这个isen.gs文件，为什么用不了呢？其他的诸如cbarn的gs都是可以直接用的，为什么isen.gs会出问题。

Joucly

是 'define part1’=isen(dvdx,tmpprs,PP,320)

平流层的萝卜

Joucly 发表于 2013-11-3 10:05

                               
登录/注册后可看大图

是 'define part1’=isen(dvdx,tmpprs,PP,320)

还是一样的