Hello,
I have tested it here on Unix and it works, so details on the kind of problem would be welcome. Do you have a crash or "simply" different results between scalar and parallel ?
I see nothing wrong in your Fortran file.
You have a minimum set of parameters and your case can probably be improved a lot. For example try :
TRAITEMENT DU SYSTEME LINEAIRE : 2
SOLVEUR : 1
and it should speed-up a lot. Are you also sure that you need so small time-steps ?
With best regards,
Jean-Michel Hervouet