The outer boundary of the plasmasphere, the plasmapause, is characterized by a sharp electron density gradient that changes under varying space weather conditions. We developed a new model, called Neustrelitz/ESOC PlasmaPause Model (NEPPM), for providing plasmapause location in terms of L-shell utilizing electron density measurements from the Van Allen probes from 2012 to 2018 and the IMAGE satellite data from 2001 to 2005. Both datasets were preprocessed, and algorithms were developed for automatic detection of plasmapause location Lpp where L denotes the McIlwain parameter. The suggested model provides a simple ellipse-based approach determined by the semi-major axis, the eccentricity, and the orientation angle of the semi-major axis. The modelled Lpp varies as a function of the Dst index and magnetic local time MLT. The NEPPM results are compared with the Global Core Plasma Model (GCPM). The plasmapause bulge in the evening hours follows the level of geomagnetic activity. The NEPPM will complete the NPSM (Neustrelitz PlasmaSphere Model), which was derived from dual-frequency GPS measurements on-board the CHAMP satellite mission