change sign of the_pol waiting for ids2cpo to change transformation

matlab/TCV_IMAS/tcv_get_ids_ec_antennas.m

@@ -34,6 +34,8 @@ end
 % $$$ ids_ec_antennas_description.time = 'time from power';
 [r_phi0,z_phi0,r1,z1,rc,zc,thetphi_L,thetphi_tor,raysbylauncher,powers,fname,z_axis,launch_params] = ...
     toray_raygeom_TCV(shot,pow.ec.t(1:10:end));
+ids_ec_antennas_description.launch_params = launch_params;
+
 
 for iant=1:nb_antennas
   ids_ec_antennas.antenna{iant}.identifier = ['L' num2str(iant)];
@@ -74,17 +76,23 @@ for iant=1:nb_antennas
       ids_ec_antennas.antenna{iant}.launching_position.r.time(i) = time_launch;
       ids_ec_antennas.antenna{iant}.launching_position.z.time(i) = time_launch;
       ids_ec_antennas.antenna{iant}.launching_position.phi.time(i) = time_launch;
-      ids_ec_antennas.antenna{iant}.launching_angle_pol.data(i) = (90.-launch_params{iant}{it_ok{iant}(i)}.theta_toray) * pi/180.;
+      ids_ec_antennas.antenna{iant}.launching_angle_pol.data(i) = - (90.-launch_params{iant}{it_ok{iant}(i)}.theta_toray) * pi/180.;
       ids_ec_antennas.antenna{iant}.launching_angle_pol.time(i) = time_launch;
       ids_ec_antennas.antenna{iant}.launching_angle_tor.data(i) = (180.-launch_params{iant}{it_ok{iant}(i)}.phi_toray) * pi/180.;
       ids_ec_antennas.antenna{iant}.launching_angle_tor.time(i) = time_launch;
-      ids_ec_antennas.antenna{iant}.beam.spot.size.data(i,1:2) = 0.05.*ones(size(time_launch));
+% $$$       ids_ec_antennas.antenna{iant}.beam.spot.size.data(i,1) = 0.023;
+% $$$       ids_ec_antennas.antenna{iant}.beam.spot.size.data(i,2) = 0.012;
+      ids_ec_antennas.antenna{iant}.beam.spot.size.data(1,i) = 0.023;
+      ids_ec_antennas.antenna{iant}.beam.spot.size.data(2,i) = 0.012;
       ids_ec_antennas.antenna{iant}.beam.spot.size.time(i) = time_launch;
-      ids_ec_antennas.antenna{iant}.beam.spot.angle.data(i,1:2) = 0.0.*ones(size(time_launch));
+      ids_ec_antennas.antenna{iant}.beam.spot.angle.data(i) = 0.0;
       ids_ec_antennas.antenna{iant}.beam.spot.angle.time(i) = time_launch;
-      ids_ec_antennas.antenna{iant}.beam.phase.angle.data(i,1:2) = 0.0.*ones(size(time_launch));
+% $$$       ids_ec_antennas.antenna{iant}.beam.phase.curvature.data(i,1) = 1./1.88;
+% $$$       ids_ec_antennas.antenna{iant}.beam.phase.curvature.data(i,2) = 1./-0.95; % negative in launch_params...
+      ids_ec_antennas.antenna{iant}.beam.phase.curvature.data(1,i) = 1./1.88;
+      ids_ec_antennas.antenna{iant}.beam.phase.curvature.data(2,i) = 1./-0.95; % negative in launch_params...
+      ids_ec_antennas.antenna{iant}.beam.phase.angle.data(i) = 0.0;
       ids_ec_antennas.antenna{iant}.beam.phase.angle.time(i) = time_launch;
-      ids_ec_antennas.antenna{iant}.beam.phase.curvature.data(i,1:2) = 1.0.*ones(size(time_launch));
       ids_ec_antennas.antenna{iant}.beam.phase.curvature.time(i) = time_launch;
     end
   end