_RL dfshdT
dfshdT = -atmcp*tau*rh(i,j)
dFlxdT(i,j) = dFlxdT(i,j) + dfshdT + dflhdT
dfshdT = -atmcp*atmrho_loc(i,j)*exf_iceCh*wsm
dfshdT = -atmcp*atmrho*exf_iceCh*wsm
dfshdT = -csha
df0dT=-dflwupdT+dfshdT+dflhdT
_RL dfshdT ! derivative of sensible heat with respect to T
df0dT = -dflwupdT+dfshdT+dflhdT