Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
G
gdat
Manage
Activity
Members
Labels
Plan
Issues
Issue boards
Milestones
Wiki
Code
Merge requests
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Model registry
Operate
Environments
Monitor
Incidents
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
Show more breadcrumbs
SPC
gdat
Commits
0f8108c4
Commit
0f8108c4
authored
7 months ago
by
Olivier Sauter
Browse files
Options
Downloads
Patches
Plain Diff
rm spaces at end of lines...
parent
2d5c99d2
No related branches found
No related tags found
1 merge request
!137
Add quantities to ids for MRE
Pipeline
#204643
passed
7 months ago
Stage: test
Stage: post-test
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
matlab/TCV_IMAS/tcv_get_ids_core_sources.m
+7
-7
7 additions, 7 deletions
matlab/TCV_IMAS/tcv_get_ids_core_sources.m
with
7 additions
and
7 deletions
matlab/TCV_IMAS/tcv_get_ids_core_sources.m
+
7
−
7
View file @
0f8108c4
...
...
@@ -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
));
...
...
This diff is collapsed.
Click to expand it.
Preview
0%
Loading
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Save comment
Cancel
Please
register
or
sign in
to comment