From 0f8108c49403a338b575ff7f99cd184abc22a96a Mon Sep 17 00:00:00 2001 From: Olivier Sauter <olivier.sauter@epfl.ch> Date: Sun, 18 Aug 2024 16:14:07 +0200 Subject: [PATCH] rm spaces at end of lines... --- matlab/TCV_IMAS/tcv_get_ids_core_sources.m | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/matlab/TCV_IMAS/tcv_get_ids_core_sources.m b/matlab/TCV_IMAS/tcv_get_ids_core_sources.m index 1a3c531d..00eb0980 100644 --- a/matlab/TCV_IMAS/tcv_get_ids_core_sources.m +++ b/matlab/TCV_IMAS/tcv_get_ids_core_sources.m @@ -103,7 +103,7 @@ vol = interp1(liuqe_time,vol_liu.data.', ohm_tgrid)'; R0 = ids_core_sources.vacuum_toroidal_field.r0; B0 = interp1(ids_core_sources.time,ids_core_sources.vacuum_toroidal_field.b0,ohm_tgrid); -jpar_tilde_to_jpar_0 = R0*T.*Rm2_fs./B0'; +jpar_tilde_to_jpar_0 = R0*T.*Rm2_fs./B0'; % map volume and conversion factor on rho_pol grid of ohm_data cd_dens jpar_tilde_to_jpar_0_mapped = ... @@ -214,15 +214,15 @@ if ~isempty(ec_gdat.ec.data) % if EC data available, fill sources T = interp1(liuqe_time,T_liu.data.', ec_tgrid_toray)'; Rm2_fs = interp1(liuqe_time,Rm2_fs_liu.data.',ec_tgrid_toray)'; vol = interp1(liuqe_time,vol_liu.data.', ec_tgrid_toray)'; - Ip = interp1(liuqe_time,Ip_liu.data, ec_tgrid_toray); - q = interp1(liuqe_time,q_liu.data.', ec_tgrid_toray)'; + Ip = interp1(liuqe_time,Ip_liu.data, ec_tgrid_toray); + q = interp1(liuqe_time,q_liu.data.', ec_tgrid_toray)'; mVprime = -2*pi*q./T./Rm2_fs; B2_fs = mu0*Ip./mVprime + T.^2.*Rm2_fs; - + % get vacuum field data from ids R0 = ids_core_sources.vacuum_toroidal_field.r0; B0 = interp1(ids_core_sources.time,ids_core_sources.vacuum_toroidal_field.b0,ec_tgrid_toray); - + jtoray_to_jpar0 = 2*pi./B0./(T.*Rm2_fs).*B2_fs; % interpolate on ec_data rho_pol grid @@ -242,7 +242,7 @@ if ~isempty(ec_gdat.ec.data) % if EC data available, fill sources p_integrated(isnan(p_integrated)) = 0; cd_dens = ec_data.cd_dens_doublewidth.data; % use double width to have realistic deposition broadening cd_dens(isnan(cd_dens)) = 0; - + cd_integrated = zeros(size(ec_data.cd_integrated.data)); for i_lau = active_launchers % convert to j// with conversion factor @@ -286,7 +286,7 @@ if ~isempty(ec_gdat.ec.data) % if EC data available, fill sources trace_cd_dens = squeeze(norm_cd_dens(irho,i_lau,:)); trace_cd_integrated = squeeze(norm_cd_integrated(irho,i_lau,:)); % interpolate on gdat powers tgrid, interpos 21 to extrapolate, - % taking constant value y_edge; needed for time window between start of EC and first TS time + % taking constant value y_edge; needed for time window between start of EC and first TS time interp_norm_p_dens( irho,i_lau,3:end-2) = interpos(21,ec_tgrid_toray,trace_p_dens, ec_tgrid_out(3:end-2)); interp_norm_p_integrated( irho,i_lau,3:end-2) = interpos(21,ec_tgrid_toray,trace_p_integrated, ec_tgrid_out(3:end-2)); interp_norm_cd_dens( irho,i_lau,3:end-2) = interpos(21,ec_tgrid_toray,trace_cd_dens, ec_tgrid_out(3:end-2)); -- GitLab