diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index f521f6fda1a608e3912cdff14b6b84245f801754..7725d0be23a2d80732e1acbd168bf205fa48a074 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -1927,6 +1927,29 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
     gdat_data.units = tracetdi.units;
     gdat_data.request_description = 'pprime=dp/dpsi';
 
+   case {'pressure', 'pressure_rho'}
+    if liuqe_matlab==0
+      nodenameeff = ['tcv_eq("' liuqefortran2liuqematlab('p_rho',1,0) '",''' psitbx_str ''')'];
+      tracetdi=tdi(nodenameeff);
+      if ~isempty(tracetdi.dim) && length(tracetdi.dim)>=1
+        tracetdi.dim{1} = sqrt(tracetdi.dim{1}); % correct x-axis psi_norm to rhopol
+      end
+      tracetdi.data = tracetdi.data ./2 ./pi; % correct node assumption (same for liuqe_fortran and fbte)
+    else
+      nodenameeff = ['tcv_eq(''p_rho'',''' psitbx_str ''')'];
+      tracetdi=tdi(nodenameeff);
+    end
+    gdat_data.data = tracetdi.data;
+    gdat_data.dim = tracetdi.dim;
+    if ~isempty(gdat_data.dim) && length(gdat_data.dim)>=2
+      gdat_data.t = gdat_data.dim{mapping_for_tcv.gdat_timedim};
+      gdat_data.x = gdat_data.dim{setdiff([1 2],mapping_for_tcv.gdat_timedim)};
+    end
+    gdat_data.data_fullpath=nodenameeff;
+    gdat_data.dimunits = tracetdi.dimunits;
+    gdat_data.units = tracetdi.units;
+    gdat_data.request_description = 'pressure';
+
     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    case {'psi_edge'}
     % psi at edge, 0 by construction in Liuqe, thus not given
diff --git a/crpptbx/TCV/tcv_requests_mapping.m b/crpptbx/TCV/tcv_requests_mapping.m
index 821097a4d053cc8c844393c9f36de4563bab2085..033c0ca62638b989b0e6a2fcac178e73d4a653f2 100644
--- a/crpptbx/TCV/tcv_requests_mapping.m
+++ b/crpptbx/TCV/tcv_requests_mapping.m
@@ -245,6 +245,10 @@ switch lower(data_request)
   mapping.timedim = 2;
   mapping.label = 'pprime';
   mapping.method = 'switchcase';
+ case {'pressure', 'p_rho'} % note: not pressure from liuqe fortran which is 2D
+  mapping.timedim = 2;
+  mapping.label = 'pressure';
+  mapping.method = 'switchcase';
  case {'psi_axis', 'psi_mag'}
   mapping.timedim = 1;
   mapping.method = 'tdiliuqe';