qwf_d -nitf1.qwfc_d*down(ht.Cp)*(ht.Td-nitf1.TdWF)/sqrt(1.0E-6+max(nitf1.Tplus_d,0)^2) "Wall function flux, downside"