CALL DIAGNOSTICS_FILL(streamice_taubx_diag,'SI_taubx',
_RL streamice_taubx_diag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
streamice_taubx_diag(i,j,bi,bj) =