用nc念资料体怎么体现具体某一天某时的水汽通量散度的

流水2012

这个是用nc资料画水汽通量散度公式中，请大家帮忙解释下怎么体现具体某一天某时的水汽通量散度的，比如我要2008年6月20日20时的水汽通量散度图？

1. 由风场、温度、相对湿度计算水汽通量散度的GS文件
   
2. 'sdfopen D:\m_d\uwnd.2008.nc'
   
3. 'sdfopen D:\m_d\vwnd.2008.nc'
   
4. 'sdfopen D:\m_d\rhum.2008.nc'
   
5. 'sdfopen D:\m_d\airtemp.2008.nc'
   
6. 'set lon 70 140'
   
7. 'set lat 15 55'
   
8. 'set lev 925'
   
9. 'set mpdset cnworld'
   
10. i=889
    
11. 
    
12. while(i<=891)
    
13. 'set t 'i
    
14. 'set grads off'
    
15. 'set grid off'
    
16. 'set gxout shaded'
    
17. 'define prs=lev'
    
18. 'define g=9.8'
    
19. 'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
    
20. 'define qs=(0.62197*es/(prs-0.378*es))'
    
21. 'define q=qs*rhum.3/100'
    
22. 'd hdivg(q*uwnd.1/g,q*vwnd.2/g)'
    
23. 'run cbarn'
    
24. 
    
25. 'printim D:\m_d\pic\rhdiv200808'i'-925hPa.png white'
    
26. 'c'
    
27. 
    
28. i=i+1
    
29. endwhile
    
30. ;
    
31. 
    
32. 
    

由风场、温度、相对湿度计算水汽通量散度的GS文件<br /> 'sdfopen D:\m_d\uwnd.2008.nc'<br /> 'sdfopen D:\m_d\vwnd.2008.nc'<br /> 'sdfopen D:\m_d\rhum.2008.nc'<br /> 'sdfopen D:\m_d\airtemp.2008.nc'<br /> 'set lon 70 140'<br /> 'set lat 15 55'<br /> 'set lev 925'<br /> 'set mpdset cnworld'<br /> i=889<br /> <br /> while(i<=891)<br /> 'set t 'i<br /> 'set grads off'<br /> 'set grid off'<br /> 'set gxout shaded'<br /> 'define prs=lev'<br /> 'define g=9.8'<br /> 'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'<br /> 'define qs=(0.62197*es/(prs-0.378*es))'<br /> 'define q=qs*rhum.3/100'<br /> 'd hdivg(q*uwnd.1/g,q*vwnd.2/g)'<br /> 'run cbarn' <br /> <br /> 'printim D:\m_d\pic\rhdiv200808'i'-925hPa.png white'<br /> 'c'<br /> <br /> i=i+1<br /> endwhile<br /> ;<br /> <br /> <br />

mofangbao

表示没看懂楼主的提问

做个霸气的木头

数据精确到什么时刻，水汽通量就可以的吧

river

这个是批量出图的一个gs，楼主想要哪个高度的那个时间的，自己改一下设置，很容易就能把所有高度上的所有时间的涂画出来。你说的公式只是这个脚本的一部分而已。楼主要是不明白我说的，那还是先看看实用手册吧，基本的东西得先学会才行

流水2012

river 发表于 2013-3-11 21:21

                               
登录/注册后可看大图

这个是批量出图的一个gs，楼主想要哪个高度的那个时间的，自己改一下设置，很容易就能把所有高度上的所有时 ...

恩，这个我已经明白了，改下时间设置就出图了，是啊，还是要把基本的弄好，哈哈，谢谢指导！

dlmuqixiang

好好学习一下，谢谢