Bugfix lambda computation for rain in case NMOM_R=1

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)

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
 !