diff --git a/src/MNH/write_lesn.f90 b/src/MNH/write_lesn.f90
index 2ef45d081b520929b6694c2ee28c0b61441f2cac..9b8a50559b7598ee1ab98b01ac4d7d7d6f967932 100644
--- a/src/MNH/write_lesn.f90
+++ b/src/MNH/write_lesn.f90
@@ -58,7 +58,7 @@ END MODULE MODI_WRITE_LES_n
 !  P. Wautelet 05/2016-04/2018: new data structures and calls for I/O
 !  C. Lac         02/2019: add rain fraction as a LES diagnostic
 !  P. Wautelet 13/10/2020: bugfix: correct some names for LES_DIACHRO_2PT diagnostics (Ri)
-!  P. Wautelet 26/10/2020: bugfix: correct some comments and conditions
+!  P. Wautelet 26/10/2020: bugfix: correct some comments and conditions + add missing RES_RTPZ
 !! --------------------------------------------------------------------------
 !       
 !*      0. DECLARATIONS
@@ -746,6 +746,10 @@ IF (LLES_RESOLVED) THEN
   CALL LES_DIACHRO_MASKS(TPDIAFILE,"RES_TLPZ",YSUBTITLE(:),  &
        "Resolved <Thldp/dz>"//YSUBTITLE(:),"K Pa m-1",XLES_RESOLVED_ThlPz,HLES_AVG)
 
+  IF (LUSERC) &
+  CALL LES_DIACHRO_MASKS(TPDIAFILE,"RES_RTPZ",YSUBTITLE(:),  &
+       "Resolved <Rtdp/dz>"//YSUBTITLE(:),"kg2 kg-2 Pa m-1",XLES_RESOLVED_RtPz,HLES_AVG)
+
   IF (LUSERV) &
   CALL LES_DIACHRO_MASKS(TPDIAFILE,"RES_RVPZ",YSUBTITLE(:),  &
        "Resolved <Rvdp/dz>"//YSUBTITLE(:),"kg2 kg-2 Pa m-1",XLES_RESOLVED_RvPz,HLES_AVG)