Hello, look this forum there is a lot information about water sat modelling
As hint use jfunction