Skip to content
Snippets Groups Projects
Commit d7b20059 authored by WAUTELET Philippe's avatar WAUTELET Philippe Committed by RODIER Quentin
Browse files

Philippe 01/06/2023: PHYEX: PWSV was not computed (most of the time) if there...

Philippe 01/06/2023: PHYEX: PWSV was not computed (most of the time) if there are flyers (PWSV is needed by flyers)
parent 2fa68610
No related branches found
No related tags found
No related merge requests found
!MNH_LIC Copyright 1994-2022 CNRS, Meteo-France and Universite Paul Sabatier !MNH_LIC Copyright 1994-2023 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 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. !MNH_LIC for details. version 1.
...@@ -204,6 +204,7 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, & ...@@ -204,6 +204,7 @@ SUBROUTINE TURB_VER_SV_FLUX(D,CST,CSTURB,TURBN,TLES,ONOMIXLG, &
!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O !! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
! Wim de Rooy 06/2019: with energycascade, 50MF nog longer necessary ! Wim de Rooy 06/2019: with energycascade, 50MF nog longer necessary
! P. Wautelet 30/11/2022: compute PWSV only when needed ! P. Wautelet 30/11/2022: compute PWSV only when needed
! P. Wautelet 01/06/2023: fix for PWSV
!!-------------------------------------------------------------------------- !!--------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -396,7 +397,7 @@ DO JSV=1,KSV ...@@ -396,7 +397,7 @@ DO JSV=1,KSV
PRHODJ(IIJB:IIJE,1:IKT)*(ZRES(IIJB:IIJE,1:IKT)-PSVM(IIJB:IIJE,1:IKT,JSV))/PTSTEP PRHODJ(IIJB:IIJE,1:IKT)*(ZRES(IIJB:IIJE,1:IKT)-PSVM(IIJB:IIJE,1:IKT,JSV))/PTSTEP
!$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT) !$mnh_end_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
! !
IF ( (TURBN%LTURB_FLX .AND. TPFILE%LOPENED) .OR. TLES%LLES_CALL ) THEN IF ( (TURBN%LTURB_FLX .AND. TPFILE%LOPENED) .OR. TLES%LLES_CALL .OR. OFLYER ) THEN
! Diagnostic of the cartesian vertical flux ! Diagnostic of the cartesian vertical flux
! !
!$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT) !$mnh_expand_array(JIJ=IIJB:IIJE,JK=1:IKT)
...@@ -444,7 +445,7 @@ DO JSV=1,KSV ...@@ -444,7 +445,7 @@ DO JSV=1,KSV
PWSV(IIJB:IIJE,IKE,JSV)=PWSV(IIJB:IIJE,IKE-IKL,JSV) PWSV(IIJB:IIJE,IKE,JSV)=PWSV(IIJB:IIJE,IKE-IKL,JSV)
!$mnh_end_expand_array(JIJ=IIJB:IIJE) !$mnh_end_expand_array(JIJ=IIJB:IIJE)
END IF END IF
END IF END IF
! !
IF (TURBN%LTURB_FLX .AND. TPFILE%LOPENED) THEN IF (TURBN%LTURB_FLX .AND. TPFILE%LOPENED) THEN
! stores the JSVth vertical flux ! stores the JSVth vertical flux
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment