;
g=def_grid(/ipcc)
;
multiplot,4
!p.multi(1)=2
if !d.name eq 'X' then begin
  window,xsize=650,ysize=950
endif else begin
  device,xsize=15,xoffset=3
  device,ysize=24,yoffset=2
endelse
;
wlims=[0.,0.,-10.,-90.,-90.,-180.,-180.,-180.,-180.,-180.,-180.]
elims=[180.,180.,170.,90.,90.,0.,0.,0.,180.,180.,180.]
oxs=0.5*(wlims+elims)
;
map=def_map(/npolar)  &  map.limit(0)=30.
coast=def_coast(/get_device)  &  coast.double=0
labels=def_labels(/off)  &  labels.gridon=1  &  labels.dlon=90.
;
restore,filename='regboxes.idlsave'
for i = 0 , nreg-1 do begin
;
pause
fd1=boxlists(*,*,i)
dummy=where(finite(fd1),nboxxx)
map.limit(1)=wlims(i)
map.limit(3)=elims(i)
map.origx=oxs(i)
inter_boxfd,fd1,g.x,g.y,map=map,coast=coast,labels=labels,$
  title=regname(i)+'  n='+string(nboxxx,format='(I3)'),$
  levels=[0.5,1.5],c_colors=[100.]
;
endfor
;
end