%diff

ANISO_A  = [0 0 0;
            0 1 0;
            0 0 0];

ANISO_Dd = [0.5 1.0 0.5;
            1.0 0.0 1.0;
            0.5 1.0 0.5];

ANISO_d  = [1 5  -2 0  -0.3 0  0 1  0.3 0  2 0  122];

% coeff should be scaled !

DIFFUSC_A = [0.10 0.15 0.10;
             0.15 0.00 0.15;
             0.10 0.15 0.10]; 

DIFFUSC_B = [0.10 0.15 0.10;
             0.15 0.00 0.15;
             0.10 0.15 0.10];