每十年平均的纬向平均气温曲线

lqouc

绯宜物语 发表于 2014-8-4 18:24
还是不行，输出没有图像。……标记怎么关掉呀？而且感觉做出的图不太对，90N要比90S高好多。

你可以用enableprint来输出
setcmark的用法官网搜索下
图像基本正确，本来两个半球就不一样，查查相关文献就知道了。

绯宜物语

lqouc 发表于 2014-8-4 22:47
你可以用enableprint来输出
setcmark的用法官网搜索下
图像基本正确，本来两个半球就不一样，查查相关 ...

还是不行……直接运行这个gs文件，窗口不显示图。一步步输进去会显示图，但输出后没有图像。而且如果我这样写，不设置X\Y\Z\T,
'sdfopen f:\nc\1.nc'
'define var=ave(t2m,lon=0,lon=360)'
'define var2=ave(var,t=1,t=1320)'
'd var2'‘printim f:\t2m1.png'
;
画出的图和之前的不一样，而且输出有图像。

lqouc

之前没看到你还输出数据了，在printim之前先disablefwrite
维度的话我也不知道你的数据是什么样，不过zdef你最好还是设置下，其他的无所谓。

绯宜物语

lqouc 发表于 2014-8-5 11:03
之前没看到你还输出数据了，在printim之前先disablefwrite
维度的话我也不知道你的数据是什么样，不过zdef ...

不输出数据也还是不行。

写了两个gs文件：
'sdfopen f:\nc\1.nc'
'define var1=ave(t2m,lon=0,lon=360)'
'define var2=ave(var1,t=1,t=1320)'
'set ccolor 2'
'set cmark 0'
'd var2'
'printim f:\t2mlon.png'
;




'sdfopen f:\nc\1.nc'
'set x 1 721'
'set y 1 361'
'set z 1 1'
'set t 1 1320'
'define var1=ave(t2m,x=1,x=721)'
'define var2=ave(var1,t=1,t=1320)'
'set ccolor 2'
'set cmark 0'
'd var2'
'printim f:\t2mx.png'
;



后面这个是截的图，输出的还是空白。两者画出的图也不一致。而且，将第一个文件中定义中lon和x的位置互换后，'define var1=ave(t2m,t=1,t=1320)'
'define var2=ave(var1,lon=0,lon=360)'做出的图和下面一样。好纠结……

lqouc

你在第一个gs里面define之前加上set t 1 1320试试，我估计图就是第二个哪样的了，这是因为grads缺省的维度设置是t=1，var1就只有一个时刻，所以你画的第一个图是1900年1月的，而不是多年平均。
还有作纬圈平均的时候，多算了一个点，x的话就是1,720，lon的话就是0,359.5

lqouc

输出图片用enableprint不行？按说你用的那个也可以啊。

绯宜物语

lqouc 发表于 2014-8-5 12:44
输出图片用enableprint不行？按说你用的那个也可以啊。

用enableprint可以输出gmf格式，输出PNG格式打不开。用printim输出都不行。

绯宜物语

lqouc 发表于 2014-8-5 12:43
你在第一个gs里面define之前加上set t 1 1320试试，我估计图就是第二个哪样的了，这是因为grads缺省的维度 ...

哦哦，这样啊。改了，现在一样了。

绯宜物语

绯宜物语 发表于 2014-8-5 14:45
用enableprint可以输出gmf格式，输出PNG格式打不开。用printim输出都不行。

不过同时输出数据的话，窗口直接没有图像。就只能分开了……

lqouc

理论上是互相不影响的，你先输出数据，然后disablefwrite之后再display一遍数据，然后画图就应该可以了。