loop_scan_3Tmodel

0001 function [norm,curr,Tph,Tph_star] = loop_scan_3Tmodel(iEcut,frac,Ecut,Tpar,Tperp,Tback,Znorm,Zcurr,ZP0,dke_out,radialDKE,equilDKE,momentumDKE,gridDKE,Zmomcoef,Zbouncecoef,hxr,hxrparam,fieldside,chordview,hxr,hxrparam,fieldside,chordview)
0002 %
0003 %
0004 %
0005 norm = zeros(length(Tpar),length(Tperp));
0006 curr = zeros(length(Tpar),length(Tperp));
0007 Tph = zeros(length(Tpar),length(Tperp));
0008 Tph_star = zeros(length(Tpar),length(Tperp));
0009 %
0010 dt = 1;%integration time step (s)
0011 %
0012 f3tempmodel(equilDKE,gridDKE,frac,Tpar,Tperp,Tback,Ecut);
0013 
0014 [Zbremchord,equilHXR] = bremchord_dke_yp(dkeparam,dkedisplay,equil,radialDKE,Zcurr,hxr,hxrparam);
0015 [Zbremplasma] = bremsstrahlung_dke_yp(dkeparam,dkedisplay,equilHXR,radialDKE,mksa,momentumDKE,dke_out,Zbremchord,hxr,hxrparam,fieldside,chordview);%No magnetic field line effect (perp view)
0016 [Zbremdiag] = bremdiag_dke_yp(Zbremplasma,hxr,hxrparam,dt);
0017 %
0018 norm(ivmin,ivmax) = Znorm.x_0(1);
0019 curr(ivmin,ivmax) = Zcurr.x_0(1);
0020 Tph(ivmin,ivmax) = Zbremdiag.tphot_plasma;
0021 Tph_star(ivmin,ivmax) = Zbremdiag.tphot_plasma_exp;
0022 end
0023 
0024