C语言翻译：BINDXAYA 子例程

void BINDXAYA(double Z00, double* XD, double* YD, double* ZD, int IIIID[][NK], int* IRRD, int NOMAD[][NK],
               int* NE34D, double XA, double YA, int* NED, int L0D, int NK) {
    double DLMIN = 1.0e50;
    double IX[4], XX[4], YY[4], ZZ[4], IBD[4], EN1[3], EN2[3];
    for (int IU = 0; IU < *NED; IU++) {
        if (IRRD[IU] == IU + 1) {
            for (int JA = 0; JA < NK; JA++) {
                IX[JA] = IIIID[IU][JA];
                XX[JA] = XD[IX[JA]];
                YY[JA] = YD[IX[JA]];
            }
            for (int JA = 0; JA < NE34D[IU]; JA++) {
                IBD[JA] = 0;
                if (NOMAD[IU][JA] == 0) {
                    IBD[JA] = JA + 1;
                }
            }
            for (int JA = 0; JA < NE34D[IU]; JA++) {
                if (IBD[JA] == JA + 1) {
                    int ID1 = JA + 1;
                    int ID2 = JA + 2;
                    if (NE34D[IU] == 3 && ID2 == 4) {
                        ID2 = 1;
                    }
                    if (NE34D[IU] == 4 && ID2 == 5) {
                        ID2 = 1;
                    }
                    double X2 = 0.5 * (XX[ID1 - 1] + XX[ID2 - 1]);
                    double Y2 = 0.5 * (YY[ID1 - 1] + YY[ID2 - 1]);
                    double DX = XA - X2;
                    double DY = YA - Y2;
                    double DL2 = DX * DX + DY * DY;
                    if (DLMIN > DL2) {
                        DLMIN = DL2;
                        IIMIN = IU + 1;
                    }
                }
            }
        }
    }
    for (int JA = 0; JA < NK; JA++) {
        XX[JA] = XD[IIIID[IIMIN - 1][JA] - 1];
        YY[JA] = YD[IIIID[IIMIN - 1][JA] - 1];
        ZZ[JA] = ZD[IIIID[IIMIN - 1][JA] - 1];
    }
    NOLXYZ(EN1, XX, YY, ZZ);
    for (int IH = 0; IH < 3; IH++) {
        EN2[IH] = EN1[IH];
    }
    double ENM = -0.707;
    if (EN1[2] < ENM) {
        if (fabs(EN1[0]) >= fabs(EN1[1])) {
            if (fabs(EN1[0]) < 1.0e-10) {
                EN1[0] = 1.0e-10;
            }
            double AB = EN1[1] / EN1[0];
            double EZ1 = sqrt((1.0 - ENM * ENM) / (1.0 + AB * AB));
            EN2[0] = EZ1;
            if (EN1[0] < 0.0) {
                EN2[0] = -EZ1;
            }
            EN2[1] = AB * EN2[0];
            EN2[2] = ENM;
        }
        if (fabs(EN1[0]) < fabs(EN1[1])) {
            if (fabs(EN2[0]) < 1.0e-10) {
                EN2[0] = 1.0e-10;
            }
            double AB = EN1[0] / EN2[0];
            double EZ1 = sqrt((1.0 - ENM * ENM) / (1.0 + AB * AB));
            EN2[1] = EZ1;
            if (EN1[1] < 0.0) {
                EN2[1] = -EZ1;
            }
            EN2[0] = AB * EN2[1];
            EN2[2] = ENM;
        }
    }
    double ANXY = EN2[0] * (XA - XX[0]) + EN2[1] * (YA - YY[0]);
    Z00 = ZZ[0] - ANXY / (EN2[2] + 1.0e-30);
    return;
}

代码解释：

函数声明：

void BINDXAYA(double Z00, double* XD, double* YD, double* ZD, int IIIID[][NK], int* IRRD, int NOMAD[][NK],
               int* NE34D, double XA, double YA, int* NED, int L0D, int NK)

该函数接受13个参数，其中包括一个输出参数Z00，代表计算得到的点到平面的距离。

局部变量定义：

    double DLMIN = 1.0e50;
    double IX[4], XX[4], YY[4], ZZ[4], IBD[4], EN1[3], EN2[3];

定义了几个局部变量，用于存储中间计算结果。

循环遍历每个点：

    for (int IU = 0; IU < *NED; IU++) {
        if (IRRD[IU] == IU + 1) {
            ...
        }
    }

循环遍历所有点，并判断当前点是否在平面内。

计算点到平面的距离：

    double ANXY = EN2[0] * (XA - XX[0]) + EN2[1] * (YA - YY[0]);
    Z00 = ZZ[0] - ANXY / (EN2[2] + 1.0e-30);

根据计算得到的平面法向量和点坐标，计算出点到平面的距离。

注意：

该代码中还包含NOLXYZ函数，该函数用于计算三维平面的法向量。需要根据具体实现情况进行补充。
代码中使用了1.0e-30来避免除零错误，可以根据实际情况进行调整。
代码中的数组下标从0开始，与FORTRAN中的下标从1开始不同，需要在翻译时进行调整。

总结：

该C语言代码实现了FORTRAN子例程BINDXAYA的功能，用于计算点到三维平面的距离。代码中包含详细的注释，可以帮助理解代码的逻辑和实现细节。在实际应用中，需要根据具体情况进行修改和补充。