WaSSI_MJ_readout
Ning Liu
16 July 2016
1 read output of WaSSI-C
- this Rscript is used for reading output result from WaSSI model of AU
- the input data for this code is generally 6 binary DAT data or TXT data
- read required parameters for R script
parameters are: [1] location of data; [2] Total grids; [3] total years ; [4] Year start for WaSSI; [5] Year End for WaSSI
if this is a Rscript it should be read by:
args <- commandArgs(TRUE)
1.1 load all required libraries and set global parameters
## Loading required package: sp##
## Attaching package: 'raster'## The following object is masked from 'package:pryr':
##
## subs## [1] "loading outputs data from R_result/RESULT_NOABNORMAL.RData"
## List of 6
## $ ANNUALCARBON:'data.frame': 1466920 obs. of 5 variables:
## ..$ ID : num [1:1466920] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR: num [1:1466920] 1982 1983 1984 1985 1986 ...
## ..$ GEP : num [1:1466920] 68.78 23.9 9.52 3.84 1.56 ...
## ..$ RECO: num [1:1466920] 214 185 176 173 171 ...
## ..$ NEE : num [1:1466920] 145 162 167 169 170 ...
## $ ANNUALFLOW :'data.frame': 1466920 obs. of 12 variables:
## ..$ ID : num [1:1466920] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR : num [1:1466920] 1982 1983 1984 1985 1986 ...
## ..$ PRE : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ PET : num [1:1466920] 350 350 351 350 350 ...
## ..$ AET : num [1:1466920] 30.567 10.621 4.233 1.705 0.691 ...
## ..$ SUN_ET : num [1:1466920] 205 205 206 205 205 ...
## ..$ RUNOFF : num [1:1466920] 2.8845 1.1807 0.4798 0.1937 0.0786 ...
## ..$ Q_P : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ ET_P : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ Q_ET : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ SNWPCKMON: num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ RFACTOR : num [1:1466920] NA NA NA NA NA NA NA NA NA NA ...
## $ HUCCARBON :'data.frame': 47320 obs. of 5 variables:
## ..$ ID : num [1:47320] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ N_Y : num [1:47320] 31 31 31 31 31 31 31 31 31 31 ...
## ..$ GEP : num [1:47320] 3.5 3.5 3.5 3.5 3.5 ...
## ..$ RECO: num [1:47320] 173 173 173 173 173 ...
## ..$ NEE : num [1:47320] 169 169 169 169 169 ...
## $ HUCFLOW :'data.frame': 47320 obs. of 9 variables:
## ..$ ID : num [1:47320] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ PRE : num [1:47320] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ PET : num [1:47320] 350 350 350 350 350 ...
## ..$ AET : num [1:47320] 1.56 1.56 1.56 1.56 1.56 ...
## ..$ RUNOFF: num [1:47320] 0.157 0.157 0.157 0.157 0.157 ...
## ..$ Q_P : num [1:47320] NA NA NA NA NA NA NA NA NA NA ...
## ..$ ET_P : num [1:47320] NA NA NA NA NA NA NA NA NA NA ...
## ..$ Q_ET_P: num [1:47320] NA NA NA NA NA NA NA NA NA NA ...
## ..$ N_Y : num [1:47320] 31 31 31 31 31 31 31 31 31 31 ...
## $ MONTHCARBON :'data.frame': 17603040 obs. of 6 variables:
## ..$ ID : num [1:17603040] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR : num [1:17603040] 1982 1982 1982 1982 1982 ...
## ..$ MONTH: num [1:17603040] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ GEP : num [1:17603040] 7.77 6.91 8.01 7.65 7.66 ...
## ..$ RECO : num [1:17603040] 19.1 18.6 19.3 19 19 ...
## ..$ NEE : num [1:17603040] 11.3 11.7 11.2 11.4 11.4 ...
## $ MONTHFLOW :'data.frame': 17603040 obs. of 13 variables:
## ..$ ID : num [1:17603040] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR : num [1:17603040] 1982 1982 1982 1982 1982 ...
## ..$ MONTH : num [1:17603040] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ PRE : num [1:17603040] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ TEMP : num [1:17603040] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ SMC : num [1:17603040] 1.59 1.64 1.36 1.28 1.12 ...
## ..$ SNWPK : num [1:17603040] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ PET : num [1:17603040] 24.8 24.2 29.3 30.9 34.1 ...
## ..$ AET : num [1:17603040] 3.45 3.07 3.56 3.4 3.41 ...
## ..$ SUN_ET : num [1:17603040] 13.8 13.3 17.2 18.4 20.8 ...
## ..$ RUNOFF : num [1:17603040] 0.356 0.304 0.316 0.284 0.27 ...
## ..$ BASEFLOW : num [1:17603040] 0.356 0.304 0.316 0.284 0.27 ...
## ..$ FLOWMCMMon: num [1:17603040] 0.000356 0.000304 0.000316 0.000284 0.00027 ...
## NULL1.2 Transfer all result data to grid raster stack and plot result
# Transfer all result data to grid raster stack and plot result
## load cellinfo data for AU
load("cellinfo.RData")
load("landcover.RData")
info<-c("latmin"=min(data_cell$LAT), "latmax"=max(data_cell$LAT), "longmin"=min(data_cell$LONG), "longmax"=max(data_cell$LONG), "ncols"=nrows, "nrows"=ncols, "year_start"=min(RESULT$ANNUALFLOW$YEAR), "year_end"=max(RESULT$ANNUALFLOW$YEAR))
##plot
outs<-c("HUCFLOW","HUCCARBON","ANNUALFLOW","ANNUALCARBON") #,"MONTHFLOW","MONTHCARBON")
for (i in 1:length(outs)){
# determ month or annual
f_grid_plot(RESULT[[outs[i]]],info=info,annual=(i %in% c(grep("ANN",outs),grep("MONTH",outs))),monthly=(i %in% grep("MONTH",outs)),plot=T)
}## Loading required package: rasterVis## Loading required package: lattice## Loading required package: latticeExtra## Loading required package: RColorBrewer##
## Attaching package: 'latticeExtra'## The following object is masked from 'package:ggplot2':
##
## layer## [1] "HUC PRE"## Saving 20 x 12.7 cm image
## [1] "HUC PET"## Saving 20 x 12.7 cm image
## [1] "HUC AET"## Saving 20 x 12.7 cm image
## [1] "HUC RUNOFF"## Saving 20 x 12.7 cm image
## [1] "HUC Q_P"## Saving 20 x 12.7 cm image
## [1] "HUC ET_P"## Saving 20 x 12.7 cm image
## [1] "HUC Q_ET_P"## Saving 20 x 12.7 cm image
## [1] "HUC N_Y"## Saving 20 x 12.7 cm image
## [1] "HUC N_Y"## Saving 20 x 12.7 cm image
## [1] "HUC GEP"## Saving 20 x 12.7 cm image
## [1] "HUC RECO"## Saving 20 x 12.7 cm image
## [1] "HUC NEE"## Saving 20 x 12.7 cm image
## [1] "Annual PRE"## Saving 20 x 12.7 cm image
## [1] "Annual PET"## Saving 20 x 12.7 cm image
## [1] "Annual AET"## Saving 20 x 12.7 cm image
## [1] "Annual SUN_ET"## Saving 20 x 12.7 cm image
## [1] "Annual RUNOFF"## Saving 20 x 12.7 cm image
## [1] "Annual Q_P"## Saving 20 x 12.7 cm image
## [1] "Annual ET_P"## Saving 20 x 12.7 cm image
## [1] "Annual Q_ET"## Saving 20 x 12.7 cm image
## [1] "Annual SNWPCKMON"## Saving 20 x 12.7 cm image
## [1] "Annual RFACTOR"## Saving 20 x 12.7 cm image
## [1] "Annual GEP"## Saving 20 x 12.7 cm image
## [1] "Annual RECO"## Saving 20 x 12.7 cm image
## [1] "Annual NEE"## Saving 20 x 12.7 cm image
f_list_summary()## [1] "print info for all list objects"
## [1] "data_cell"
## 'data.frame': 47320 obs. of 7 variables:
## $ ID : int 1 2 3 4 5 6 7 8 9 10 ...
## $ LAT : num 33.2 33.2 33.2 33.2 33.1 ...
## $ LONG : num 102 102 102 102 102 ...
## $ VEG : num 10 10 10 10 10 10 10 10 10 10 ...
## $ ALT : num 3725 3795 3795 3800 3830 ...
## $ FOREST_LS: num 255 255 255 255 255 255 255 255 255 255 ...
## $ BASIN : num NA NA NA NA NA NA NA NA NA NA ...
## [1] "data_LC"
## 'data.frame': 567840 obs. of 3 variables:
## $ ID : int 1 1 1 1 1 1 1 1 1 1 ...
## $ YEAR: int 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 ...
## $ VEG : num 10 10 10 10 10 10 10 10 10 10 ...
## [1] "RESULT"
## List of 6
## $ ANNUALCARBON:'data.frame': 1466920 obs. of 5 variables:
## ..$ ID : num [1:1466920] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR: num [1:1466920] 1982 1983 1984 1985 1986 ...
## ..$ GEP : num [1:1466920] 68.78 23.9 9.52 3.84 1.56 ...
## ..$ RECO: num [1:1466920] 214 185 176 173 171 ...
## ..$ NEE : num [1:1466920] 145 162 167 169 170 ...
## $ ANNUALFLOW :'data.frame': 1466920 obs. of 12 variables:
## ..$ ID : num [1:1466920] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR : num [1:1466920] 1982 1983 1984 1985 1986 ...
## ..$ PRE : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ PET : num [1:1466920] 350 350 351 350 350 ...
## ..$ AET : num [1:1466920] 30.567 10.621 4.233 1.705 0.691 ...
## ..$ SUN_ET : num [1:1466920] 205 205 206 205 205 ...
## ..$ RUNOFF : num [1:1466920] 2.8845 1.1807 0.4798 0.1937 0.0786 ...
## ..$ Q_P : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ ET_P : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ Q_ET : num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ SNWPCKMON: num [1:1466920] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ RFACTOR : num [1:1466920] NA NA NA NA NA NA NA NA NA NA ...
## $ HUCCARBON :'data.frame': 47320 obs. of 5 variables:
## ..$ ID : num [1:47320] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ N_Y : num [1:47320] 31 31 31 31 31 31 31 31 31 31 ...
## ..$ GEP : num [1:47320] 3.5 3.5 3.5 3.5 3.5 ...
## ..$ RECO: num [1:47320] 173 173 173 173 173 ...
## ..$ NEE : num [1:47320] 169 169 169 169 169 ...
## $ HUCFLOW :'data.frame': 47320 obs. of 9 variables:
## ..$ ID : num [1:47320] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ PRE : num [1:47320] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ PET : num [1:47320] 350 350 350 350 350 ...
## ..$ AET : num [1:47320] 1.56 1.56 1.56 1.56 1.56 ...
## ..$ RUNOFF: num [1:47320] 0.157 0.157 0.157 0.157 0.157 ...
## ..$ Q_P : num [1:47320] NA NA NA NA NA NA NA NA NA NA ...
## ..$ ET_P : num [1:47320] NA NA NA NA NA NA NA NA NA NA ...
## ..$ Q_ET_P: num [1:47320] NA NA NA NA NA NA NA NA NA NA ...
## ..$ N_Y : num [1:47320] 31 31 31 31 31 31 31 31 31 31 ...
## $ MONTHCARBON :'data.frame': 17603040 obs. of 6 variables:
## ..$ ID : num [1:17603040] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR : num [1:17603040] 1982 1982 1982 1982 1982 ...
## ..$ MONTH: num [1:17603040] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ GEP : num [1:17603040] 7.77 6.91 8.01 7.65 7.66 ...
## ..$ RECO : num [1:17603040] 19.1 18.6 19.3 19 19 ...
## ..$ NEE : num [1:17603040] 11.3 11.7 11.2 11.4 11.4 ...
## $ MONTHFLOW :'data.frame': 17603040 obs. of 13 variables:
## ..$ ID : num [1:17603040] 1 1 1 1 1 1 1 1 1 1 ...
## ..$ YEAR : num [1:17603040] 1982 1982 1982 1982 1982 ...
## ..$ MONTH : num [1:17603040] 1 2 3 4 5 6 7 8 9 10 ...
## ..$ PRE : num [1:17603040] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ TEMP : num [1:17603040] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ SMC : num [1:17603040] 1.59 1.64 1.36 1.28 1.12 ...
## ..$ SNWPK : num [1:17603040] 0 0 0 0 0 0 0 0 0 0 ...
## ..$ PET : num [1:17603040] 24.8 24.2 29.3 30.9 34.1 ...
## ..$ AET : num [1:17603040] 3.45 3.07 3.56 3.4 3.41 ...
## ..$ SUN_ET : num [1:17603040] 13.8 13.3 17.2 18.4 20.8 ...
## ..$ RUNOFF : num [1:17603040] 0.356 0.304 0.316 0.284 0.27 ...
## ..$ BASEFLOW : num [1:17603040] 0.356 0.304 0.316 0.284 0.27 ...
## ..$ FLOWMCMMon: num [1:17603040] 0.000356 0.000304 0.000316 0.000284 0.00027 ...
## [1] "ANNUALCARBON"
## ID YEAR GEP RECO
## Min. : 1 Min. :1982 Min. : 0.0 Min. : 0.0
## 1st Qu.:11831 1st Qu.:1989 1st Qu.: 0.0 1st Qu.: 176.4
## Median :23661 Median :1997 Median : 462.9 Median : 475.9
## Mean :23661 Mean :1997 Mean : 447.7 Mean : 743.0
## 3rd Qu.:35490 3rd Qu.:2005 3rd Qu.: 714.5 3rd Qu.: 952.7
## Max. :47320 Max. :2012 Max. :3617.2 Max. :6828.6
## NA's :43626 NA's :43626
## NEE
## Min. :-559.68
## 1st Qu.: -7.23
## Median : 170.36
## Mean : 295.24
## 3rd Qu.: 340.77
## Max. :3211.43
## NA's :43626
## [1] "ANNUALFLOW"
## ID YEAR PRE PET
## Min. : 1 Min. :1982 Min. : 0.0 Min. : 188.4
## 1st Qu.:11831 1st Qu.:1989 1st Qu.: 0.0 1st Qu.: 350.2
## Median :23661 Median :1997 Median : 827.3 Median : 364.1
## Mean :23661 Mean :1997 Mean : 581.3 Mean : 433.4
## 3rd Qu.:35490 3rd Qu.:2005 3rd Qu.: 982.4 3rd Qu.: 482.5
## Max. :47320 Max. :2012 Max. :1768.0 Max. :1193.3
##
## AET SUN_ET RUNOFF Q_P
## Min. : 0.0 Min. :136.2 Min. : 0.0 Min. :0.000
## 1st Qu.: 0.0 1st Qu.:205.1 1st Qu.: 0.0 1st Qu.:0.000
## Median :206.7 Median :255.8 Median : 478.8 Median :0.536
## Mean :190.6 Mean :290.8 Mean : 406.0 Mean :0.400
## 3rd Qu.:308.7 3rd Qu.:341.3 3rd Qu.: 692.4 3rd Qu.:0.715
## Max. :803.8 Max. :871.0 Max. :9972.0 Max. :8.206
## NA's :64842 NA's :64879 NA's :14793
## ET_P Q_ET SNWPCKMON RFACTOR
## Min. :0.000 Min. :0.000 Min. : 0.000 Min. : 112.2
## 1st Qu.:0.000 1st Qu.:0.000 1st Qu.: 0.000 1st Qu.: 373.9
## Median :0.196 Median :0.987 Median :12.000 Median : 479.7
## Mean :0.197 Mean :0.597 Mean : 6.942 Mean : 517.7
## 3rd Qu.:0.333 3rd Qu.:1.002 3rd Qu.:12.000 3rd Qu.: 619.7
## Max. :0.993 Max. :8.257 Max. :12.000 Max. :1972.9
## NA's :14756 NA's :14793 NA's :582649
## [1] "HUCCARBON"
## ID N_Y GEP RECO
## Min. : 1 Min. :31 Min. : 0.000 Min. : 0.0
## 1st Qu.:11831 1st Qu.:31 1st Qu.: 3.537 1st Qu.: 173.7
## Median :23661 Median :31 Median : 474.441 Median : 480.0
## Mean :23661 Mean :31 Mean : 447.734 Mean : 743.0
## 3rd Qu.:35490 3rd Qu.:31 3rd Qu.: 710.911 3rd Qu.: 949.5
## Max. :47320 Max. :31 Max. :3344.972 Max. :6325.0
## NA's :1408 NA's :1408
## NEE
## Min. :-474.901
## 1st Qu.: -8.372
## Median : 169.061
## Mean : 295.240
## 3rd Qu.: 388.504
## Max. :2980.027
## NA's :1408
## [1] "HUCFLOW"
## ID PRE PET AET
## Min. : 1 Min. : 0.0 Min. : 227.6 Min. : 0.3166
## 1st Qu.:11831 1st Qu.: 0.0 1st Qu.: 350.4 1st Qu.: 1.4623
## Median :23661 Median : 845.5 Median : 372.5 Median :212.0230
## Mean :23661 Mean : 581.3 Mean : 433.4 Mean :190.5666
## 3rd Qu.:35490 3rd Qu.: 978.3 3rd Qu.: 479.9 3rd Qu.:306.1187
## Max. :47320 Max. :1342.7 Max. :1140.5 Max. :743.6083
## NA's :2094
## RUNOFF Q_P ET_P Q_ET_P
## Min. : 0.0001 Min. :20.89 Min. :0.030 Min. :0.950
## 1st Qu.: 0.1571 1st Qu.:60.11 1st Qu.:0.230 1st Qu.:0.996
## Median : 502.5417 Median :68.93 Median :0.306 Median :0.997
## Mean : 405.9980 Mean :66.99 Mean :0.326 Mean :0.996
## 3rd Qu.: 687.2831 3rd Qu.:76.55 3rd Qu.:0.394 3rd Qu.:0.997
## Max. :1237.7963 Max. :96.64 Max. :0.792 Max. :1.003
## NA's :2094 NA's :18405 NA's :18405 NA's :18405
## N_Y
## Min. :31
## 1st Qu.:31
## Median :31
## Mean :31
## 3rd Qu.:31
## Max. :31
##
## [1] "MONTHCARBON"
## ID YEAR MONTH GEP
## Min. : 1 Min. :1982 Min. : 1.00 Min. : 0.0
## 1st Qu.:11831 1st Qu.:1989 1st Qu.: 3.75 1st Qu.: 0.0
## Median :23661 Median :1997 Median : 6.50 Median : 0.0
## Mean :23661 Mean :1997 Mean : 6.50 Mean : 37.2
## 3rd Qu.:35490 3rd Qu.:2005 3rd Qu.: 9.25 3rd Qu.: 71.5
## Max. :47320 Max. :2012 Max. :12.00 Max. :635.6
## NA's :462420
## RECO NEE
## Min. : 0.0 Min. :-117.8
## 1st Qu.: 14.2 1st Qu.: 2.7
## Median : 24.4 Median : 14.2
## Mean : 61.8 Mean : 24.6
## 3rd Qu.: 75.3 3rd Qu.: 24.4
## Max. :1187.2 Max. : 551.6
## NA's :462420 NA's :462420
## [1] "MONTHFLOW"
## ID YEAR MONTH PRE
## Min. : 1 Min. :1982 Min. : 1.00 Min. : 0.000
## 1st Qu.:11831 1st Qu.:1989 1st Qu.: 3.75 1st Qu.: 0.000
## Median :23661 Median :1997 Median : 6.50 Median : 9.506
## Mean :23661 Mean :1997 Mean : 6.50 Mean : 48.440
## 3rd Qu.:35490 3rd Qu.:2005 3rd Qu.: 9.25 3rd Qu.: 84.863
## Max. :47320 Max. :2012 Max. :12.00 Max. :433.022
##
## TEMP SMC SNWPK PET
## Min. :-19.244 Min. :0.0 Min. : 0.000 Min. : 5.905
## 1st Qu.: 0.000 1st Qu.:0.0 1st Qu.: 0.000 1st Qu.: 24.646
## Median : 0.000 Median :3.1 Median : 0.099 Median : 30.733
## Mean : 1.697 Mean :2.6 Mean : 16.613 Mean : 36.114
## 3rd Qu.: 4.564 3rd Qu.:4.6 3rd Qu.: 11.461 3rd Qu.: 41.573
## Max. : 27.542 Max. :8.6 Max. :13116.099 Max. :186.719
## NA's :777027
## AET SUN_ET RUNOFF BASEFLOW
## Min. : 0.0 Min. : 0.5221 Min. : 0.0 Min. : 0.0
## 1st Qu.: 0.0 1st Qu.: 13.8128 1st Qu.: 0.0 1st Qu.: 0.0
## Median : 0.0 Median : 18.4690 Median : 0.4 Median : 0.3
## Mean : 15.8 Mean : 24.2328 Mean : 33.8 Mean : 1.6
## 3rd Qu.: 31.2 3rd Qu.: 30.7131 3rd Qu.: 55.9 3rd Qu.: 1.9
## Max. :141.7 Max. :142.5124 Max. :6536.0 Max. :49.9
## NA's :710773 NA's :777027 NA's :777027
## FLOWMCMMon
## Min. :0.0
## 1st Qu.:0.0
## Median :0.0
## Mean :0.0
## 3rd Qu.:0.1
## Max. :6.5
## NA's :777027