Hello Thong,
Basically it is very simple to do a quick and dirty solution : in subroutine prosou.f, you need to have the array SMH%R set to the rain expressed in m/s (so a little number even in monsoon season). In prosou look for the line YASMH=.FALSE. and the previous line CALL OS('X=0 ',X=SMH) must be changed into CALL OS('X=C ',X=SMH,C=your rain in m/s).
However this will not work really well if you have already point sources (they will interfere) and if you have tracers (mass balance, etc.). The correct treatment of rain will be in version 6.2 (end of august, so very close), including the temperature of the rain !
With best regards and souvenir,
Jean-Michel Hervouet