From 5349abacb6459e807e9d92b38253824c5d7a2833 Mon Sep 17 00:00:00 2001
From: Antonia Frank <antonia.frank@epfl.ch>
Date: Tue, 15 Oct 2024 13:55:20 +0200
Subject: [PATCH] Add total magnetic field to have a default output of gdat
bfields
---
matlab/TCV/gdat_tcv.m | 19 +++++++++++++++++++
1 file changed, 19 insertions(+)
diff --git a/matlab/TCV/gdat_tcv.m b/matlab/TCV/gdat_tcv.m
index 4885f998..564eebba 100644
--- a/matlab/TCV/gdat_tcv.m
+++ b/matlab/TCV/gdat_tcv.m
@@ -907,6 +907,24 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
gdat_data.Btor.x = {L.rx,L.zx};
gdat_data.Btor.data_fullpath = datapath;
gdat_data.Btor.label = 'Toroidal magnetic field map in (R,Z)';
+
+ % total magnetic field
+ gdat_data.Btot.data = sqrt(gdat_data.Btor.data.^2 + gdat_data.Bz.data.^2 + gdat_data.Br.data.^2);
+ gdat_data.Btot.units = 'T';
+ gdat_data.Btot.dim = {L.rx,L.zx,LY.t};
+ gdat_data.Btot.dimunits = {'m','m','s'};
+ gdat_data.Btot.t = LY.t;
+ gdat_data.Btot.x = {L.rx,L.zx};
+ gdat_data.Btot.data_fullpath = datapath;
+ gdat_data.Btot.label = 'Total magnetic field map in (R,Z)';
+
+ gdat_data.data = gdat_data.Btot.data;
+ gdat_data.units = 'T';
+ gdat_data.dim = {L.rx,L.zx,LY.t};
+ gdat_data.dimunits = {'m','m','s'};
+ gdat_data.t = gdat_data.Btot.t;
+ gdat_data.x = {L.rx,L.zx};
+ gdat_data.label = 'Total magnetic field map in (R,Z)';
else
@@ -914,6 +932,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
gdat_data.Btor = [];
gdat_data.Br = [];
gdat_data.Bz = [];
+ gdat_data.Btot = [];
end
--
GitLab