grads把多个nc文件写入一个grb文件中，为什么写不进去，也没显示什么错误。

LemonTree007

想把时间连续的三个nc文件写入一个二进制文件，（1901-2014）但运行之后没什么反应，不知道错在哪里？
代码如下，数据文件太大传不上来。

1. 'reinit'
   
2. 'set gxout fwrite'
   
3. 'set fwrite D:/matlabht/pdsi/pdsi1991-2014.grb'
   
4. 'sdfopen D:/matlabht/remap05x05.pdsi.mon.mean.selfcalibrated_190001-194912.nc'
   
5. 'sdfopen D:/matlabht/remap05x05.pdsi.mon.mean.selfcalibrated_195001-199912.nc'
   
6. 'sdfopen D:/matlabht/remap05x05.pdsi.mon.mean.selfcalibrated_200001-201412.nc'
   
7. 'set x 1 720'
   
8. 'set y 1 360'
   
9. 'set lev 0'
   
10. 'set dfile 1'
    
11. 'set t 13'
    
12. while('t'<601)
    
13.     'd ave(pdsi,t='t',t='t'+11)'
    
14.     'set t 't+1
    
15. endwhile
    
16. 
    
17. 'set dfile 2'
    
18. 'set t 1'
    
19. while('t'<601)
    
20.     if('t'=1)
    
21.         'd ave(pdsi,t='t',t='t'+12)'
    
22.         'set t 't+1
    
23.     else if(t>12)
    
24.         'd ave(pdsi,t='t',t='t'+11)'
    
25.         'set t 't+1
    
26.     endif
    
27. endwhile
    
28. 
    
29. 'set dfile 3'
    
30. 'set t 1'
    
31. while('t'<181)
    
32.     if('t'=1)
    
33.         'd ave(pdsi,t='t',t='t'+12)'
    
34.         'set t 't+1
    
35.     else if('t'>12)
    
36.         'd ave(pdsi,t='t',t='t'+11)'
    
37.         'set t 't+1
    
38.     endif
    
39. endwhile
    
40. 'disable fwrite'
    
41. ;
    

复制代码

LemonTree007

问题已经解决，之前因为自己没有认真看兰溪的分享，有些细节错误，今天终于弄懂了。
code如下：

1. 'reinit'
   
2. 'set gxout fwrite'
   
3. 'set fwrite d:/matlabht/pdsi/pdsi1991-2014.grd'
   
4. 'sdfopen d:/matlabht/remap05x05.pdsi.mon.mean.selfcalibrated_190001-194912.nc'
   
5. 'sdfopen d:/matlabht/remap05x05.pdsi.mon.mean.selfcalibrated_195001-199912.nc'
   
6. 'sdfopen d:/matlabht/remap05x05.pdsi.mon.mean.selfcalibrated_200001-201412.nc'
   
7. 
   
8. 'set dfile 1'
   
9. 'set x 1 720'
   
10. 'set y 1 360'
    
11. 'set lev 0'
    
12. tt=13
    
13. while(tt<=600)
    
14.     'd ave(pdsi,t='tt',t='tt'+11)'
    
15.     tt=tt+12
    
16. endwhile
    
17. 
    
18. 'set dfile 2'
    
19. 'set x 1 720'
    
20. 'set y 1 360'
    
21. 'set lev 0'
    
22. tt=1
    
23. while(tt<=600)
    
24.     if(tt=1)
    
25.         'd ave(pdsi,t='tt',t='tt'+11)'
    
26.         tt=tt+12
    
27.     else if(tt>12)
    
28.         'd ave(pdsi,t='tt',t='tt'+11)'
    
29.         tt=tt+12
    
30.     endif
    
31. endwhile
    
32. 
    
33. 'set dfile 3'
    
34. 'set x 1 720'
    
35. 'set y 1 360'
    
36. 'set lev 0'
    
37. tt=1
    
38. while(tt<=180)
    
39.     if(tt=1)
    
40.         'd ave(pdsi,t='tt',t='tt'+11)'
    
41.         tt=tt+12
    
42.     else if(tt>12)
    
43.         'd ave(pdsi,t='tt',t='tt'+11)'
    
44.         tt=tt+12
    
45.     endif
    
46. endwhile
    
47. 'disable fwrite'
    
48. ;
    

复制代码

LemonTree007

LemonTree007 发表于 2018-4-23 16:04
问题已经解决，之前因为自己没有认真看兰溪的分享，有些细节错误，今天终于弄懂了。
code如下：

但是出来的数据似乎有错误，因为画图时等值线数值有问题，比如-1e+..