diff --git a/docs/hegyhatsal/daily_1.png b/docs/hegyhatsal/daily_1.png index f1624e6..a8be5c7 100644 Binary files a/docs/hegyhatsal/daily_1.png and b/docs/hegyhatsal/daily_1.png differ diff --git a/docs/hegyhatsal/daily_2.png b/docs/hegyhatsal/daily_2.png index 794ab82..97849b9 100644 Binary files a/docs/hegyhatsal/daily_2.png and b/docs/hegyhatsal/daily_2.png differ diff --git a/docs/hegyhatsal/daily_3.png b/docs/hegyhatsal/daily_3.png new file mode 100644 index 0000000..1cb175f Binary files /dev/null and b/docs/hegyhatsal/daily_3.png differ diff --git a/docs/hegyhatsal/daily_4.png b/docs/hegyhatsal/daily_4.png new file mode 100644 index 0000000..0519481 Binary files /dev/null and b/docs/hegyhatsal/daily_4.png differ diff --git a/docs/hegyhatsal/daily_5.png b/docs/hegyhatsal/daily_5.png new file mode 100644 index 0000000..001b252 Binary files /dev/null and b/docs/hegyhatsal/daily_5.png differ diff --git a/docs/hegyhatsal/lai_base.png b/docs/hegyhatsal/lai_base.png new file mode 100644 index 0000000..c316a4e Binary files /dev/null and b/docs/hegyhatsal/lai_base.png differ diff --git a/docs/hegyhatsal/yield_base.png b/docs/hegyhatsal/yield_base.png new file mode 100644 index 0000000..001b252 Binary files /dev/null and b/docs/hegyhatsal/yield_base.png differ diff --git a/docs/index.html b/docs/index.html index 1d470cd..a3f21be 100644 --- a/docs/index.html +++ b/docs/index.html @@ -26,38 +26,39 @@
-

Szimuláció RBBGCMuso-val

Hollós Roland

Created: 2018-09-23 Sun 19:21

+

Szimuláció RBBGCMuso-val

Hollós Roland

2018. 09. 27

+
-

1 Bevezetés

+

1 Bevezetés

-

2 Kérdésfeltevés

+

2 Kérdésfeltevés

-

3 RBBGCMuso használata

+

3 RBBGCMuso használata

-

3.1 devtools telepítése

+

3.1 devtools telepítése

install.packages("devtools")
@@ -66,7 +67,7 @@
 
 
-

3.2 RBBGCMuso telepítése

+

3.2 RBBGCMuso telepítése

devtools::install_github("hollorol/RBBGCMuso/RGBBCMuso")
@@ -74,7 +75,7 @@
 
-

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

+

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

library(RBBGCMuso)
@@ -84,7 +85,7 @@
 
 
-

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

+

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

beallitasok <- setupMuso()
@@ -93,11 +94,97 @@
 
 
-

3.5 modellfuttatás: calibMuso()

+

3.5 modellfuttatás: calibMuso()

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

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

+
+ +
plot(kimenet[,"daily_gpp"])
+
+
+ + +
+

daily_1.png +

+
+ +
+
+

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")
+
+
+ + +
+

daily_2.png +

+
+ +
+
+

3.8 Mindez sokkal egyszerűbben

+
+ +
plotMuso(variable="daily_gpp")
+
+
+ + +
+

daily_3.png +

+
+ +
+
+

3.9 Egy apró javítás

+
+ +
plotMuso(variable="daily_gpp")+
+    labs(y=expression(paste("daily GPP [",kg[C],m^-2,"]")))
+
+
+ + +
+

daily_4.png +

+
+ +
+
+

3.10 musoMappingFind()

+ +
+
+

3.11 musoMappingFind()

+
+ +
musoMappingFind("yield")
 
@@ -107,102 +194,61 @@ kimenet[1:6,1:3] - - - + - - -tsoil0 -evapotransp -soilwevap - - -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 +3037 +cum_yieldC_HRV
-
-

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

+
+

3.12 Nézzük meg, hogy hogyan változott a termés 2013 és 2017 között!

-
plot(kimenet[,"daily_gpp"])
+
plotMuso(variable="cum_yieldC_HRV")+
+    labs(y=expression(paste("yield [",kg[C],m^-2,"]")))
 
-

daily_1.png +

yield_base.png

+
-
-

3.7 egy kicsit cicomázva:

+
+

3.13 Mentsük el az alapábránkat!

-
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")
+
termesAlap <- plotMuso(variable="cum_yieldC_HRV")+
+    labs(y=expression(paste("yield [",kg[C],m^-2,"]")))
+
+
+ + + + +
+
+

3.14 Nézzük meg, majd mentsük le a LAI-t is!

+
+ +
LAIAlap <- plotMuso(variable="proj_lai")+
+    labs(y=expression(
+             paste("Levélfelületi index ",over(m^2,m^2))))
+LAIAlap
 
-

daily_2.png -

-
- -
-
-

3.8 Mindez sokkal egyszerűbben

-
- -
plot(kimenet[,"daily_gpp"])
-
-
- - -
-

daily_1.png +

lai_base.png