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.