请问GrADS可以对nc数据中的温度进行设置吗？

yangchao888

首先向关注我提问的各位大神表示感谢。


我遇到的问题是这样的，我想利用GrADS画一个积冰指数图，用到的公式是Ic=[(RH-50 )×2 ]×[T×(T+14)/(-49)],RH-相对湿度，T-温度。
我编了一个gs文件进行绘图，所用的数据是ECMWF-Interim的nc数据，最后得到的积冰指示图在下面可以看到。我发现红圈区域其实并没有发生积冰的可能性，但计算出来的结果却是正值。主要是因为该区域温度大于零，但是相对湿度小于零，所以计算出的Ic为正值，其实并没有发生积冰的可能。
所以，我想在850hPa上画图时，只使用那部分相对湿度大于50%，温度在-14到0摄氏度的数据进行计算，并画图。请问各位大神，我的gs程序该怎么修改？
另外如果我画积冰指数的垂直剖面图，同样可以只使用那部分相对湿度大于50%，温度在-14到0摄氏度的数据进行计算吗？该怎么办

1. 'reinit'
   
2. ***打开文件
   
3. 'sdfopen d:/data/ecmwf25/ice-2.nc'
   
4. 'set lon 90 120'
   
5. 'set lat 10 50'
   
6. 'set lev 850'
   
7. 'set t 1'
   
8. ***积冰指数计算
   
9. 'define a=(r-50)*2*(t-273)*(t-273+14)/(-49)'
   
10. ***经纬度、高度、时间设置
    
11. 'set lon 95 115'
    
12. 'set lat 15 35'
    
13. 'set lev 850'
    
14. 'set t 1'
    
15. *绘图
    
16. 'set ccolor 4'
    
17. 'set cint 10'
    
18. 'set cmin 0'
    
19. 'set cmax 100'
    
20. 'd a'
    
21. 'draw mark 3 4.925 6.175 0.1'
    
22. 'enable print d:/pic/ICcontrast/ice-2-e-1-850.gmf'
    
23. 'print'
    
24. 'disable print'
    
25. 'reinit'

复制代码

问题问的有些小白，请海涵，谢谢。悬赏帖子的电梯请问GrADS可以对nc数据中的温度进行设置吗？http://bbs.06climate.com/forum.php?mod=viewthread&tid=51105&fromuid=61531
(出处: 气象家园)
还请各位大神帮忙出出主意啊