diff --git a/src/MNH/lima_warm.f90 b/src/MNH/lima_warm.f90 index 58098af61d7995fab66b8062db5cac79e2dd5174..653ae66e31cf3005735db69228fa9d279cef8fec 100644 --- a/src/MNH/lima_warm.f90 +++ b/src/MNH/lima_warm.f90 @@ -120,6 +120,7 @@ END MODULE MODI_LIMA_WARM !! ------------- !! Original ??/??/13 !! C. Barthe * LACy * jan. 2014 add budgets +!! J. Escobar : for real*4 , use XMNH_HUGE !! !------------------------------------------------------------------------------- ! @@ -283,7 +284,7 @@ END IF ! ---------------------------------------- ! ! -ZWLBDC3(:,:,:) = 1.E45 +ZWLBDC3(:,:,:) = XMNH_HUGE ZWLBDC(:,:,:) = 1.E15 ! WHERE (PRCT(:,:,:)>XRTMIN(2) .AND. PCCT(:,:,:)>XCTMIN(2)) diff --git a/src/MNH/lima_warm_nucl.f90 b/src/MNH/lima_warm_nucl.f90 index af3de3889d5e2f0e00ab6422f305b2ca8ab29301..ea473aa913b81b35f13c2d75b75ea03afa1d1e43 100644 --- a/src/MNH/lima_warm_nucl.f90 +++ b/src/MNH/lima_warm_nucl.f90 @@ -95,6 +95,7 @@ END MODULE MODI_LIMA_WARM_NUCL !! MODIFICATIONS !! ------------- !! Original ??/??/13 +!! J. Escobar : 10/2017 , for real*4 use XMNH_EPSILON !! !------------------------------------------------------------------------------- ! @@ -739,7 +740,7 @@ INTEGER :: PIVEC1 ALLOCATE(PFUNCSMAX(NPTS)) ! PFUNCSMAX(:) = 0. -PZVEC1 = MAX( 1.00001,MIN( FLOAT(NHYP)-0.00001, & +PZVEC1 = MAX( ( 1.0 + 10.0 * XMNH_EPSILON ) ,MIN( FLOAT(NHYP)*( 1.0 - 10.0 * XMNH_EPSILON ) , & XHYPINTP1*LOG(PPZSMAX)+XHYPINTP2 ) ) PIVEC1 = INT( PZVEC1 ) PZVEC1 = PZVEC1 - FLOAT( PIVEC1 )