function psi = initmom(s)
psi=0.*s;