From 4a0421531a4a8a41552dad7ea80f04b25eb7afec Mon Sep 17 00:00:00 2001
From: Juan Escobar <juan.escobar@aero.obs-mip.fr>
Date: Thu, 1 Oct 2015 14:46:34 +0000
Subject: [PATCH] Juan 1/10/2015: add computation of CFL for L1D case
---
src/MNH/advection_metsv.f90 | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/src/MNH/advection_metsv.f90 b/src/MNH/advection_metsv.f90
index 01f9a90dd..5cee9d2d8 100644
--- a/src/MNH/advection_metsv.f90
+++ b/src/MNH/advection_metsv.f90
@@ -127,6 +127,7 @@ END MODULE MODI_ADVECTION_METSV
!! 09/2014 (G.Delautier) close OUTPUT_LISTING before STOP
!! 04/2015 (J.Escobar) remove/commente some NHALO=1 test
!! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1
+!! J.Escobar : 01/10/2015 : add computation of CFL for L1D case
!!
!-------------------------------------------------------------------------------
!
@@ -275,6 +276,7 @@ END IF
!
!* 2.2 computes CFL numbers
!
+
IF (.NOT. L1D) THEN
ZCFLU = 0.0 ; ZCFLV = 0.0 ; ZCFLW = 0.0
ZCFLU(IIB:IIE,IJB:IJE,:) = ABS(ZRUCPPM(IIB:IIE,IJB:IJE,:) * PTSTEP)
@@ -285,6 +287,10 @@ IF (.NOT. L1D) THEN
ELSE
ZCFL = SQRT(ZCFLU**2+ZCFLW**2)
END IF
+ELSE
+ ZCFLU = 0.0 ; ZCFLV = 0.0 ; ZCFLW = 0.0
+ ZCFLW(IIB:IIE,IJB:IJE,:) = ABS(ZRWCPPM(IIB:IIE,IJB:IJE,:) * PTSTEP)
+ ZCFL = SQRT(ZCFLW**2)
END IF
!
!* prints in the file the 3D Courant numbers (one should flag this)
--
GitLab