From 003521e97cb4690ace5350bed8e6b2b3c711b18a Mon Sep 17 00:00:00 2001
From: VIE Benoit <vie@belenoslogin3.belenoshpc.meteo.fr>
Date: Mon, 20 Feb 2023 15:05:36 +0000
Subject: [PATCH] bugfix timestep in aro_adjust_lima
---
src/arome/ext/aro_adjust_lima.F90 | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/src/arome/ext/aro_adjust_lima.F90 b/src/arome/ext/aro_adjust_lima.F90
index 018e5aa73..0e5c9cdc3 100644
--- a/src/arome/ext/aro_adjust_lima.F90
+++ b/src/arome/ext/aro_adjust_lima.F90
@@ -225,8 +225,8 @@ ZLS(:,:,:)=XLSTT +(XCPV-XCI) *(ZT(:,:,:)-XTT)
ZCPH(:,:,:)=XCPD +XCPV*2.*PTSTEP*PRS(:,:,:,1)
!set concentration for LIMA
-PRS = PRS * PTSTEP
-PSVS = PSVS * PTSTEP
+PRS = PRS * 2.*PTSTEP
+PSVS = PSVS * 2.*PTSTEP
IF (LLIMAINIT) THEN
CALL SET_CONC_LIMA (1,'ICE3',PRHODREF,PRT,PSVT)
CALL SET_CONC_LIMA (1,'ICE3',PRHODREF,PRS,PSVS)
@@ -234,8 +234,8 @@ ELSE
CALL SET_CONC_LIMA_LBC (1,'ICE3',PRHODREF,PRT,PSVT)
CALL SET_CONC_LIMA_LBC (1,'ICE3',PRHODREF,PRS,PSVS)
END IF
-PRS = PRS / PTSTEP
-PSVS = PSVS / PTSTEP
+PRS = PRS / (2.*PTSTEP)
+PSVS = PSVS / (2.*PTSTEP)
!print *, "aro_adjust_lima 2"
!
--
GitLab