{"id":381991,"date":"2024-09-09T13:42:33","date_gmt":"2024-09-09T13:42:33","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/reply\/381991\/"},"modified":"2024-09-09T13:42:33","modified_gmt":"2024-09-09T13:42:33","slug":"381991","status":"publish","type":"reply","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/reply\/381991\/","title":{"rendered":"Reply To: code won’t print all the temperature values"},"content":{"rendered":"

<p>Hi,<\/p><p> <\/p><p>Thank you for your reply. Yes I have this code and it runs well in serial. I’ve tried to paralellize the code and it prints me only values for two boundary faces ID’s. This is the code using macros for parallel and debug messages.<\/p><div>#include “udf.h”<\/div><div> <\/div><div>DEFINE_EXECUTE_AT_END(get_left_boundary_temperatures)<\/div><div>{<\/div><div>#if !RP_HOST <\/div><div>    Domain* d = Get_Domain(1);  \/\/ Get the domain<\/div><div>    Thread* t;<\/div><div>    face_t f;<\/div><div>#endif<\/div><div>    FILE* fp;<\/div><div>    int i, ID;<\/div><div>    int zone[3] = { 5, 6, 7 };<\/div><div> <\/div><div>    \/\/ Open file to write temperatures<\/div><div>    fp = fopen(“left_boundary_temperatures2.txt”, “w”);<\/div><div>    if (fp == NULL)<\/div><div>    {<\/div><div>        Message(“Error opening file for writing temperatures.\\n”);<\/div><div>        return;<\/div><div>    }<\/div><div>    Message(“File opened successfully.\\n”);<\/div><div> <\/div><div> <\/div><div>    for (i = 0; i < 3; i++) <\/div><div>    {  \/\/ Loop over relevant zones (left side)<\/div><div>        ID = zone[i];<\/div><div>#if !RP_HOST<\/div><div>        t = Lookup_Thread(d, ID);<\/div><div> <\/div><div>        \/\/ Debug message for thread ID<\/div><div>        Message(“Processing Thread ID: %d\\n”, ID);<\/div><div> <\/div><div>        \/\/ Check if the thread is a boundary thread<\/div><div>        if (BOUNDARY_FACE_THREAD_P(t))<\/div><div>        {<\/div><div>            Message(“Thread ID %d is a boundary thread.\\n”, ID);  \/\/ Debug message<\/div><div>            \/\/ Loop through faces in the thread<\/div><div>            begin_f_loop(f, t)<\/div><div>            {<\/div><div>                \/\/ Get the temperature at the face<\/div><div>                real temperature = F_T(f, t);<\/div><div> <\/div><div>                \/\/ Write the temperature to the file<\/div><div>                fprintf(fp, “Thread ID: %d, Face ID: %d, Temperature: %f\\n”, ID, f, temperature);<\/div><div> <\/div><div>                \/\/ Debug message for face ID and temperature<\/div><div>                Message(“Thread ID: %d, Face ID: %d, Temperature: %f\\n”, ID, f, temperature);<\/div><div>            }<\/div><div>            end_f_loop(f, t);<\/div><div>        }<\/div><div> <\/div><div>        else<\/div><div>        {<\/div><div>            Message(“Thread ID %d is NOT a boundary thread.\\n”, ID);  \/\/ Debug message<\/div><div>        }<\/div><div>#endif   <\/div><div>        PRF_GSYNC();<\/div><div> <\/div><div>    }<\/div><div> <\/div><div>   <\/div><div>    \/\/ Close the file<\/div><div>    fclose(fp);<\/div><div>    Message(“File closed successfully.\\n”);<\/div><div>}<\/div><\/p>\n","protected":false},"template":"","class_list":["post-381991","reply","type-reply","status-publish","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t