grads极地投影等值线断节问题

我是彤-Tina

在用GRADS画北半球极地投影的时候，发现图中红色框内的等值线出现断节的问题。搜了家园的很多帖子，感觉应该是2.5*2.5格点数据第144个格点是357.5，不是360？但是不知道怎么办了。。
然后现在大概可能因为这个问题所以画出来的图也不是正圆形。
求助各位大佬，怎么样才能解决图中数据等值线断节的问题？
不画shaded的时候，等值线是正常可以画出来的！


以下是CTL：

1. dset E:\hecheng\chun\result\hgt1000-chun.grd
   
2. title 1000 hpa      
   
3. undef -9.9621e+33
   
4. xdef  144 linear 0 2.5
   
5. ydef 73 linear -90 2.5
   
6. zdef 1 levels 1000
   
7. tdef 1 linear 1sep1961 1yr
   
8. vars 3
   
9. strong 1 99
   
10. jpq 1 99
    
11. stq 1 99
    
12. endvars

复制代码

然后是GS

1. 'reinit'
   
2. 'open E:\hecheng\chun\result\HECHENG.ctl'
   
3. 
   
4. 
   
5. 
   
6. *'set grid on 1'
   
7. 'set grads off'
   
8. 'set map 1 1 5'
   
9. 
   
10. 'set mproj nps'
    
11. 'set frame circle'
    
12. 'set lon -90 270'
    
13. 'set lat 0 90'
    
14. 
    
15. 'set mpdset  cnriver lowres'
    
16. 'run D:\tools\colors.gs'
    
17. 'set gxout shaded'
    
18. 'set clevs -8 -4 -2 -1 0 1 2 4 8'
    
19. 'd smth9(jpq)'
    
20. 
    
21. 'set gxout contour '
    
22. 'set ccolor 1'
    
23. 'd smth9(strong)'
    
24. 
    
25. 'xcbar 2 9 0.2 0.4 -line on'
    
26. 
    
27. 
    
28. 'set strsiz 0.15'
    
29. 'q w2xy 90 -3'
    
30. x1=subwrd(result,3);y1=subwrd(result,6)
    
31. 'set string 1 c 8 0'
    
32. 'draw string 'x1' 'y1' 90E'
    
33. 
    
34. 'q w2xy 180 -3'
    
35. x2=subwrd(result,3);y2=subwrd(result,6)
    
36. 'set string 1 l 8 0'
    
37. 'draw string 'x2' 'y2' 180'
    
38. 
    
39. 'q w2xy 270 -3'
    
40. x3=subwrd(result,3);y3=subwrd(result,6)
    
41. 'set string 1 c 8 0'
    
42. 'draw string 'x3' 'y3' 90W'
    
43. 
    
44. 'q w2xy 360 -3'
    
45. x4=subwrd(result,3);y4=subwrd(result,6)
    
46. 'set string 1 r 8 0'
    
47. 'draw string 'x4' 'y4' 0'
    
48. 'gxprint E:\hecheng\chun\result\p\chun-1000.png x1000 y800 white'
    
49. 
    
50. ;

复制代码

wwwwxxxx

楼主求指导
我的ctl:
DSET SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-%y4%m2%d2%h2.grd

UNDEF -999.0


OPTIONS  little_endian template

TITLE China Hourly Merged Precipitation Analysis

XDEF 700 LINEAR 70.05 0.10


YDEF 440 LINEAR 15.05 0.10

TDEF  24 LINEAR 00Z08JUL2010 1hr
ZDEF   1 LEVELS 1
VARS 2                           

crain 1 00 CH01 combined analysis (mm/Hour)

gsamp 1 00 CH02 gauge numbers

ENDVARS
在grads中打开这个ctl，时间范围是这个
Time values set:2010:7:8:0 2010:7:8:0
是不是打开了一个时间的数据？
求指导！！！感谢感谢！！！

我是彤-Tina

wwwwxxxx 发表于 2018-2-1 11:37
楼主求指导
我的ctl:
DSET SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10-%y4%m2%d2%h2.grd

啊？你要用 q ctl命令才是查看时间段啊。而且你要打开某个时次的话要 set t 吧

wwwwxxxx

我是彤-Tina 发表于 2018-2-1 15:42
啊？你要用 q ctl命令才是查看时间段啊。而且你要打开某个时次的话要 set t 吧

能加个qq吗？我实在着急，不知道怎么弄

看云淡风清

smth9这是什么命令啊？直接画看看正不正常呢

我是彤-Tina

看云淡风清 发表于 2018-2-1 19:34
smth9这是什么命令啊？直接画看看正不正常呢

smth9 是九点平滑  
我去掉之后真的正常了诶。。
谢谢！