grads求假相当位温任意两点剖面

茄子

'reinit' 'sdfopen d:\by\xt.nc' 'enable print d:\by\wwpm-6-27-20.gmf' 'set grads off' 'set grid off' 'set mpdset hires' 'set mpdset  cnworld' 'set lon 105 115' 'set lat 28 34' 'set lev 1000 150' 'set t 13' 'define prs=lev' 'define mm=r' 'define e1=17.27*(t-273.16)/(t-35.86)' 'define es=6.112*exp(e1)' 'define qs=mm*(0.62197*es/(prs-0.378*es))/100' 'define e=prs*q/(0.62197+q)+1e-10' 'define tk=55.0+2840.0/(3.5*log(t)-log(e)-4.805)' 'define pot=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))' 'define eqt=pot*exp(((3376./tk)-2.54)*q*(1.0+0.81*q))' lon1 = 105.0 lon2 = 115.0 lat1 = 40.0 lat2 = 30.0 lon = lon1 'collect 1 free' while(lon<=lon2) lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1) 'collect 1 gr2stn(eqt,'lon','lat')' lon=lon+1 endwhile 'set lon 105 115' 'set lev 1000 150' *'set xlabs 34N,108E|33N,109E|32N,110E|31N,111E|30N,112E|29N,113E|28N,114E|' 'set gxout contour' 'd coll2gr(1,-u)' 'print' 'disable print' ; 上面是求假相当位温，任意两点剖面，运行后出现问题，请大家帮忙看看。

茄子

运行时，出现collect error Z or T must be the only varying dimension