水汽通量画图出错，资料问题吗？

★冷霜子★

1. 'set t 1'
   
2. 'set lon 0 180'
   
3. 'set lat -90 90'
   
4. 'set lev 850'
   
5. 'define prs=lev'
   
6. 'define g=9.8'
   
7. 'define es=(6.112*exp((17.67*(air.1-273.15))/(air.1-29.65)))'
   
8. 'define qs=(0.62197*es/(prs-0.378*es))'
   
9. 'define q=qs*rhum.2/100'
   
10. 'define qu1=(uwnd.3*q/g)'
    
11. 'define qv1=(vwnd.4*q/g)'

'set t 1'<br /> 'set lon 0 180'<br /> 'set lat -90 90'<br /> 'set lev 850'<br /> 'define prs=lev'<br /> 'define g=9.8'<br /> 'define es=(6.112*exp((17.67*(air.1-273.15))/(air.1-29.65)))'<br /> 'define qs=(0.62197*es/(prs-0.378*es))'<br /> 'define q=qs*rhum.2/100'<br /> 'define qu1=(uwnd.3*q/g)'<br /> 'define qv1=(vwnd.4*q/g)'

同样的GS程序，我用月平均的NC资料计算水汽通量，画阴影图就提示：logic error 8 in gxshad，error in gxshad，但是画contour就没问题，很正常；后来我换了2005年的日资料后，两种图都可以正常画出，为什么呢？实在想不通啊。

chenhongfengye

http://bbs.06climate.com/forum.p ... &extra=page%3D1

★冷霜子★

chenhongfengye 发表于 2014-5-5 22:40
http://bbs.06climate.com/forum.php?mod=viewthread&tid=22268&extra=page%3D1

家园里的我都搜过了，计算部分基本差不太多，可是我想不明白为什么月平均资料画图会出现这种异常

星魂¤爱雪

过来学习当中

赤橙黄绿青蓝

在用fortran对变量进行平均计算时，若分母为0，则计算结果为NAN（nan在一个数被零除的时候会出现，即为无效数字）。故出错。

qlguai

楼主你好，你用这个gs画阴影图，色标的值大概是什么范围啊，我用这个画500hpa的图，出来色标的值是10以内的，我感觉不太对

ynccjy

我也是出现了同样的问题，不知道咋解决呢！