Nothing
#
# vim:set ff=unix expandtab ts=2 sw=2:
# This test checks the code published along with Sierra et al. (2012, Geos. Model Devel. 5: 1045)
test.SoilRPaper1=function(){
#DEACTIVATED("The examples have to be sourced out, linkeddocs had to learn this")
#require(RUnit)
source(sprintf('%s/SoilR/examples/exICBMModel.R',.libPaths()[1]))
exICBMModel_paper() #Runs the example and shows Figure 2
#
# #Figure 3
# attr(TwopFeedbackModel, "ex") #Shows the code
# attr(TwopFeedbackModel, "ex")() #Runs the example and shows Figure 3
#
# #Figure 4. Implementing the RothC model from scratch
# attr(RothCModel, "ex") #Shows the code
# attr(RothCModel, "ex")() #Runs the example and shows Figure 4
#
# #Figure 5
# #Note: Since this part requires a large dataset, it is not included yet in the test.
## library(sp)
## library(raster)
## library(ncdf)
## library(colorspace)
## mypal=rev(sequential_hcl(10))
##
## setwd(".../Data/") #Set to appropriate directory where accompanying data is stored
##
## PETb=(brick("PET_PT.WATCH.MonthlyMean.1980.2001.nc"))
## Pb=(brick("Precip.WATCH.MonthlyMean.1950.2001.nc"))
## Tbair=(brick("Tair.WATCH.MonthlyMean.1950.2001.nc"))
##
## PET=as.array(PETb)
## P=as.array(Pb)
## Tair=as.array(Tbair)
##
## TCent1=fT.Century1(Temp=Tair-273)
## WCent=fW.Century(PPT=P,PET=PET)
##
## CDIm=TCent1*WCent
## CDI=brick(CDIm)
##
## CDIa=calc(CDI,mean)
##
##
## plot(CDIa,ylab="Latitude",xlab="Longitude",col=mypal,xaxt="n",yaxt="n")
#
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.