ctl一个描述多个dat

四月雨之音

这个照理不是大问题，但是今个遇到了，也解决不了，就发个帖子希望遇到高手指点一下迷津。首先是我单个去写ctl描述单个dat，没问题，接着我改成批量描述的ctl，也即是加个options template,吧路径对应的日期改为%d2，最后改改t列的数据从单个45变成多个乘以45的值，但是呢一画图就发现t=45起都是缺测；我把ctl放上来
dset E:/hw/test/cfsv2.daily.pratesfc.2011.01.%d2.00.dat
undef 9.999E+20
options template
title E:\hw\PRATEsfc.2011010100.6hr.grb2
*  produced by g2ctl v0.0.3d
* griddef=1:0:(384 x 190):grid_template=40:winds(N/S): Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN number of latitudes between pole-equator=95 #points=72960 lat 89.277000 to -89.277000 lon 0.000000 to 359.062000 这两个不重要因为是我把grb2的数据fwrite成dat用的还是grb2的ctl改改来用，单个可以用，批量不行
xdef 384 linear 0.000000 0.9375
ydef 190 levels
-89.277 -88.340 -87.397 -86.454 -85.509 -84.565 -83.620 -82.676 -81.731 -此处省略多个具体纬度
zdef 1 linear 1 1
tdef 4500 linear 06Z01jan2011 1dy
vars 1
PRATEsfc  0,1,0   0,1,7 ** surface none Precipitation Rate [kg/m^2/s]
ENDVARS

四月雨之音

郁闷了，此前写批处理也就这三个步骤也都行的通，这次那出问题了呢？？

四月雨之音

我单独写dat的ctl是没问题的，t都有45层，而且和fwrite前对应的数据画出来的图一模一样，说明dat是没有错误的，而且单独为每个dat去写描述文件的话也是可以打开的，该批处理的形式以后呢，只能画到t=44…… 小伙伴们，有没有遇到过类似的情况呢~是不是因为我的数据原本的是grb2然后fwrite成dat，而不是一开始就是dat写成dat,因此才会出错呢