From f90db7229df528c41f33fb5af44b7987d9fe1768 Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Thu, 26 Nov 2015 15:04:08 +0000
Subject: [PATCH] C.Lac 11/2015 : Correction on aerosols for radiation + New
 LES diagnostics

---
 src/MNH/ini_modeln.f90 | 18 +++++++++++++-----
 1 file changed, 13 insertions(+), 5 deletions(-)

diff --git a/src/MNH/ini_modeln.f90 b/src/MNH/ini_modeln.f90
index e931f1404..ebb581f09 100644
--- a/src/MNH/ini_modeln.f90
+++ b/src/MNH/ini_modeln.f90
@@ -308,7 +308,7 @@ USE MODD_DIM_n
 USE MODD_BUDGET
 USE MODD_RADIATIONS_n
 USE MODD_SHADOWS_n
-USE MODD_PARAM_RAD_n,   ONLY : CLW, CAER, LRAD_DUST
+USE MODD_PARAM_RAD_n,   ONLY : CLW, CAER, CAOP     
 USE MODD_VAR_ll,        ONLY : IP
 !
 USE MODD_STAND_ATM,     ONLY : XSTROATM, XSMLSATM, XSMLWATM, XSPOSATM, XSPOWATM
@@ -700,11 +700,21 @@ IF (CTURB /= 'NONE') THEN
   ALLOCATE(XRTKES(IIU,IJU,IKU))
   ALLOCATE(XRTKEMS(IIU,IJU,IKU)); XRTKEMS = 0.0
   ALLOCATE(XWTHVMF(IIU,IJU,IKU))
+  ALLOCATE(XDYP(IIU,IJU,IKU))
+  ALLOCATE(XTHP(IIU,IJU,IKU))
+  ALLOCATE(XTR(IIU,IJU,IKU))
+  ALLOCATE(XDISS(IIU,IJU,IKU))
+  ALLOCATE(XLEM(IIU,IJU,IKU))
   XTKEMIN=XKEMIN
 ELSE
   ALLOCATE(XTKET(0,0,0))
   ALLOCATE(XRTKES(0,0,0))
   ALLOCATE(XWTHVMF(0,0,0))
+  ALLOCATE(XDYP(0,0,0))
+  ALLOCATE(XTHP(0,0,0))
+  ALLOCATE(XTR(0,0,0))
+  ALLOCATE(XDISS(0,0,0))
+  ALLOCATE(XLEM(0,0,0))
 END IF
 IF (CTOM == 'TM06') THEN
   ALLOCATE(XBL_DEPTH(IIU,IJU))
@@ -1967,10 +1977,8 @@ IF (CRAD   == 'ECMW') THEN
                                XSTATM,ZCOVER,XOZON, XAER,XDST_WL, LSUBG_COND              )
 !
     DEALLOCATE(ZCOVER)
-    IF (LRAD_DUST) THEN
-       ALLOCATE (XAER_CLIM(SIZE(XAER,1),SIZE(XAER,2),SIZE(XAER,3),SIZE(XAER,4)))
-        XAER_CLIM(:,:,:,:) =XAER(:,:,:,:)
-    END IF    
+    ALLOCATE (XAER_CLIM(SIZE(XAER,1),SIZE(XAER,2),SIZE(XAER,3),SIZE(XAER,4)))
+    XAER_CLIM(:,:,:,:) =XAER(:,:,:,:)
 !
   END IF
 ELSE
-- 
GitLab