Skip to content
Snippets Groups Projects
Commit 9dc96e65 authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Quentin 29/05/2024: add kernels on turb_ver_dyn_flux and complete testing of... 

Quentin 29/05/2024: add kernels on turb_ver_dyn_flux and complete testing of shuman/gradients conversion function to routines with the last version of pyft (of 29/05/2024)
parent 348bf320
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
src/PHYEX/aux/shuman_phy.f90
+ 203
− 1
View file @ 9dc96e65
......@@ -1397,9 +1397,110 @@ DO JI=1,JPHEXT
END DO
END DO
!
END SUBROUTINE DXM_PHY
!-------------------------------------------------------------------------------
!
END SUBROUTINE DXM_PHY
SUBROUTINE DXM2D_PHY(D,PA,PDXM)
! ###############################
!
!!**** *DXM* - Shuman operator : finite difference operator in x direction
!! for a variable at a mass localization
!!
!! PURPOSE
!! -------
! The purpose of this function is to compute a finite difference
! along the x direction (I index) for a field PA localized at a mass
! point. The result is localized at a x-flux point (u point).
!
!!** METHOD
!! ------
!! The result PDXM(i,:,:) is defined by (PA(i,:,:)-PA(i-1,:,:))
!! At i=1, PDXM(1,:,:) are replaced by the values of PDXM,
!! which are the right values in the x-cyclic case.
!!
!!
!! EXTERNAL
!! --------
!! NONE
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!! Module MODD_PARAMETERS: declaration of parameter variables
!! JPHEXT: define the number of marginal points out of the
!! physical domain along the horizontal directions.
!!
!! REFERENCE
!! ---------
!! Book2 of documentation of Meso-NH (SHUMAN operators)
!! Technical specifications Report of The Meso-NH (chapters 3)
!!
!!
!! AUTHOR
!! ------
!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
!! Original 05/07/94
!! Modification to include the periodic case 13/10/94 J.Stein
!! optimisation 20/08/00 J. Escobar
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
!
USE MODD_PARAMETERS, ONLY: JPHEXT
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
!
!* 0.1 Declarations of argument and result
! ------------------------------------
!
TYPE(DIMPHYEX_t), INTENT(IN) :: D
REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PA ! variable at mass
! localization
REAL, DIMENSION(D%NIT,D%NJT), INTENT(OUT) :: PDXM ! result at flux
! side
!
!* 0.2 Declarations of local variables
! -------------------------------
!
!
INTEGER :: JI ! Loop index in x direction
INTEGER :: IIU ! upper bound in x direction of PA
!
INTEGER :: JJK,IJU
INTEGER :: JIJK,JIJKOR,JIJKEND
!
!-------------------------------------------------------------------------------
!
!* 1. DEFINITION OF DXM
! ------------------
!
IIU = SIZE(PA,1)
IJU = SIZE(PA,2)
!
JIJKOR = 1 + 1
JIJKEND = IIU*IJU
!
!CDIR NODEP
!OCL NOVREC
DO JIJK=JIJKOR , JIJKEND
PDXM(JIJK,1) = PA(JIJK,1) - PA(JIJK-1,1)
END DO
!
!CDIR NODEP
!OCL NOVREC
DO JI=1,JPHEXT
DO JJK=1,IJU
PDXM(JI,JJK) = PDXM(IIU-2*JPHEXT+JI,JJK) ! for reprod JPHEXT <> 1
END DO
END DO
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE DXM2D_PHY
!
! ###############################
SUBROUTINE DYM_PHY(D,PA,PDYM)
......@@ -1503,6 +1604,107 @@ END DO
!-------------------------------------------------------------------------------
!
END SUBROUTINE DYM_PHY
! ###############################
SUBROUTINE DYM2D_PHY(D,PA,PDYM)
! ###############################
!
!!**** *DYM* - Shuman operator : finite difference operator in y direction
!! for a variable at a mass localization
!!
!! PURPOSE
!! -------
! The purpose of this function is to compute a finite difference
! along the y direction (J index) for a field PA localized at a mass
! point. The result is localized at a y-flux point (v point).
!
!!** METHOD
!! ------
!! The result PDYM(:,j,:) is defined by (PA(:,j,:)-PA(:,j-1,:))
!! At j=1, PDYM(:,1,:) are replaced by the values of PDYM,
!! which are the right values in the y-cyclic case.
!!
!!
!! EXTERNAL
!! --------
!! NONE
!!
!! IMPLICIT ARGUMENTS
!! ------------------
!! Module MODD_PARAMETERS: declaration of parameter variables
!! JPHEXT: define the number of marginal points out of the
!! physical domain along the horizontal directions.
!!
!! REFERENCE
!! ---------
!! Book2 of documentation of Meso-NH (SHUMAN operators)
!! Technical specifications Report of The Meso-NH (chapters 3)
!!
!!
!! AUTHOR
!! ------
!! V. Ducrocq * Meteo France *
!!
!! MODIFICATIONS
!! -------------
!! Original 05/07/94
!! Modification to include the periodic case 13/10/94 J.Stein
!! optimisation 20/08/00 J. Escobar
!! correction of in halo/pseudo-cyclic calculation for JPHEXT<> 1
!! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1
!-------------------------------------------------------------------------------
!
!* 0. DECLARATIONS
! ------------
!
USE MODD_PARAMETERS, ONLY: JPHEXT
USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
!
IMPLICIT NONE
!
!* 0.1 Declarations of argument and result
! ------------------------------------
!
TYPE(DIMPHYEX_t), INTENT(IN) :: D
REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PA ! variable at mass
! localization
REAL, DIMENSION(D%NIT,D%NJT), INTENT(OUT) :: PDYM ! result at flux
! side
!
!* 0.2 Declarations of local variables
! -------------------------------
!
INTEGER :: JJ ! Loop index in y direction
INTEGER :: IJU ! Size of the array in the y direction
!
!
INTEGER :: IIU,IKU
INTEGER :: JIJK,JIJKOR,JIJKEND
!
!-------------------------------------------------------------------------------
!
!* 1. DEFINITION OF DYM
! ------------------
!
IIU=SIZE(PA,1)
IJU=SIZE(PA,2)
!
JIJKOR = 1 + IIU
JIJKEND = IIU*IJU
!
!CDIR NODEP
!OCL NOVREC
DO JIJK=JIJKOR , JIJKEND
PDYM(JIJK,1) = PA(JIJK,1) - PA(JIJK-IIU,1)
END DO
!
DO JJ=1,JPHEXT
PDYM(:,JJ) = PDYM(:,IJU-2*JPHEXT+JJ) ! for reprod JPHEXT <> 1
END DO
!
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE DYM2D_PHY
! ###############################
SUBROUTINE DYF_PHY(D,PA,PDYF)
! ###############################
......
src/PHYEX/turb/mode_turb_ver_dyn_flux.f90
+ 228
− 414
View file @ 9dc96e65
This diff is collapsed.
src/Rules.LXnvhpc2202.mk
+ 2
− 2
View file @ 9dc96e65
......@@ -340,6 +340,6 @@ PHYEX_OPTDEFAULT = --addMPPDB_CHECKS --addStack MESONH --stopScopes toto --mathF
PHYEX_LIST = $(notdir $(shell find PHYEX/micro PHYEX/turb -follow -type f -name "*.f*" -not -name "minpack.f90" | sed -e 's/\(.*\)\(\.\).*/\1.D/g' ))
$(PHYEX_LIST) : PYFT = pyft_tool.py $(PHYEX_OPTDEFAULT)
PHYEX_SHUMAN = mode_tke_eps_sources.D mode_turb_ver_thermo_flux.D mode_turb_ver_thermo_corr.D
$(PHYEX_SHUMAN) : PYFT = pyft_tool.py --shumanFUNCtoCALL $(PHYEX_OPTDEFAULT)
PHYEX_SHUMAN = mode_tke_eps_sources.D mode_turb_ver_thermo_flux.D mode_turb_ver_thermo_corr.D mode_turb_ver_dyn_flux.D
$(PHYEX_SHUMAN) : PYFT = pyft_tool.py --shumanFUNCtoCALL --expandAllArraysConcurrent $(PHYEX_OPTDEFAULT)
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme