Replication climate R script