From 75709e37253dcc5e0c0a599cf3d603ad41433ddc Mon Sep 17 00:00:00 2001
From: Jean-Pierre Chaboureau <jean-pierre.chaboureau@aero.obs-mip.fr>
Date: Mon, 18 Sep 2017 09:35:32 +0200
Subject: [PATCH] JP Chaboureau 30/05/2017 exclude the first layer when
considering clouds
---
src/MNH/call_rttov11.f90 | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/src/MNH/call_rttov11.f90 b/src/MNH/call_rttov11.f90
index 220a7f2fe..69bcd9fbe 100644
--- a/src/MNH/call_rttov11.f90
+++ b/src/MNH/call_rttov11.f90
@@ -78,6 +78,7 @@ SUBROUTINE CALL_RTTOV11(KDLON, KFLEV, PEMIS, PTSRAD, &
!! JP Chaboureau 27/03/2008 Vectorization
!! JP Chaboureau 02/11/2009 move GANGL deallocation outside the sensor loop
!! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1
+!! JP Chaboureau 30/05/2017 exclude the first layer when considering clouds
!!----------------------------------------------------------------------------
!!
!!* 0. DECLARATIONS
@@ -318,6 +319,9 @@ DO JSAT=1,IJSAT ! loop over sensors
ELSE
opts % rt_ir % addclouds = .FALSE. ! Include cloud effects
END IF
+ opts % config % verbose = .FALSE. ! Enable printing of warnings
+ opts % config % do_checkinput = .FALSE.
+
! Read and initialise coefficients
! -----------------------------------------------------------------------------
@@ -461,7 +465,7 @@ DO JSAT=1,IJSAT ! loop over sensors
IF( coef_rttov%coef% id_sensor /= sensor_id_mw) THEN
profiles(1)%ish = 2 ! Aggregates
profiles(1)%idg = 4 ! McFarquar et al (2003)
- DO JK=IKB,IKE-1 ! nlayers
+ DO JK=IKB+1,IKE-1 ! nlayers
JKRAD = nlev-JK+1 !INVERSION OF VERTICAL LEVELS!
profiles(1) %cfrac(JKRAD) = PCLDFR(JI,JJ,JK)
profiles(1) %cloud(1,JKRAD) = PRT(JI,JJ,JK,2)*XRHODREF(JI,JJ,JK)*1.0E03
@@ -571,7 +575,7 @@ DO JSAT=1,IJSAT ! loop over sensors
! ENDIF
IGRID =1
ILENCH =LEN(YCOMMENT)
- PRINT *,'YRECFM='//TRIM(YRECFM)
+! PRINT *,'YRECFM='//TRIM(YRECFM)
CALL FMWRIT(HFMFILE,YRECFM,CLUOUT,'XY',ZBT(:,:,JCH), &
IGRID,ILENCH,YCOMMENT,IRESP)
END DO
--
GitLab