add filling at T=0 to calculation step