This is the source code used in Landis et al., [2017] on the vapor production from various configurations of water ice on Ceres (DOI:10.1002/2017JE005335)
Buried ice table sublimation
"Ceres_annual_average_temperatures.m" should be run first. It calculates the annual average temperatures with latitude over Ceres, assuming it is a smooth sphere. This must be run first to give temperature files used later. The default is to calculate annual average temperatures for 2-20 obliquity.
"Ceres_obliquity_variation_regolith_buildup.m" calculates depth of sublimation lag for a pore-filling water ice table using the varying obliquity described in the literature (Ermakov et al., [2017]).
"Ceres_constant_obliquity_vapor_production.m" loads one of the annual average temperature files generated from "Ceres_annual_average_temperatures.m" and calculates the vapor and depth of sublimation lag produced on Ceres over solar system history. Can handle pore-filling and excess water ice cases.
Exposed surface ice sublimation
"Flat_and_intra_crater_terrain_vapor_production.m" calculates sublimation for exposed surface ice (e.g. ice without a dessicated regolith overlayer). The code can be set to model flat terrain or terrain within an impact crater (impact crater parameters currently set to those of Oxo crater, Combe et al., [2016])
"Ceres_brightspot_fading_code.m" determines the fading lifetime of water ice patch on Ceres (depending on latitude and assumed dust content) by having a monolayer of regolith particles build up over the surface.
Particle lofting
"finite_area_particle_lofting.m" plots the particle diameter lofted at certain heights at a given temperature above a water ice patch on Ceres.