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)
b=gm('xbxrb/x.01m/',m=m,sc=sc,/sub)
;c=gm('ctl3/x.01m',m=m,sc=sc)
c=gm('seaice-motion-fowler/x.01m',m=m,sc=sc)
d=gm('seaice-motion-eis/x.01m',m=m,sc=sc,/sub)

d=wpp_ff(d,get_mask(/sea,/rmg),1,/nom)
d.data(where(d.data eq d.bmdi))=0

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

pp_plot,pp_mag(d),xmarg=xm,ymarg=ym,abc='Obs: EIS',lev=lev,nh=keyword_set(nh),fil=1,/top
pp_plot,pp_mag(c),xmarg=xm,ymarg=ym,abc='Obs: Fowler',lev=lev,nh=keyword_set(nh),fil=1,/top
pp_plot,pp_mag(b),xmarg=xm,ymarg=ym,abc='HiGEM',lev=lev,nh=keyword_set(nh),fil=1,/top
pp_plot,pp_mag(a),xmarg=xm,ymarg=ym,abc='HadGEM',lev=lev,nh=keyword_set(nh),fil=1,/top

gettwogifs,out=file

end

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

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

sc=[148,149]

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

end