calc_dtmax_pd function

Calculate dTmax by the daily predawn method