Glaciers in Peruvian Andes play an important role for local hydrology. To quantify the availability of glacier meltwater in a climate change scenario, surface energy (SEB) and mass balance (MB) models are required. However, few studies of SEB and MB modeling using coupled process-based models were carried out for this region. In this study, we use the newly updated COupled Snowpack and Ice surface energy and mass balance model in PYthon (COSIPY) in its 2-dimesnional spatially distributed version to simulate energy and mass fluxes over Artesonraju Glacier, Peruvian Andes, from 2016 to 2018. A glacier mass balance program was started at this glacier in 1995, with installation of a network of 22 ablation stakes and two automatic weather stations were started operating since 2005. COSIPY is forced by weather parameters like solar radiation, air temperature, relativity humidity, air pressure,...