GRADS输出两个变量到.GRD文件中出现错误——文件小了1KB？

1434272731

    今天想把计算得到的两个水汽通量的分量vpfU、vpfV写入到一个GRD文件里，但是发现结果文件的大小少了1KB，有人知道是什么原因吗？
    单独输出一个变量vpfU到一个GRD文件里大小是27102KB，按理来说将两个变量vpfU和vpfV都写入到一个GRD文件的话大小是54204KB才正确，但是实际的文件大小只有54203KB，画图的时候画到后面的几个时次就报错了，不知道是哪里出了问题，望知道的大神解答一下，谢谢！
    附上.gs脚本文件：

1. 'reinit'
   
2. 'sdfopen D:/cryosphere_environment/data/uwnd.mon.mean.nc'
   
3. 'sdfopen D:/cryosphere_environment/data/vwnd.mon.mean.nc'
   
4. 'sdfopen D:/cryosphere_environment/data/shum.mon.mean.nc'
   
5. 'sdfopen D:/cryosphere_environment/data/pres.mon.mean.nc'
   
6. 'set gxout fwrite'
   
7. 'set z 1 8'
   
8. 'set lat -90 90'
   
9. 'set lon 0 357.5'
   
10. 'set fwrite D:/cryosphere_environment/circulation/quqv/vpf.grd'
    
11. yr0=1948
    
12. yr1=1961
    
13. yr2=2015
    
14. while yr1<=yr2
    
15. imon=1
    
16. while imon<=12
    
17. it=(yr1-yr0)*12+imon
    
18. 'set t 'it''
    
19. 'define qu=uwnd*shum.3/9.8'
    
20. 'define qv=vwnd.2*shum.3/9.8'
    
21. 'set z 1'
    
22. *计算整层水汽通量define vpfU=vint(pres.4(z=1),qu,300)/1000'
    
23. 'define vpfV=vint(pres.4(z=1),qv,300)/1000'
    
24. *'define vpf=mag(vpfU,vpfV)'
    
25. 'd vpfU'
    
26. 'd vpfU'
    
27. imon=imon+1
    
28. endwhile
    
29. yr1=yr1+1
    
30. endwhile
    
31. 'disable fwrite'
    
32. ;

复制代码

1434272731

22到26行打错了，修改一下：

1. 22.*计算整层水汽通量
   
2. define vpfU=vint(pres.4(z=1),qu,300)/1000'
   
3. 
   
4. 23.'define vpfV=vint(pres.4(z=1),qv,300)/1000'
   
5. 
   
6. 25.'d vpfU'
   
7. 
   
8. 26.'d vpfV'

复制代码

maiyatang_xihui

学习学习

river

最前面的'set z 1 8'有什么作用吗？
你把两个变量分别在循环里计算输出，不要放在一个循环里试试。而且你这个循环写的有些复杂，我没理解错的话，你是把所有时次的整层水汽通量都输出了，那么直接使用第一个时次到最后一个时次循环就可以了啊，没必要区分月份

旺旺小小酥

我也遇到这个问题，再写两个变量到grd的时候就出现错误了