Cannot output rotational granular temperature

I added the code to calculate the rotational granular temperature in vtk_out.f file, but it failed to output successfully. Can you help me find out what the problem is? Thank you!
rotating_drum_5.26_2022-05-30T163205.041000.zip (27.7 KB)

Please attach your code (all needed files) and explain what you expect to get and what you mean exactly by “it failed to output successfully” (no output in the vtk file, incorrect output etc.). Are you running in serial or parallel?

Thank you for your reply!@jeff.dietiker

Hello,@jeff.dietiker!
I added the code to calculate the rotational granular temperature in vtk_out.f file, I tried to output the calculated data of the rotational granular temperature by vtk_u_g(2)= True, but failed.Please focus on the vtk_out.f file and help me see where the problem is. Thank you!
rotating_drum1.571rad_2022-06-01T091354.373000.zip (6.8 MB)

Your conditional is based on AU_S, AV_S and AW_S being allocated and they are not allocated. You need to allocate AU_S, AV_S and AW_S before the conditional and deallocate after writing the data.

Thank you for your help!@jeff.dietiker
I’ll give it a try!

Success,thank you very much!@jeff.dietiker