Added the vessel annular description. Added the vessel resistivity as an...

Added the vessel annular description. Added the vessel resistivity as an average value of the poloidally space varying resistivity given in Jean-Marc LIUQE paper. Added type 0 to wall description

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@12013 d63d8f72-b253-0410-a779-e742ad2e26cf

crpptbx/TCV_IMAS/tcv_get_ids_wall.m

@@ -14,7 +14,7 @@ function [ids_wall,ids_wall_description,varargout] = tcv_get_ids_wall(shot, ids_
 ids_wall.global_quantities.neutral = {}; % wall does not work with Antoine's addpath
 ids_wall.description_ggd = {};
 ids_wall.description_2d{1}.mobile.unit = {};
-ids_wall.description_2d{1}.vessel.unit = {};
+
 
 %% TODO 
 ids_description = [];

crpptbx/TCV_IMAS/tcv_ids_wall_description_2d.m

 function [ids_struct_out,ids_struct_out_description] =  tcv_ids_wall_description_2d(shot, ids_structures)
 
-%% 
-r_l = mdsdata('STATIC("R_T")');
-z_l = mdsdata('STATIC("Z_T")');
+ids_struct_out =  ids_structures;
 
 % There might be many 2d descriptions
-ids_struct_out =  ids_structures;
+%% Type 0 "index = 0 for equilibrium codes (single closed limiter and vessel)"
+
+ids_struct_out{1}.type.index = 0;
+% limiter
+r_l = mdsdata('STATIC("R_T")');
+z_l = mdsdata('STATIC("Z_T")');
 ids_struct_out{1}.limiter.unit{1}.outline.r = r_l;
 ids_struct_out_description{1}.limiter.unit{1}.outline_r = ['from ''STATIC("R_T")'];
 ids_struct_out{1}.limiter.unit{1}.outline.z = z_l;
 ids_struct_out_description{1}.limiter.unit{1}.outline_z = ['from ''STATIC("Z_T")'];
 
+% vessel
+r_v_in = mdsdata('STATIC("R_V:IN")');
+z_v_in = mdsdata('STATIC("Z_V:IN")');
+r_v_out = mdsdata('STATIC("R_V:OUT")');
+z_v_out = mdsdata('STATIC("Z_V:OUT")');
+ids_struct_out{1}.vessel.unit{1}.annular.outline_inner.r = r_v_in ;
+ids_struct_out_description{1}.vessel.unit{1}.annular.outline_inner.r = ['from ''STATIC("R_V:IN")'];
+ids_struct_out{1}.vessel.unit{1}.annular.outline_inner.z = z_v_in;
+ids_struct_out_description{1}.vessel.unit{1}.annular.outline_inner.z = ['from ''STATIC("Z_V:IN")'];
+ids_struct_out{1}.vessel.unit{1}.annular.outline_outer.r = r_v_out;
+ids_struct_out_description{1}.vessel.unit{1}.annular.outline_outer.r = ['from ''STATIC("R_V:OUT")'];
+ids_struct_out{1}.vessel.unit{1}.annular.outline_outer.z = z_v_out;
+ids_struct_out_description{1}.vessel.unit{1}.annular.outline_outer.z = ['from ''STATIC("Z_V:OUT")'];
+
+% The value of the resistivity is taken as an average of the value in Fig.3
+% of Jean-Marc LIUQE paper.
+ids_struct_out_description{1}.vessel.unit{1}.annular.resistivity = 0.25*1e-3; %[ohm m]
+
 %% TODO Add the other description of the wall and the input pharser