; IDL procedure to make plots of Cyclone chi vs. R/LT ; pro cyclone ; !p.multi=[0,1,2] ; 1*2 plots on each page RLT_gk = [5.3, 6, 6.9, 8.97, 12, 14, 16, 18] chi_gk = [0.0, 1.32, 2.4, 5, 8, 9, 9.5, 10] RLT_gkfull = [5.3, 6.9, 8.97] chi_gkfull = [0.4, 1.0, 2.2] RLT_gf = [4.83, 6.9, 8.97] chi_gf = [5.0, 8.0, 11.0] ; select the first line for postscript, next line for X-windows: plot_setup_ps, 'cyclone.ps' ;plot_setup_x8 @color_names.h !y.title=TextoIDL("\chi_i (\rho^2v_t /L_n)") !x.title=TextoIDL("R/L_{Ti}") !x.range=[0,20] !y.range=[0,12] !p.charsize=1.5 plot, RLT_gk, chi_gk, /nodata oplot, RLT_gk, chi_gk, psym=2, color=green oplot, RLT_gf, chi_gf, psym=1 oplot, RLT_gkfull, chi_gkfull, psym=4 xyouts, 13, 10.5, "GF fit" xyouts, 15, 8., "GK fit", color=green xyouts, 1., 10.8, "GyroFluid flux-tube" xyouts, 10., 4.7, "GyroKinetic flux-tube", color=green xyouts, 10., 2, "GyroKinetic full-torus" RLT_ifs96 = 4.0 + findgen(100)/5.0 chi_ifs96 = 3.735*( (rlt_ifs96-4.0) < (rlt_ifs96-4.0)^0.5) chi_ifsgk98 = 1.2*( (rlt_ifs96-5.0) < 2.3*( (rlt_ifs96-5.0 > 0) )^0.5) oplot, RLT_ifs96, chi_ifs96 oplot, RLT_ifs96, chi_ifsgk98, color=green ; xyouts, 3.5, 5, TextoIDL("IFS-PPPL(96, w/ trap. elec.)") ; tmp1='' ; read,tmp1,Prompt='Return when ready for next plot' ; if tmp1 eq 'q' then return ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; IF !D.name EQ 'PS' THEN BEGIN device,/close_file ; close postscript graphics file ENDIF ; set_plot, 'X' ; return to X-windows graphics end