;;;;;;;;;;;装换时间
t = f->time
YYYY =cd_calendar(t, -1)
yrstrt =201901 ；起始年、月
yrlast =201912 ；结束年、月
ist =ind(yrstrt.eq.YYYY)
iet =ind(yrlast.eq.YYYY)
pre19 =pre0(ist:iet,: