plot_statefile

fsize=20
figure()
st=OMFIT['ASIPP_IM']['HISTORY']['evo_iterations']['it_613']['ONETWO_ohmic']['FILES']['statefile']
rho=linspace(0,1,201)

Jrf=st['currf']['data']
J_t=st['curden']['data']
J_ohm=st['curohm']['data']
J_beam=st['curbeam']['data']
J_boot=st['curboot']['data']

totboot_cur=str(st['totboot_cur']['data']*1e-3)
totbeam_cur=str(st['totbeam_cur']['data']*1e-3)
totrf_cur=str(st['totrf_cur']['data']*1e-3)
totohm_cur=str(st['totohm_cur']['data']*1e-3)
totcur_bc=str(st['totcur_bc']['data']*1e-3)

qerf=st['qrfe']['data']
qbeame=st['qbeame']['data']
qbeami=st['qbeami']['data']
subplot(1,2,1)
plot(rho,J_t,'k-',lw=4,label='J_total,'+totcur_bc[:5]+'kA')
plot(rho,Jrf,'r-',lw=2,label='J_rf,'+totrf_cur[:5]+'kA')
plot(rho,J_boot,'k--',lw=2,label='J_boot,'+totboot_cur[:5]+'kA')
plot(rho,J_beam,'b-',lw=2,label='J_beam,'+totbeam_cur[:5]+'kA')
plot(rho,J_ohm,'g-',lw=2,label='J_ohm,'+totohm_cur[:5]+'kA')
xlabel(r'$\rho$',fontsize=fsize)
ylabel(r'J(A/m^2)',fontsize=fsize)
tick_params(labelsize=18,direction='in')
legend(loc='best',fontsize=20,frameon=0).draggable=True

subplot(1,2,2)
plot(rho,qerf,'r-',lw=4,label='qe_rf')
plot(rho,qbeame,'k-',lw=2,label='qbeame')
plot(rho,qbeami,'k--',lw=2,label='qbeami')
xlabel(r'$\rho$',fontsize=fsize)
ylabel(r'q(W/m^3)',fontsize=fsize)
tick_params(labelsize=18,direction='in')
legend(loc='best',fontsize=20,frameon=0).draggable=True