real(kind=dp), dimension(6) :: external_pressure

This stores the external pressure on the unit cell. The elements of this array map on to the pressure tensor as follows:

$$\begin{eqnarray*} \verb ...$$