Skip to content
Snippets Groups Projects
Commit d4341889 authored by WAUTELET Philippe's avatar WAUTELET Philippe
Browse files

Naima / Philippe 14/01/2022: OpenACC: PRINT_MSG not allowed in kernels region

parent 3fb48933
No related branches found
No related tags found
No related merge requests found
!MNH_LIC Copyright 2006-2019 CNRS, Meteo-France and Universite Paul Sabatier
!MNH_LIC Copyright 2006-2022 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 version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
!MNH_LIC for details. version 1.
......@@ -358,20 +358,26 @@ REAL, DIMENSION(:,:,:), INTENT(INOUT) :: PFRAC_ICE ! Ice fraction (1 for ice on
!------------------------------------------------------------------------
! 1. Compute FRAC_ICE
!
!$acc kernels
SELECT CASE(HFRAC_ICE)
CASE ('T') !using Temperature
!$acc kernels
PFRAC_ICE(:,:,:) = ( XTMAXMIX - PT(:,:,:) ) / ( XTMAXMIX - XTMINMIX ) ! freezing interval
!$acc end kernels
CASE ('O') !using Temperature with old formulae
!$acc kernels
PFRAC_ICE(:,:,:) = ( XTT - PT(:,:,:) ) / 40. ! freezing interval
!$acc end kernels
CASE ('N') !No ice
!$acc kernels
PFRAC_ICE(:,:,:) = 0.
!$acc end kernels
CASE ('S') !Same as previous
!nothing to do
CASE DEFAULT
CALL PRINT_MSG(NVERB_FATAL,'GEN','COMPUTE_FRAC_ICE3D_DEVICE','invalid option')
END SELECT
!
!$acc kernels
PFRAC_ICE(:,:,:) = MAX( 0., MIN(1., PFRAC_ICE(:,:,:) ) )
!$acc end kernels
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment