%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];

FLOW_A = [0 1 0;
          1 -2.9 1;
          0 1 0];        
FLOW_B = [0 0 0;
          0 1 0;
          0 0 0];
FLOW_I =1;
         
HEATING_A = 10*[0.5 1 0.5;
          1 -6 1;
          0.5 1 0.5];        
HEATING_B = [0 0 0;
          0 0.1 0;
          0 0 0];
HEATING_I =0;    
         
HEAT2D_A = [0.5 1 0.5;
          1 -6 1;
          0.5 1 0.5];        
HEAT2D_B = [0 0 0;
          0 0 0;
          0 0 0];
HEAT2D_I =0;         
         
HEAT_A = [0 1 0;
          0 -2 0;
          0 1 0];        
HEAT_B = [0 0 0;
          0 0 0;
          0 0 0];
HEAT_I =0;
         
         
WAVE_A = [0 1 0;
          0 0 0;
          0 0 0];        
WAVE_B = [0 0 0;
          0 0 0;
          0 0 0];
WAVE_I =0;