From ac34c31796e02bc778294f0d2c54c88d6b7aedd6 Mon Sep 17 00:00:00 2001
From: Antonia Frank <antonia.frank@epfl.ch>
Date: Thu, 8 Aug 2024 18:04:18 +0200
Subject: [PATCH] Correct units for current density in ec_data call of gdat
---
matlab/TCV/gdat_tcv.m | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/matlab/TCV/gdat_tcv.m b/matlab/TCV/gdat_tcv.m
index 7bf32e2f..adfcc621 100644
--- a/matlab/TCV/gdat_tcv.m
+++ b/matlab/TCV/gdat_tcv.m
@@ -1704,7 +1704,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
ec_data.cd_dens.data = currentdrive_dens.data * 1e6;
ec_data.cd_dens.data(:,end+1,:) = sum(ec_data.cd_dens.data,2,'omitnan'); % add total
ec_data.cd_dens.label = [strrep(currentdrive_dens.comment,'MA','A') ' ; last index is total'];
- ec_data.cd_dens.units = 'A/m^2';
+ ec_data.cd_dens.units = 'A/m^3';
ec_data.cd_dens.x = currentdrive_dens.rgrid;
ec_data.cd_dens.grids = currentdrive_dens.grids;
ec_data.cd_dens.t = currentdrive_dens.tgrid;
@@ -1723,7 +1723,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
%
ec_data.max_cd_dens.data = icdmax.data * 1e6;
ec_data.max_cd_dens.label = strrep(icdmax.comment,'MA','A');
- ec_data.max_cd_dens.units = 'A/m^2';
+ ec_data.max_cd_dens.units = 'A/m^3';
ec_data.max_cd_dens.x = [];
ec_data.max_cd_dens.t = icdmax.tgrid;
ec_data.max_cd_dens.dim = {ec_data.max_cd_dens.t};
--
GitLab