* this is where I've been making the comparisons:
*
* cd /aosc/horse/carton/fluxes/
*
* compute monthly mean fields
*
!cdo monmean /aosc/greenland/forcing/era_i/era_i_q_6hr_19790101_20181231.nc era_i_q_1x1_1mn.nc
!cdo monmean /aosc/greenland/forcing/era_i/era_i_t2m_6hr_19790101_20181231.nc era_i_t2m_1x1_mn.nc
!cdo monmean /aosc/greenland/forcing/CERA-20C/combined/CERA20C_q_1901_2010.nc CERA20C_q_1901_2010_mn.nc
!cdo monmean /aosc/greenland/forcing/CERA-20C/combined/CERA20C_t2m_1901_2010.nc CERA20C_t2m_1901_2010_mn.nc 
*
* Compare downwelling shortwave and longwave from a new reanalysis with ERA-I.  
*
* The upper two panels compare the 
* LEFT:  difference in record means and 
* RIGHT:  difference in monthly variances.
* The lower panel compares a one-year long time series at a point
* (42N, 163E) in the central North Pacific
reinit
set display color white
c
* load the color table
/homes/metofac/carton/colors
*
* plot q2m
*
sdfopen  CERA20C_q_1901_2010_mn.nc
mean = ave(q,time=jan1980,time=dec2017)
var = ave((q-mean)*(q-mean),time=jan1980,time=dec2017))
set lon 162
set lat 42
set t 1 last
ts = q
close 1
sdfopen /aosc/horse/carton/fluxes/era_i/era_i_q_1x1_1mn.nc
set lon 162
set lat 42
set t 1 last
tss = q
set  lon 0 360
set lat -90 90
set t 1
means = ave(q,time=jan1980,time=dec2017)
vars = ave((q-means)*(q-means),time=jan1980,time=dec2017))
means2 = (lterp(mean,means)-means)
vars2 = (lterp((var),vars)-vars)
set vpage 0 4.5 8.5 10.5
set xlopts 1 3 .15
set ylopts 1 3 .15
set grads off
set grid off
set cmark 0
set gxout shaded
set clevs -2.0 -1. -.5 -.2 .2 .5 1. 2.0 
set ccols 47 45 43 42 0 61 63 65 67
d means2*1.e3
basemap L 15 15
cbarn
draw title M_1-M_era (x10^3)
set ccolor 1
set vpage 4 8.5 8.5 10.5
set xlopts 1 3 .15
set ylopts 1 3 .15
set grads off
set grid off
set cmark 0
set clevs  -2. -1 -.5 -.1 .1 .5 1 2  
set ccols 47 45 43 42 0 61 63 65 67
d vars2*1.e6
basemap L 15 15
cbarn
draw title VAR_1-VAR_era (x10^6)
set vpage off
set strsiz .14
set string 1 c 4 0
draw string 4.25 10.8  SPECIFIC HUMIDITY (kg/kg)
set x 1
set y 1
set time jan1980 dec2019
set vpage 0 4.25 6.5 8.75
set cmark 0
set ccolor 1
set cthick 8
set grads off
set grid off
set axlim 00 15
set yaxis 00 15 5
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab on
set ylpos 0 l 
d ts*1.e3
set axlim 00 15
set yaxis 00 15 5
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab on
set ylpos 0 l 
set cthick 2
set ccolor 2
set cmark 0
d tss*1.e3
set vpage off
set strsiz .13
set string 1 c 4 0
draw string 2.45606 8.55556 MONTHLY 42N,162E
set string 2 l 3 0
set strsiz .09
draw string 2.4787 8.21605 -- ERA-I
set line 1 1 1
draw line 2 6.3 7 6.3
close 1
*
* plot T2m
*
sdfopen  CERA20C_t2m_1901_2010_mn.nc 
mean = ave(t2m,time=jan1980,time=dec2017)
var = ave((t2m-mean)*(t2m-mean),time=jan1980,time=dec2017))
set lon 162
set lat 42
set t 1 last
ts = t2m
close 1
sdfopen /aosc/horse/carton/fluxes/era_i/era_i_t_2m_1x1_mn.nc
set lon 162
set lat 42
set t 1 last
tss = t2m
set  lon 0 360
set lat -90 90
set t 1
means = ave(t2m,time=jan1980,time=dec2017)
vars = ave((t2m-means)*(t2m-means),time=jan1980,time=dec2017))
means2 = (lterp(mean,means)-means)
vars2 = (lterp((var),vars)-vars)
set vpage 0 4.5 3.5 5.5
set xlopts 1 3 .15
set ylopts 1 3 .15
set grads off
set grid off
set cmark 0
set gxout shaded
set clevs  -2  -1 -.5 .5 1 2  
set ccols  47 45 43  0  63 65 67
d means2
basemap L 15 15
cbarn
draw title M_1-M_era
set ccolor 1
set vpage 4 8.5 3.5 5.5
set xlopts 1 3 .15
set ylopts 1 3 .15
set grads off
set grid off
set cmark 0
set clevs  -4  -2 -1 1 2 4  
set ccols  47 45 43  0  63 65 67
d vars2
basemap L 15 15
cbarn
draw title VAR_1-VAR_era
set vpage off
set strsiz .14
set string 1 c 4 0
draw string 4.25 5.8  T2m
set x 1
set y 1
set time jan1980 dec2019
set vpage 0 4.25 1.5 3.75
set cmark 0
set cthick 8
set grads off
set grid off
set axlim 270 300 
set yaxis 270 300 10
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab on
set ylpos 0 l 
d ts
*set vpage 3.75 8 1.5 3.75
set cmark 0
set axlim 270 300 
set yaxis 270 300 10
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab off
set ccolor 2
set cthick 2
set grads off
set grid off
d tss
set vpage off
set strsiz .13
set string 1 c 4 0
draw string  2.45606 3.4856 42N,162E
set string 2 l 3 0
set strsiz .09
draw string  1.08655 3.27058 -- ERA-I
gxprint /aosc/greenland/forcing/CERA-20C/combined/q2_T2_flux-comparison.pdf