ctl 和nc文件同时打开

杨雨轩

我想画温度距平和风矢量的叠加场，但这两个物理量在不同的文件中，画了以后总是报错，只能画出温度距平，不知道错误出在哪

杨雨轩

grads报错

这是grads的报错
下面的是gs：
'reinit'
file='f:\yangyuxuan\theory\201512\sst\sst_anomalies.gmf'
'enable print 'file
'set vpage 0 11 0 8.5'
'set mpdset lowres'
'set mpdset cnworld'
'set poli on'
'draw map'
'set map 15 1 8'
'sdfopen f:\yangyuxuan\theory\201512\sst\hadisst_sst.nc'
'open f:\yangyuxuan\theory\201512\fnl_20151205_00_00.ctl'
'set grads off'
'set grid off'
'set lat 10 45'
'set lon 110 150'
'set xlint 5';'set ylint 5'
'define asst=ave(sst.1,t=1356,t=1704,12)'
'define ansst=sst.1(t=1716)-asst'
'd ansst'
'set gxout contour'
'set lev 500'
'set t 18'
'set gxout vector'
'd ugrdprs.2;vgrdprs.2'
'draw title 'sst_ave_12
'print'
'clear'
'set grads off'
'disable print'
'reinit'

jolincai

分辨率不一样

杨雨轩

jolincai 发表于 2015-12-23 21:12
分辨率不一样

两个文件的资料都是1度x1度的

杨雨轩

jolincai 发表于 2015-12-23 21:12
分辨率不一样

我知道了，是命令的位置放的不对，关闭nc文件后再开ctl文件就可以了

jolincai

杨雨轩 发表于 2015-12-25 18:24
我知道了，是命令的位置放的不对，关闭nc文件后再开ctl文件就可以了

你应该把正确的再次放上来 看看改动哪里才能正确画图

杨雨轩

jolincai 发表于 2015-12-26 15:25
你应该把正确的再次放上来 看看改动哪里才能正确画图

'reinit'
file='f:\yangyuxuan\theory\201512\sst\sst_anomalies.gmf'
'enable print 'file
'set vpage 0 11 0 8.5'
'set mpdset lowres'
'set mpdset cnworld'
'set poli on'
'draw map'
'set map 15 1 8'
'sdfopen f:\yangyuxuan\theory\201512\sst\hadisst_sst.nc'

'set grads off'
'set grid off'
'set lat 10 45'
'set lon 110 150'
'set xlint 5';'set ylint 5'

'define asst=ave(sst.1,t=1356,t=1704,12)'
'define ansst=sst.1(t=1716)-asst'
'close 1'

i=13
while(i<=29)
'd ansst'
'open f:\yangyuxuan\theory\201512\fnl_20151205_00_00.ctl'
'set lat 10 45'
'set lon 110 150'
'set lev 500'
'set t 'i
'set gxout vector'
'd skip(ugrdprs,2);skip(vgrdprs,2)'
'draw title 'sst_ave_12
'print'
'clear'
'set grads off'
i=i+1
endwhile

'disable print'
'reinit'

其实改动不大

IAMQWW1

杨雨轩 发表于 2015-12-26 19:03
'reinit'
file='f:\yangyuxuan\theory\201512\sst\sst_anomalies.gmf'
'enable print 'file

你好，我现在也在尝试这么做，楼主可不可以把你的ctl贴出来给我看一下，因为不是要求时间、经纬度都要一致才能画吗

杨雨轩

IAMQWW1 发表于 2016-1-20 14:24
你好，我现在也在尝试这么做，楼主可不可以把你的ctl贴出来给我看一下，因为不是要求时间、经纬度都要一 ...

dset f:/yangyuxuan/theory/201512/fnl_201512%d2_%h2_00.grib2
options template
index f:/yangyuxuan/theory/201512/fnl_20151205_00_00.grib2.idx
undef 9.999E+20
title f:/yangyuxuan/theory/201512/fnl_20151205_00_00.grib2
*  produced by g2ctl v0.0.3g
* griddef=1:0:(360 x 181):grid_template=0: lat-lon grid:(360 x 181) units 1e-006 input WE:NS output WE:SN res 48 lat 90.000000 to -

90.000000 by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160

dtype grib2
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 29 linear 00Z05dec2015 6hr
* PROFILE hPa
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000

15000 10000 7000 5000 3000 2000 1000
options pascals
vars 121
...
变量就不列出来了，希望能对你有帮助。

cch513

问什么非得关闭第一个，才能打开第二个呢？grads不是支持打开多个文件吗？我也是同时打开三个，前两个是为了站点资料插值，可以使用。第三个是fnl数据，想把高度场和降水场叠加，结果fnl资料死活不出来。刚才试了一下楼主的方法，先把前两个关闭，再打开第三个绘图，也不行呀。这两天和grads较上劲了。有谁指点迷津呀