_RL dflwupdT
dflwupdT = emiss*stefanBoltzmann*Ts2*Tsf * 4. _d 0
dFlxdT(i,j) = -dflwupdT
dflwupdT=4. _d 0*ocean_emissivity*stefan*tsf**3
dflwupdT=4. _d 0*snow_emissivity*stefan*tsf**3
dflwupdT=4. _d 0*ice_emissivity*stefan*tsf**3
df0dT=-dflwupdT+dfshdT+dflhdT
_RL dflwupdT ! derivative of long wave with respect to T
dflwupdT= ocean_emissivity*stefan*Ts2*Tsf*4. _d 0
dflwupdT = snow_emissivity*stefan*Ts2*Tsf*4. _d 0
dflwupdT = ice_emissivity*stefan*Ts2*Tsf*4. _d 0
df0dT = -dflwupdT+dfshdT+dflhdT