From 4e88cbacd58db32b8eda4b4a2113e43cb481c53b Mon Sep 17 00:00:00 2001
From: Olivier Sauter <Olivier.Sauter@epfl.ch>
Date: Fri, 8 Jul 2022 18:07:55 +0200
Subject: [PATCH] add li3 and conversion factor with li efit in AUG requests
---
matlab/AUG/aug_requests_mapping.m | 11 +++++++++++
1 file changed, 11 insertions(+)
diff --git a/matlab/AUG/aug_requests_mapping.m b/matlab/AUG/aug_requests_mapping.m
index 4112db7b..c755d4ba 100644
--- a/matlab/AUG/aug_requests_mapping.m
+++ b/matlab/AUG/aug_requests_mapping.m
@@ -182,6 +182,17 @@ switch lower(data_request)
mapping.label = 'l_i';
mapping.method = 'signal';
mapping.expression = [{'FPG'},{'li'}];
+ case 'li3'
+ mapping.timedim = 1;
+ mapping.label = 'li3';
+ mapping.method = 'expression';
+ mapping.expression = ['params_eff = gdat_data.gdat_params;params_eff.data_request=''li'';' ...
+ 'gdat_tmp=gdat_aug(shot,params_eff);params_eff.data_request=''volume'';r0exp=1.65;gdat_tmp.r0=r0exp;' ...
+ 'gdat_tmp2=gdat_aug(shot,params_eff);ij=find(gdat_tmp2.data==0);gdat_tmp2.data(ij)=NaN;' ...
+ 'params_eff.data_request=[{''FPG''},{''Circumf''}];' ...
+ 'gdat_tmp3=gdat_aug(shot,params_eff);ij=find(gdat_tmp3.data==0);gdat_tmp3.data(ij)=NaN;' ...
+ 'gdat_tmp.conversion_factor_liefit_li3=gdat_tmp.r0./2.*gdat_tmp3.data.^2./gdat_tmp2.data;' ...
+ 'gdat_tmp.data = gdat_tmp.data./gdat_tmp.conversion_factor_liefit_li3;'];
case 'mhd'
mapping.timedim = 1;
mapping.label = {'n\_odd','n\_even'};% this superseeds the expression
--
GitLab