INTEGER ind2
atm_tair(ind2)= (a1*atm_tair(ind1) + a2*atm_tair(ind2))*rsuma
atm_precip(ind2)= (a1*atm_precip(ind1) +
& a2*atm_precip(ind2))*rsuma
atm_runoff(ind2)= atm_runoff(ind1)+ atm_runoff(ind2)
atm_evap_ice(ind2)= (a1*atm_evap_ice(ind1) +
& a2*atm_evap_ice(ind2))*rsuma
atm_evap_ocn(ind2)= (a1*atm_evap_ocn(ind1) +
& a2*atm_evap_ocn(ind2))*rsuma
atm_qnet_ice(ind2)= (a1*atm_qnet_ice(ind1)+
& a2*atm_qnet_ice(ind2))*rsuma
atm_qnet_ocn(ind2)= (a1*atm_qnet_ocn(ind1) +
& a2*atm_qnet_ocn(ind2))*rsuma
atm_dFdt_ice(ind2)= (a1*atm_dFdt_ice(ind1)+
& a2*atm_dFdt_ice(ind2))*rsuma
atm_dFdt_ocn(ind2)= (a1*atm_dFdt_ocn(ind1)+
& a2*atm_dFdt_ocn(ind2))*rsuma
atm_dLdt_ice(ind2)= (a1*atm_dLdt_ice(ind1)+
& a2*atm_dLdt_ice(ind2))*rsuma
atm_dLdt_ocn(ind2)= (a1*atm_dLdt_ocn(ind1)+
& a2*atm_dLdt_ocn(ind2))*rsuma
atm_dFdt_iceq(ind2)= (a1*atm_dFdt_iceq(ind1)+
& a2*atm_dFdt_iceq(ind2))*rsuma
atm_dFdt_ocnq(ind2)= (a1*atm_dFdt_ocnq(ind1)+
& a2*atm_dFdt_ocnq(ind2))*rsuma
atm_dLdt_iceq(ind2)= (a1*atm_dLdt_iceq(ind1)+
& a2*atm_dLdt_iceq(ind2))*rsuma
atm_dLdt_ocnq(ind2)= (a1*atm_dLdt_ocnq(ind1)+
& a2*atm_dLdt_ocnq(ind2))*rsuma
atm_solarinc(ind2)= (a1*atm_solarinc(ind1) +
& a2*atm_solarinc(ind2))*rsuma
atm_solar_ocn(ind2)= (a1*atm_solar_ocn(ind1)+
& a2*atm_solar_ocn(ind2))*rsuma
atm_solar_ice(ind2)= (a1*atm_solar_ice(ind1)+
& a2*atm_solar_ice(ind2))*rsuma
atm_windspeed(ind2)= (a1*atm_windspeed(ind1) +
& a2*atm_windspeed(ind2))*rsuma
atm_slp(ind2)= (a1*atm_slp(ind1) + a2*atm_slp(ind2))*rsuma
atm_pco2(ind2)= (a1*atm_pco2(ind1)+a2*atm_pco2(ind2))*rsuma
SUBROUTINE COMBINE_ENDS(a1,a2,ind1,ind2,rsuma)