Update

Not sure whether the way I compute R_vir is fine….