diff --git a/src/MNH/budget.f90 b/src/MNH/budget.f90 index bed79a905fbfa201cad576828aa8f4ae23e4c5b2..e5bc6aa85657b2d33c8e3610083d5f122551d583 100644 --- a/src/MNH/budget.f90 +++ b/src/MNH/budget.f90 @@ -129,10 +129,10 @@ REAL :: XPRECISION ! for reproductibility checks ! Warning: requires an adaptation of the makefile in order to run two runs in ! parallel for comparison ! -XPRECISION = 1E-10 +!XPRECISION = 1E-10 IF (LCHECK) THEN print*,'BUDGET :',HBUVAR - CALL MPPDB_CHECK3D(PVARS,HBUVAR,XPRECISION) + CALL MPPDB_CHECK3D(PVARS,HBUVAR,PRECISION) END IF ! ! diff --git a/src/MNH/modeln.f90 b/src/MNH/modeln.f90 index 8b07c2ce1b3fb29c37b4abbbfd3b8cde02b976a6..c2e05af8a444e7358afa855ce842dc796ca996e6 100644 --- a/src/MNH/modeln.f90 +++ b/src/MNH/modeln.f90 @@ -893,7 +893,7 @@ ZTIME1=ZTIME2 IF( LLG .AND. IMI==1 ) CALL SETLB_LG ! IF (CCONF == "START" .OR. (CCONF == "RESTA" .AND. KTCOUNT /= 1 )) THEN -CALL MPPDB_CHECK3DM("before BOUNDARIES:XUT, XVT, XWT, XTHT, XTKET",PRECISION,& +CALL MPPDB_CHECK3DM("before BOUNDARIES:XU/V/W/TH/TKET",PRECISION,& & XUT, XVT, XWT, XTHT, XTKET) CALL BOUNDARIES ( & XTSTEP,CLBCX,CLBCY,NRR,NSV,KTCOUNT, & @@ -903,7 +903,7 @@ CALL BOUNDARIES ( & XLBYUS,XLBYVS,XLBYWS,XLBYTHS,XLBYTKES,XLBYRS,XLBYSVS, & XRHODJ, & XUT, XVT, XWT, XTHT, XTKET, XRT, XSVT, XSRCT ) -CALL MPPDB_CHECK3DM("after BOUNDARIES:XUT, XVT, XWT, XTHT, XTKET",PRECISION,& +CALL MPPDB_CHECK3DM("after BOUNDARIES:XU/V/W/TH/TKET",PRECISION,& & XUT, XVT, XWT, XTHT, XTKET) END IF ! diff --git a/src/MNH/tridiag_tke.f90 b/src/MNH/tridiag_tke.f90 index 389a1a7a7517810122f5de4a16816a0c88510fdd..a91fb9d1552b6b6bedf711e31f1b1d3c57c53a1e 100644 --- a/src/MNH/tridiag_tke.f90 +++ b/src/MNH/tridiag_tke.f90 @@ -147,6 +147,7 @@ END MODULE MODI_TRIDIAG_TKE !* 0. DECLARATIONS ! USE MODD_PARAMETERS +USE MODE_MPPDB ! IMPLICIT NONE ! @@ -185,6 +186,8 @@ INTEGER :: IKTB,IKTE ! start, end of k loops in physical domain ! !* 1. COMPUTE THE RIGHT HAND SIDE ! --------------------------- +CALL MPPDB_CHECK3DM("tridiag_tke::PVARM,PA,PRHODJ,PSOURCE,PDIAG",PRECISION,& + & PVARM,PA,PRHODJ,PSOURCE,PDIAG) ! !$acc kernels IKT=SIZE(PVARM,3) diff --git a/src/MNH/turb_ver_thermo_flux.f90 b/src/MNH/turb_ver_thermo_flux.f90 index 9c66901d6203d3ac9c0b1dd319989b8e6950d024..b5d05b69df3c397357034bebeb5d351b9c261686 100644 --- a/src/MNH/turb_ver_thermo_flux.f90 +++ b/src/MNH/turb_ver_thermo_flux.f90 @@ -370,6 +370,8 @@ USE MODI_TM06_H USE MODE_PRANDTL ! USE MODI_SECOND_MNH +USE MODE_MPPDB + ! IMPLICIT NONE ! @@ -723,6 +725,8 @@ ELSE !$acc end kernels END IF ! +CALL MPPDB_CHECK3DM("turb_ver_thermo_flux::ZF,PRHODJ,PREDTH1,PREDR1,PRED2TH3,PRED2THR3",PRECISION,& + & ZF,PRHODJ,PREDTH1,PREDR1,PRED2TH3,PRED2THR3 ) ! Compute the splitted conservative potential temperature at t+deltat CALL TRIDIAG_THERMO(KKA,KKU,KKL,PTHLM,ZF,ZDFDDTDZ,PTSTEP,PIMPL,PDZZ,& PRHODJ,PTHLP) @@ -797,6 +801,8 @@ ELSE !$acc end kernels END IF #endif +CALL MPPDB_CHECK3DM("turb_ver_thermo_flux::PTP,PETHETA,ZFLXZ,PBETA",PRECISION,& + & PTP,PETHETA,ZFLXZ,PBETA) ! ! Buoyancy flux at flux points !