温度对数压力图中等饱和比湿线的绘制

shuuichi

各位大神好：
小弟最近在尝试使用Matlab绘制温度对数压力图，其中等饱和比湿线的绘制遇到了点困难，就是那个绿色实线？

等饱和比湿线的计算是通过固定比湿，计算不同温度处的气压值，
公式是这个：q=0.622e/(P-0.378e)
但是对于饱和水气压的计算遇到了困难，目前利用的Magnus公式，就是那个
e=e0 exp(at/(273.16+t-b))
对于冰面和水面a和b取不同常数，当t>-15°当水面处理，当t<-40是当冰面处理，当-40到-15°间时
按e=0.022((80+2t)*e水-(30+2t)*e冰)处理
但是当绘制0.01,0.05,0.1,0.2时按冰面处理得到的等饱和比湿线与标准图较为一致
当绘制10，12，14,16，。。。30,35,40按水面处理得到的等饱和比湿线与标准图较为一致
但是0.2-10之间的等饱和比湿线绘制的并不准确，请问这一部分用什么公式计算能够绘出与标准图一样的等饱和比湿线呢？
或者说标准艾玛图上的等饱和比湿线是按什么标准绘制的呢》？？？
感谢！！！！

shuuichi

求教啊，谢谢谢谢