visu_luke

0001 h = findobj(0,'type','figure','tag','luke');
0002 if isempty(h)
0003        h=figure('tag','luke');
0004 else
0005        figure(h);
0006 end   
0007 clf
0008 set(h,'defaultaxesfontsize',10,'defaultaxesfontweight','bold','defaultaxesfontname','times', ...
0009     'defaultlinelinewidth',2,'color',[1 1 1])
0010    
0011 subplot(2,2,1)        
0012 ,
0013       graph1D_jd(equilDKE.xrho,j_OHM_c,0,0,'\rho','J_{ohm} (MA/m^2)',[par.conf,',I_{ohm} = ',num2str(I_OHM_c*1000,3),' (kA)'],NaN,'','','-','none','b',2);hold off,drawnow
0014 %
0015 %title('Ohmic current profile')
0016 subplot(2,2,2)   
0017 
0018 
0019       graph1D_jd(equilDKE.xrho,p_OHM_c,0,0,'\rho','P_{ohm} (MW/m^3)',['P_{ohm} = ',num2str(P_OHM_c*1000,4),' (kW) / P_{synohm}= ',num2str(P_OHMsyn_c*1000,4),' kw'],NaN,'','','-','none','r',2);drawnow
0020 %title('Ohmic power')
0021 subplot(2,2,3)   
0022         nolh = 1;
0023         nofce = 0;
0024         if isstruct(cons)
0025            if isfield(cons,'lh_cons')
0026              nolh = 1;
0027            else
0028              nolh = 0;
0029            end
0030            if isfield(cons,'ecrh_cons')
0031              nofce = 1;
0032            else
0033              nofce = 0;
0034            end
0035         end
0036     if nolh == 1 & nofce == 0 
0037       graph1D_jd(equilDKE.xrho,j_RF_c,0,0,'\rho','J_{LH} (MA/m^2)',['I_{LH} = ',num2str(I_RF_c,3),' (MA) / Ip =',num2str(equi.ip/1e6,3),' (MA)'],NaN,'','','-','none','b',2);hold off,drawnow
0038         elseif nolh==0 & nofce == 1
0039       graph1D_jd(equilDKE.xrho,j_RF_c,0,0,'\rho','J_{EC} (MA/m^2)',['I_{EC} = ',num2str(I_RF_c,3),' (MA) / Ip =',num2str(equi.ip/1e6,3),' (MA)'],NaN,'','','-','none','b',2);hold off,drawnow
0040         else
0041       graph1D_jd(equilDKE.xrho,j_RF_c,0,0,'\rho','J_{tot} (MA/m^2)',['I_{tot} = ',num2str(I_RF_c,3),' (MA) / Ip =',num2str(equi.ip/1e6,3),' (MA)'],NaN,'','','-','none','b',2);hold off,drawnow
0042         end
0043 %title('LH Current profile (consistent with E-field)')
0044 
0045 subplot(2,2,4)
0046         if nolh == 1 & nofce == 0 
0047           if strcmp(par.lob,'No')
0048             if isstruct(cons)
0049               Pinj = sum(abs(cons.lh_cons))/1e3*0.7;
0050             else
0051               Pinj = sum(abs(cons))/1e3*0.7;
0052             end
0053           else
0054             if isstruct(cons)
0055               Pinj = sum(abs(cons.lh_cons))/1e3;
0056             else
0057               Pinj = sum(abs(cons))/1e3;
0058             end
0059           end
0060       graph1D_jd(equilDKE.xrho,p_RF_c,0,0,'\rho','P_{LH} (MW/m^3)',['P_{LH} = ',num2str(P_RF_c*1000,4),' (kW) / Pinj=',num2str(Pinj,4),' kW'],NaN,'','','-','none','r',2);drawnow
0061         elseif nolh==0 & nofce == 1
0062           Pinj = sum(abs(cons))/1e3;
0063       graph1D_jd(equilDKE.xrho,p_RF_c,0,0,'\rho','P_{EC} (MW/m^3)',['P_{EC} = ',num2str(P_RF_c*1000,4),' (kW) / Pinj=',num2str(Pinj,4),' kW'],NaN,'','','-','none','r',2);drawnow
0064         else
0065           if strcmp(par.lob,'No')
0066             Plh  = sum(abs(cons.lh_cons))*0.7;
0067           else
0068             Plh  = sum(abs(cons.lh_cons));
0069           end
0070           Pec  = sum(abs(cons.ecrh_cons));
0071           Pinj = Pec+Plh;
0072        graph1D_jd(equilDKE.xrho,p_RF_c,0,0,'\rho','P_{tot} (MW/m^3)',['P_{tot} = ',num2str(P_RF_c*1000,4),' (kW) / Pinj=',num2str(Pinj,4),' kW'],NaN,'','','-','none','r',2);drawnow
0073        end
0074 
0075 %
0076 %title('LH power')