YasserSelima
Subscriber
Hmmm ... if it doesn't work and you get the same error , try to include all of this inside the if statementnnnn  n        ID = 10000;n    }nn    if (ID != 10000)n    {n        Thread* thread_upstream = Lookup_Thread(domain_upstream, ID);nn        mdot_sum = 0.0;n        A_sum = 0.0;n        n        begin_f_loop(f_upstream, thread_upstream)n            if PRINCIPAL_FACE_P(f_upstream, thread_upstream)n            {n                F_AREA(area, f_upstream, thread_upstream);n                Amagnitude = NV_MAG(area);n                Flux = F_FLUX(f_upstream, thread_upstream);nn                mdot_sum += Flux;n                A_sum += Amagnitude;n            }n        end_f_loop(f_upstream, thread_upstream)nn        mdot_glob = PRF_GRSUM1(mdot_sum);n        A_glob = PRF_GRSUM1(A_sum);n        mflux_glob = mdot_glob / A_glob;n        printf(mflux_glob %f
, mflux_glob);n        printf(mdot_glob %f
, mdot_glob);n        printf(A_glob %f
, A_glob);n    }n    elsen    {n        mflux_glob = 0.0;n    }nnnn