From 9cd59c0ec2a3d0d5f2b09f1e1236e9e21f928c00 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Mon, 26 Oct 2020 16:14:30 +0100
Subject: [PATCH] Philippe 26/10/2020: add missing RES_RTPZ variable

---
 src/MNH/write_lesn.f90 | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/src/MNH/write_lesn.f90 b/src/MNH/write_lesn.f90
index 2ef45d081..9b8a50559 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)
-- 
GitLab