diff --git a/src/LIB/SURCOUCHE/src/mode_toolsz_ll.f90 b/src/LIB/SURCOUCHE/src/mode_toolsz_ll.f90 index 67ec731d0d40e551b5b8390ef9b79236ac55420f..8794b9e36973903f350540e66c5b0062c57f1d8f 100644 --- a/src/LIB/SURCOUCHE/src/mode_toolsz_ll.f90 +++ b/src/LIB/SURCOUCHE/src/mode_toolsz_ll.f90 @@ -339,11 +339,10 @@ END IF IF ( ( (1+TPROC(IP)%NXEND-TPROC(IP)%NXOR) == 0 ) & .OR. ( (1+TPROC(IP)%NYEND-TPROC(IP)%NYOR) == 0 ) ) THEN - PRINT*, "/!\ SPLITZ: some proc have 0 size local domaine , to much processors used for domaine size ", & - " IP=" ,IP , & - " DIMX=",1+TPROC(IP)%NXEND-TPROC(IP)%NXOR, & - " DIMY=",1+TPROC(IP)%NYEND-TPROC(IP)%NYOR - CALL ABORT() + write( yval1, '( I10 )' ) 1+TPROC(IP)%NXEND-TPROC(IP)%NXOR + write( yval2, '( I10 )' ) 1+TPROC(IP)%NYEND-TPROC(IP)%NYOR + call Print_msg( NVERB_FATAL, 'GEN', 'SPLITZ', 'zero-size local-domain dimensions: DIMX='// & + trim(yval1)//', DIMY='//trim(yval2)//'. Too many processes for this domain size.') END IF ! !* 3. shift from physical to extended domain diff --git a/src/MNH/compute_updraft_hrio.f90 b/src/MNH/compute_updraft_hrio.f90 index 6ec162b20d6bd4c14587325a2e5e4c266ca4dc19..8086333465e5cd97bfc66dacda454c363555e725 100644 --- a/src/MNH/compute_updraft_hrio.f90 +++ b/src/MNH/compute_updraft_hrio.f90 @@ -1,7 +1,8 @@ -!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier +!MNH_LIC Copyright 1994-2019 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 version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt !MNH_LIC for details. version 1. +!----------------------------------------------------------------- ! ######spl MODULE MODI_COMPUTE_UPDRAFT_HRIO ! ########################### @@ -135,6 +136,7 @@ END MODULE MODI_COMPUTE_UPDRAFT_HRIO! ######spl !! S. Riette Jan 2012: support for both order of vertical levels !! V.Masson, C.Lac : 02/2011 : SV_UP initialized by a non-zero value !! Q.Rodier 01/2019 : support RM17 mixing length +! P. Wautelet 12/04/2019: replace ABORT and STOP calls by Print_msg !! -------------------------------------------------------------------------- ! !* 0. DECLARATIONS @@ -150,6 +152,8 @@ USE MODD_GRID_n, ONLY : XDXHAT, XDYHAT USE MODD_BLANK USE MODD_TURB_n, ONLY :CTURBLEN +use mode_msg + !USE MODI_COMPUTE_ENTR_DETR USE MODI_TH_R_FROM_THL_RT_1D USE MODI_SHUMAN_MF @@ -394,7 +398,7 @@ PSV_DO(:,:,:)=0. PTHL_UP(:,KKB)= ZTHLM_F(:,KKB)+MAX(0.,MIN(ZTMAX,(PSFTH(:)/SQRT(ZTKEM_F(:,KKB)))*XALP_PERT)) PRT_UP(:,KKB) = ZRTM_F(:,KKB)+MAX(0.,MIN(ZRMAX,(PSFRV(:)/SQRT(ZTKEM_F(:,KKB)))*XALP_PERT)) !------------------------ -print*,OENTR_DETR +! print*,OENTR_DETR !------------------------ IF (OENTR_DETR) THEN ZTHM_F (:,:) = MZM_MF(KKA,KKU,KKL,PTHM (:,:)) @@ -767,7 +771,8 @@ DO JK=KKB,KKE-KKL,KKL ! on cherche à savoir s'il y a des vitesses verticales non définies ! je n'utilise que ZW_UP2 pour pouvoir avoir une valeur si ZW_UP ! n'est pas défini -IF (maxval(ZW_UP2(:,JK+KKL)) .NE. maxval(ZW_UP2(:,JK+KKL))) STOP 'probleme ici' +IF (maxval(ZW_UP2(:,JK+KKL)) .NE. maxval(ZW_UP2(:,JK+KKL))) & + call Print_msg( NVERB_FATAL, 'GEN', 'COMPUTE_UPDRAFT_HRIO', 'maxval(ZW_UP2(:,JK+KKL)) /= maxval(ZW_UP2(:,JK+KKL))' ) ! si on est dans la zone grise la définition du flux de masse change ! donc celle de alpha aussi WHERE(GTEST) @@ -858,6 +863,6 @@ ENDDO ! boucle JK GWORK1(:)= (GTESTLCL(:) .AND. (PDEPTH(:) > ZDEPTH_MAX1) ) GWORK2(:,:) = SPREAD( GWORK1(:), DIM=2, NCOPIES=MAX(KKU,KKA) ) ZCOEF(:,:) = SPREAD( (1.-(PDEPTH(:)-ZDEPTH_MAX1)/(ZDEPTH_MAX2-ZDEPTH_MAX1)), DIM=2, NCOPIES=SIZE(ZCOEF,2)) -print*,"je sors de compute_updraft" +! print*,"je sors de compute_updraft" END SUBROUTINE COMPUTE_UPDRAFT_HRIO