A general 2D-$hp$-adaptive Finite Element (FE) implementation inFortran 90 is described. The implementation is based on an abstractdata structure, which allows to incorporate the full $hp$-adaptivityof triangular and quadrilateral finite elements.The $h$-refinement strategies are based on $h2$-refinement ofquadrilaterals and $h4$-refinement of triangles. For $p$-refinementwe allow the approximation order to vary within any element. The ...