From 46f37a39cca34c9a1280a8fbae0dcea70f2b3620 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Wed, 8 Mar 2023 13:26:41 +0100
Subject: [PATCH] Pierre T. 08/03/2023: bugfix LIMA with aerosols PHYEX part
---
src/PHYEX/micro/init_aerosol_properties.f90 | 10 +++-------
1 file changed, 3 insertions(+), 7 deletions(-)
diff --git a/src/PHYEX/micro/init_aerosol_properties.f90 b/src/PHYEX/micro/init_aerosol_properties.f90
index eeb747265..f7ea0bad7 100644
--- a/src/PHYEX/micro/init_aerosol_properties.f90
+++ b/src/PHYEX/micro/init_aerosol_properties.f90
@@ -44,7 +44,6 @@ END MODULE MODI_INIT_AEROSOL_PROPERTIES
!* 0. DECLARATIONS
! ------------
!
-!USE MODD_LUNIT, ONLY : TLUOUT0
USE MODD_PARAM_LIMA, ONLY : NMOD_CCN, HINI_CCN, HTYPE_CCN, &
XR_MEAN_CCN, XLOGSIG_CCN, XRHO_CCN, &
XKHEN_MULTI, XMUHEN_MULTI, XBETAHEN_MULTI, &
@@ -90,8 +89,6 @@ REAL, DIMENSION(3) :: RHOCCN
!
INTEGER :: I,J,JMOD
!
-!INTEGER :: ILUOUT0 ! Logical unit number for output-listing
-!INTEGER :: IRESP ! Return code of FM-routines
!
REAL :: X1, X2, X3, X4, X5
! REAL, DIMENSION(7) :: diameters=(/ 0.01E-6, 0.05E-6, 0.1E-6, 0.2E-6, 0.5E-6, 1.E-6, 2.E-6 /)
@@ -103,7 +100,6 @@ INTEGER :: II, IJ, IK
!
!-------------------------------------------------------------------------------
!
-!ILUOUT0 = TLUOUT0%NLU
!
!!!!!!!!!!!!!!!!
! CCN properties
@@ -410,10 +406,10 @@ IF (LORILAM) THEN
XMDIAM_IFN(3) = 2 * XINIRADIUSI * 1E-6
XMDIAM_IFN(4) = 2 * XINIRADIUSI * 1E-6
END IF
- XSIGMA_IFN(3) = LOG(XINISIGJ)
- XSIGMA_IFN(4) = LOG(XINISIGJ)
XRHO_IFN(3) = XRHOI(JP_AER_BC)
XRHO_IFN(4) = XRHOI(JP_AER_OC)
+ XSIGMA_IFN(3) = XINISIGI
+ XSIGMA_IFN(4) = XINISIGI
END IF
IF (LDUST) THEN
@@ -424,7 +420,7 @@ IF (LDUST) THEN
ELSE
XMDIAM_IFN(JMOD) = 2 * XINIRADIUS(JPDUSTORDER(JMOD)) * 1E-6
END IF
- XSIGMA_IFN(JMOD) = LOG(XINISIG(JPDUSTORDER(JMOD)))
+ XSIGMA_IFN(JMOD) = XINISIG(JPDUSTORDER(JMOD))
XRHO_IFN(JMOD) = XDENSITY_DUST
ENDDO
END IF
--
GitLab