Quick Hacks - Revision history

Root: /* Save to LaTex */

2011-02-23T08:18:37Z

Save to LaTex

← Older revision		Revision as of 08:18, 23 February 2011
Line 111:		Line 111:
	set terminal push		set terminal push
	set terminal tikz standalone header "\\usepackage{txfonts}"		set terminal tikz standalone header "\\usepackage{txfonts}"
	set output "$0.tex		set output "$0.tex"
	replot		replot
	set output		set output
Line 121:		Line 121:
	set terminal push		set terminal push
	set terminal tikz standalone header "\\usepackage{txfonts}"		set terminal tikz standalone header "\\usepackage{txfonts}"
	set output "$0.tex		set output "$0.tex"
	call "$0.gp"		call "$0.gp"
	set output		set output

Root: /* Multiple Plots */

2011-02-22T22:00:58Z

Multiple Plots

@@ Line 137: / Line 137: @@
 ===Multiple Plots===
 <source lang=gnuplot>
-call "latex1.gp" "<name>" # where <name>.tex is the output file.
+set xdata time
 </source>

Root: /* Multiple Plots */

2011-02-22T22:00:27Z

Multiple Plots

← Older revision		Revision as of 22:00, 22 February 2011
Line 137:		Line 137:
	===Multiple Plots===		===Multiple Plots===
	<source lang=gnuplot>		<source lang=gnuplot>
	~~<pre>~~
	call "latex1.gp" "<name>" # where <name>.tex is the output file.		call "latex1.gp" "<name>" # where <name>.tex is the output file.
	~~</pre>~~
	</source>		</source>

Root: /* Multiple Plots with the script */

2011-02-22T21:59:53Z

Multiple Plots with the script

@@ Line 104: / Line 104: @@
 plot "<awk -v intv=168 -f tmp.awk cntr_hly_99.csv.dat" using 1:2 title "weekly" with boxes fs solid
 unset multiplot
 </source>

Root at 21:34, 22 February 2011

2011-02-22T21:34:48Z

@@ Line 82: / Line 82: @@
 EOF
 rm -f $$.awk
 </source>

Root at 16:16, 22 February 2011

2011-02-22T16:16:28Z

← Older revision		Revision as of 16:16, 22 February 2011
Line 35:		Line 35:
	awk -v intv=24 -f accum.awk cntr_hly_1.csv.dat		awk -v intv=24 -f accum.awk cntr_hly_1.csv.dat
	</pre>		</pre>

			===Plotting those data with GNUplot===
			<source lang=bash>
			#!/bin/bash
			if [ $# -lt 2 ]; then
			echo "Usage: $0 <timeseriesfile> <accumulation> [<startdate>] [<enddate>]"
			exit 5
			fi
			file=$1
			intv=$2
			startd=`head -n1 $file\|awk '{print $1}'`
			endd=`tail -n1 $file\|awk '{print $1}'`
			if [ $# -ge 3 ]; then
			startd=$3
			endd=$4
			fi
			cat << EOF > $$.awk
			#accum.awk
			BEGIN{
			ct=0
			avg=0.0
			}
			{
			avg+=\$2
			if(++ct%intv==0){
			x=\$1
			print x,avg/ct
			ct=0
			avg=0.0
			}

			}
			EOF
			gnuplot -persist << EOF
			set xdata time
			set timefmt x "%m/%d/%Y-%H:%M"
			set key title ""
			set xlabel "Time"
			set xlabel offset character 0, 0, 0 font "" textcolor lt -1 norotate
			set xrange [ "$startd" : "$endd" ] noreverse nowriteback
			set ylabel "Rainfall (mm)"
			set ylabel offset character 0, 0, 0 font "" textcolor lt -1 rotate by -270
			#set boxwidth 1
			plot "<awk -v intv=$intv -f $$.awk $file" using 1:2 with boxes fs solid
			# EOF
			EOF
			rm -f $$.awk
			</source>

Root: /* Awk script for accumulating time series data */

2011-02-22T12:52:15Z

Awk script for accumulating time series data

← Older revision		Revision as of 12:52, 22 February 2011
Line 3:		Line 3:

	==Awk script for accumulating time series data==		==Awk script for accumulating time series data==
	<source lang=~~awk~~>		<source lang=bash>
	#accum.awk		#accum.awk
	BEGIN{		BEGIN{

Root: Created page with '==What is this== These are primarily notes for myself. Things that I often need to do in computing. Instead of keeping them in a private location I list them here. Perhaps someon…'

2011-02-22T12:51:13Z

Created page with '==What is this== These are primarily notes for myself. Things that I often need to do in computing. Instead of keeping them in a private location I list them here. Perhaps someon…'

New page

==What is this==
These are primarily notes for myself. Things that I often need to do in computing. Instead of keeping them in a private location I list them here. Perhaps someone else will also find a use..

==Awk script for accumulating time series data==
<source lang=awk>
#accum.awk
BEGIN{
ct=0
avg=0.0
}
{
if(++ct%intv==0){
print x,avg/ct
ct=0
avg=0.0
}else{
if(ct%intv==1){
x=$1
}
avg+=$2
}
}
</source>
Data:
<source lang=text>
01/01/1945-00:00 0.0
01/01/1945-01:00 0.0
01/01/1945-02:00 0.0
01/01/1945-03:00 0.0
...
</source>

Running:
<pre>
awk -v intv=24 -f accum.awk cntr_hly_1.csv.dat
</pre>