Unable to output viscosity under GD model in vtk_out.f

I added the code to calculate the eta_k_star, eta_star, eta0 belonging to the GD model in vtk_out.f file, but it failed to output successfully, that is, the output is 0. Can you help me find out what the problem is? Thank you!
taghipour3d_gd.zip (27.8 MB)

You declare the array eta_k_star and allocate it but you never assign a value to it. The variable eta_k_star in calc_mu_s is a local scalar variable and its value is not saved to an array. You need to store the value into an array (use an different name).

Thank you, Jeff. I’ll try.