still an issue with units - qevap was being computed in SI units but wasn't being converted to cal/min/cm2, but this is now corrected, and have also corrected the calculation of HC in dsub prior to HD being calculated, because this value of HC must be in SI units for EVAP but is being derived from qconv which is in non-SI units