救命！！温度平流显示error from cdiff

弄晴

计算温度平流显示cdiff错误
求助是什么原因
报错截图已上传！
gs文件如下：
'sdfopen f:\shuj\2021040105.nc'
'set t 1'
'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 6'
'set parea 0.9 9.8 1 8'
'set clopts -1 -1 0.12'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.16'
'set font 0'
'set z 1 9'
'set lon 81.53'
'set lat 41.99'
'set zlog on'
'set ylevs 1000 900 850 700 600 500 300 200'
'set xlabs 08|20|08|20|08|20|08|20|08|20'

****************************** 温度平流 ********************************
'define pi=3.1415726'
'define r=6.356e6'
'define dtx=cdiff(t,x)'
'define dty=cdiff(t,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi/180'
'define dy=cdiff(lat,y)*pi'
'define advT=-(u*dtx/dx+v*dty/dy)/r'

*设置出图
'set gxout shaded'
'set gxout contour'
'run C:\OpenGrADS\colors\10colors.gs'
'set clevs  2 4 6 8 10'
'set ccols   0  30 31 32 33 34 35 '
'set black -0.001 0.001'
'd smth9(advT*100000)'
'cbarn 1.3 1'

'gxprint F:\shuj\Tse-4.png x1500 y1000 white'
*'gxprint F:\shuj\Tse-4.eps x1500 y1000 white'
'c'
;

river

这不是显而易见的问题吗，你一上来啥都固定好了他怎么计算？给你一个点你能弄出一个面吗？

弄晴

river 发表于 2022-10-9 12:21
这不是显而易见的问题吗，你一上来啥都固定好了他怎么计算？给你一个点你能弄出一个面吗？

谢谢大佬，还有一个疑问，就是温度变量定义的是t，会和时间t计算温度平流的时候混淆吗

wode_q_x

楼主解决了吗？同样的问题

river

弄晴 发表于 2022-10-9 14:54
谢谢大佬，还有一个疑问，就是温度变量定义的是t，会和时间t计算温度平流的时候混淆吗

这种情况下应该不会，不过以防万一最好不要使用GrADS默认的来定义变量。你的NC资料是自己输出的吗，变量名怎么会使 t 呢

bbhofseven

{:5_275:}{:5_275:}