! cd /aosc/horse/carton/fluxes/
*
* compute monthly mean fields
*
!cdo monmean /aosc/greenland/forcing/era_i/era_i_ssrd_1dy_19790101_20181231.nc era_i_ssrd_1x1_1mn.nc
!cdo monmean /aosc/greenland/forcing/era_i/era_i_strd_1dy_19790101_20181231.nc era_i_strd_1x1_1mn.nc
!cdo seldate,1980-01-01,2010-12-31 /aosc/greenland/forcing/CERA-20C/combined/CERA20C_ssrd_1901_2010.nc CERA20C_ssrd_1980_2010_dy.nc
!cdo monmean CERA20C_ssrd_1980_2010_dy.nc  CERA20C_ssrd_1980_2010_mn.nc
!cdo seldate,1980-01-01,2010-12-31 /aosc/greenland/forcing/CERA-20C/combined/CERA20C_strd_1901_2010.nc CERA20C_strd_1980_2010_dy.nc
!cdo monmean CERA20C_strd_1980_2010_dy.nc CERA20C_strd_1980_2010_mn.nc 
!rm CERA20C_ssrd_1980_2010_dy.nc CERA20C_strd_1980_2010_dy.nc
*
* Compare downwelling shortwave and longwave from a new reanalysis with ERA-I.  
*
* The upper two panels compare the 
* LEFT:  normalized difference in record means and 
* RIGHT: normalized difference in standard deviations.
* 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 shortwave
*
sdfopen   CERA20C_ssrd_1980_2010_mn.nc
mean = ave( ssrd,time=jan1980,time=dec2018)
var = ave(( ssrd-mean)*( ssrd-mean),time=jan1980,time=dec2018))
set lon 162
set lat 42
set t 1 last
ts =  ssrd
close 1
sdfopen /aosc/horse/carton/fluxes/era_i/era_i_ssrd_1x1_1mn.nc
set lon 162
set lat 42
set t 1 last
tss = ssrd
set  lon 0 360
set lat -90 90
set t 1
means = ave(ssrd,time=jan1980,time=dec2018)
vars = ave((ssrd-means)*(ssrd-means),time=jan1980,time=dec2018))
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 -20 -15 -10 -5 5 10 15 20 
set ccols 47 45 43 42 0 61 63 65 67
d means2
basemap L 15 15
cbarn
draw title M_1-M_era
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 -500 -250 -100  100 250 500
set ccols 45 43 42 0 61 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 10.8  SHORTWAVE DOWN
set x 1
set y 1
set time jan1980 dec2018
set vpage 0 4.25 6.5 8.75
set cmark 0
set ccolor 1
set grads off
set grid off
set axlim 00 300
set yaxis 00 300 100
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab on
set ylpos 0 l 
set cthick 7
d ts
set axlim 00 300
set yaxis 00 300 100
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab on
set ylpos 0 l 
set ccolor 2
set cmark 0
set cthick 2
d tss
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 longwave
*
sdfopen  CERA20C_strd_1980_2010_mn.nc  
mean = ave( strd,time=jan1980,time=dec2018)
var = ave(( strd-mean)*( strd-mean),time=jan1980,time=dec2018))
set lon 162
set lat 42
set t 1 last
ts =  strd
close 1
sdfopen /aosc/horse/carton/fluxes/era_i/era_i_strd_1x1_1mn.nc
set lon 162
set lat 42
set t 1 last
tss = strd
set  lon 0 360
set lat -90 90
set t 1
means = ave(strd,time=jan1980,time=dec2018)
vars = ave((strd-means)*(strd-means),time=jan1980,time=dec2018))
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 -20 -15 -10 -5 5 10 15 20 
set ccols 47 45 43 42 0 61 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 -300 -200 -100 -50 50 100 200 300
set ccols 47 45 43 42 0 61 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  LONGWAVE DOWN
set x 1
set y 1
set time jan1980 dec2018
set vpage 0 4.25 1.5 3.75
set cmark 0
set ccolor 1
set grads off
set grid off
set axlim 220 420
set yaxis 220 420 100
set xlopts 1 3 .15
set ylopts 1 3 .15
set ylab on
set ylpos 0 l 
set cthick 6
d ts
*set vpage 3.75 8 1.5 3.75
set cmark 0
set axlim 220 420
set yaxis 220 420 100
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/radiative_flux-comparison.pdf