cfile='climate/hadat2_monthly_zonal_and_global_lr.nc' o=ncdf_open(cfile,/nowrite) ncdf_varget,o,'lat',lat ncdf_varget,o,'time',t ncdf_varget,o,'pressure',p ncdf_varget,o,'zonal',z ncdf_close,o nlat=36 nt=569 np=9 t=1958+findgen(nt)/12 str=strarr(576) openr,1,'solar/Moscow_1958-2005.dat' readf,1,str close,1 n0=where(strmid(str,0,7) eq '1958.01') n1=where(strmid(str,0,7) eq '2005.05') str=str(n0:n1) gcr=intarr(n1-n0+1) reads,str,gcr,format='(7x,f)' gcr=-gcr z_corr=fltarr(nlat,np) for i=0,nlat-1 do for j=0,np-1 do begin zz=z(i,j,*) ok=where(zz gt -1.e10 and gcr lt 10000,nok) if(nok ge 50) then z_corr(i,j)=correlate(zz(ok),gcr(ok)) else z_corr(i,j)=-10 endfor set_plot,'ps' device,/color,bits_per_pixel=8,xsize=16,ysize=22,filename='GCR_HAD_zonal.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) plot,indgen(nlat)*5,indgen(np),xrange=[-90,90],yrange=[0,9],xstyle=1,ystyle=1,title='Radiosonde Zonal Temp - GCR',pos=[.02,.5,.98,.9],/nodata,yticks=9,ytickname=['850','700','500','300','200','150','100','50','30',' '] for i=0,nlat-1 do for j=0,np-1 do if(z_corr(i,j) ge -1) then begin lats=[lat(i)-2.49,lat(i)+2.49] ps=[j,j+.99] val(*,*)=z_corr(i,j)>(-.4) contour,val,lats,ps,/fill,/overplot,levels=-.4+findgen(40)/50,max_value=1,min_value=-1 endif 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.38,.8,.42], $ xrange=[-.4,.4],xstyle=1,/noerase,ystyle=1,yrange=[0,1],yticks=1, $ ytickname=[' ',' '],title='Corr.: GCR - HAD zonal temp. anom.',charsize=1 device,/close_file end