G.Delautier 19/09/2017 : bug with CLOUD_ONLY
authorGaelle DELAUTIER <gaelle.delautier@meteo.fr>
Tue, 19 Sep 2017 08:47:20 +0000 (10:47 +0200)
committerGaelle DELAUTIER <gaelle.delautier@meteo.fr>
Tue, 19 Sep 2017 08:47:20 +0000 (10:47 +0200)
src/MNH/radiations.f90

index 880f903..b2de10d 100644 (file)
@@ -208,6 +208,7 @@ END MODULE MODI_RADIATIONS
 !!      C.Lac       11/2015   Correction on aerosols
 !!      J.Escobar 30/03/2017  : Management of compilation of ECMWF_RAD in REAL*8 with MNH_REAL=R4
 !!      J.Escobar 29/06/2017  : Check if Pressure Decreasing with height <-> elsif PB & STOP 
+!!      Q.LIBOIS  06/2017     : correction on CLOUD_ONLY
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -2393,8 +2394,17 @@ END DO
 !final  THETA_radiative tendency and surface fluxes 
 !
 IF(OCLOUD_ONLY) THEN
-  !
-  ZWORKL(:,:) = SUM(PCLDFR(:,:,:),DIM=3) > 0.0
+  !! Q.LIBOIS 06/2017
+  !ZWORKL(:,:) = SUM(PCLDFR(:,:,:),DIM=3) > 0.0
+  DO JJ=IJB,IJE
+    DO JI=IIB,IIE
+        IIJ = 1 + (JI-IIB) + (IIE-IIB+1)*(JJ-IJB)
+        GCLOUD_SURF(JI,JJ) = GCLOUD(IIJ,1)
+    END DO
+  END DO
+  ZWORKL(:,:) = GCLOUD_SURF(:,:) ! nouvelle condition
+  !! Q.LIBOIS 06/2017
   DO JK = IKB,IKE
     WHERE( ZWORKL(:,:) )
       PDTHRAD(:,:,JK) = (ZDTRAD_LW(:,:,JK)+ZDTRAD_SW(:,:,JK))/ZEXNT(:,:,JK)