From 93b1cd1c17e88c452a0707e10e7fd29b09d408fc Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Thu, 26 Nov 2015 14:49:00 +0000
Subject: [PATCH] C.Lac 11/2015 : correction for ZP_QS
---
src/SURFEX/diag_surf_atmn.F90 | 9 ++++++++-
1 file changed, 8 insertions(+), 1 deletion(-)
diff --git a/src/SURFEX/diag_surf_atmn.F90 b/src/SURFEX/diag_surf_atmn.F90
index 213f4199e..53224105b 100644
--- a/src/SURFEX/diag_surf_atmn.F90
+++ b/src/SURFEX/diag_surf_atmn.F90
@@ -27,6 +27,7 @@ SUBROUTINE DIAG_SURF_ATM_n(HPROGRAM)
!! Original 01/2004
!! Modified 01/2006 : sea flux parameterization.
!! Modified 08/2008 : cumulated fluxes
+!! C.Lac 11/2015 : correction for ZP_QS
!!------------------------------------------------------------------
!
@@ -465,11 +466,17 @@ IF (LCOEF) THEN
XCD_TILE (KMASK(JJ),KTILE) = ZP_CD (JJ)
XCH_TILE (KMASK(JJ),KTILE) = ZP_CH (JJ)
XCE_TILE (KMASK(JJ),KTILE) = ZP_CE (JJ)
- XQS_TILE (KMASK(JJ),KTILE) = ZP_QS (JJ)
+! XQS_TILE (KMASK(JJ),KTILE) = ZP_QS (JJ)
XZ0_TILE (KMASK(JJ),KTILE) = ZP_Z0 (JJ)
XZ0H_TILE (KMASK(JJ),KTILE) = ZP_Z0H (JJ)
ENDDO
END IF
+!
+IF (LSURF_VARS) THEN
+ DO JJ=1,KSIZE
+ XQS_TILE (KMASK(JJ),KTILE) = ZP_QS (JJ)
+ ENDDO
+END IF
!----------------------------------------------------------------------
IF (LHOOK) CALL DR_HOOK('DIAG_SURF_ATM_n:TREAT_SURF',1,ZHOOK_HANDLE)
!
--
GitLab