Problem with negative precipitation

[maignan]

It happens that some forcing files have some slightly negative precipitation. In this case the program stops in hydrol_soil_infilt with the message:

Error in the calculation of infilt tot -4.176114610371373E-008
 k, ji, jst, mc   1.17784270216726        2.16374183793683             2394
           1  0.277241865204730

FATAL ERROR FROM ROUTINE hydrol_soil_infilt
 --> We will STOP after hydrol_soil_infilt.
 -->
 -->

Fatal error from IOIPSL. STOP in ipslerr with code

We have to decide if we correct the forcing files or if we modify the code.
Such a modification could be made in intersurf_initialize_2d/intersurf_main_2d:

!@ PROPOSED MODIFICATION
    WHERE(zprecip_rain(:) .LT. 0.)
	zprecip_rain(:)=0.		
    ENDWHERE				
!@ END PROPOSED MODIFICATION
 !
    IF (check_INPUTS) THEN
...

[jpolcher]

For me this should clearly be solved in the forcing.
How can a forcing file provide negative rainfall ... this makes no physical sense.
We can put some test into the driver, but that is the most I would concede.

Jan