CRU的降水PRE数据单位为mm/month,PET数据的单位则是mm/day
为了方便比较,需要将其转换为同一个量级
首先判断每月的具体日数
建立月份数据,判断如果当年是闰年,则2月应为29天,不是28天
然后将每月的数据除以对应月的日数
得到mm/day
ntim=dimsizes(pre&time)
nyear=ntim/12
md=(/31,28,31,30,31,30,31,31,30,31,30,31/)
do i=0,11
nn=0
do j=0,ntim-1,12
c=i+nn
isleap=isleapyear(YYYYMM(j)/100)
md(1)=where(i .eq. 1 .and. isleap.eq. "True" , 29, 28)
pre2(c,:,:)=pre(c,:,:)/md(i)
nn=nn+12
end do
end do