<div dir="ltr">Can the documentation for the relhum() function be more explicit that it is calculating with respect to ice below -20C and a mix phase between 0 and -20C.   ( I checked the values in the lookup table against a textbook. )<div><br></div><div>Also that the other moisture functions mixhum... dewtemp... are not using the same tables/formula.  Even relhum_ttd() isn&#39;t.  Therefore values from relhum are not entirely compatible with the other functions.  Given that many of those functions require rh as an input I imagine a lot of people are first using relhum to then use a follow up function.   This could lead to unexpected results given the lack of details in the documentation at the moment. </div><div><br></div><div>To clarify this isn&#39;t an issue above 0C only once below ~-10C do the differences start to become noticeable. </div><div><br></div><div>Thanks, </div><div>Alan. </div><div><br></div><div>Some code below to exemplify how using the functions back to back can add moisture to a system. </div><div><br></div><div><div>tempc = -31.3</div><div>q_gkg = 0.6</div><div>hpa = 250</div><div>dewp_c= -37.7   ; calculated from <a href="http://www.rotronic.com/humidity_measurement-feuchtemessung-mesure_de_l_humidite/humidity-calculator-feuchterechner-mr">http://www.rotronic.com/humidity_measurement-feuchtemessung-mesure_de_l_humidite/humidity-calculator-feuchterechner-mr</a></div><div><br></div><div>tempk = tempc+273.15</div><div>q_kgkg = q_gkg/1000.</div><div>pa = hpa * 100.</div><div>dewp_k = dewp_c + 273.15</div><div><br></div><div>w_gkg = q_gkg / (1 - q_gkg)</div><div>w_kgkg = q_kgkg / (1 - q_kgkg)</div><div><br></div><div>print(&quot; Differences between Relative Humidity with respect to ice and water &quot; )</div><div><br></div><div>print(&quot; Starting Conditions&quot;)</div><div>print(&quot;   Air Temperature   : &quot;+ tempc+&quot; C&quot;   )</div><div>print(&quot;   Specific Humidity : &quot;+ q_gkg+&quot; g/kg&quot;)</div><div>print(&quot;   Pressure Level    : &quot;+  hpa +&quot; hPa&quot; )</div><div><br></div><div>print(&quot;  &quot;)</div><div>print(&quot; Using relhum() function; which calculates wrt ice below -20C&quot;)</div><div>r_hum = relhum( tempk, w_kgkg, pa)   ;;; wrt ice RH = ~73%</div><div>print(&quot;   Relative Humidity  : &quot;+ sprintf(&quot;%2.1f&quot;,r_hum) +&quot; % &quot;)</div><div><br></div><div>print(&quot;  &quot;)</div><div>print(&quot; Using that RH what is our specific humidity and dewpoint? &quot;)</div><div>q_gkg_ptrh = mixhum_ptrh(hpa, tempk, r_hum, -2)   ;;  answer is now 0.8g/kg not 0.6g/kg</div><div>dew_k = dewtemp_trh( tempk, r_hum)</div><div>dew_c = dew_k - 273.15</div><div>print(&quot;   Specific Humidity  : &quot;+sprintf(&quot;%2.1f&quot;, q_gkg_ptrh )+&quot; g/kg   -- 0.2g/kg higher than we started with!!&quot;)</div><div>print(&quot;   Dew Point Temp     : &quot;+sprintf(&quot;%2.1f&quot;, dew_c )+&quot; C&quot;)</div><div><br></div><div>print(&quot;  &quot;)</div><div>print(&quot; Using an online calculator a dew point of -37.7 was calculated&quot;)</div><div>print(&quot; So what is the dew point? Using this dew point how does our Specific humidiy look&quot;)</div><div>q_gkg_ptd = mixhum_ptd( pa, dewp_k, -2)          ;; now wrt to water so ~0.6g/kg</div><div>print(&quot;   Specific Humidity  : &quot;+sprintf(&quot;%2.1f&quot;, q_gkg_ptd )+&quot; g/kg&quot;)</div></div><div><br></div><div><br></div></div>