In the plots comparing data and models, can you use the legend to indicate which is which? e.g. use dots for data and solid lines for models, and then label them as that in the legend. Also nice to include error bars on the temperautre measurements; I think there's a python way to plot this as a shaded region as well.