& streamice_kx_b_init, streamice_ky_b_init,
_RL streamice_kx_b_init, streamice_ky_b_init
& sin(2*streamice_ky_b_init*PI*y/leny)))
streamice_ky_b_init = 1.