#!/bin/tcsh
# plots the routed monthly flows for the Colorado basin against naturalized flows

set PROJ = -JX5.5/2
set PROJ2 = -JX5.5/1.5
set COORD = -R1917/2002/0/12000
set OUTF = figureExtraHydroandClimate.ps
set YSHFT = -2.5
set YSHFT2 = -2

gmtset BASEMAP_TYPE plain

set temp = ./DATA/ExtraEnergy.txt
tail +0 $temp >! STRG
gmtset LABEL_FONT_SIZE 8 ANOT_FONT_SIZE 8 ANOT_FONT Helvetica  \
LABEL_FONT Helvetica ANOT_FONT_SIZE 10 LABEL_FONT_SIZE 10 TICK_LENGTH 0.025

#open file
pstext -N -K -Y9.4 -P -X1.2 $COORD $PROJ << END >! $OUTF
#1950 50000 10 0 5 2 "Simulated Hydropower Generation in MegaW-hr"
END

##################################################################

# First plot, extra energy versus years
awk '{print $1, $2}' STRG | \
psxy -P -O -K -Y$YSHFT -B10f500/2000\:"Extra hydropower, AMW - hr":SWne $COORD $PROJ  -W3 >> $OUTF

# Second plot, ENSO signals
##### plot positive anomaly in red
awk '{print $1, 0, 2-$3}' STRG | \
pswiggle -O -K -P -V -Y$YSHFT2 -Z1.5 -G250/0/0 -T  -B10f10/1\:"ENSO":SWne -R1917/2002/-1.5/1.5 $PROJ2 >> $OUTF
##### plot negative anomalies in blue
awk '{print $1, 0, 2-$3}' STRG | \
pswiggle -O -K -P -V -Z1.5 -N -G0/0/250  -T -W0.25 -R1917/2002/-1.5/1.5 $PROJ2 >> $OUTF

# Third plot, idem as above for PDO
awk '{print $1, 0, 2-$4}' STRG | \
pswiggle -O -K -P -V -Y$YSHFT2 -Z1.5 -G250/0/0 -T -B10f10/1\:"PDO":SWne -R1917/2002/-1.5/1.5 $PROJ2 >> $OUTF
awk '{print $1, 0, 2-$4}' STRG | \
pswiggle -O -K -P -V -Z1.5 -N -G0/0/250 -T -W0.25  -R1917/2002/-1.5/1.5 $PROJ2 >> $OUTF

pstext -N -O  $COORD $PROJ << END >> $OUTF
#1950 35000  14 0 5 6 PNW Simulated Extra Hydropower 
END


gs $OUTF
rm STRG
rm temp