Fixed a few bugs: a flag compiling with gfortran, thetae calculation in flag = 2