Szimuláció RBBGCMuso-val

Hollós Roland

Created: 2018-09-23 Sun 19:21

Tartalomjegyzék

1. Bevezetés
2. Kérdésfeltevés
3. RBBGCMuso használata

1 Bevezetés

2 Kérdésfeltevés

3 RBBGCMuso használata

3.1 devtools telepítése

install.packages("devtools")

3.2 RBBGCMuso telepítése

devtools::install_github("hollorol/RBBGCMuso/RGBBCMuso")

3.3 csomagok betöltése, munkakönyvtár beállítása:

library(RBBGCMuso)
library(ggplot2)

3.4 futtatás előkészítése: setupMuso()

beallitasok <- setupMuso()

3.5 modellfuttatás: calibMuso()

kimenet <- calibMuso(settings = beallitasok,skipSpinup=FALSE)
kimenet[1:6,1:3]

tsoil₀	evapotransp	soilw_evap
3.09415192745925	0.214592813613201	0.214592813613201
2.5353521572297	0.0115324270940953	0.0115324270940953
1.62868961612864	0.0222151624832023	0.0222151624832023
2.34137184558436	0.1	0.1
3.66348482398438	0.454946076512924	0.454946076512924
5.0076149044176	0.14425963841612	0.14425963841612

3.6 kimeneti adatok ábrázolása (első verzió)

plot(kimenet[,"daily_gpp"])

3.7 egy kicsit cicomázva:

par(mar=c(5,5,4,4))
plot(y=kimenet[,"daily_gpp"],
     x=as.Date(rownames(kimenet),"%d.%m.%Y"), xlab="Days",
     ylab = expression(paste("daily GPP[",kg[C],m^-2,"]")),
     type="l", col="blue", lwd="1.8")

3.8 Mindez sokkal egyszerűbben

plot(kimenet[,"daily_gpp"])