The absolute unrelaxed surface energy and its full orientation-dependent behaviour of 13 HCP metals are studied via a broken-bond base geometric model. The model is integrated with the Rose-Vinet universal potential to investigate arbitrary orientations which are not assessable by other methods. Using only three materials constants, the calculated results show only marginal discrepancies with reported experimental values, except for divalent sp metals Mg, Zn and Cd where the calculated values are lower by a factor of 2. Stereographic projections of all 13 metals show global minimum on (0001) pole with an overall anisotropy of 15% to 21%. The equilibrium crystal shape of HCP is found to be truncated hexagonal bi-prismatic, with (0001) always favoured but the bi-prismatic facets vary from one metal to another. All projection patterns show strong six-fold symmetries but are unique for every element. The patterns are found to be largely determined by an anharmonicity factor η. Best agreement with experimental findings are found for Be, Sc, Ti,Y, Zr and Hf which possess comparatively low η. We believe the stereographic projections of these elements are the more representative for HCP metals.