The Kardar-Parisi-Zhang equation (KPZ equation) is solved via Cole-Hopf transformation h=logu, where u is the solution of the multiplicative stochastic heat equation(SHE). In [CD20, CSZ20, G20], they consider the solution of two dimensional KPZ equation via the solution uε of SHE with the flat initial condition and with noise which is mollified in space on scale in ε and its strength is weakened as βε=βˆ 2π −logε, and they prove that when βˆ∈(0,1), 1 βε(loguε−E[loguε]) converges in distribution as a random field to a solution of Edwards-Wilkinson equation. In this paper, we consider a stochastic heat equation uε with a general initial condition u0 and its transformation F(uε) for F in a class of functions F, which contains F(x)=xp (0<p≤1) and F(x)=logx. Then, we prove that 1 βε(F(uε(t,⋅))−E[F(uε(t,⋅))]) converges in distribution as a random field to a centered Gaussian field jointly in finitely many F∈F, t, and u0. In particular, we show the fluctuations of solutions of stochastic heat equations and KPZ equations jointly converge to solutions of SPDEs which depend on u0. Our main tools are Itô’s formula, the martingale central limit theorem, and the homogenization argument as in [CNN22]. To this end, we also prove a local limit theorem for the partition function of intermediate disorder 2d directed polymers.