grads色彩设置出图错误，不明白为什么？

LemonTree007

用'run define_colors.gs''set clevs -4 -3 -2 -1 1 2 3 4''set ccols 26 25 24 22 41 43 44 45 46'和不用，出图差别很大，前者明显有错误，但不知道哪里错了?第一张图代码：

1. 'reinit'
   
2. 'open D:\matlabht\pdsi\pdsi1991-2014.ctl'
   
3. 'set mpdset cnbasemap'
   
4. 'set lon 70 140'
   
5. 'set lat 10 60'
   
6. 'set lev 0'
   
7. 'set grads off'
   
8. 'set gxout shaded'
   
9. 'set grid off'
   
10. 'set xlint -10'
    
11. 'set ylint -10'
    
12. 'set xlopts 1 4 0.2'
    
13. 'set ylopts 1 4 0.2'
    
14. 'run define_colors.gs'
    
15. 'set clevs -4 -3 -2 -1 1 2 3 4'
    
16. 'set ccols 26 25 24 22 41 43 44 45 46'
    
17. 'cnbasemap ave(pdsi,t=1,t=114)'
    
18. 'cbarn'
    
19. 'run southsea.gs'
    
20. 'printim D:\matlabht\tp\pdsi000.png white'
    
21. pull dummy
    
22. 'c'

复制代码

第二张图代码：

1. 'reinit'
   
2. 'open D:\matlabht\pdsi\pdsi1991-2014.ctl'
   
3. 'set mpdset cnbasemap'
   
4. 'set lon 70 140'
   
5. 'set lat 10 60'
   
6. 'set lev 0'
   
7. 'set grads off'
   
8. 'set gxout shaded'
   
9. 'set grid off'
   
10. 'set xlint -10'
    
11. 'set ylint -10'
    
12. 'set xlopts 1 4 0.2'
    
13. 'set ylopts 1 4 0.2'
    
14. 'cnbasemap ave(pdsi,t=1,t=114)'
    
15. 'cbarn'
    
16. 'run southsea.gs'
    
17. 'printim D:\matlabht\tp\pdsi000.png white'
    
18. pull dummy
    
19. 'c'
    

复制代码

LemonTree007

应该是grd中生成的数据有问题。所以最后等值线的大小不对。

柚子皮

'set clevs -4 -3 -2 -1 1 2 3 4'  设定的数值与图2中数值数量级不一致，检查数据

LemonTree007

柚子皮 发表于 2018-4-24 09:49
'set clevs -4 -3 -2 -1 1 2 3 4'  设定的数值与图2中数值数量级不一致，检查数据

对的，生成的grd文件数据有问题，但不知道错误出现在哪里？这是我的脚本，麻烦您帮我看看数据生成过程有什么问题。脚本的目的是将三个月平均数据的nc文件先求出年平均，再综合写入一个grd文件中。逻辑清晰，不长，急求指导。

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. '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. tt=13
   
10. while(tt<600)
    
11.     'set lon 70 140'
    
12.     'set lat 10 60'
    
13.     'set lev 0'
    
14.     'set t 'tt''
    
15.     'd ave(pdsi,t='tt',t='tt'+11)'
    
16.     tt=tt+12
    
17. endwhile
    
18. 
    
19. 'set dfile 2'
    
20. tt=1
    
21. while(tt<600)
    
22.     'set lon 70 140'
    
23.     'set lat 10 60'
    
24.     'set lev 0'
    
25.     if(tt=1)
    
26.         'set t 'tt''
    
27.         'd ave(pdsi,t='tt',t='tt'+11)'
    
28.         tt=tt+12
    
29.     else if(tt>12)
    
30.         'set t 'tt''
    
31.         'd ave(pdsi,t='tt',t='tt'+11)'
    
32.         tt=tt+12
    
33.     endif
    
34. endwhile
    
35. 
    
36. 'set dfile 3'
    
37. tt=1
    
38. while(tt<180)
    
39.     'set lon 70 140'
    
40.     'set lat 10 60'
    
41.     'set lev 0'
    
42.     if(tt=1)
    
43.         'set t 'tt''
    
44.         'd ave(pdsi,t='tt',t='tt'+11)'
    
45.         tt=tt+12
    
46.     else if(tt>12)
    
47.         'set t 'tt''
    
48.         'd ave(pdsi,t='tt',t='tt'+11)'
    
49.         tt=tt+12
    
50.     endif
    
51. endwhile
    
52. 'disable fwrite'
    
53. ;
    

复制代码

LemonTree007

柚子皮 发表于 2018-4-24 09:49
'set clevs -4 -3 -2 -1 1 2 3 4'  设定的数值与图2中数值数量级不一致，检查数据

描述文件如下：打开什么的都没问题。
dset D:\matlabht\pdsi\pdsi1901-2014.grd
undef -99999
title  1991-2014 annual average sc_PDSI
xdef 141 linear 70 0.5
ydef 101 linear 10 0.5
zdef 1 linear 0 1
tdef 114 linear 00Z01JAN1901 1yr
vars  1
pdsi=>pdsi  0 t,y,x Monthly Self-calibrated Palmer Drought Severity Index using Peman-Monteith PE
ENDVARS

柚子皮

LemonTree007 发表于 2018-4-24 10:00
对的，生成的grd文件数据有问题，但不知道错误出现在哪里？这是我的脚本，麻烦您帮我看看数据生成过程有 ...

建议现将pdsi的三个时间段写入一个数组，构成连续的月平均资料，然后直接求一次年平均。

LemonTree007

柚子皮 发表于 2018-4-25 11:00
建议现将pdsi的三个时间段写入一个数组，构成连续的月平均资料，然后直接求一次年平均。

但是我不会定义多维数组，您能给个简单的类似的例子吗？z主要是不知道如何连续读入，分三个nc文件的话，数组会自动接起来吗？
‘set lon 70 140’
'set lat 10 60'
'set lev 0'
'set t 1 1368'
‘set dfile 1’
'define sc_pdsi=pdsi'
‘set dfile 2’
'define sc_pdsi=pdsi'
‘set dfile 3’
'define sc_pdsi=pdsi'

LemonTree007

柚子皮 发表于 2018-4-25 11:00
建议现将pdsi的三个时间段写入一个数组，构成连续的月平均资料，然后直接求一次年平均。

set t 1 1380'

柚子皮

搜论坛看看，我对grads的语法已经模糊了

LemonTree007

柚子皮 发表于 2018-4-25 13:34
搜论坛看看，我对grads的语法已经模糊了

好的，我再看看