; ; Plots maps of the change in correlation between MXD and precip or temperature ; from the early to the late period ; ; Now prepare for plotting ; loadct,39 multi_plot,nrow=1,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=14 endelse def_1color,20,color='red' def_1color,21,color='blue' def_1color,22,color='green' def_1color,23,color='black' def_1color,24,color='orange' def_1color,25,color='purple' cpl_usersym,/circle,/fill ; map=def_map(/npolar) & map.limit(0)=25. coast=def_coast(/get_device) & coast.double=1 labels=def_labels(/off) ; ; Restore all correlations ; restore,filename='mxd_moncorr_early.idlsave' ; allr,ncorr,nvar,nchron,varname,corrname,statlat,statlon,allp,moir,moir2 tearly=reform(allr(*,18,1)) pearly=reform(allp(*,19,0)) restore,filename='mxd_moncorr_late.idlsave' ; allr,ncorr,nvar,nchron,varname,corrname,statlat,statlon,allp,moir,moir2 ; pause inter_boxfd,/nodata,map=map,coast=coast,labels=labels,$ title='Change in r(MXD,April-September temperature)' ; oner=reform(allr(*,18,1))-tearly ; if !d.name eq 'PS' then device,font_size=11 ; ol=where(oner lt -0.15,norange) if norange gt 0 then plots,statlon(ol),statlat(ol),psym=8,color=24 ; rl=where((oner ge -0.15) and (oner lt -0.05),nred) if nred gt 0 then plots,statlon(rl),statlat(rl),psym=8,color=20 ; bl=where((oner ge -0.05) and (oner lt 0.),nblue) if nblue gt 0 then plots,statlon(bl),statlat(bl),psym=8,color=21 ; gl=where((oner ge 0.) and (oner lt 0.05),ngreen) if ngreen gt 0 then plots,statlon(gl),statlat(gl),psym=8,color=22 ; bl=where((oner ge 0.05) and (oner lt 0.15),nblack) if nblack gt 0 then plots,statlon(bl),statlat(bl),psym=8,color=23 ; pl=where(oner ge 0.15,npurple) if npurple gt 0 then plots,statlon(pl),statlat(pl),psym=8,color=25 ; if !d.name eq 'PS' then device,font_size=14 ; legpos=convert_coord([-160.],[35.],/data,/to_normal) if !d.name eq 'PS' then device,/courier,/bold legend,$ ['del-r <= -0.15 (N='+string(npurple,format='(I3)')+')',$ 'del-r = -0.15 to -0.05 (N='+string(nblack,format='(I3)')+')',$ 'del-r = -0.05 to +0.00 (N='+string(ngreen,format='(I3)')+')',$ 'del-r = +0.00 to +0.05 (N='+string(nblue,format='(I3)')+')',$ 'del-r = +0.05 to +0.15 (N='+string(nred,format='(I3)')+')',$ 'del-r > +0.15 (N='+string(norange,format='(I3)')+')'],$ psym=[8,8,8,8,8,8],color=[25,23,22,21,20,24],position=legpos if !d.name eq 'PS' then device,/helvetica,/bold ; oner=reform(allp(*,19,0))-pearly ; pause inter_boxfd,/nodata,map=map,coast=coast,labels=labels,$ title='Change in r(MXD,May-Aug precip,holding May-Aug temp constant)' ; if !d.name eq 'PS' then device,font_size=11 ; ol=where(oner lt -0.15,norange) if norange gt 0 then plots,statlon(ol),statlat(ol),psym=8,color=24 ; rl=where((oner ge -0.15) and (oner lt -0.05),nred) if nred gt 0 then plots,statlon(rl),statlat(rl),psym=8,color=20 ; bl=where((oner ge -0.05) and (oner lt 0.),nblue) if nblue gt 0 then plots,statlon(bl),statlat(bl),psym=8,color=21 ; gl=where((oner ge 0.) and (oner lt 0.05),ngreen) if ngreen gt 0 then plots,statlon(gl),statlat(gl),psym=8,color=22 ; bl=where((oner ge 0.05) and (oner lt 0.15),nblack) if nblack gt 0 then plots,statlon(bl),statlat(bl),psym=8,color=23 ; pl=where(oner ge 0.15,npurple) if npurple gt 0 then plots,statlon(pl),statlat(pl),psym=8,color=25 ; if !d.name eq 'PS' then device,font_size=14 ; legpos=convert_coord([-160.],[35.],/data,/to_normal) if !d.name eq 'PS' then device,/courier,/bold legend,$ ['del-r <= -0.15 (N='+string(npurple,format='(I3)')+')',$ 'del-r = -0.15 to -0.05 (N='+string(nblack,format='(I3)')+')',$ 'del-r = -0.05 to +0.00 (N='+string(ngreen,format='(I3)')+')',$ 'del-r = +0.00 to +0.05 (N='+string(nblue,format='(I3)')+')',$ 'del-r = +0.05 to +0.15 (N='+string(nred,format='(I3)')+')',$ 'del-r > +0.15 (N='+string(norange,format='(I3)')+')'],$ psym=[8,8,8,8,8,8],color=[25,23,22,21,20,24],position=legpos if !d.name eq 'PS' then device,/helvetica,/bold ; end