change time axis name for annual mean + test for density up to 31