tab_coeffs1(1,1)=theta(prof_i,prof_j,k,bi,bj) !SO

               tab_coeffs1(1,2)=theta(prof_i+1,prof_j,k,bi,bj) !SE

               tab_coeffs1(2,1)=theta(prof_i,prof_j+1,k,bi,bj) !NO

               tab_coeffs1(2,2)=theta(prof_i+1,prof_j+1,k,bi,bj) !NZ              

               tab_coeffs1(1,1)=salt(prof_i,prof_j,k,bi,bj) !SO

               tab_coeffs1(1,2)=salt(prof_i+1,prof_j,k,bi,bj) !SE

               tab_coeffs1(2,1)=salt(prof_i,prof_j+1,k,bi,bj) !NO

               tab_coeffs1(2,2)=salt(prof_i+1,prof_j+1,k,bi,bj) !NZ         

               tab_coeffs1(1,1)=uVel(prof_i,prof_j,k,bi,bj) !SO

               tab_coeffs1(1,2)=uVel(prof_i+1,prof_j,k,bi,bj) !SE

               tab_coeffs1(2,1)=uVel(prof_i,prof_j+1,k,bi,bj) !NO

               tab_coeffs1(2,2)=uVel(prof_i+1,prof_j+1,k,bi,bj) !NZ

               tab_coeffs1(1,1)=vVel(prof_i,prof_j,k,bi,bj) !SO

               tab_coeffs1(1,2)=vVel(prof_i+1,prof_j,k,bi,bj) !SE

               tab_coeffs1(2,1)=vVel(prof_i,prof_j+1,k,bi,bj) !NO

               tab_coeffs1(2,2)=vVel(prof_i+1,prof_j+1,k,bi,bj) !NZ

               tab_coeffs1(1,1)=pTracer(prof_i,prof_j,k,bi,bj,1) !SO

               tab_coeffs1(1,2)=pTracer(prof_i+1,prof_j,k,bi,bj,1) !SE

               tab_coeffs1(2,1)=pTracer(prof_i,prof_j+1,k,bi,bj,1) !NO

               tab_coeffs1(2,2)=pTracer(prof_i+1,prof_j+1,k,bi,bj,1) !NZ

               tab_coeffs1(1,1)=0 !SO

               tab_coeffs1(1,2)=0 !SE

               tab_coeffs1(2,1)=0 !NO

               tab_coeffs1(2,2)=0 !NZ

               tab_coeffs1(1,1)=etan(prof_i,prof_j,bi,bj) !SO

               tab_coeffs1(1,2)=etan(prof_i+1,prof_j,bi,bj) !SE

               tab_coeffs1(2,1)=etan(prof_i,prof_j+1,bi,bj) !NO

               tab_coeffs1(2,2)=etan(prof_i+1,prof_j+1,bi,bj) !NZ

               tab_coeffs1(1,1)=0.

               tab_coeffs1(2,1)=0.

               tab_coeffs1(1,2)=0.

               tab_coeffs1(2,2)=0.

         tab_coeffs1(1,1)=tab_coeffs1(1,1)*ponderations(1,1)/pondsSUM

         tab_coeffs1(1,2)=tab_coeffs1(1,2)*ponderations(1,2)/pondsSUM

         tab_coeffs1(2,1)=tab_coeffs1(2,1)*ponderations(2,1)/pondsSUM

         tab_coeffs1(2,2)=tab_coeffs1(2,2)*ponderations(2,2)/pondsSUM

         traj_cur(k)=tab_coeffs1(1,1)+tab_coeffs1(2,1)+

     &   tab_coeffs1(1,2)+tab_coeffs1(2,2)   

      _RL tab_coeffs1(2,2),tab_coeffs2(2,2),tab_coeffs3(2,2)