grads脚本画等熵位涡图总出错，哪位大神帮解决一下，给发红包

xingbin0804

'reinit''open d:\python\fnl_20190511_00_00.ctl''set mpdset cnworld''set lon 60 140''set lat 10 60''set lev 700 100'"define t=TMPprs""define u=UGRDprs""define v=VGRDprs""define pp=lev""define coriol=2*7.29e-5*sin(lat*3.1415/180)""define dudy=cdiff(u,y)/(111177*cdiff(lat,y))""define dvdx=cdiff(v,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))""define dt=t(z-1)*pow(1000/PP(z-1),0.286)-t(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,t,PP,320)"define part2="isen(dudy,t,PP,320)"define part3="isen(dtdp,t,PP,320)"define pv320=-9.8*(coriol+part1-part2)*part3""set z 1"'set lon 60 140''set lat 10 60'"d pv320*1000000"'printim d:/python/PVU20190511_00_00.jpg x1024 y768 white'function isen(field,tgrid,pgrid,tlev)
能解决的问题的 加微信给发红包
来自群组: 成于大气 信达天下

tulalang

你是不是直接在gseditor上跑的程序？？你可以试试打开grads后把脚本拖进grads窗口去跑跑

bincol

tulalang：你是不是直接在gseditor上跑的程序？？你可以试试打开grads后把脚本拖进grads窗口去跑跑
是的，，可以在opengrads中run 等熵位涡.gs就可以了！

mumululu

function isen(field,tgrid,pgrid,tlev)这句是什么意思啊

mumululu

提示isen这个公式不存在 这个要怎么解决？

zaqbob

请问楼主，问题解决了吗？

xingbin0804

没有人知道怎么改吗？

xingbin0804

求大神帮解决 加本人微信号Icewordforever

xingbin0804

求大神帮解决 加本人微信号Iceworldforever

xingbin0804

'reinit'
'open d:\python\fnl_20190511_00_00.ctl'
'set mpdset cnworld'
'set lon 60 140'
'set lat 10 60'
'set lev 700 100'
"define t=TMPprs"
"define u=UGRDprs"
"define v=VGRDprs"
"define pp=lev"
"define coriol=2*7.29e-5*sin(lat*3.1415/180)"
"define dudy=cdiff(u,y)/(111177*cdiff(lat,y))"
"define dvdx=cdiff(v,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))"
"define dt=t(z-1)*pow(1000/PP(z-1),0.286)-t(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,t,PP,320)
"define part2="isen(dudy,t,PP,320)
"define part3="isen(dtdp,t,PP,320)
"define pv320=-9.8*(coriol+part1-part2)*part3"
"set z 1"
'set lon 60 140'
'set lat 10 60'
"d pv320*1000000"

'printim d:/python/PVU20190511_00_00.jpg x1024 y768 white'

function isen(field,tgrid,pgrid,tlev)