Simulation model cells of a reservoir having gas condensate production from layers of a subsurface hydrocarbon reservoir are determined by the Total Dynamic Productivity Index (TDPI) method. The determinations are based on rock and fluid physics, including non-Darcy flow and pseudo-pressure integral, which are responsible for pressure drop near wellbores in the gas condensate wells. A three-dimensional (3D) grid matrix is linked with a well placement optimization algorithm to target high gas energy spots and increase productivity, efficiency and recovery from the gas condensate wells.