$$$$取第一个曲面的点生成测量点
$$$--------如果是三十个点把下面的4都换成30
DECL/REAL,XA[4],YA[4],ZA[4],IA[4],JA[4],KA[4]
$$_______________________________________

DECL/INTGR,K1
DECL/CHAR,20,PONAME

$$$--------如果是三十个点把下面的4换成30
DO/K1,1,4,1
$$$----------------------------


XA[K1]=OBTAIN/F(GSU_1)[K1],4
YA[K1]=OBTAIN/F(GSU_1)[K1],5
ZA[K1]=OBTAIN/F(GSU_1)[K1],6

JA[K1]=OBTAIN/F(GSU_1)[K1],8
KA[K1]=OBTAIN/F(GSU_1)[K1],9
IA[K1]=OBTAIN/F(GSU_1)[K1],10
PONAME=ASSIGN/CONCAT('PI_',STR(K1))
F(@PONAME)=FEAT/POINT,CART,XA[K1],YA[K1],ZA[K1],JA[K1],KA[K1],IA[K1]
XA[K1]=OBTAIN/FA(GSU_1)[K1],4
YA[K1]=OBTAIN/FA(GSU_1)[K1],5
ZA[K1]=OBTAIN/FA(GSU_1)[K1],6

JA[K1]=OBTAIN/FA(GSU_1)[K1],8
KA[K1]=OBTAIN/FA(GSU_1)[K1],9
IA[K1]=OBTAIN/FA(GSU_1)[K1],10
PONAME=ASSIGN/CONCAT('PI_',STR(K1))
FA(@PONAME)=FEAT/POINT,CART,XA[K1],YA[K1],ZA[K1],JA[K1],KA[K1],IA[K1]
ENDDO



$$$$取第二个曲面的点生成测量点
$$$--------如果是三十个点把下面的4都换成30
DECL/REAL,XA[4],YA[4],ZA[4],IA[4],JA[4],KA[4]
$$_______________________________________
DECL/INTGR,K1
DECL/CHAR,20,PONAME
$$$--------如果是三十个点把下面的4换成30
DO/K1,1,4,1
$$$----------------------------


XA[K1]=OBTAIN/F(GSU_2)[K1],4
YA[K1]=OBTAIN/F(GSU_2)[K1],5
ZA[K1]=OBTAIN/F(GSU_2)[K1],6

JA[K1]=OBTAIN/F(GSU_2)[K1],8
KA[K1]=OBTAIN/F(GSU_2)[K1],9
IA[K1]=OBTAIN/F(GSU_2)[K1],10
PONAME=ASSIGN/CONCAT('PII_',STR(K1))

F(@PONAME)=FEAT/POINT,CART,XA[K1],YA[K1],ZA[K1],JA[K1],KA[K1],IA[K1]


XA[K1]=OBTAIN/FA(GSU_2)[K1],4
YA[K1]=OBTAIN/FA(GSU_2)[K1],5
ZA[K1]=OBTAIN/FA(GSU_2)[K1],6

JA[K1]=OBTAIN/FA(GSU_2)[K1],8
KA[K1]=OBTAIN/FA(GSU_2)[K1],9
IA[K1]=OBTAIN/FA(GSU_2)[K1],10
PONAME=ASSIGN/CONCAT('PII_',STR(K1))
FA(@PONAME)=FEAT/POINT,CART,XA[K1],YA[K1],ZA[K1],JA[K1],KA[K1],IA[K1]

ENDDO
$$$$----------求两个对应点的距离

V(TEXT)=VFORM/ALL
DISPLY/TERM,V(TEXT)
DECL/INTGR,i
DECL/CHAR,20,POINAME1,POINAME2

$$$--------如果是三十个点把下面的4换成30
DO/i,1,4,1
POINAME1=ASSIGN/CONCAT('PI_',STR(i))
POINAME2=ASSIGN/CONCAT('PII_',STR(i))
T(DISTB_1)=TOL/DISTB,NOMINL,13.8771,0.0000,0.0000,PT2PT
OUTPUT/FA(@POINAME1),FA(@POINAME2),TA(DISTB_1)
ENDDO