Hybrid discontinuous Galerkin method for Stokes-Darcy problem

We used the hybrid discontinuous Galerkin method to solve the Stokes - Darcy problem and the code in the attachment is mainly from FGS. Why can I calculate the L2 error of speed well, but not the L2 error of pressure well? Does anyone know what’s wrong ?

Attachment: general-SD-HDG.py