37 lines
764 B
R
37 lines
764 B
R
#' ===================================================
|
|
#' Analysis the output results from SHUD model
|
|
#' Author: Lele Shu <shulele@lzb.ac.cn>
|
|
#' Website: www.shud.xyz
|
|
#' Date: 2022-10-10
|
|
#' Copyright: MIT-Licence
|
|
#' ===================================================
|
|
#' Function:
|
|
#'
|
|
#'
|
|
|
|
|
|
rm(list=ls())
|
|
library(rSHUD)
|
|
library(xts)
|
|
library(raster)
|
|
library(rgdal)
|
|
library(sp)
|
|
library(ggplot2)
|
|
|
|
|
|
shud.env(prjname = 'ccw', inpath = '../input/ccw/', outpath = '../output/ccw.out/')
|
|
AA = getArea()
|
|
oid = getOutlets()
|
|
qo = ts2Daily(read.tsd('../input/ccw/ccw.tsd.obs')[[1]])
|
|
qdown = readout('rivqdown')
|
|
qs = ts2Daily(qdown[, oid])
|
|
ts = time(qs)
|
|
to = time(qo)
|
|
ct= ts[to %in% ts]
|
|
qq = cbind('SIM' = qs, 'OBS' = qo)[ts,]
|
|
|
|
autoplot(qq[,])
|
|
|
|
xl = loaddata()
|
|
wb = wb.all(xl=xl)
|