Add taur

taur(rho=0.4) used as averaged proxy in B. Labit 2023 IBL paper Should add also volume averaged tau_R, as used in S. Van Mulders DEMO ramp-down paper

matlab/TCV/tcv_requests_mapping.m

@@ -446,6 +446,17 @@ switch lower(data_request)
   mapping.timedim = 1;
   mapping.gdat_timedim = 2;
   mapping.method = 'switchcase';
+ case 'tau_r'
+  mapping.timedim = 1;
+  mapping.gdat_timedim = 1;
+  mapping.method = 'expression';
+  mapping.expression = ['params_eff = gdat_data.gdat_params;params_eff.data_request=''a_minor''; ' ...
+                    'gdat_tmp=gdat_tcv([],params_eff); ' ...
+                    'params_eff = gdat_data.gdat_params;params_eff.data_request=''\results::ibs:signeo''; ' ...
+                    'aa=gdat_tcv([],params_eff);ix=iround_os(aa.x,0.4); ' ...
+                    'bb=interpos(aa.t,aa.data(ix,:),gdat_tmp.t,-100);' ...
+                    'gdat_tmp.data = 4e-7*pi.*gdat_tmp.data.^2.*bb./1.22; gdat_tmp.label=''\tau_R'';gdat_tmp.gdat_request=''tau_r'';' ...
+                    'gdat_tmp.units=''s'';gdat_tmp.help=''tauR=mu0 a^2 / 1.22 etaneo(rhopol=0.4)'''];
  case 'te'
   mapping.timedim = 2;
   mapping.label = 'Te';