In this research, a method for considering the spatial variation of deal layer (DL) thickness at the lateral and top surfaces of HPGe detectors is presented to achieve a more precise simulation model for efficiency calculation especially in marinelli beaker geometries. As a good approximation, the lateral surface of the detector is divided into 12 segments, assuming each segment covers 30 degrees and the top surface of detector is divided to 25 sections. The detector efficiencies are experimentally measured in each lateral segment and also in 25 points at the top of detector, using Am-241 source at each mentioned positions. The nearest DL thickness for each segment is selected through MCNP calculations of Full Energy Peak Efficiency (FEPE) with different DL thicknesses. Then a proposed detector model with a lateral DL thickness equal to average of selected DL thicknesses on lateral surface and a top DL thickness equal to average of selected DL thicknesses on top surface is used for FEPE calculation of bulk sample geometries such as Marinelli beaker containers. The resulting efficiencies are compared with the experimentally measured efficiencies of a Marinelli beaker containing a set of standard radiation source solution with specified activities. The experimental and simulation results of FEPE show a very good agreement.