CALL DIAGNOSTICS_FILL(trac2d_diag,'SI2Dtrac',
_RL trac2d_diag
trac2d_diag(i,j,bi,bj)=trac2d(i,j,bi,bj)v