pro make_rad_plot,data ; ; Author: Mike Hapgood ; ; Revision ; 23 Feb 2000 original version ; ; Specify name of CDF file to read and of the variables to be read day_string='20001001' base_name='C1_JP_PMP_' + day_string + '_V01.cdf' cdf_name='/home/mah/data/jsoc_cats/jpsp/C1/' + base_name vnames='Epoch%C1_JP_PMP,INVAR_LAT%C1_JP_PMP,MAG_LOCAL_TIME%C1_JP_PMP' + $ 'L_VALUE%C1_JP_PMP,PRED_B_MAG%C1_JP_PMP' ; ; and then read them into data structure 'data' print,'Reading ',cdf_name data=read_mycdf(vnames,cdf_name,/all) ; close the channels left open by CDAWlib close,/all ; ; Check for valid data, i.e. exclude null/fill values ind=where(data.(4).dat gt 1.0, count) ; all if we have at least two valid data points, make a plot if count gt 1 then begin l_value=data.(4).dat(ind) b_value=data.(5).dat(ind) ; set black on white plot s1=!p.background s2=!p.color !p.background=7 !p.color=0 ptitle='Cluster - ' + base_name plot_oo,l_value,b_value,charsize=1.5, $ xtitle='L value',ytitle='B (nT)', $ title=ptitle ; restore settings !p.background=s1 !p.color=s2 endif ; return end