有没有人用过TRMM 3A12的nc资料，不能批量打开

银酱赛高

前几天下载了TRMM 3A12的对流降水率资料，为nc格式。

单个文件都能正常打开，而且能正常画图，如下：


现在我有201001-201412的逐月平均资料，按照“q ctlinfo”出来的信息，编写了一个批量读取的ctl文件如下：

1. dset d:\ACADEMY\DATA\heatisland\convp\%y4%m2.nc
   
2. options template
   
3. undef -9999.9
   
4. dtype netcdf
   
5. xdef 60 linear 105.25 0.5
   
6. ydef 60 linear 0.25 0.5
   
7. zdef 1  linear 0 1
   
8. tdef 1  linear 00z1jan2010 1mn
   
9. edef 1
   
10. 1 1 00z1jan2010
    
11. endedef
    
12. vars 2
    
13. time_bnds 0 -103,7.20215e+252,1.35681e-153,7.09859e-154 time_bnds
    
14. convp     0 -103,-101,-100,1.30867e-47 surface instantaneous convective rain rate
    
15. endvars

复制代码

但是读出来的出现以下错误：


附上3个nc文件和ctl，请问哪位大神能解答一下这个问题？？？

银酱赛高

{:eb303:}   单个都能读出来，为什么批量就读不出来呢

银酱赛高

查了一下ctlinfo里面的edef到endedef，可能批量读取的ctl需要修改这一项，但是我对TRMM数据不太了解，这一项改成：
edef 60
m01 1 00z1jan2010
m02 2 00z1feb2010
……
m60 60 00z1dec2014
endedef

好像也不行，有了解这方面的亲吗？

MeteoInfo

把edef那几行去掉

银酱赛高

MeteoInfo 发表于 2017-6-5 16:05
把edef那几行去掉

老师，我一开始的ctl就是去掉edf那几行的，读不出来，所以才试着加回去再修改一下
请问一下这几行是什么意思呢？

银酱赛高

MeteoInfo 发表于 2017-6-5 16:05
把edef那几行去掉

老师，原来我犯了一个低级错误，tdef忘记改了
不过文件能打开了，却还是画不出图。。。
文件单独打开并画变量convp是正常的——

但批量打开显示没有convp这个变量——


以下是我的ctl文件：

1. dset d:\ACADEMY\DATA\heatisland\convp\%y4%m2.nc
   
2. options template
   
3. undef -9999.9
   
4. dtype netcdf
   
5. xdef 60 linear 105.25 0.5
   
6. ydef 60 linear 0.25 0.5
   
7. zdef 1  linear 0 1
   
8. tdef 60 linear 00z01jan2010 1mn
   
9. vars 2
   
10. time_bnds 0 99 time_bnds
    
11. convp     0 99 surface instantaneous convective rain rate
    
12. endvars

复制代码

对照单个文件的ctl信息file:///C:/Users/panqiaoy/Desktop/3.png



主要差别是：1）批量读取表述的不同；2）删除了edef三行信息；3）两个变量描述中，单个文件ctl的数值每个nc文件都有差别，而批量读取ctl中改为了99（试过改成0）。

为了排除数据源的问题，我把前10个文件（201001-201010）逐一打开，查询ctl信息都是一致的，然后ctl的tdef数量改成10，还是画不出图。。。

银酱赛高

这图片怎么放上来啊，难道要上传文件吗？