PPT=seq(0,1500,by=10)
PET=rep(1500,length(PPT))
PPT.PET=fW.Century(PPT,PET)
plot(PPT/PET,PPT.PET,
ylab="f(PPT, PET) (unitless)",
main="Effects of precipitation and potential evapotranspiration on decomposition rates")
Run the code above in your browser using DataLab