I figured out why res = res*0.5 was in femsim.f90 finally. It was because of the difference between setting dr in fem_initialize when square_pixel was true vs when it was false. I corrected things now so that this version matches the original. Also, it looks like I am placing pixels exactly like Paul did! That is very good since it means I am understanding.
Merge branch 'rmc_dev' of https://github.com/refreshx2/OdieCode into rmc_dev