;
; Plots maps of correlations between MXD and Apr-Sep temperature and
; May-Aug precipitation.
;
dotem=0         ; 0=use CRUTEM1, 1=use CRUTEM2v
if dotem eq 0 then begin
  savtem=''
endif else begin
  savtem='ct2v_'
endelse
;
; 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='vlgrey'
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
coast.fill=1 & coast.fillcolor=30
labels=def_labels(/off)
;
;trv=0           ; selects tree-ring-variable: 0=MXD 1=TRW
for trv = 0 , 0 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='datastore/'+savtem+fnadd+'_moncorr.idlsave'
    1: restore,filename='datastore/'+savtem+fnadd+'_moncorr_early.idlsave'
    2: restore,filename='datastore/'+savtem+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(10)
  ;
  bl=where(oner ge 0.22,nblack)
  if nblack gt 0 then plots,statlon(bl),statlat(bl),psym=def_sym(15)
  ;
  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