
set grads off

set parea 1 10 0.5 8

set clopts 1 4 0.12
set cthick 3
set clab off
set xlopts 1 4 0.18
set ylopts 1 4 0.18

set gxout shaded

set lat -90 90
set lon -180 180
!set mproj robinson

set display color white

psl1c = ave(psl.1, t=4, t=6)/100.
psl2c = ave(psl.1, t=16, t=18)/100.
psl3c = ave(psl.1, t=28, t=30)/100.
psl4c = ave(psl.1, t=40, t=42)/100.
psl5c = ave(psl.1, t=52, t=54)/100.
psl6c = ave(psl.1, t=64, t=66)/100.
psl7c = ave(psl.1, t=76, t=78)/100.
psl8c = ave(psl.1, t=88, t=90)/100.
psl9c = ave(psl.1, t=100, t=102)/100.
psl10c = ave(psl.1, t=112, t=114)/100.

psl1ng = ave(psl.2, t=4, t=6)/100.
psl2ng = ave(psl.2, t=16, t=18)/100.
psl3ng = ave(psl.2, t=28, t=30)/100.
psl4ng = ave(psl.2, t=40, t=42)/100.
psl5ng = ave(psl.2, t=52, t=54)/100.
psl6ng = ave(psl.2, t=64, t=66)/100.
psl7ng = ave(psl.2, t=76, t=78)/100.
psl8ng = ave(psl.2, t=88, t=90)/100.
psl9ng = ave(psl.2, t=100, t=102)/100.
psl10ng = ave(psl.2, t=112, t=114)/100.

diff1 = psl1c - psl1ng
diff2 = psl2c - psl2ng
diff3 = psl3c - psl3ng
diff4 = psl4c - psl3ng
diff5 = psl5c - psl5ng
diff6 = psl6c - psl6ng
diff7 = psl7c - psl7ng
diff8 = psl8c - psl8ng
diff9 = psl9c - psl9ng
diff10 = psl10c - psl10ng

sqrtn = pow(10,0.5)

aver = diff1+diff2+diff3+diff4+diff5+diff6+diff7+diff8+diff9+diff10
aver = aver/10.

numer = pow(diff1-aver,2) + pow(diff2-aver,2) + pow(diff3-aver,2)
numer = numer + pow(diff4-aver,2) + pow(diff5-aver,2) + pow(diff6-aver,2)
numer = numer + pow(diff7-aver,2) + pow(diff8-aver,2) + pow(diff9-aver,2)
numer = numer + pow(diff10-aver,2)
stdev = numer/(9.)
stdev = pow(stdev,0.5)

tdist = aver*sqrtn/stdev

set gxout shaded
set clevs -3.250 -2.262 -1.833 -1.383 1.383 1.833 2.262 3.250
d tdist
run cbarn

set gxout contour
set clevs -3.250 -2.262 -1.833 -1.383 1.383 1.833 2.262 3.250
d tdist


