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).

Bắt đầu trong cpltur.f90 và tái hiện segfault với CS_TURB_NONE và iturt(iscalt)=10, bao gồm cả thiết lập hệ số khuếch tán biến thiên theo không gian và thời gian. Theo dõi nhánh GGDH cho biến nhiệt và so sánh giả định tensor đối xứng của nhánh này với trường hợp vô hướng dòng tầng. Được xem là hoàn thành khi các thiết lập không tương thích không còn gây ra segfault và либо được bỏ qua an toàn hoặc được báo cáo kèm cảnh báo.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.