_RL dflhdT
dflhdT = -lath*dEvdT(i,j)
dFlxdT(i,j) = dFlxdT(i,j) + dfshdT + dflhdT
dflhdT = -lath*dEvdT
df0dT=-dflwupdT+dfshdT+dflhdT
_RL dflhdT ! derivative of latent heat with respect to T
df0dT = -dflwupdT+dfshdT+dflhdT