diff --git a/src/MNH/ini_lima_warm.f90 b/src/MNH/ini_lima_warm.f90
index 5ef1977aad916d21d9d50b68048b60f85c1051c4..38664ac88dae17bb5eaf2c2c87d83d7325326cf6 100644
--- a/src/MNH/ini_lima_warm.f90
+++ b/src/MNH/ini_lima_warm.f90
@@ -166,7 +166,7 @@ XLBEXC = 1.0/XBC
 IF (NMOM_R.EQ.1) THEN
    XCCR   = 8.E6
    XCXR   = -1.
-   XLBEXR = 1.0/XBR
+   XLBEXR = 1.0/(XCXR-XBR)
    XLBR   = ( XAR*XCCR*MOMG(XALPHAR,XNUR,XBR) )**(-XLBEXR)
 ELSE
    XLBR   = XAR*ZGAMR(2)
diff --git a/src/MNH/lima_rain_evaporation.f90 b/src/MNH/lima_rain_evaporation.f90
index f4c0df7d033880a071937fc0c696bacc2899879a..c7211f2fcd62960a43373cb63f19c38314ea6d0f 100644
--- a/src/MNH/lima_rain_evaporation.f90
+++ b/src/MNH/lima_rain_evaporation.f90
@@ -150,8 +150,7 @@ ELSE
 !
 ! Compute the function G(T)
 !
-      ZZW2(:) = 1. / ( XRHOLW*((((PLV(:)/PT(:))**2)/(XTHCO*XRV)) +          & ! G
-           (XRV*PT(:))/(XDIVA*PEVSAT(:))))
+      ZZW2(:) = 1. / ( XRHOLW*((((PLV(:)/PT(:))**2)/(XTHCO*XRV)) + (XRV*PT(:))/(XDIVA*PEVSAT(:)))) !G
 !
 ! Compute the evaporation tendency
 !