pro test_fucking_mrr_read mrr_file='0319_short.ave' mrr_path='/Users/u0029340/Documents/data/Capricorn/mrr/' ;;;;;;;;;;; get the MRR dbz read_mrr_ave_txt_sub, mrr_file, mrr_path, mrr_dbz,mrr_jday,mrr_height ;mrr_height=rotate(reform(mrr_height),5) ;mrr_dbz=rotate(mrr_dbz, 7) caldat,mrr_jday,mm,dd,yy,hh,mi,ss ; Calculate fractional hour mrr_hrfrac=hh+mi/60.0+ss/(60.0*60.0) ;;;;;;;;;;;;;;;;;;;;;;; common colors,r_orig,g_orig,b_orig,r_curr,g_curr,b_curr loadct, 11 ; This is the stern color scheme r_curr[250:255]=255 & r_orig=r_curr g_curr[250:255]=255 & g_orig=g_curr b_curr[250:255]=255 & b_orig=b_curr gamma_ct, 1.0 set_plot, 'Z' device, set_resolution=[3000,2000] !p.multi = [0,1,1] !p.charsize=3. !p.font=-1 !p.background=!d.n_colors plot, mrr_hrfrac, mrr_dbz[*,3], color=255,xtitle='Time (UTC)',ytitle='gate value',$ thick=1, charsize=3.0, charthick=3, xthick=5, ythick=5, psym=3, $ xstyle=1 write_gif, mrr_path+'test_shit2.gif', tvrd() stop end