求教：知道温度和相对湿度，怎么求出露点

兰北

求教：知道温度和相对湿度（比湿），如何求出露点温度？有公式吗？非常感谢！

喜欢伽菲

北大版大气物理学书19-23页

兰北

在线等啊，先谢谢啦

兰北

感谢，已看到了……有个问题是计算的时候怎么会出现负值呢？我写的是
'define EW=(log(es*RHprs/6.1078))/17.2693'
'define dp=(35.86*EW-273.16)/(EW-1.0)'
es是饱和水汽压，RHprs是相对湿度，dp是露点

兰北

估计是数据的问题，我再看看。

西风江南

以下两个子程序在vc环境下运行，可以改编成适合自己的语言
已经知道气温t,可以求出此时饱和水气压esw, 乘以相对湿度得到此时实际水气压e，
double esw(double t){
          double pol,pol1;
          pol=0.99999683     +t*(-0.90826951e-02+
           t*(0.78736169e-04 +t*(-0.61117958e-06+
           t*(0.43884187e-08 +t*(-0.29883885e-10+
           t*(0.21874425e-12 +t*(-0.17892321e-14+
           t*(0.11112018e-16 +t*(-0.30994571e-19)))))))));
      pol1=6.1078/pow(pol,8);
          return pol1;
        };


//由实际水气压e求露点dpt
        double dpt(double e)
        {   double lu,x,dnm,fac,t,edp,dtdew,dt;
                if(e<=0.06||e>=1013.)
                {lu=9999.;
                return lu;}
        x=log(e/6.1078);
                dnm=17.269388-x;
        t=237.3*x/dnm;
        fac=1./(e*dnm);
        edp=esw(t);
        dtdew=(t+237.3)*fac;
                dt=dtdew*(e-edp);
                t=t+dt;
                while(fabs(dt)>=1.e-4)
                { edp=esw(t);
        dtdew=(t+237.3)*fac;
                dt=dtdew*(e-edp);
                t=t+dt;
                }
                return t;
        };

兰北

谢谢，已解决了！非常支持气象家园！

曾记否

兰北 发表于 2013-6-18 20:40

                               
登录/注册后可看大图

谢谢，已解决了！非常支持气象家园！

能否说一下怎么解决的？我也遇到同样的问题。谢谢！

兰北

曾记否 发表于 2013-9-23 15:55

                               
登录/注册后可看大图

能否说一下怎么解决的？我也遇到同样的问题。谢谢！

参见大气物理书（北大版）21页

chuanhao

兰北 发表于 2013-9-26 08:44
参见大气物理书（北大版）21页

你好，能把你的gs拿出来分享一下么？