In this paper, we consider robust channel estimation for a millimeter wave (mmWave) massive MIMO system with uniform planar arrays (UPA). For many gridless angle estimation methods of mmWave channels, the channel gains needs to be time-invariant during training. We propose a gridless method that is applicable to time-invariant and time-varying channels, and the proposed method is robust to channel variations. This is done by employing time-varying precoders and combiners that are banded block Toeplitz (BBT) matrices. The advantage of time-varying BBT precoders is two-fold. First, with time-varying BBT precoding we can fabricate time-varying channels. Thus the channel can be estimated as if it is a time-varying channel whether the actual channel is time-varying or not. We will show that, thanks to fabricated time variation in the channel, a small training overhead can be used. Secondly, BBT precoding and combining preserve the shift invariance property of the received data matrix. Due to the preservation of shift invariance property, we can employ most existing direction finding methods for estimating the angles and delays, including low complexity gridless subspace based algorithms. Simulation results are given to show that the proposed time-varying BBT precoding provides robust and accurate channel estimates with a low complexity.