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

yangchao888

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


我遇到的问题是这样的，我想利用GrADS画一个积冰指数图，用到的公式是Ic=[(RH-50 )×2 ]×[T×(T+14)/(-49)],RH-相对湿度，T-温度。
我编了一个gs文件进行绘图，所用的数据是ECMWF-Interim的nc数据，最后得到的积冰指示图在下面可以看到。我发现红圈区域其实并没有发生积冰的可能性，但计算出来的结果却是正值。主要是因为该区域温度大于零，但是相对湿度小于50%，所以计算出的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'

复制代码

yangchao888

{:5_201:}{:5_201:}第一次发这样的帖子，正在修炼毕业论文。谢谢大家

yangchao888

顶顶呢{:eb513:}{:eb513:}

yangchao888

{:eb328:}{:eb328:}{:eb328:}{:eb328:}没人吗？

淡水河

相对湿度小于零?

yangchao888

淡水河 发表于 2017-3-12 07:53
相对湿度小于零?

小于50%，是我写错了。请问你知道怎么写这个gs吗？就是只是用温度小于零，相对湿度大于50%的那部分数据。

yangchao888

顶顶吧，大神们，救救我呀{:eb303:}{:eb303:}

yangchao888

帖子就这样沉了？不科学啊

yangchao888

就这样沉了？

river

使用判断语句来辨别数据的大小啊，if  ····  endif