The diag variable in the source code

[Hejitian Pan]

Hello all,

In the source code file there is a vadiable int* diag, it is often modified to 0 or 1 like in line 490: diag[ index(c_offset_T,j) ] = 0 ,line 453: diag[ index(c_offset_Y + k,j) ] = 1. The line number may not be accurate because i have modified some of the lines. What does this vadiable diag do? From the name i guess it's the diagonal of Jacobian or some matrix. When should i set it to 0 and when to 1? And I would be more than grateful if somebody can tell me how does it take effect in the numerical solving process.

Regards,

Pan

[Hejitian Pan]

i miss out the file name. this variable appears in StFlow.h and .cpp file