Philippe 28/06/2021: small optimisation in the call to ice4_rainfr_vert (create arrays with zeros instead of multiplications by zero)