分享将多个nc文件写到一个grd文件的脚本

hello呀喵

学习了

haiyanzhang

希望有所帮助，我遇到的nc文件更麻烦，不能用sdfopen打开

haiyanzhang

楼主，貌似你说的做循环批处理不可以实现呢

LemonTree007

楼主你好，谢谢你的分享，但是我的这个运行出来没问题，但结果得出的数据明显是错的，您能帮我看看有什么不合理的code吗？ctl也能正确写出来，图也能出来，但等值线的数值明显是有错误的，不知道哪里计算错误了，导致出来的不合理。

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

复制代码