grads读取ctl全球数据，怎么跨越0度经度画图？

维克多·蓝晶

dset K:\ecdaily\z,u,v\grd\12.1.grd
title geopotential
undef -9.99e+8
xdef 721 linear 0 0.5
ydef 361 linear -90 0.5
zdef   1 linear 0 1
tdef 1116 linear 01jan1979 1dy
vars 1
z 1 99 the geopotential
ENDVARS

ctl的经度设置从0度开始到360度，想画跨0度线的图应该怎么修改ctl或设置输出方式？
直接用
set lon -40  40
只能画出0到40的值

四叶草

建议楼主去好好看看基础知识，比如：1、ctl文件扮演什么角色；2、数据文件是哪一个；3、ctl的经度设置从0度开始到360度，那楼主你的 set  lon -40 40 中的-40 从何而来？？？

兰溪之水

一般设置set lon -40  40 应该可以实现的，你这问题估计是你的0和360两经线重复的原因，你可以测试下～
如将你的数据存成xdef 720 linear 0 0.5，这样就不会多一根360和0度重复，试试是不是就能实现你要的效果～
比如grads官网提供的测试数据里面的：
title 5 Days of Sample Model Output
undef -2.56e+33
xdef 72 linear 0 5
ydef 46 linear -90 4
zdef 7 levels 1000 850 700 500 300 200 100
tdef 5 linear 00Z02JAN1987 1440mn
vars 8
ps  0  99  Surface Pressure
u  7  99  U Winds
v  7  99  V Winds
hgt  7  99  Geopotential Heights
tair  7  99  Air Temperature
q  5  99  Specific Humidity
tsfc  0  99  Surface Temperature
p  0  99  Precipitation
endvars
它的最后一根经线是355～
如果是这个问题的话，那恭喜你，你又得到个小经验了～

维克多·蓝晶

兰溪之水 发表于 2015-4-9 23:26
一般设置set lon -40  40 应该可以实现的，你这问题估计是你的0和360两经线重复的原因，你可以测试下～
如 ...

谢谢，应该就是这个问题，我的这组数据的NC文件里面的CTL的x是720个，用grads转成GRD时设置的  
lon 0 360，结果多出来一条经度线，所以对应的ctl我写成xdef 721才行。

兰溪之水

维克多·蓝晶 发表于 2015-4-10 08:55
谢谢，应该就是这个问题，我的这组数据的NC文件里面的CTL的x是720个，用grads转成GRD时设置的  
lon 0 3 ...

用fwrite时直接设置set x 1 720