function [i,j]=ij(ii,nx);

i=rem(ii-1,nx)+1;
j=floor((ii-1)/nx)+1;
