Utilities - Revision history

Jchoate: /* RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA */

2012-07-12T05:33:06Z

RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA

2012-07-12T05:29:53Z

RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA

2012-07-12T05:27:08Z

RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA

2012-07-12T05:26:32Z

RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA

2012-07-12T05:20:23Z

PATCH DAILY WATER BALANCE FUNCTIONS

2012-07-12T05:18:48Z

2012-03-28T18:03:17Z

BASIN DAILY WATER BALANCE FUNCTION

2011-01-19T17:55:09Z

PATCH DAILY WATER BALANCE FUNCTIONS

2011-01-19T17:52:33Z

UPDATE A WORLDFILE FOR COMPATIBILITY WITH RHESSYS VERSION 5.14

2011-01-19T17:51:46Z

BASIN DAILY WATER BALANCE FUNCTION

@@ Line 245: / Line 245: @@
 cyr = new[i]                            # checks what yr for var seq(?)
 tmpn = subset(clim,clim$wy==cyr)        # extracts that var/new yr f/orig clim
 lnn = nrow(tmpn)                        # n days in extracted yr (leap?)
 oyr = orig[i]                           # checks what year for orig seq
 tmpo = subset(clim, clim$wy==oyr)       # extracts that "in seq" yr f/orig clim
 lno = nrow(tmpo)                        # n days in orig yr (leap?)

@@ Line 205: / Line 205: @@
 orig = unique(clim$wy)
 nyrs = length(orig)
 # determine when leap years are
 dyrs = aggregate(clim$wyd, by=list(clim$wy), max)
 colnames(dyrs) = c("wyd","nday")
 leaps = subset(dyrs, dyrs$nday>365)
+# initialize a vector called "new" of years that you want climate for
 # you can make this however you want
 # here we show how to randomly sample from original years
 new = sample(orig, replace=T, size=nyrs)

← Older revision		Revision as of 05:27, 12 July 2012
Line 197:		Line 197:

	clim = read_rhessys_met("../../clim/sage")		clim = read_rhessys_met("../../clim/sage")
		+
	clim=mkdate(clim)		clim=mkdate(clim)

@@ Line 193: / Line 193: @@
 == RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA ==
-# read rhesus climate data into a data file, for this example we
+# read RHESSys climate data (you can use the read_rhessys_met script) into an R data file, for this example we call this new R object "clim"
-# call this clim
+# you will also need to run the R script mkdate on the R object "clim" so all of the needed columns are included in the R object
 clim = read_rhessys_met("../../clim/sage")
 clim=mkdate(clim)
-#  note you must have complete years or water years (depending on which you will use) for
+#  note: you must have complete calendar years or water years (depending on which you will use) for selection
-# selection, if not get rid of incomplete years, using subset eg. clim = subset(clim, clim$wy > 1953),
+# if years are not complete - get rid of incomplete years (you can use subset, i.e.: clim = subset(clim, clim$wy > 1953),  this would get rid of first year 1953 if it didn't include a full wateryear
-# would get rid of first year 1953 if it didn't have a full year
+# this may occur at beginning and end of the record
 orig = unique(clim$wy)

@@ Line 191: / Line 191: @@
+}
-== PATCH DAILY WATER BALANCE FUNCTIONS ==
+== RESAMPLE EXISTING RHESSys METEOROLOGICAL DATA ==
 # read rhesus climate data into a data file, for this example we

← Older revision		Revision as of 05:18, 12 July 2012
Line 1:		Line 1:
−

	== UPDATE A WORLDFILE FOR COMPATIBILITY WITH RHESSYS VERSION 5.14 ==		== UPDATE A WORLDFILE FOR COMPATIBILITY WITH RHESSYS VERSION 5.14 ==
Line 191:		Line 190:

	}		}
		+
		+	== PATCH DAILY WATER BALANCE FUNCTIONS ==
		+
		+	# read rhesus climate data into a data file, for this example we
		+	# call this clim
		+	# you will also need to run mkdate on the clim input so we have all of the columns that we need
		+	# "la" was the prefix for the RHESSys input files
		+
		+
		+	clim = read_rhessys_met("../../clim/sage")
		+	clim=mkdate(clim)
		+
		+	# note you must have complete years or water years (depending on which you will use) for
		+	# selection, if not get rid of incomplete years, using subset eg. clim = subset(clim, clim$wy > 1953),
		+	# would get rid of first year 1953 if it didn't have a full year
		+	# this occurs at beginning and end of the record
		+
		+	orig = unique(clim$wy)
		+	nyrs = length(orig)
		+	# determine when leap years are
		+	dyrs = aggregate(clim$wyd, by=list(clim$wy), max)
		+	colnames(dyrs) = c("wyd","nday")
		+	leaps = subset(dyrs, dyrs$nday>365)
		+
		+
		+	# initialize a vector new of years that you want clim for
		+	# you can make this however you want
		+	# here we show how to randomly sample from original years
		+
		+	new = sample(orig, replace=T, size=nyrs)

@@ Line 32: / Line 32: @@
 rb$sddiff=tmp
-tmp=diff(rb$snow)
+tmp=diff(rb$snowpack)
 tmp=c(0,tmp)

@@ Line 158: / Line 158: @@
+}
-MULTIPATCH RUN FUNCTION
+PASTE THE FOLLOWING COMMANDS INTO R FOR THE MULTIPATCH RUN FUNCTION
 pwatbalmult = function(rp,rc) {

@@ Line 12: / Line 12: @@
-To apply to an existing worldfile and create a new worldfile compatible with rhessys5.14
+To apply to an existing worldfile and create a new worldfile compatible with rhessys5.14, type the following command:
 awk -f add_litfrac_lines.awk < worldfile_to_change > new_worldfile_name
 == BASIN DAILY WATER BALANCE FUNCTION ==

@@ Line 1: / Line 1: @@
 == BASIN DAILY WATER BALANCE FUNCTION ==
@@ Line 62: / Line 80: @@
 (on the order of 10^-6).  If your watbal values are negative then
 water inputs are less than water outputs, and vice versa.
 == PATCH DAILY WATER BALANCE FUNCTIONS ==