求助，用GRADS画水汽通量差值检验图出错

默戏

我先分析了降水资料在总体年份里分出了多雨年和少雨年，然后进行多雨年和少雨年水汽通量的差值检验，程序可以运行，但是图画出来以后明显不对。
以下是我用的gs代码：

1. 'reinit'
   
2. 'sdfopen d:\data\micaps\shum.mon.mean.nc'
   
3. 'sdfopen d:\data\micaps\uwnd.mon.mean.nc'
   
4. 'sdfopen d:\data\micaps\vwnd.mon.mean.nc'
   
5. 'sdfopen d:\data\micaps\pres.sfc.mon.mean.nc'
   
6. 'set gxout stream'
   
7. 'enable print d:\data\data8\jianyan/chazhijianyan2.gmf white'
   
8. 
   
9. *************************多雨年*****************
   
10. yr.1=1952; yr.2 = 1957; yr.3=1963; yr.4=1964;yr.5=1967;yr.6=1970;yr.7=1983;
    
11. yr.8=1985;yr.9=1990;yr.10=1991;yr.11=1998;yr.12=2002
    
12. mo.3='Mar'; mo.4='Apr'; mo.5='May'
    
13. i=1
    
14. while(i<=12)
    
15. j=3
    
16. while(j<=5)
    
17. 'set x 1 144'
    
18. 'set y 1 73'
    
19. 'set z 1 8'
    
20. 'set time 'mo.j%yr.i
    
21. 'define duq=uwnd.2*shum/9.8'
    
22. 'define dvq=vwnd.3*shum/9.8'
    
23. 'set z 1'
    
24. 'define dqu=vint(pres.4(z=1),duq,300)*9.8/1000'
    
25. 'define dqv=vint(pres.4(z=1),dvq,300)*9.8/1000'
    
26. 'define dsd=hdivg(dqu,dqv)'
    
27. j=j+1
    
28. endwhile
    
29. i=i+1
    
30. endwhile
    
31. ***********************少雨年***************************
    
32. yr.1=1953; yr.2 = 1955; yr.3=1962; yr.4 = 1979;yr.5=1981;yr.6=1982;yr.7=1994;
    
33. yr.8=1995;yr.9=1996;yr.10=2000;yr.11=2001;yr.12=2004
    
34. mo.3='Mar'; mo.4='Apr'; mo.5='May'
    
35. n=1
    
36. while(n<=12)
    
37. m=3
    
38. while(m<=5)
    
39. 'set x 1 144'
    
40. 'set y 1 73'
    
41. 'set z 1 8'
    
42. 'set time 'mo.m%yr.n
    
43. 'define suq=uwnd.2*shum/9.8'
    
44. 'define svq=vwnd.3*shum/9.8'
    
45. 'set z 1'
    
46. 'define squ=vint(pres.4(z=1),suq,300)*9.8/1000'
    
47. 'define sqv=vint(pres.4(z=1),svq,300)*9.8/1000'
    
48. 'define ssd=hdivg(squ,sqv)'
    
49. m=m+1
    
50. endwhile
    
51. n=n+1
    
52. endwhile
    
53. **********************************差值的均值***************************************
    
54. 'define qua=ave((dqu-squ),t=1,t=12)'
    
55. 'define qva=ave((dqv-sqv),t=1,t=12)'
    
56. 'define sd=ave((dsd-ssd),t=1,t=12)'
    
57. ************************************总体****
    
58. 
    
59. yr.1=1951; yr.2 = 1952; yr.3=1953; yr.4 = 1954;yr.5=1955;yr.6=1956;yr.7=1957;
    
60. yr.8=1958;yr.9=1959;yr.10=1960;yr.11=1961;yr.12=1962;yr.13=1963;yr.14=1964;
    
61. yr.15=1965; yr.16 = 1966; yr.17=1967; yr.18 = 1968;yr.19=1969;yr.20=1970;yr.21=1971;
    
62. yr.22=1972;yr.23=1973;yr.24=1974;yr.25=1975;yr.26=1976;yr.27=1977;yr.28=1978;
    
63. yr.29=1979; yr.30 = 1980; yr.31=1981; yr.32 = 1982;yr.33=1983;yr.34=1984;yr.35=1985;
    
64. yr.36=1986;yr.37=1987;yr.38=1988;yr.39=1989;yr.40=1990;yr.41=1991;yr.42=1992;
    
