In this paper, we present the first results of the ionospheric potential (IP) calculations with the chemistry-climate model (CCM) SOCOL (Solar Climate Ozone Links). For the study, we exploit a parameterization of the difference in electric potential between Earth's surface and lower boundary of the ionosphere as a function of thunderstorm and electrified cloud properties. The model shows a good enough agreement with the IP obtained by balloon soundings. The simulated UT variation of IP exhibits a maximum at 20 Universal time (UT) and minimum at about 2 UT which agree with the UT cycle of the lightning activity. The obtained results allow understanding of IP variability pattern at diurnal, seasonal and annual timescales. We also compare our results with the IP simulated with the climate model INMCM4 using similar IP parameterization. The comparison shows a good agreement of UT cycles especially before 12 UT. Simulated IP annual cycle reaches its maximum in late spring in both models. However, the comparison also reveals some differences in amplitudes of IP variability on different time scales. The large deviations occur after 12 UT for all seasons except summer where the maximum of both results happens before 12 UT. The UT cycle of IP simulated with CCM SOCOL is in a better agreement with observations after 12 UT in terms of phase with similar timing of maximum values. The calculation of IP using climate models can help to fill up the gaps when the observed IP is not available. The interactive calculation of IP is also a step forward in coupling atmospheric and ionospheric processes.