;
; Plots maps of correlations between MXD and Apr-Sep temperature and
; May-Aug precipitation.
;
; Now prepare for plotting
;
loadct,39
multi_plot,nrow=2,layout='large'
if !d.name eq 'X' then begin
  window,ysize=850
  !p.font=-1
endif else begin
  !p.font=0
  device,/helvetica,/bold,font_size=11
endelse
def_1color,19,color='white'
def_1color,20,color='green'
def_1color,21,color='deepblue'
def_1color,22,color='purple'
def_1color,23,color='red'
def_1color,24,color='black'
def_1color,30,color='mgrey'
cpl_usersym,/circle,/fill
;
map=def_map(/npolar)  &  map.limit(0)=25.
map.xmargin=[0,0] & map.ymargin=[0,0]
coast=def_coast(/get_device) & coast.double=0
labels=def_labels(/off)
;
;trv=0           ; selects tree-ring-variable: 0=MXD 1=TRW
for trv = 0 , 1 do begin
case trv of
  0: begin
    fnadd='mxd'
    iseas=18
    end
  1: begin
    fnadd='trw'
    iseas=20
    end
endcase
titadd=strupcase(fnadd)
;
; Repeat for full, early and late periods
;
for iper = 0 , 0 do begin
  ;
  ; Restore all correlations
  ;
  case iper of
    0: restore,filename=fnadd+'_moncorr.idlsave'
    1: restore,filename=fnadd+'_moncorr_early.idlsave'
    2: restore,filename=fnadd+'_moncorr_late.idlsave'
  endcase
  ;  allr,ncorr,nvar,nchron,varname,corrname,statlat,statlon
  ;
  mtit='   (1881-1984)'
  ;
  ; Plot temperature correlations
  ;
  inter_boxfd,/nodata,map=map,coast=coast,labels=labels
;    title='r( MXD , April-September temp.)'+mtit
  ;
  oner=reform(allr(*,iseas,1))
  ;
  if !d.name eq 'PS' then device,font_size=10
  ;
  cpl_usersym,/circle,/fill
  plots,statlon,statlat,psym=8,color=19
  ;
  wl=where(oner lt 0.22,nwhite)
  if nwhite gt 0 then plots,statlon(wl),statlat(wl),psym=def_sym(15)
  ;
  gl=where((oner ge 0.22) and (oner lt 0.4),ngreen)
  if ngreen gt 0 then plots,statlon(gl),statlat(gl),psym=def_sym(1)
  ;
  pl=where((oner ge 0.4) and (oner lt 0.5),npurple)
  if npurple gt 0 then plots,statlon(pl),statlat(pl),psym=def_sym(2)
  ;
  rl=where((oner ge 0.5) and (oner lt 0.6),nred)
  if nred gt 0 then plots,statlon(rl),statlat(rl),psym=def_sym(16)
  ;
  bl=where(oner ge 0.6,nblack)
  if nblack gt 0 then plots,statlon(bl),statlat(bl),psym=def_sym(10)
  ;
  if !d.name eq 'PS' then device,font_size=11
;
legpos=convert_coord([-143.],[26.],/data,/to_normal)
  ;
if !d.name eq 'PS' then device,/courier,/bold,font_size=10
legend,$
  ['r('+titadd+','+corrname(iseas)+' temp.) 1881-1984',$
   'r <  +0.22         (N='+string(nwhite,format='(I3)')+')',$
   'r =  +0.22 to +0.4 (N='+string(ngreen,format='(I3)')+')',$
   'r =  +0.4  to +0.5 (N='+string(npurple,format='(I3)')+')',$
   'r =  +0.5  to +0.6 (N='+string(nred,format='(I3)')+')',$
   'r >  +0.6          (N='+string(nblack,format='(I3)')+')'],$
  linestyle=[-1,0,0,0,0,0],psym=[-1,15,1,2,16,10],position=legpos,$
  charsize=0.85,spacing=1.25,margin=[0.45]
if !d.name eq 'PS' then device,/helvetica,/bold,font_size=11
;
endfor
;
endfor
;
end