位涡诊断

杨雨轩

我想画700hPa高度上，干位涡的纬度—时间剖面图（经度固定在东经116度），所用的资料是ECMWF的逐日再分析（0.5°×0.5°）。尝试了很多次，但总是报错：

希望各位能帮我找出错误。
gs文件如下：
'reinit'
file='f:\yangyuxuan\theory\201512\mpv\pv_lat_t_0.5.gmf'
'enable print 'file
'set vpage 0 11 0 8.5'
'set parea 2 7 2 7'
'set mpdset hires'
'set map 15 1 8'

'sdfopen f:\yangyuxuan\theory\201512\201512.nc'
'set grads off'
'set grid off'

'set t 1 60'
'set lat 0 60'
'set lon 110 120'
'set z 1 11'

if(t>258.0)
'define a=17.2694'
'define b=35.86'
else
'define a=21.8764'
'define b=7.66'
endif

'define es=6.112*exp(a*(t-273.15)/(t-b))'
'define rs=0.622*es/(lev-es)'
'define theta=t*pow(1000/(lev-es),(287/1004))'
'define lv=2.5e6-2323*(t-273.15)'
*'define thetase=theta*exp(rs*lv/1004/t)'

'define alpha=287/lev*theta*pow((lev-es)/1000,(287/1004))'
'define g=9.8'
'define wz=-alpha*w/g'
'define dx=cdiff(lon,x)*cos(lat*3.1415926/180)*3.1415926/180*6371000'
'define dy=cdiff(lat,y)*3.1415926/180*6371000'
'define dz=z(z+1)-z(z-1)'
'define xu=cdiff(u,x)/dx'
'define yu=cdiff(u,y)/dy'
'define zu=(u(z+1)-u(z-1))/dz'
'define xv=cdiff(v,x)/dx'
'define yv=cdiff(v,y)/dy'
'define zv=(v(z+1)-v(z-1))/dz'
'define xw=cdiff(wz,x)/dx'
'define yw=cdiff(wz,y)/dy'
'define zw=(wz(z+1)-wz(z-1))/dz'
'define xtheta=cdiff(theta,x)/dx'
'define ytheta=cdiff(theta,y)/dy'
'define ztheta=(theta(z+1)-theta(z-1))/dz'

'define zizhuan=7.292e-5'
'define ee=2*zizhuan*cos(lat*3.14/180)'
'define f=2*zizhuan*sin(lat*3.14/180)'

'define PV1=alpha*(yw-zv)*xtheta'
'define PV2=alpha*(ee+zu-xw)*ytheta'
'define PV3=alpha*(f+xv-yu)*ztheta'
'define PV=PV1+PV2+PV3'

*'set gxout shaded'
*'set cmin 30'
*'set xyrev on'
*'d PV*1e6'

'set gxout contour'
'set ccolor 1'
'set lat 10 50';'set lon 116.5'
'set lev 700'
'set t 1 60'
'd tloop(ave(PV,t,t+3,4))*1e6'
*'drawfigstr TL (a) 1.5'

'print'
'clear'
'disable print'
'close 1'

Ashe

检查一下你的设置维度设定和数据的维度是否吻合，还可以先不去设定具体的lat lon 整体计算完毕后再设置绘图区域

杨雨轩

Ashe 发表于 2016-3-18 09:26
检查一下你的设置维度设定和数据的维度是否吻合，还可以先不去设定具体的lat lon 整体计算完毕后再设置绘图 ...

红色的报错我知道问题在哪了，是时间设错了，应该设成
'set t 3 59'
'd tloop(ave(PV,t-2,t+1,4))*1e6'
之前那么设就溢出了
但是request is……还不知道问题在哪，不过现在可以出图

Ashe

我觉得还是设置的问题，你可以尝试着把程序分段运行来检查

tyq嗯

楼主，默默问下，数据哪里找

杨雨轩

tyq嗯 发表于 2016-12-19 15:48
楼主，默默问下，数据哪里找

可以到欧洲中心的官网上下载

tyq嗯

杨雨轩 发表于 2016-12-19 16:08
可以到欧洲中心的官网上下载

下下来的数据是nc格式的么

杨雨轩

tyq嗯 发表于 2016-12-19 16:22
下下来的数据是nc格式的么

是nc格式的，可以直接画

May-xiaobai

弱弱问一句，欧洲中心官网下数据好像需要账号？我点注册，但是一直不成功。。。怎么破

杨雨轩

May-xiaobai 发表于 2017-5-13 21:15
弱弱问一句，欧洲中心官网下数据好像需要账号？我点注册，但是一直不成功。。。怎么破

现在应该已经无法注册了，因为验证码被屏蔽了