Problème d'implémentation de conditions initiales

Bonjour,
Lors de la construction d'un intégrateur pour simuler l'évolution temporelle d'un glacier monodimensionnel, nous devons implémenter une certaine équation contenant un terme de source. Cependant, nous rencontrons un problème pour ce qui est d'implémenter ce terme de source. Voici la définition de celui-ci :
Pièce jointe 388091

Ainsi que son implémentation dans notre code :

Pièce jointe 388095

Dans notre code, seule l'implémentation de S0 est réalisée à cet endroit.
Pour implémenter la partie concernant S(x,t,h), nous avons utiliser l'option nonnegative que nous avons ajoutée à notre intégrateur ode45.

Pièce jointe 388096

Apparemment, nous avons un problème avec ce terme de source mais nous ne voyons pas d'autre façon d'implémenter celui-ci. Quelqu'un pourrait-il nous aider ?
Si besoin, je peux envoyer l'entièreté de notre énoncé et de notre code par mail.

Merci beaucoup,
Terry