翻译成C语言SUBROUTINE PREZD0Z0X0Y0XDYDZDIIIIDIIDEIRRDNOMADNE34D & IQDL0L0DNEDNK IMPLICIT REAL8A-HO-Z DIMENSION Z0L0X0L0Y0L0XDL0DYDL0DZDL0D & IIIIDNEDNKIIDEL0

void PREZD0(double *Z0, double *X0, double *Y0, double *XD, double *YD, double *ZD, int *IIIID, int *IIDE, int *IRRD, double *NOMAD, double *NE34D, int IQD, int L0, int L0D, int NED, int NK) { int IXYD; double XX[4], YY[4], ZZ[4]; for (int I0=0; I0<L0; I0++) { double XA=X0[I0]; double YA=Y0[I0]; XX[0]=XD[IIIID[IIDE[I0]NK]]; XX[1]=XD[IIIID[IIDE[I0]NK+1]]; XX[2]=XD[IIIID[IIDE[I0]NK+2]]; XX[3]=XD[IIIID[IIDE[I0]NK+3]]; YY[0]=YD[IIIID[IIDE[I0]NK]]; YY[1]=YD[IIIID[IIDE[I0]NK+1]]; YY[2]=YD[IIIID[IIDE[I0]NK+2]]; YY[3]=YD[IIIID[IIDE[I0]NK+3]]; PREXYA(&IXYD,XX,YY,&XA,&YA,NK); if (IXYD==1) { ZZ[0]=ZD[IIIID[IIDE[I0]NK]]; ZZ[1]=ZD[IIIID[IIDE[I0]NK+1]]; ZZ[2]=ZD[IIIID[IIDE[I0]NK+2]]; ZZ[3]=ZD[IIIID[IIDE[I0]NK+3]]; PREXYZ1(&Z0[I0],XX,YY,ZZ,&XA,&YA,NK); continue; } if (IXYD==2) { ZZ[0]=ZD[IIIID[IIDE[I0]NK]]; ZZ[1]=ZD[IIIID[IIDE[I0]NK+1]]; ZZ[2]=ZD[IIIID[IIDE[I0]NK+2]]; ZZ[3]=ZD[IIIID[IIDE[I0]NK+3]]; PREXYZ2(&Z0[I0],XX,YY,ZZ,&XA,&YA,NK); continue; } for (int II=0; II<NED; II++) { XX[0]=XD[IIIID[IINK]]; XX[1]=XD[IIIID[IINK+1]]; XX[2]=XD[IIIID[IINK+2]]; XX[3]=XD[IIIID[IINK+3]]; YY[0]=YD[IIIID[IINK]]; YY[1]=YD[IIIID[IINK+1]]; YY[2]=YD[IIIID[IINK+2]]; YY[3]=YD[IIIID[IINK+3]]; double EPS=1.0E-8; IXYD=0; double XMAX=XX[0]; double XMIN=XX[0]; double YMAX=YY[0]; double YMIN=YY[0]; for (int I=1; I<NK; I++) { if (XMAX<XX[I]) XMAX=XX[I]; if (YMAX<YY[I]) YMAX=YY[I]; if (XMIN>XX[I]) XMIN=XX[I]; if (YMIN>YY[I]) YMIN=YY[I]; } if (XA<XMAX+EPS && XA>XMIN-EPS) { if (YA<YMAX+EPS && YA>YMIN-EPS) { IXYD=1; } } if (IXYD==1) PREXYA(&IXYD,XX,YY,&XA,&YA,NK); if (IXYD==1) { ZZ[0]=ZD[IIIID[IINK]]; ZZ[1]=ZD[IIIID[IINK+1]]; ZZ[2]=ZD[IIIID[IINK+2]]; ZZ[3]=ZD[IIIID[IINK+3]]; PREXYZ1(&Z0[I0],XX,YY,ZZ,&XA,&YA,NK); IIDE[I0]=II; break; } if (IXYD==2) { ZZ[0]=ZD[IIIID[IINK]]; ZZ[1]=ZD[IIIID[IINK+1]]; ZZ[2]=ZD[IIIID[IINK+2]]; ZZ[3]=ZD[IIIID[IINK+3]]; PREXYZ2(&Z0[I0],XX,YY,ZZ,&XA,&YA,NK); IIDE[I0]=II; break; } } if (IXYD==0 && IQD==1) { BINDXAYA(&Z0[I0],XD,YD,ZD,IIIID,IRRD,NOMAD,NE34D,XA,YA,NED,L0D,NK); } }