; IDL procedure to make plots of Cyclone chi vs. R/LT ; Modified to show improved Gyrofluid closures and new GK results. ; pro bolton ; !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, 0.0, 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, 5.3, 6.9, 8.97] chi_gf = [ 5.0, 6.0, 8.0, 11.0] RLT_gf98 = [4.83, 5.3, 6.9, 9.0] chi_gf98 = [ 0.0, 1.5, 5.0, 7.5] ; select the first line for postscript, next line for X-windows: plot_setup_ps, 'bolton.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=2.0 !p.thick=4.0 !x.thick=4.0 !y.thick=4.0 symsize=2.0 plot, RLT_gk, chi_gk, /nodata oplot, RLT_gk, chi_gk, psym=2, symsize=symsize, color=green oplot, RLT_gf, chi_gf, psym=1, symsize=symsize oplot, RLT_gf98, chi_gf98, psym=4, symsize=symsize, color=red ; oplot, RLT_gkfull, chi_gkfull, psym=4, symsize=symsize 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) chi_ifsgk98b = 0.8*3.735*( (rlt_ifs96-6.0) < (rlt_ifs96-6.0 > 0)^0.5) chi_ifsgk98c = 15*(rlt_ifs96-6.0 >0)/rlt_ifs96 chi_bolton1 = 2.4*( (rlt_ifs96-6.0 > 0)^0.5) chi_bolton2 = 0.5*3.735*( (rlt_ifs96-6.0) < (rlt_ifs96-6.0 > 0)^0.5) oplot, RLT_ifs96, chi_ifs96 oplot, RLT_ifs96, chi_ifsgk98c, color=green oplot, RLT_ifs96, chi_ifsgk98b, color=blue oplot, RLT_ifs96, chi_bolton1, color=black oplot, RLT_ifs96, chi_bolton2, color=black xyouts, 2.8, 9.2, "IFS-PPPL94" ;xyouts, 15, 8., "GK fit", color=green xyouts, 5.3, 10.6, "GF94" xyouts, 7.5, 6.5, "GF98", color=red ; xyouts, 2.5, 9.8, "flux-tube" xyouts, 11.5, 8.5, "GK", color=green ;xyouts, 11., 5.7, "GyroKinetic", color=green ;xyouts, 11., 4.7, "flux-tube", color=green ; xyouts, 10., 2, "GyroKinetic full-torus" xyouts, 16., 8.3, "B1" xyouts, 18., 5.5, "B2" ; 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