求助！关于画温度距平曲线的问题

呆萌萱

我想画如图一的温度距平曲线，我用的是ncep月平均资料，我的提取文件是


'reinit'
'sdfopen F:\benke\juping\1979_2015_3.nc'
'set fwrite F:\benke\juping\juping_100.grd'
'set gxout fwrite'
'set lev 100'
it=13
while(it<=444)
'set t 'it''
'define  qave=aave(t,lon=75,lon=105,lat=25,lat=40)'
'd qave'
it=it+1
endwhile
'define tave=ave(qave,t=13,t=444)'
'd tave'
it=13
while(it<=444)
'define ttave=ave(qave,t='it',t='it+11')'
'define jp=ttave-tave'
'd jp'
it=it+12
endwhile
'disable fwrite'
'reinit'
;

描述文件是：

dset F:\benke\juping\juping_100.grd
title t data
undef -9.99e+33
xdef 13 linear 75 2.5
ydef 2 linear 25 2.5
zdef 1 levels 100
tdef 36 linear 00Z1JAN1980 1yr
vars 1
jp 0 999 temps
endvars


画图文件是：

'reinit'
'open F:\benke\juping\juping_100.ctl'
'enable print F:\benke\juping\juping_100.gmf'

'set t 36'
*'set zlog on'
'set grads off'
'set grid off'
'set xlopts 1 1 0.15'
'set ylopts 1 2 0.15'
*'set yaxis 1000 100 25'

*'set xaxis -60 30 5'
'set gxout line'
'set ccolor 1'
'set cmark 2'
'set cthick 6'
'd jp'
*'printim e:\tdata\2015jan\line1.png white'
'print'
'disable print'

;


但是出来的结果如图二所示，这个问题怎么解决啊？求助大神们

呆萌萱

{:eb302:}{:eb302:}大神呢

liyf

fwrite如下:
'reinit'
'sdfopen F:\benke\juping\1979_2015_3.nc'
'set fwrite F:\benke\juping\juping_100.grd'
'set gxout fwrite'

it=13
while(it<=444)
'set lat 25 40'
'set lon 75 105'
'set lev 100'
'd ave(t,t='it',t='it+11')'*我不知道你的nc里面的变量名是什么，自己改
it=it+12
endwhile
'disable fwrite'
;
对应ctl:
dset F:\benke\juping\juping_100.grd
title t data
undef -9.99e+33
xdef 13 linear 75 2.5
ydef 2 linear 25 2.5
zdef 1 levels 100
tdef 36 linear 00Z1JAN1980 1yr
vars 1
jp 0 999 temps
endvars
相应gs:
'reinit'
'open F:\benke\juping\juping_100.ctl'
'enable print F:\benke\juping\juping_100.gmf'
'set t 1 36'
'set lev 100'
'define a=aave(ip,lon=75,lon=105,lat=25,lat=40)'
'set lat 25'
'set lon 75'
'set gxout line'
'd a'
'print'
'disable print'
;
至于图形的美观自己调。。。。。

要吃早饭

{:eb513:}{:eb513:}

要吃早饭

{:eb513:}{:eb513:}

呆萌萱

liyf 发表于 2016-3-31 10:50
fwrite如下:
'reinit'
'sdfopen F:\benke\juping\1979_2015_3.nc'

谢谢~~~太感激了

CRYJKW

看看，逛逛