ripformfactor.py seems to work. ripintensity.py needs to be fixed related to common_scale variable. Can reproduce Wenjun's S1G result. See reproduce_WenjunSun_PNAS.py