From bf10d66036c39c3791f13e93b185fdaca52b0bb4 Mon Sep 17 00:00:00 2001
From: Antonia Frank <antonia.frank@epfl.ch>
Date: Wed, 14 Aug 2024 10:57:45 +0200
Subject: [PATCH] start changing to nodes for LIUQE calls, in progress
---
matlab/TCV_IMAS/tcv_get_ids_core_sources.m | 50 +++++++++++++++++++---
1 file changed, 45 insertions(+), 5 deletions(-)
diff --git a/matlab/TCV_IMAS/tcv_get_ids_core_sources.m b/matlab/TCV_IMAS/tcv_get_ids_core_sources.m
index 96e8def3..a025571d 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_core_sources.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_core_sources.m
@@ -48,6 +48,29 @@ params_eff = params_eff_ref; params_eff.data_request='ohm_data';
ohm_gdat = gdat(params_core_sources.shot,params_eff);
params_eff = params_eff_ref; params_eff.data_request='powers';
powers_gdat = gdat(params_core_sources.shot,params_eff);
+
+%% load liuqe data from nodes for conversions
+params_eff = params_eff_ref;
+% normalized sqrt poloidal flux including axis
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:rho';
+rho_pol_norm_liu = gdat(shot,params_eff);
+% <1/R^2>
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:gpsi0_r2_fsa';
+Rm2_fs = gdat(shot,params_eff);
+% Volume
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:vol';
+vol= gdat(shot,params_eff);
+% T(\psi) = R*B_tor
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:rbtor_rho';
+T = gdat(shot,params_eff);
+% Q1Q: -dpsi/dV
+%mVprime = 1./LY.Q1Q; % Q1Q: -dpsi/dV
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:DPSI_VOL';
+Vprime = gdat(shot,params_eff);
+% I_tor = I_p, total plasma current
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:i_pl';
+Ip = gdat(shot,params_eff);
+
%% initialize source from template
% ohm
ohm_t_grid = powers_gdat.ohm.t; ohm_n_t = numel(ohm_t_grid);
@@ -64,13 +87,26 @@ ohm_data = ohm_gdat.ohm.ohm_data;
% load LIUQE data to convert
% \tilde{j}_// = <jdotB>/(R0<Bphi/R>) to j_//0 = <jdotB>/B0; Bphi=T(psi)/R
[L,~,LY] = liuqe(shot,ohm_t_grid,'iterq',50,'ilim',3,'icsint',true);
-rho_pol_norm_liu = L.pQ; % normalized sqrt poloidal flux including axis
-T = LY.TQ;
-Rm2_fs = LY.Q2Q;
-vol = LY.VQ;
+params_eff = params_eff_ref;
+% normalized sqrt poloidal flux including axis
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:rho';
+rho_pol_norm_liu = gdat(shot,params_eff);
+% <1/R^2>
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:gpsi0_r2_fsa';
+Rm2_fs = gdat(shot,params_eff);
+% Volume
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:vol';
+vol= gdat(shot,params_eff);
+% T(\psi) = R*B_tor
+params_eff.data_request='\tcv_shot::top.results.equil_1.results:rbtor_rho';
+T = gdat(shot,params_eff);
+%rho_pol_norm_liu = L.pQ; % normalized sqrt poloidal flux including axis
+%T = LY.TQ;
+%Rm2_fs = LY.Q2Q;
+%vol = LY.VQ;
R0 = ids_core_sources.vacuum_toroidal_field.r0;
B0 = interp1(ids_core_sources.time,ids_core_sources.vacuum_toroidal_field.b0,ohm_t_grid);
-jpar_tilde_to_jpar_0 = R0*T.*Rm2_fs./B0';
+jpar_tilde_to_jpar_0 = R0*T.data.*Rm2_fs.data./B0';
% map volume and conversion factor on rho_pol grid of ohm_data cd_dens
jpar_tilde_to_jpar_0_mapped = interp1(rho_pol_norm_liu,jpar_tilde_to_jpar_0,ohm_data.cd_dens.rhopol_norm');
@@ -180,6 +216,10 @@ if ~isempty(ec_gdat.ec.data) % if EC data available, fill sources
rho_pol_norm_liu = L.pQ; % normalized sqrt poloidal flux including axis
mu0 = 4*pi*10^-7;
Rm2_fs = LY.Q2Q;
+ params_eff = params_eff_ref;
+ params_eff.data_request='\tcv_shot::top.results.equil_1.results:gpsi0_r2_fsa';
+ Rm2_fs = gdat(shot,params_eff);
+
mVprime = 1./LY.Q1Q; % Q1Q: -dpsi/dV
Ip_psi = LY.ItQ; % Ip(psi) = I_tor(psi)
T = LY.TQ;
--
GitLab