Create surfacewater.m