After the group meeting on Oct. 3, I need to do some modification on the code and plot new figures