From a9227562be0bd48ff790b8f7ec2459b4af3dd99b Mon Sep 17 00:00:00 2001 From: Olivier Sauter <Olivier.Sauter@epfl.ch> Date: Thu, 17 Oct 2019 13:32:17 +0200 Subject: [PATCH] add error_bar tests in commecnted part TCV_IMAS/tcv_get_ids_equilibrium.m --- matlab/TCV_IMAS/tcv_get_ids_equilibrium.m | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) diff --git a/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m b/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m index 54b00bb6..54dbc48c 100644 --- a/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m +++ b/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m @@ -385,15 +385,22 @@ end % $$$ ids_equilibrium.time_slice{it}.coordinate_system.grid.dim2 = linspace(0,2*pi,ldim2); % $$$ ids_equilibrium.time_slice{it}.coordinate_system.tensor_contravariant = 2.*ones(ldim1,ldim2,3,3); % $$$ ids_equilibrium.time_slice{it}.coordinate_system.tensor_covariant = 0.5*ones(ldim1,ldim2,3,3); -% $$$ ids_equilibrium.time_slice{it}.coordinate_system.g11_contravariant = 11.*ones(ldim1,ldim2,3,3); % $$$ ids_equilibrium.time_slice{it}.coordinate_system.g13_contravariant = 13.*ones(ldim1,ldim2,3,3); -% $$$ for it=3:10 +% $$$ ids_equilibrium.time_slice{it}.coordinate_system.g13_contravariant_error_upper = 14.*ones(ldim1,ldim2,3,3); +% $$$ ids_equilibrium.time_slice{it}.coordinate_system.g13_contravariant_error_lower = 12.*ones(ldim1,ldim2,3,3); +% $$$ for it=1:2100 +% $$$ ids_equilibrium.time_slice{it}.coordinate_system.g11_contravariant = 11.*ones(ldim1,ldim2,3,3); +% $$$ ids_equilibrium.time_slice{it}.coordinate_system.tensor_covariant = 0.5*ones(ldim1,ldim2,3,3); % $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid_type.name = profiles_2d.grid_type.name; % $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid_type.index = 11; % $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid_type.description = profiles_2d.grid_type.description; % $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid.dim1 = linspace(0,1,ldim1)'; +% $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid.dim1_error_upper = 1.2.*linspace(0,1,ldim1)'; +% $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid.dim1_error_lower = 0.8.*linspace(0,1,ldim1)'; % $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.grid.dim2 = linspace(0,2*pi,ldim2); % $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.psi(:,:) = 11.*ones(ldim1,ldim2); +% $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.psi_error_upper(:,:) = 12.*ones(ldim1,ldim2); +% $$$ ids_equilibrium.time_slice{it}.profiles_2d{2}.psi_error_lower(:,:) = 10.*ones(ldim1,ldim2); % $$$ end % cocos automatic transform -- GitLab