UDF граничные условия

[djeak11 0]

Работаю с 3d моделью. Написал код для граничных условий(https://disk.yandex.ru/i/i63pturKk-pJZA), но не запускается, пишет ошибки: 

error C2198: pow: ??? ???? ?????? ??? ????

..\..\src\udf_CH4.c(45): warning C4020: sqrt: ??? ????? ??????? ??????

..\..\src\udf_CH4.c(45): error C2143: ???????? ????: ???????? ")" ??? ";"

..\..\src\udf_CH4.c(52): error C2143:

 

Вот код: 

#include "udf.h"
DEFINE_PROFILE(inlet_mf,th,i)
{

face_t f;
cell_t c;
Domain *d = Get_Domain(1);
Thread *t = Lookup_Thread(d,7);

real pw;
real pin;
real pcr;
real Tin;
real gamma;
real g;
real gmax;
real R;
real rho1;
real rho2;
real area;
real A[ND_ND];

R = 368.9; // J/kg k
gamma = 1.29;

pin = 1209000;
pcr = pin * pow( (2/(gamma+1)), (gamma/(gamma-1)) );
Tin = 665;

begin_f_loop(f,th)
    {

    pw = F_P(f,th); //pressure on one face

    F_AREA(A,f,th);
    area = NV_MAG(A); // Find area of one face

    if(pw >= pin)
    {
        F_PROFILE(f,th,i) = 0.0000001;
    }

    else if(pw > pcr && pw < pin)
    {
        g = sqrt( ((2*gamma)/(gamma-1)) * ((pin*pin)/(R*Tin))) * sqrt( pow(pw/pin), (2/gamma) - pow( (pw/pin), ((gamma+1)/gamma) );
        rho1 = F_R(f,th);
        F_PROFILE(f,th,i) = rho1*g*area;
    }

    else if(pw <= pcr)
    {
        gmax = sqrt( ((gamma*pin*pin)/(R*Tin)) * pow( (pcr/pin), ((gamma+1)/gamma) );
    rho2 = F_R(f,th);
        F_PROFILE(f,th,i) = rho2*gmax*area;
    }

    }
end_f_loop(f,th)
Возможно, я со скобками перемудрил, или для 3d по другому написать нужно?