& bdot_streamice1
_RL bdot_streamice1
& bdot_streamice1,
& bdot_streamice1(i,j,bi,bj)
& + (1-fac) * bdot_streamice1(i,j,bi,bj)