pro p1a,a,b,c,d,file,lev,nh=nh,sc=sc,m=m

a=gm('adzbf/x.01m',m=m,sc=sc)
b=gm('eacew/x.01m',m=m,sc=sc)
c=gm('ctl3/x.01m',m=m,sc=sc)
d=gm('seaice-comiso/x.01m',m=m,sc=sc)

xm=[0,0]
ym=[3,2]

pp_plot,d,xmarg=xm,ymarg=ym,abc='Obs',lev=lev,nh=keyword_set(nh)
pp_plot,c,xmarg=xm,ymarg=ym,abc='CM3',lev=lev,nh=keyword_set(nh)
pp_plot,b,xmarg=xm,ymarg=ym,abc='HiGEM',lev=lev,nh=keyword_set(nh)
pp_plot,a,xmarg=xm,ymarg=ym,abc='HadGEM',lev=lev,nh=keyword_set(nh)

gettwogifs,out=file

end

; -------------------------

!p.multi=[0,2,2]

sc=146
lev=[0.15,.5,.8,.9,.95]

p1a,a,b,c,d,'all4.sep.sh',lev,m=9,sc=sc
p1a,a,b,c,d,'all4.sep.nh',lev,/nh,m=9,sc=sc
p1a,a,b,c,d,'all4.mar.sh',lev,m=3,sc=sc
p1a,a,b,c,d,'all4.mar.nh',lev,/nh,m=3,sc=sc

end