From 0a5ebbb7a10e04f926237ecdbc642ae8af8f6ca6 Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Wed, 14 Dec 2022 17:00:46 +0100
Subject: [PATCH] Quentin 14/12/2022: correction of delete of modd_ref :
ini_snow use XTHVREF as 300.0 hard coded
---
src/arome/ext/aroini_micro.F90 | 3 +--
src/common/micro/ini_snow.F90 | 9 ++-------
2 files changed, 3 insertions(+), 9 deletions(-)
diff --git a/src/arome/ext/aroini_micro.F90 b/src/arome/ext/aroini_micro.F90
index 7bc27187a..347033c9c 100644
--- a/src/arome/ext/aroini_micro.F90
+++ b/src/arome/ext/aroini_micro.F90
@@ -64,7 +64,6 @@ USE MODD_RAIN_ICE_DESCR
USE MODD_RAIN_ICE_PARAM
USE MODD_PARAM_ICE
USE MODD_PARAM_C1R3
-USE MODD_REF, ONLY: XTHVREFZ
USE MODI_INI_RAIN_ICE
USE MODI_INI_TIWMX
@@ -156,7 +155,7 @@ IF(PRESENT(PFRMIN))THEN
WRITE(UNIT=KULOUT,FMT='('' UPDATED VALUES OF XFRMIN FROM NAMPARAR :'')')
WRITE(UNIT=KULOUT,FMT='('' XFRMIN = '',40E10.3)') XFRMIN
IF(XFRMIN(16) > 0.) THEN
- CALL INI_SNOW(KULOUT,XTHVREFZ) ! Recalculate snow parameters : XCCS = XFRMIN(16),XCXS = XFRMIN(17)
+ CALL INI_SNOW(KULOUT) ! Recalculate snow parameters : XCCS = XFRMIN(16),XCXS = XFRMIN(17)
ENDIF
ENDIF
diff --git a/src/common/micro/ini_snow.F90 b/src/common/micro/ini_snow.F90
index 8ff219640..c61f54a0d 100644
--- a/src/common/micro/ini_snow.F90
+++ b/src/common/micro/ini_snow.F90
@@ -1,5 +1,5 @@
! ######spl
- SUBROUTINE INI_SNOW ( KLUOUT, PTHVREFZ )
+ SUBROUTINE INI_SNOW ( KLUOUT )
USE PARKIND1, ONLY : JPRB
USE YOMHOOK , ONLY : LHOOK, DR_HOOK
! ###########################################################
@@ -69,12 +69,8 @@ IMPLICIT NONE
!
!
INTEGER, INTENT(IN) :: KLUOUT ! Logical unit number for prints
-REAL, DIMENSION(:) :: PTHVREFZ
!* 0.2 Declarations of local variables :
!
-INTEGER :: IKB ! Coordinates of the first physical
- ! points along z
-
REAL :: ZRHO00 ! Surface reference air density
REAL :: ZCONC_MAX ! Maximal concentration for snow
@@ -86,8 +82,7 @@ IF (LHOOK) CALL DR_HOOK('INI_RAIN_ICE',0,ZHOOK_HANDLE)
XCCS = XFRMIN(16)
XCXS = XFRMIN(17)
-IKB = 1 + JPVEXT
-ZRHO00 = XP00/(XRD*PTHVREFZ(IKB))
+ZRHO00 = XP00/(XRD*300.0)
! recalculate ini_rain_ice stuff:
! 3.4 Constants for shape parameter
--
GitLab