add read_write_data.f90, self-consistent loop is running
run self-consistent loop and get some results for Chi and Sigma