Why is the calculated Power Spectrum peak at $l=302$ different from the chart peak at $l=220$? This is a follow-up to a really basic question I had: How is the first acoustic peak calculated in CMB? Plank quotes an angular size of the sound horizon at 0.0104, but that corresponds to an multipole of $l=\frac{\pi}{0.0104}=302$, but the first peak on the chart is clearly around $l=220$.
Edit: Apparently there's a "Shift Parameter" that changes the angular distance, $D_A(z*)$.  This parameter appears to be linked to the amount of Dark Energy, however, it appears that this parameter has already been factored in by the Plank Collaboration in their calculation of the angular scale.
 A: This question is partially answered by Relation between multipole moment and angular scale of CMB .
The relationship that $\theta \sim \pi/l$ is just an approximate relation and there is no direct one-to-one mapping between angular scale and multipole moment. Angular variations with a set value of $\theta$ contribute power over a range of spherical harmonics defined in terms of $l$ and $m$.
The approximate relationship is derived by just noting that there are $2l$ nodes around the equator. The angular separation of these nodes is $\pi/l$. But if you go to a "latitude" $\delta$ on your sphere and draw a line of constant $\delta$, then you also intercept $2l$ nodes, but the angular separation of each node is now $\pi \sin \delta/l$.
The solid angle of sky visible at any $\delta$ is $2\pi \cos \delta\ d\delta$, so
an "average" value of $\delta$ is $2/\pi$ and $\sin \bar{\delta} = 0.594$, so an alternative "average" value for angular scale (along circumferences of constant $\delta$) is 
$$\Delta \theta \simeq 0.594 \frac{\pi}{l}$$
However at the same time you still have $2l$ nodes along a full $2\pi$ angle along a line of longitude, each node separated by an angle $\pi/l$.
The two approximations bracket the actual value of the peak $l$ and so perhaps a better (best?) hand-waving correspondence between $\theta$ and $l$ would be to take the geometric mean:
$$ \bar{\Delta \theta} \simeq \sqrt{0.594}\, \frac{\pi}{l} = 0.77\, \frac{\pi}{l}$$
This factor of 0.77 is indeed quite close to $220/302 = 0.73$.
Your thoughts about the "shift parameter" are a red herring. The shift parameter is just something that can be measured from the CMB. It is closely related to where the first acoustic peak occurs and is theoretically defined in terms of a combination of cosmological parameters - it is strongly dependent on the curvature of the universe but only weakly dependent on the cosmological constant or equation of state for dark energy.
A: This is an interesting question, I hadn't come across the shift parameter, its conspicuous absence from the Planck 2015 paper (https://arxiv.org/abs/1502.01589) suggests that this effect is taken into account in a different way by the Planck collaboration (I would be interested to hear suggestions about this).
In the (somewhat limited) sources I managed to find, these slides, http://bccp.berkeley.edu/o/beach_program-oops/presentations11/COTB11Melchiorri2.pdf, suggest that the shift parameter has a simple geometric interpretation,
$R=\frac{H_0}{\sqrt{\Omega_m}}\int_0^{z_{dec}}\, \frac{dz}{H(z)}$
(from slide 21), but there is still no physical intuition for what is causing this shift. 
