...
Date | Position | Tare (g) (g) | Wet (g) (g) | Dry (g) (g) | rho (g/cm^3) | Moist (%) | Moist (%) |
|---|---|---|---|---|---|---|---|
6/15 | soil.a 0-3cm | 7.981 | 74.357 | 48.354 | 0.61 | 39.2 |
|
6/15 | soil.a 3-6cm | 7.940 | 97.473 | 72.210 | 0.97 | 38.0 | 30 |
6/15 | soil.b 0-3cm | 7.945 | 90.836 | 67.464 | 0.90 | 35.2 |
|
6/15 | soil.b 3-6cm | 7.991 | 98.365 | 78.623 | 1.06 | 29.7 | 23 |
6/21 | soil.a 0-3cm | 7.972 | 51.536 | 29.353 | 0.48 | 50.1 |
|
6/21 | soil.a 3-6cm | 7.940 | 102.324 | 68.270 | 0.91 | 51.3 | 27 |
6/21 | soil.b 0-3cm | 7.927 | 59.973 | 42.641 | 0.78 | 39.1 |
|
6/21 | soil.b 3-6cm | 7.987 | 74.350 | 55.189 | 0.71 | 28.8 | 25 |
6/26 | soil.a 0-3cm | 7.972 | 116.089 | 79.391 | 1.08 | 55.2 |
|
6/26 | soil.a 3-6cm | 7.939 | 120.337 | 92.536 | 1.27 | 41.8 | 34 |
6/26 | soil.b 0-3cm | 7.926 | 78.270 | 59.924 | 0.78 | 27.6 |
|
6/26 | soil.b 3-6cm | 7.984 | 103.500 | 84.423 | 1.15 | 28.7 | 21 |
7/3 | soila 0-3cm | 7.990 | 73.581 | 46.360 | 0.58 | 41.0 |
|
7/3 | soila 3-6cm | 7.945 | 96.626 | 69.548 | 0.93 | 40.8 | 34 |
7/3 | soilb 0-3cm | 7.913 | 105.195 | 77.396 | 1.05 | 41.8 |
|
7/3 | soilb 3-6cm | 7.986 | 103.252 | 83.068 | 1.13 | 30.4 | 22 |
tare = c(c(7.981,7.940,7.945,7.991),c(7.972,7.940,7.927,7.987),c(7.972,7.939,7.926,7.984),c(7.990,7.945,7.913,7.985))
...
rho = dry/vol[c(2,2,2,2,1,2,1,2,2,2,2,2,2,2,2,2)]
m = array(moist,c(2,2,4))
grav.moist = t(0.5*(m[1,,]+m[2,,]))
grav.comp = t(m[2,,])
ec5 = matrix(c(30,23,27,25,34,21,34,22),ncol=2,byrow=T)
matplot(grav.comp,ec5,xlim=c(0,60),ylim=c(0,60)); abline(0,1,lty=2); abline(-8,1,col=3,lty=2)