位势高度叠水汽输送场图

彭小民

1. 'reinit'
   
2. 'sdfopen e:\CNRM_CM5\paleo\hus_lgm.nc'
   
3. 'sdfopen e:\CNRM_CM5\paleo\ps_lgm.nc'
   
4. 'sdfopen e:\CNRM_CM5\paleo\ua_lgm.nc'
   
5. 'sdfopen e:\CNRM_CM5\paleo\va_lgm.nc'
   
6. 'sdfopen e:\CNRM_CM5\paleo\zg_lgm.nc'
   
7. 'sdfopen e:\CNRM_CM5\paleo\hus_pi.nc'
   
8. 'sdfopen e:\CNRM_CM5\paleo\ps_pi.nc'
   
9. 'sdfopen e:\CNRM_CM5\paleo\ua_pi.nc'
   
10. 'sdfopen e:\CNRM_CM5\paleo\va_pi.nc'
    
11. 'sdfopen e:\CNRM_CM5\paleo\zg_pi.nc'
    
12. 'set lon -180 180'
    
13. 'set lat -60 80'
    
14. 'set grads off'
    
15. 'set grid off'
    
16. 'set parea 1 10 2 7'
    
17. 'set xlint 60'
    
18. 'set ylint 30'
    
19. 'set xlopts 1 9 0.2'
    
20. 'set ylopts 1 9 0.2'
    
21. 'set z 1 17'
    
22. 'define uql=ave(ua.3,t=1,t=12)*ave(hus,t=1,t=12)/9.8'
    
23. 'define vql=ave(va.4,t=1,t=12)*ave(hus,t=1,t=12)/9.8'
    
24. 'define uqp=ave(ua.8,t=1,t=12)*ave(hus.6,t=1,t=12)/9.8'
    
25. 'define vqp=ave(va.9,t=1,t=12)*ave(hus.6,t=1,t=12)/9.8'
    
26. 'set z 1'
    
27. 'define uqsl=vint(ps.2(z=1)/100,uql,300)*9.8'
    
28. 'define vqsl=vint(ps.2(z=1)/100,vql,300)*9.8'
    
29. 'define uqsp=vint(ps.7(z=1)/100,uqp,300)*9.8'
    
30. 'define vqsp=vint(ps.7(z=1)/100,vqp,300)*9.8'
    
31. 
    
32. 
    
33. *~~~~~~~~~~~~~~hgt
    
34. 
    
35. 'set lev 500'
    
36. 'define hl=ave(zg.5,t=1,t=12)'
    
37. 'define hp=ave(zg.10,t=1,t=12)'
    
38. 'set cthick 9'
    
39. 'set clskip 3'
    
40. 'set cint 100'
    
41. 'set gxout shaded'
    
42. 'C:\OpenGrADS\Classic\scripts\128colors.gs'
    
43. 'd hl-hp'
    
44. 
    
45. *~~~~~~~~~~~~~~~vapor
    
46. 
    
47. 'd uqsl-uqsp;vqsl-vqsp'
    
48. 'printim E:\CNRM_CM5\paleo\l_p.png x1800 y1200 white'
    
49. ;
    
50. 
    

复制代码

想要画lgm时段减去pi时段的位势高度叠水汽输送的图，但是画出来的图一直显示Entire Grid Undefined，求大神解答一下