standalone/equilibrium_moisture.m

   1 function [d,w]=equilibrium_moisture(RH,T)
   2 RH = RH * 100;
   3 d=0.942*RH.^0.679 + 0.000499*exp(0.1*RH) + 0.18*(21.1+273.15-T).*(1-exp(-0.115*RH)); % equilibrium moisture for drying
   4 w=0.618*RH.^0.753 + 0.000454*exp(0.1*RH) + 0.18*(21.1+273.15-T).*(1-exp(-0.115*RH)); % equilibrium moisture for wetting
   5 d = d/100;
   6 w = w/100;
   7 end