diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index cba8b5ff64f757c7bbdc9956c11eec7916984172..63fb9ad663d33d2adee89c859040dcd40f75fb49 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -1560,6 +1560,42 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
     % add grids_1d to have rhotor, etc
     gdat_data = get_grids_1d(gdat_data,1,1,gdat_params.nverbose);
 
+    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+   case {'rbphi_rho', 'rbtor_rho'}
+    % R*Bphi(rho,t) from F from FFprime
+    if liuqe_matlab==0
+      disp('not yet implemented for liuqe fortran')
+      return
+    else
+      nodenameeff=['tcv_eq(''rbtor_rho'',''LIUQE.M' substr_liuqe_tcv_eq ''')'];
+    end
+    if liuqe_version_eff==-1
+      disp('not yet implemented for fbte')
+      return
+    end
+    tracetdi=tdi(nodenameeff);
+    if isempty(tracetdi.data) || isempty(tracetdi.dim)  % || ischar(tracetdi.data) (to add?)
+      if (gdat_params.nverbose>=1); warning(['problems loading data for ' nodenameeff ' for data_request= ' data_request_eff]); end
+      if (gdat_params.nverbose>=3); disp(['rerun LIUQE?']); end
+      return
+    end   
+    gdat_data.data = tracetdi.data;
+    gdat_data.dim = tracetdi.dim;
+    gdat_data.t = gdat_data.dim{2};
+    gdat_data.data_fullpath=[nodenameeff ' on rhopol'];
+    if liuqe_matlab==0
+      rhopol_eff = ones(size(tracetdi.dim{1}));
+      rhopol_eff(:) = sqrt(linspace(0,1,length(tracetdi.dim{1})));
+      gdat_data.dim{1} = rhopol_eff;
+    end      
+    gdat_data.x = gdat_data.dim{1};
+    gdat_data.dimunits{1} = 'rho_pol~sqrt(\psi_norm)';
+    gdat_data.dimunits{2} = 's';
+    gdat_data.units = '';
+    gdat_data.request_description = nodenameeff;
+    % add grids_1d to have rhotor, etc
+    gdat_data = get_grids_1d(gdat_data,1,1,gdat_params.nverbose);
+
    case {'phi_tor', 'phitor', 'toroidal_flux'}
     % Phi(LCFS) = int(Bphi dSphi), can use Eq.(11) of "Tokamak coordinate conventions: COCOS" paper:
     % O. Sauter, S.Yu. Medvedev, Comput. Phys. Commun. 184 (2013) 293–302
diff --git a/crpptbx/TCV/tcv_requests_mapping.m b/crpptbx/TCV/tcv_requests_mapping.m
index ecc6774e1b32e95ab030aba635089e41d8ce72a9..51505bb27c3ffb4960d285484037bcbb96cb111e 100644
--- a/crpptbx/TCV/tcv_requests_mapping.m
+++ b/crpptbx/TCV/tcv_requests_mapping.m
@@ -269,6 +269,10 @@ switch lower(data_request)
   mapping.timedim = 2;
   mapping.label = 'q';
   mapping.method = 'switchcase';
+ case {'rbphi_rho', 'rbtor_rho'}
+  mapping.timedim = 2;
+  mapping.label = 'R B_\phi';
+  mapping.method = 'switchcase';
  case 'r_contour'
   mapping.timedim = 2;
   mapping.method = 'tdiliuqe';