mp=matrix(rnorm(20,0,1),ncol=1) # 20-year precipitation of a specific month (e.g., August)
mt=matrix(rnorm(20,0,1),ncol=1) # 20-year temperature of a specific month (e.g., August)
threp=20
thret=80
DH<-GetDH(mp,mt,threp,thret)
Run the code above in your browser using DataLab