The geothermal heat flux determines how much heat enters from beneath the ice sheet, and thus impacts the temperature and the flow of the ice sheet. In this study we investigate how much geothermal heat flux impacts the initialization of the Greenland ice sheet. We use the Community Ice Sheet Model with two different initialization methods. We find a non-trivial influence of the choice of heat flow boundary conditions on the ice sheet initializations for further designs of ice sheet modeling.