; IDL procedure to make plots of Q_fusion vs. T_pedestal. ; using runs fit to to Dimit's IAEA 98 GK runs. ; pro ped_plot98c ; !p.multi=[0,1,2] ; 1*2 plots on each page Tped = [0.7, 1.0, 1.5, 2.0, 3.0, 4.0, 5.0, 7.0] Qfus_gf96 = [127., 145., 193., 270., 555., 950., 1377., 2179.]/100.0 plot_setup_ps, 'ped_plot98c.ps' ;plot_setup_x8 @color_names.h !y.title="Q (Paux = 100 MW)" !x.title=TextoIDL("Pedestal temperature (keV)") !x.range=[0,8] !y.range=[0,20] !p.charsize=1.5 plot, Tped, Qfus_gf96 ;xyouts, 3.0, 3.5, TextoIDL("IFS-PPPL, n_e=1.5 n_{GW}") xyouts, 3.0, 4.5, TextoIDL("IFS-PPPL") Qfus_gk98c = [439, 500, 658, 867, 1351, 1830, 2256, 2878]/100.0 ; Qfus_gk98c_low_n = [386, 488, 628, 968, 1320, 1641, 2133, 2878]/100.0 Qfus_gk98c2 = [1151., 1191., 1292., 1430., 1768., 2124., 2454., 2950.]/100.0 Qfus_gk98c2_lown = [636., 656., 707., 780., 967., 1174., 1371., 1684.]/100.0 oplot, Tped, Qfus_gk98c, color=red ;xyouts, 3.0, 12.0, TextoIDL("GK fit, n_e=1.5 n_{GW}"), color=red xyouts, 3.0, 12.5, TextoIDL("GK fit"), color=red oplot, Tped, Qfus_gk98c2, color=green ;xyouts, 2.5, 14.0, TextoIDL("GK fit/2, n_e=1.5 n_{GW}"), color=green xyouts, 1.0, 15.0, TextoIDL("GK fit/2"), color=green oplot, Tped, Qfus_gk98c2_lown, color=green ;xyouts, 2.5, 7.5, TextoIDL("GK fit/2, n_e=1.3 n_{GW}"), color=green ;xyouts, 2.5, 7.5, TextoIDL("GK fit/2, n_e=1.15 n_{GW}"), color=green xyouts, 5.5, 13.5, TextoIDL("GK fit/2, low n_e"), color=green ; 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