Bug: Symmetric tensor diffusivity

The combination of `turb_model->iturb = CS_TURB_NONE;` and `iturt(iscalt) = 10` (i.e. GGDH) leads to a segfault. I would expect the code to ignore `iturt` for laminar flow so the user can simply change one parameter when switching on the turbulence model. IMO, we should at least throw a warning that the settings are incompatible, rather than segfault.

It seems the issue stems from `cpltur.f90` where a test for GGDH leads to an assumption of a symmetric tensor diffusivity for the thermal variable, when in reality the diffusivity is a scalar (my settings have diffusivity as variable in space and time, if it makes a difference).

Inizia in cpltur.f90 e riproduci il segfault con CS_TURB_NONE e iturt(iscalt)=10, includendo l'impostazione della diffusività variabile nello spazio e nel tempo. Segui il ramo GGDH per la variabile termica e confronta la sua assunzione di tensore simmetrico con il caso scalare laminare. Il lavoro è completato quando le impostazioni incompatibili non causano più un segfault e vengono ignorate in modo sicuro oppure segnalate con un avviso.

Scritto dal modello di indicizzazione a partire dal testo della issue.