modified:   time-variation-average-N1.py