calc_dtmax_sp function

Calculate dTmax by the successive predawn method