hflux: SELECT CASE( config_flags%isfflx )
CASE (0,2) ! with fixed surface heat flux given in the namelist
heat_flux = config_flags%tke_heat_flux ! constant heat flux value
DO j = j_start, j_end
DO i = i_start, i_end
cpm = cp * (1. + 0.8 * moist(i,kts,j,P_QV))
hfx(i,j)= heat_flux*cpm*rho(i,1,j)
ENDDO
ENDDO