65. yr.43=1993; yr.44 = 1994; yr.45=1995; yr.46 = 1996;yr.47=1997;yr.48=1998;yr.49=1999;
    
66. yr.50=2000;yr.51=2001;yr.52=2002;yr.53=2003;yr.54=2004;yr.55=2005;yr.56=2006;
    
67. yr.57=2007; yr.58 = 2008; yr.59=2009; yr.60 = 2010;yr.61=2011;yr.62=2012;yr.63=2013;
    
68. yr.64=2014;yr.65=2015;yr.66=2016
    
69. mo.3='Mar'; mo.4='Apr'; mo.5='May'
    
70. n1=1
    
71. while(n1<=66)
    
72. m1=3
    
73. while(m1<=5)
    
74. 'set x 1 144'
    
75. 'set y 1 73'
    
76. 'set z 1 8'      
    
77. 'set time 'mo.m1%yr.n1
    
78. 'define aluq=uwnd.2*shum/9.8'
    
79. 'define alvq=vwnd.3*shum/9.8'
    
80. 'set z 1'
    
81. 'define alqu=vint(pres.4(z=1),aluq,300)*9.8/1000'
    
82. 'define alqv=vint(pres.4(z=1),alvq,300)*9.8/1000'
    
83. 'define alsd=hdivg(alqu,alqv)'
    
84. m1=m1+1
    
85. endwhile
    
86. n1=n1+1
    
87. endwhile
    
88. ******************************总体的均值*****************************
    
89. 
    
90. 'define allqu=ave(alqu,t=1,t=66)'
    
91. 'define allqv=ave(alqv,t=1,t=66)'
    
92. 'define allsd=ave(alsd,t=1,t=66)'
    
93. ********************************求方差**********************************
    
94. 
    
95. 'define fcu=sqrt(ave(pow((dqu-squ),2),t=1,t=12))'
    
96. 'define fcv=sqrt(ave(pow((dqv-sqv),2),t=1,t=12))'
    
97. 'define fcs=sqrt(ave(pow((dsd-ssd),2),t=1,t=12))'
    
98. ***********************样本量开根号*********************************
    
99. 
    
100. 'define nn=sqrt(12)'
     
101. **************************计算t******************************************
     
102. 'define tu=((qua-allqu)/fcu)/nn'
     
103. 'define tv=((qva-allqv)/fcv)/nn'
     
104. 'define ts=((sd-allsd)/fcs)/nn'
     
105. 
     
106. ******显著性检验0.01、0.05、0.1 0.5 t=2.718  2.201  1.796 0.697******************
     
107. 'define ttu=maskout(abs(tu),abs(tu)-0.697)'
     
108. 'define tttu=smth9(ttu)'      
     
109. 'define ttv=maskout(abs(tv),abs(tv)-0.697)'
     
110. 'define tttv=smth9(ttv)'  
     
111. 'define tts=maskout(abs(ts),abs(ts)-0.697)'
     
112. 'define ttts=smth9(tts)'  
     
113. 
     
114. ***************************
     
115. 
     
116. 'set x 1 144'
     
117. 'set y 1 73'
     
118. 'set grid off'
     
119. 'set grads off'
     
120. 'set map 15'
     
121. 'draw map'
     
122. 'set rgb 19 255 194 194'
     
123. 'set xlint 40'
     
124. 'set ylint 20'
     
125. 'set xlopts 1 4 0.2'
     
126. 'set ylopts 1 4 0.2'
     
127. 'set clopts 1 1 0.2'
     
128. 'set gxout shaded'
     
129. 'set clevs 0'
     
130. 'set ccols 19'
     
131. 'set cterp off'
     
132. 'set clab on'
     
133. 'set csmooth on'
     
134. 'd tttu;tttv'
     
135. 'set gxout contour'
     
136. 'd qua;qva'
     
137. 'set ccolor 1'
     
138. 'set annot 1 3'                 
     
139. 'print'
     
140. 'disable print'
     
141. ;

复制代码

其中qu、qv、sd分别是纬向、经向水汽通量和散度，但是我画出来的图是这样的：

不知道那里出错了，还请各位指点迷津。

river

你所有的都是 'd tttu;tttv'  
这能出来阴影和等值线吗？

默戏

不好意思回晚了，这个gs的毛病找到了，是t的循环问题，应该把t改成36才对。