Desculpe-me, o problema é:
prect <- matrix(nrow= 15000,ncol=i)
prect[[,t]] <- prec
Erro em prect[[, t]] <- prec :
mais elementos fornecidos do que podem ser substituídos.
Eu penso que é por conta do número de observações ser diferente nas estações (por conta dos dados faltantes), será?