ax[0, 0].plot(data1[:, ny_loc], 'k', label='Simulation')