AutoSHUD/Step5_ResultVisualization.R at master · ynzys/AutoSHUD · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
# Task:
# 1. Configure the Model input/output path and project names.
# 2. The values of your interest.
# 3. Load the time-series (TS) data
# 4. Do the TS-plot
# 5. 2D spatial plot.
# 6. Water balance calculation.
# 7.
# 8.
source('GetReady.R')

pp = shud.env(prjname = xfg$prjname, inpath = xfg$dir$modelin, outpath = xfg$dir$modelout)
ia=getArea(); AA = sum(ia)
res=round(sqrt(mean(ia)), -2)
rmask=shud.mask(cellsize = res)
spr=readriv.sp()
oid=getOutlets()
vns= c("eleysurf","eleyunsat","eleygw",
       "elevprcp","elevetp",
       "elevinfil","elevrech",
       # "eleqsurf", "eleqsub",
       "elevetic", "elevettr", "elevetev",'elevetp',
       "rivqdown","rivqsub", "rivqsurf","rivystage")
# undebug(BasicPlot)
xl=loaddata(varname = vns)
stop()
png.control('WaterBalance.png', path=pp$anapath)
wb=wb.all(xl=xl, apply.weekly, plot = T)
dev.off()

# gw.mon =apply.monthly(
#   xl$eleygw, mean)
# gw.r.mon = MeshData2Raster(gw.mon, stack = T)
# animate(gw.r.mon)