From b7ee659d533060aaef333a505fa55a7e7994f601 Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Tue, 22 Jan 2019 17:02:04 +0100
Subject: [PATCH] Philippe 22/01/2019: correct kind of complex datatype
---
src/MNH/bhmie.f90 | 13 +++++++------
1 file changed, 7 insertions(+), 6 deletions(-)
diff --git a/src/MNH/bhmie.f90 b/src/MNH/bhmie.f90
index 60614aefb..5e1d0e340 100644
--- a/src/MNH/bhmie.f90
+++ b/src/MNH/bhmie.f90
@@ -1,6 +1,6 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
! #################
MODULE MODI_BHMIE
@@ -63,6 +63,7 @@ END MODULE MODI_BHMIE
!! portable. In event that portable version is
!! needed, use src/bhmie_f77.f
!! 93/06/01 (BTD): Changed AMAX1 to generic function MAX
+!! 22/01/2019 (P.Wautelet): correct kind of complex datatype
!!***********************************************************************
!
!* 0. DECLARATIONS
@@ -187,8 +188,8 @@ DO J = 1,ISTOP
DO JJ = 1,KNANG
ZPI(JJ) = ZPI1(JJ)
ZTAU(JJ) = ZEN*ZAMU(JJ)*ZPI(JJ) - (ZEN+1.)*ZPI0(JJ)
- PPS1(JJ) = PPS1(JJ) + ZFN*(ZZAN*ZPI(JJ)+ZZBN*ZTAU(JJ))
- PPS2(JJ) = PPS2(JJ) + ZFN*(ZZAN*ZTAU(JJ)+ZZBN*ZPI(JJ))
+ PPS1(JJ) = PPS1(JJ) + CMPLX(ZFN*(ZZAN*ZPI(JJ)+ZZBN*ZTAU(JJ)),kind=kind(PPS1(1)))
+ PPS2(JJ) = PPS2(JJ) + CMPLX(ZFN*(ZZAN*ZTAU(JJ)+ZZBN*ZPI(JJ)),kind=kind(PPS2(1)))
ENDDO
!
!*** Now do angles greater than 90 using PI and TAU from
@@ -198,8 +199,8 @@ DO J = 1,ISTOP
ZONE = -ZONE
DO JJ = 1,KNANG-1
JJJ = 2*KNANG-JJ
- PPS1(JJJ) = PPS1(JJJ) + ZFN*ZONE*(ZZAN*ZPI(JJ)-ZZBN*ZTAU(JJ))
- PPS2(JJJ) = PPS2(JJJ) + ZFN*ZONE*(ZZBN*ZPI(JJ)-ZZAN*ZTAU(JJ))
+ PPS1(JJJ) = PPS1(JJJ) + CMPLX(ZFN*ZONE*(ZZAN*ZPI(JJ)-ZZBN*ZTAU(JJ)),kind=kind(PPS1(1)))
+ PPS2(JJJ) = PPS2(JJJ) + CMPLX(ZFN*ZONE*(ZZBN*ZPI(JJ)-ZZAN*ZTAU(JJ)),kind=kind(PPS2(1)))
ENDDO
ZPSI0 = ZPSI1
ZPSI1 = ZPSI
--
GitLab