A few generic Splus scripts have been developed to allow the user to quickly plot VIC model output. The complete library of scripts is available via ftp. When untarred all of the Splus functions are installed in the directory .Src (the default directory for Splus scripts), while the example script plot.hourly.wb.s is put into the current directory. The example script compares the flux and snow outputs from two model runs, with time series plots of the data, time series plots of the differences, and scatter plots.
Descriptions of the Splus functions are included in the table below. Hopefully they will provide the user with a frame work with which they can build their own routines.
| Function Name | Description |
|---|---|
| compare.daily.variables.s | This function reads results from two daily water balance model simulations and produces a postscript file with time series of simulated values, time series of the differences between simulated values, and scatter plots of the simulated values. |
| compare.hourly.variables.s | This function reads results from two sub-daily energy balance model simulations and produces a postscript file with time series of simulated values, time series of the differences between simulated values, and scatter plots of the simulated values. |
| func.plot.compare.nottime.s | This function takes two time series of data and produces two plots: the first is a time series of both sets of data, the second is a time series of their differences (first - second). |
| func.plot.scatter.compare.s | This function takes two time series of data and produces a scatter plot of the data values. |
| func.read.flux.files.s | This function reads a VIC model output flux file, and returns the time series of each variable as components in a list. e.g. the first 10 evaporation values for data set tmp.data can be indexed as tmp.data$evap[1:10]. |
| func.read.snow.files.s | This function reads a VIC model output snow file, and returns the time series of each variable as components in a list. e.g. the first 10 snow depth values for data set tmp.data can be indexed as tmp.data$depth[1:10]. |
