def CalcTimeDer(q2mean, nt, delt):
  
  """Calculate time derivative using centered differences scheme """
  dqdt = (q2mean[2:]-q2mean[:-2])/(2.*delt) # approximation of time derivative of v-w mean q^2
    
  return dqdt