用Fortran实现《湿度查算表》电子化.doc

用Fortran实现《湿度查算表》电子化 　　（民航华东空管局气象中心 上海 200335） 摘 要：空气湿度湿表征大气物理状态的一个要素，是民航气象观测基本观测项目之一。在应急人工观测时，通过人工查阅《湿度查算表》的方式得到相对湿度和露点温度值。人工查表的过程比较繁琐，如遇复杂天气，值班人员还要忙于通报天气，很容易造成迟发报。本文将以中国气象局编写的《湿度查算表》为依据，介绍用Fortran语言作为工具进行程序设计，实现电子化《湿度查算表》的思路与方法。经检验，程序计算湿度结果与人工查算结果相吻合 关键词：湿度查算表，Fortran，二分法 一、引言 空气湿度湿表征大气物理状态的一个要素，是民航气象观测基本观测项目之一。当自动观测设备发生故障时，气象观测员需要在室外进行应急人工观测。其中湿度项目的观测采用传统的干湿表法，将百叶箱内读取到的干湿球温度，通过人工查阅《湿度查算表》的方式得到相对湿度和露点温度值。在查表过程中，由于需要对湿球温度进行订正，因此需要来回翻阅多次，不仅费时，而且很容易出错，如遇复杂天气，值班人员还要忙于通报天气，很容易造成迟发报。如能由计算机计算代替人工查算，将会大大提高观测员的工作效率，帮助观测员节约时间，从而有更多精力进行报文的编制与检查等工作 本文将以中国气象局编写的《湿度查算表》为依据，介绍用Fortran语言作为工具进行程序设计，实现电子化《湿度查算表》的思路与方法 二、《湿度查算表》的查算方法简介 《湿度查算表》主要由表1湿球结冰部分，表2湿球未结冰部分以及附表的湿球温度订正值组成。表1和表2中，每栏居中的数值为干球温度，其余每行分别给出了订正参数（n），湿球温度（tw），水汽压（e），相对湿度（U）和露点温度（td）的一组数值。附表为不同型号干湿表的湿球温度订正值，不同干湿表经过各自的湿球温度订正值订正后，就可以从表1或表2查取湿度。比如，浦东机场气象台使用通风速度为0.4m/s的球状干湿表，应使用附表3进行湿球温度值的订正 查表时，根据湿球结冰与否，决定使用表1或表2。用读取到的干湿球温度从表2（或表1）中查得订正参数n值，再用n值和本站气压（个位数四舍五入）查附表3中相应的湿球温度订正值dtw，将此订正值加在tw上，然后用干球温度和经订正后的湿球温度再查表2（或表1），即可查出相对湿度和露点温度 由此可见，要严格按照查表方法进行查算，要求观测员做十分细致的工作，才能保证工作质量。如能将这部分工作交给计算机去完成，将起到事半功倍的效果 三、开发工具简介 Fortran作为一种擅长科学计算的语言，其语法简明易读，非常接近数学公式的自然描述。在数值计算上，有强大的IMSL函数库的支持，使用Fortran语言进行科学计算能够保证很好的精度。Fortran源代码编译成.exe文件后可以直接运行。在对界面要求不高的情况下，通过读写语句即可与用户进行交互。基于以上优点，笔者选择使用Fortran进行开发 四、电子《湿度查算表》的主要功能 电子《湿度查算表》对人工查算过程做了简化，用户只需输入读取到的干湿球温度和本站气压，即可由程序计算出相应的相对湿度和露点温度，而无需关注中间的计算细节。程序以《湿度查算表》提供的有关公式为依据，计算结果与人工查算结果是一致的 五、湿度计算相关公式 1、水汽压e（hPa） （1） 式中， etw（hPa）为湿球温度所对应的纯水平液（冰）面饱和水汽压，A为干湿表系数（），P为本站气压（hPa），这里取1000hPa，t为干球温度（℃），tw为湿球温度（℃） 2、饱和水汽压 纯水平液（冰）面饱和水汽压根据Goff-Gratch公式计算得出 （2） （3） 式中， ew为纯水平液面饱和水汽压（hPa）， ei为纯水平冰面饱和水汽压（hPa）， T1=273.16K，T°K=273.15+t℃（绝对温度）。湿球结冰且湿球温度小于等于0℃时，才使用纯水平冰面饱和水汽压公式 3、相对湿度U（%） （4） 4、露点温度 根据露点温度的定义，与当时空气中水汽压e相等的饱和水汽压值所对应的温度值，即为该水汽压e对应的露点温度。因此实际水汽压e与露点温度的关系可以用Goff-Gratch公式来描述。由于露点温度无法用水汽压e直接表示，《湿度差算表》中没有给出露点温度的直接计算公式。但根据实际水汽压可以反查露点温度，为解决该问题可用数值迭代的方法求出露?c温度的数值解 露点温度的初值由马格拉斯公式给出： （5） 式中，为0℃时的饱和水汽压，为6.1078hPa，a为系数，取7.69，b也为系数，取243.92。本程序采用最简单的二分法进行数值