求助：如何在场计算中解决零对数的问题

liusx2007

    各位高手，利用NCEP（FNL1°X1°grib1资料，求饱和水汽压和实际水汽压湿完全正确，
为何在计算露点温度时出现如下一大堆错误？
    露点温度公式中 log(Ee/6.112))为自然对数，Ee为实际水汽压。
    在实数范围内，负数和零没有对数。
    也就是说，计算的实际水汽压Ee有可能为零。
    最有可能实际水汽压Ee在对流层上层（300hPa以上）的计算中出现零。
    我在计算实际水汽压Ee时，后面加0.1, 0.5, 1等数据后，计算仍出现类似的问题。
    请高手指点，如何在grads场计算中解决零对数的问题？

Warning from LOG:  Data has 16 values <= zero
                   These were set to the undefined value
Warning from LOG:  Data has 16 values <= zero
                   These were set to the undefined value
Warning from LOG:  Data has 221 values <= zero
                   These were set to the undefined value
Warning from LOG:  Data has 221 values <= zero
                   These were set to the undefined value
。。。。。。。。。。。。

警告日志：数据有16个值≤0

这些被设置为未定义的值

。。。。。。。。。。。。
-------------------------------------------------------------------------------------------------------
  'define Es=6.112*exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))';*饱和水汽压
  'define Ee=RHprs*Es/100+0.1'   
  'define Td=(243.5*log(Ee/6.112))/(17.62-log(Ee/6.112))+273.15'

kkloutf

{:loveliness:}{:loveliness:}{:loveliness:}{:loveliness:}{:loveliness:}

传说中的谁

set clevs之后就不要再set cint了，同样是设置间隔的，后面一个把前一个覆盖掉了，也就是说起作用的是后一个

chengwei06am

传说中的谁 发表于 2014-3-24 19:32

                               
登录/注册后可看大图

set clevs之后就不要再set cint了，同样是设置间隔的，后面一个把前一个覆盖掉了，也就是说起作用的是后一个 ...

麻烦问一下，像这张降水量间隔不同的怎么填色？是分别设置cmin和cmax限制出图数据填色吗，还是其他方法？谢谢~

river

这个看着没问题啊，这个出错提示一般是你的标点写成中文状态下的了，或者少了 '  之类的

liusx2007

  各位高手，求饱和水汽压和实际水汽压湿完全正确，为何在计算露点温度时
出现如下一大堆错误？
    露点温度公式中 log(Ee/6.112))为自然对数，Ee为实际水汽压。    在实数范围内，负数和零没有对数。
    也就是说，计算的实际水汽压Ee有可能为零，最有可能在对流层上层的计算中出现零。
    我在计算的实际水汽压Ee后面加0.1, 0.5, 1等数据后，计算仍出现类似的问题。
    请高手指点，如何在grads场计算中解决零对数的问题？

Warning from LOG:  Data has 16 values <= zero
                   These were set to the undefined value
Warning from LOG:  Data has 16 values <= zero
                   These were set to the undefined value
Warning from LOG:  Data has 221 values <= zero
                   These were set to the undefined value
Warning from LOG:  Data has 221 values <= zero
                   These were set to the undefined value
。。。。。。。。。。。。


警告日志：数据有16个值≤0
这些被设置为未定义的值
。。。。。。。。。。。。
-------------------------------------------------------------------------------------------------------
  'define Es=6.112*exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))';*饱和水汽压
  'define Ee=RHprs*Es/100+0.1'   
  'define Td=(243.5*log(Ee/6.112))/(17.62-log(Ee/6.112))+273.15'

爽歪歪

请问楼主这个问题解决了吗，我也出现了相同的情况

sunflower

感谢楼主分享，我也尝试一下看看

气候-谭志强

sunflower 发表于 2015-9-23 10:51
感谢楼主分享，我也尝试一下看看

楼主是不是同时设置了set lev 和set z 了

xgm324

同样求助该问题