cfile='climate/channel_TLS_Tb_anom_v02_1.nc' o=ncdf_open(cfile,/nowrite) ncdf_varget,o,'longitude_edges',lon ncdf_varget,o,'latitude_edges',lat ncdf_varget,o,'years',y ncdf_varget,o,'months',m ncdf_varget,o,'monthly_maps',c ncdf_close,o nlon=144 nlat=72 nt=12*28 lon=(lon(1:nlon)+lon(0:nlon-1))/2 lat=(lat(1:nlat)+lat(0:nlat-1))/2 t=fltarr(nt) for i=0,27 do for j=0,11 do t(i*12+j)=y(i)+(m(j)-1)/12. c=reform(c,nlon,nlat,nt) c=c(*,*,10:nt-3) t=t(10:nt-3) nt=nt-12 c=c(*,3:nlat-4,*) lat=lat(3:nlat-4) nlat=nlat-6 zonal=fltarr(nlat,nt) for i=0,nlat-1 do for j=0,nt-1 do begin cc=c(*,i,j) ok=where(finite(cc),nok) zonal(i,j)=total(cc(ok))/nok endfor w=cos(lat*!dtor) zw=zonal for i=0,nlat-1 do zw(i,*)=zw(i,*)*w(i) global=total(zonal,1)/total(w) str=strarr(576) openr,1,'solar/Moscow_1958-2005.dat' readf,1,str close,1 n0=where(strmid(str,0,7) eq '1978.11') n1=where(strmid(str,0,7) eq '2005.10') str=str(n0:n1) gcr=intarr(n1-n0+1) reads,str,gcr,format='(7x,f)' gcr=-gcr c_corr=fltarr(nlon,nlat) z_corr=fltarr(nlat) for i=0,nlon-1 do for j=0,nlat-1 do begin cc=c(i,j,*) ok=where(finite(cc),nok) if(nok ge 50) then c_corr(i,j)=correlate(cc(ok),gcr(ok)) else c_corr(i,j)=-10 endfor for i=0,nlat-1 do begin zz=zonal(i,*) ok=where(finite(zz),nok) if(nok ge 50) then z_corr(i)=correlate(zz(ok),gcr(ok)) else z_corr(i)=-10 endfor lag=indgen(25)-6 g_corr=c_correlate(gcr,global,lag) set_plot,'ps' device,/color,bits_per_pixel=8,xsize=16,ysize=22,filename='GCR_MSU_strat.ps',xoffset=2,yoffset=5 !p.charsize=1.3 !p.charthick=4.0 !x.thick=4.0 !y.thick=4.0 !p.thick=4.0 !x.margin=[8,10] loadct,40 tvlct,red,green,blue,/get red(64:127)=4*indgen(64) blue(100:127)=255 blue(128:191)=255-4*indgen(64) green(192:255)=255-4*indgen(64) red(128:255)=255 red(128:151)=255-3*indgen(24) red(152:175)=192+3*indgen(24) red(0)=0 green(0)=0 blue(0)=0 ;red(255)=255 ;green(255)=255 ;blue(255)=255 tvlct,red,green,blue val=fltarr(2,2) map_set,/cyl,limit=[-89,-180,89,180],title=' MSU Strat Temp - GCR',pos=[.02,.5,.98,.9] for i=0,nlon-1 do for j=0,nlat-1 do if(c_corr(i,j) ge -1) then begin lats=[lat(j)-1.24,lat(j)+1.24] lons=[lon(i)-1.24,lon(i)+1.24] val(*,*)=c_corr(i,j)>(-.4) contour,val,lons,lats,/fill,/overplot,levels=-.4+findgen(40)/50,max_value=1,min_value=-1 endif map_set,/continent,/cyl,limit=[-89,-180,89,180],/noerase,pos=[.02,.5,.98,.9] z=fltarr(39,2) z(*,0)=-.38+findgen(39)/39*.76 z(*,1)=z(*,0) x=-.4+findgen(39)/38*.8 y=[0,1] contour,z,x,y,levels=-.4+findgen(40)/50,/fill,pos=[0.13,0.43,.8,.47], $ xrange=[-.4,.4],xstyle=1,/noerase,ystyle=1,yrange=[0,1],yticks=1, $ ytickname=[' ',' '],title='Corr.: GCR - MSU stratos. temp. anom.',charsize=1 plot,lat,z_corr,xrange=[-90,90],pos=[.02,.22,.98,.37],/noerase,xstyle=1,title='zonal correlation' plot,lag,g_corr,pos=[.02,.02,.98,.15],/noerase,xrange=[-6,18],xstyle=1,title='Lagged global correlation',xtitle='Lag (months)' device,/close_file end