Skip to content
Snippets Groups Projects
Commit 4e7b2f93 authored by Gaelle Tanguy's avatar Gaelle Tanguy Committed by WAUTELET Philippe
Browse files

(Redelsperger & Pianezze)11/2015 Introduction of ponderation coefficient for LBC

parent 1bd90442
No related branches found
No related tags found
No related merge requests found
...@@ -170,6 +170,7 @@ END MODULE MODI_BOUNDARIES ...@@ -170,6 +170,7 @@ END MODULE MODI_BOUNDARIES
!! Modification 01/2015 (JL Redelsperger) Introduction of ponderation !! Modification 01/2015 (JL Redelsperger) Introduction of ponderation
!! for non normal velocity and potential temp !! for non normal velocity and potential temp
!! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1 !! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1
!! Redelsperger & Pianezze : 08/2015 : add XPOND coefficient
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -188,9 +189,11 @@ USE MODD_ELEC_DESCR ...@@ -188,9 +189,11 @@ USE MODD_ELEC_DESCR
USE MODD_ELEC_n USE MODD_ELEC_n
USE MODD_REF_n USE MODD_REF_n
USE MODD_PARAM_n, ONLY : CELEC USE MODD_PARAM_n, ONLY : CELEC
USE MODD_LBC_n, ONLY : XPOND
#ifdef MNH_FOREFIRE #ifdef MNH_FOREFIRE
USE MODD_FOREFIRE, ONLY : LFOREFIRE USE MODD_FOREFIRE, ONLY : LFOREFIRE
#endif #endif
USE MODD_NESTING, ONLY : NDAD
! !
USE MODE_MODELN_HANDLER USE MODE_MODELN_HANDLER
! !
...@@ -439,7 +442,7 @@ END IF ...@@ -439,7 +442,7 @@ END IF
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! PONDERATION COEFF for Non-Normal velocities and pot temperature ! PONDERATION COEFF for Non-Normal velocities and pot temperature
! !
ZPOND = 0.2 ZPOND = XPOND
! !
!* 4. LBC FILLING IN THE X DIRECTION (LEFT WEST SIDE): !* 4. LBC FILLING IN THE X DIRECTION (LEFT WEST SIDE):
! ------------------------------------------------ ! ------------------------------------------------
...@@ -1039,7 +1042,7 @@ IF ( LFOREFIRE .AND. IMI == 1) THEN ...@@ -1039,7 +1042,7 @@ IF ( LFOREFIRE .AND. IMI == 1) THEN
ENDIF ENDIF
#endif #endif
! !
IF ( CELEC /= 'NONE' .AND. IMI == 1) THEN IF ( CELEC /= 'NONE' .AND. (NSV_ELEC_A(NDAD(IMI)) == 0 .OR. IMI == 1)) THEN
CALL ION_BOUNDARIES (HLBCX,HLBCY,PUT,PVT,PSVT) CALL ION_BOUNDARIES (HLBCX,HLBCY,PUT,PVT,PSVT)
ENDIF ENDIF
! !
......
...@@ -41,6 +41,7 @@ ...@@ -41,6 +41,7 @@
!! 26/06/13 (C.Lac) Introduction of CPHASE_PBL !! 26/06/13 (C.Lac) Introduction of CPHASE_PBL
!! 01/04/14 (C.Lac) Introduction of CARPKMAX different to !! 01/04/14 (C.Lac) Introduction of CARPKMAX different to
!! RIMKMAX !! RIMKMAX
!! 08/2015 (Redelsperger & Pianezze) add XPOND coefficient for LBC
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -71,6 +72,7 @@ TYPE LBC_t ...@@ -71,6 +72,7 @@ TYPE LBC_t
! velocity in the PBL if constant ! velocity in the PBL if constant
REAL :: XCARPKMAX ! relaxation coefficient for REAL :: XCARPKMAX ! relaxation coefficient for
!Carpenter, different to XRIMKMAX !Carpenter, different to XRIMKMAX
REAL :: XPOND ! relaxation coefficient for LBC
END TYPE LBC_t END TYPE LBC_t
TYPE(LBC_t), DIMENSION(JPMODELMAX), TARGET, SAVE :: LBC_MODEL TYPE(LBC_t), DIMENSION(JPMODELMAX), TARGET, SAVE :: LBC_MODEL
...@@ -83,6 +85,7 @@ INTEGER, DIMENSION(:), POINTER :: NLBLY=>NULL() ...@@ -83,6 +85,7 @@ INTEGER, DIMENSION(:), POINTER :: NLBLY=>NULL()
REAL, POINTER :: XCPHASE=>NULL() REAL, POINTER :: XCPHASE=>NULL()
REAL, POINTER :: XCPHASE_PBL=>NULL() REAL, POINTER :: XCPHASE_PBL=>NULL()
REAL, POINTER :: XCARPKMAX=>NULL() REAL, POINTER :: XCARPKMAX=>NULL()
REAL, POINTER :: XPOND=>NULL()
CONTAINS CONTAINS
...@@ -107,6 +110,7 @@ NLBLY=>LBC_MODEL(KTO)%NLBLY ...@@ -107,6 +110,7 @@ NLBLY=>LBC_MODEL(KTO)%NLBLY
XCPHASE=>LBC_MODEL(KTO)%XCPHASE XCPHASE=>LBC_MODEL(KTO)%XCPHASE
XCPHASE_PBL=>LBC_MODEL(KTO)%XCPHASE_PBL XCPHASE_PBL=>LBC_MODEL(KTO)%XCPHASE_PBL
XCARPKMAX=>LBC_MODEL(KTO)%XCARPKMAX XCARPKMAX=>LBC_MODEL(KTO)%XCARPKMAX
XPOND=>LBC_MODEL(KTO)%XPOND
END SUBROUTINE LBC_GOTO_MODEL END SUBROUTINE LBC_GOTO_MODEL
......
...@@ -31,6 +31,7 @@ ...@@ -31,6 +31,7 @@
!! right(2) boundaries ( number of delta y) !! right(2) boundaries ( number of delta y)
!! XCPHASE : prescribed value of the phase velocity if constant !! XCPHASE : prescribed value of the phase velocity if constant
!! XCPHASE_PBL : prescribed value of the phase velocity in the PBL if constant !! XCPHASE_PBL : prescribed value of the phase velocity in the PBL if constant
!! XPOND : prescribed value for relaxation of LBC
!! !!
!! REFERENCE !! REFERENCE
!! --------- !! ---------
...@@ -44,6 +45,7 @@ ...@@ -44,6 +45,7 @@
!! ------------- !! -------------
!! Original 13/09/94 !! Original 13/09/94
!! Modification : 26/06/13 (C.Lac) Introduction of CPHASE_PBL !! Modification : 26/06/13 (C.Lac) Introduction of CPHASE_PBL
!! 25/08/2015 (Redelsperger & Pianezze) Introduction of ponderation coefficient for LBC
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
!* 0. DECLARATIONS !* 0. DECLARATIONS
...@@ -56,7 +58,8 @@ USE MODD_LBC_n, ONLY: & ...@@ -56,7 +58,8 @@ USE MODD_LBC_n, ONLY: &
NLBLY_n => NLBLY, & NLBLY_n => NLBLY, &
XCPHASE_n => XCPHASE, & XCPHASE_n => XCPHASE, &
XCPHASE_PBL_n => XCPHASE_PBL, & XCPHASE_PBL_n => XCPHASE_PBL, &
XCARPKMAX_n => XCARPKMAX XCARPKMAX_n => XCARPKMAX, &
XPOND_n => XPOND
! !
IMPLICIT NONE IMPLICIT NONE
! !
...@@ -67,8 +70,9 @@ INTEGER,SAVE, DIMENSION(2) :: NLBLY ...@@ -67,8 +70,9 @@ INTEGER,SAVE, DIMENSION(2) :: NLBLY
REAL,SAVE :: XCPHASE REAL,SAVE :: XCPHASE
REAL,SAVE :: XCPHASE_PBL REAL,SAVE :: XCPHASE_PBL
REAL,SAVE :: XCARPKMAX REAL,SAVE :: XCARPKMAX
REAL,SAVE :: XPOND
! !
NAMELIST/NAM_LBCn/CLBCX,CLBCY,NLBLX,NLBLY,XCPHASE,XCPHASE_PBL,XCARPKMAX NAMELIST/NAM_LBCn/CLBCX,CLBCY,NLBLX,NLBLY,XCPHASE,XCPHASE_PBL,XCARPKMAX,XPOND
! !
CONTAINS CONTAINS
! !
...@@ -80,6 +84,7 @@ SUBROUTINE INIT_NAM_LBCn ...@@ -80,6 +84,7 @@ SUBROUTINE INIT_NAM_LBCn
XCPHASE = XCPHASE_n XCPHASE = XCPHASE_n
XCPHASE_PBL = XCPHASE_PBL_n XCPHASE_PBL = XCPHASE_PBL_n
XCARPKMAX = XCARPKMAX_n XCARPKMAX = XCARPKMAX_n
XPOND = XPOND_n
END SUBROUTINE INIT_NAM_LBCn END SUBROUTINE INIT_NAM_LBCn
SUBROUTINE UPDATE_NAM_LBCn SUBROUTINE UPDATE_NAM_LBCn
...@@ -90,6 +95,7 @@ SUBROUTINE UPDATE_NAM_LBCn ...@@ -90,6 +95,7 @@ SUBROUTINE UPDATE_NAM_LBCn
XCPHASE_n = XCPHASE XCPHASE_n = XCPHASE
XCPHASE_PBL_n = XCPHASE_PBL XCPHASE_PBL_n = XCPHASE_PBL
XCARPKMAX_n = XCARPKMAX XCARPKMAX_n = XCARPKMAX
XPOND_n = XPOND
END SUBROUTINE UPDATE_NAM_LBCn END SUBROUTINE UPDATE_NAM_LBCn
END MODULE MODN_LBC_n END MODULE MODN_LBC_n
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment