added calc_crustcor_plot.m