; ; Makes a movie of yearly maps of calibrated (PCR-infilled or not) MXD ; reconstructions ; of growing season temperatures. Uses "corrected" MXD - but shouldn't usually ; plot past 1960 because these will be artificially adjusted to look closer to ; the real temperatures. ; doinfill=1 dofilt=1 thalf=5. for imov = 0 , 5 do begin dost=imov*100+1400 if imov ne 5 then ndo=100 else ndo=51 moviefn='recon'+string([dost,dost+ndo-1],format='(2I4)') if dofilt ne 0 then moviefn=moviefn+'sm'+string(thalf,format='(I2.2)') print,moviefn ; ; Now prepare for plotting ; loadct,39 multi_plot,nrow=1,layout='centred' if !d.name eq 'X' then begin window,ysize=300,xsize=300 !p.font=-1 !p.charsize=!p.charsize*1.5 endif else begin !p.font=0 device,/helvetica,/bold,font_size=11 endelse def_1color,20,color='palepurple' def_1color,21,color='lpurple' def_1color,22,color='deepblue' def_1color,23,color='mlblue' def_1color,24,color='vlblue' def_1color,25,color='vvlgreen' def_1color,26,color='lsand' def_1color,27,color='orange' def_1color,28,color='red' def_1color,29,color='dred' ; ; Define map parameters ; map=def_map(/npolar) & map.limit(0)=25 & map.xmargin=[0,0] & map.ymargin=[0,0] labels=def_labels(/off) coast=def_coast(/get_device) & coast.thick=[0.5,2.5] levs=[-100,-2,-1.2,-0.8,-0.4,-0.2,0,0.3,0.6,1,100] cols=indgen(10)+20 ; ; Get the calibrated data ; if doinfill eq 0 then begin restore,'calibmxd5.idlsave' ; Gets: g,mxdyear,mxdnyr,fdcalibu,fdcalibc,mxdfd2,timey,fdseas endif else begin restore,'calibmxd5_pcr.idlsave' ; Gets: g,mxdyear,mxdnyr,fdcalibc,timey,fdseas endelse ; ; Filter the data in time if required! ; if dofilt ne 0 then begin print,'Time-filtering of the data; for ix = 0 , g.nx-1 do begin print,g.x(ix),format='($,F8.1)' for iy = 0 , g.ny-1 do begin onets=reform(fdcalibc(ix,iy,*)) filter_cru,thalf,/nan,tsin=onets,tslow=lowts fdcalibc(ix,iy,*)=lowts(*) endfor endfor print endif ; for iyr = dost, dost+ndo-1 do begin iii=where(mxdyear eq iyr) iii=iii(0) fd=reform(fdcalibc(*,*,iii)) ; inter_boxfd,fd,g.x,g.y,$ coast=coast,map=map,labels=labels,$ levels=levs,c_colors=cols xyouts,-47,28,string(iyr,format='(I4)') lon_polar,map=map,[-135.,-45.,45.,135.] ; tvlct,/get,colr,colg,colb img1=tvrd() wr_movie,moviefn,img1,iyr-dost,ndo,/gif,/noscale,/pseudo,$ colr=colr,colg=colg,colb=colb,framedelay=25 ; endfor ; endfor ; ;!p.multi(0)=1 ;!p.position=[0.3,0.015,0.7,ykey-0.015] ;scale_horiz,levels=levs,c_colors=cols,noextremes=['Below','Above'],$ ; title='!Uo!NC wrt 1961-90' ; end