def compute_RG(snapshot, atom='CA'):
    '\n    Radius of gyration...\n    '
    (coords, resis) = read_PDB(snapshot, atom)
    R_cm = np.mean(coords, axis=0)
    dR = (coords - R_cm)
    mag_R = np.sum((dR * dR), axis=1)
    RG = np.sqrt(np.mean(mag_R))
    return RG
Radius of gyration...
def compute_RG(snapshot, atom='CA'):
    '\n    \n    '
    (coords, resis) = read_PDB(snapshot, atom)
    R_cm = np.mean(coords, axis=0)
    dR = (coords - R_cm)
    mag_R = np.sum((dR * dR), axis=1)
    RG = np.sqrt(np.mean(mag_R))
    return RG