SUBROUTINE FDEQCT (K,INDG)
COMMON /CIWZP1/ IWZP, JD, ID, IF, JB1, IB1
COMMON /CJ2/ JTM, JPL, JVO, JTAUO, IJB1, JRZ1, JRZ2, JRZ3
COMMON /CVP/ A(21,21,9), G1(21,21), G2(21,21), R0(21,21), ZMU(21,21)
COMMON /CCOEF/ CE, CW,CN,CS,C(5)
COMMON /CCHECK/ RSDU(9), RP(9), CC, DC, NMAX, NPRINT
COMMON /CGRID/ IMIN(21), IMAX(21), IN, JN, INM, JNM, IE
COMMON /CNUMBR/ NW, NF, NZK, NHS, NVT, NZML, NZMM, NZMN, NP
COMMON /CDI/ DIFMAX(9)
COMMON /CYLI/ JA, JA1, JB, JC, IAB, IC, RADA, RADA1, RADB, RADC, RADN,
1 DA, DA1,DB,DCN,DN,X2AXIS

COMMON /CX1X1/ X1(21), X2(21), DELX1(21), DELX2(21)
COMMON /CZW/ ZW, ZWHALF, IL,IH
COMMON /CRAD/ R(21)
WC1H (WK, DELF, ROP, DELRO, RAD)= - (DELF/(DX2*DX2))
4 /(RAD*(ROP*(RAD/3.+ 5.* WK * DX2/24.)
5 + DELRO*(5.*RAD/24.+3/ *WK*DX2/20.)))
7 -GGI /(ROP+DX2*(RAD/3.+5.*WK*DX2/24.))
WC2H(WK,ROP, DELRO, RAD)= -(RAD+WK*DX2)*
2 (ROP*(RAD/6.+WK*DX2/8.)
3 +DELRO+(RAD/8.+WK*DXP/10.))
4 /(RAD*(ROP*(RAD/3.+5.*WK*DX2/24.)
5 +DELRO*(5.*RAD/24.+3.*WK*DX2/20.)))
WC1V(DELF, ROP, DELRO,ETA2)= - DELF/(((ETA2*DX1)**2)*
1 (ROP/3.+3.*DELRO/24.))
7 = 3.0*GG1/(ROP*RAD*DX1)
WC2V(ROP, DELRO)= - (ROP/6.+DELRO/8.)/(ROP/3.+5/*DELRO/24.)
D0 10 J=2, JNM
IL=IMIN(J)+1
IH=IMAX(J)-1
IF(JA.EQ.JA1) GO TO 20
IF (J.EQ. JA. OR. J. EQ. JA1) IL=IAB+1
20 CONTINUE
IF (J.EQ.JB) IL=IAB+1
IF (J.EQ.JC) IH=IC-1
DO 10 I=IL, IH
IF (JRZ1. EQ. 1) GO TO 110
IF (K. EQ. NW. AND. (J. EQ. JA. AND. I.EQ. (IAB+1))) GO TO 100
IF (K. EQ. NW. AND. (J. EQ. JA1. AND. I.EQ. (IAB+1))) GO TO 100
IF (K. EQ. NW. AND. (J. EQ. JB. AND. I.EQ. (IAB+1))) GO TO 100
IF (K. EQ. NW. AND. (J. EQ. JD. AND. I.EQ. (IC+1)) .AND.JRZ3.EQ.1)
4 GO TO 100
5577 GO T0 110
CALL COFFCT (I,J,NF,INDG)
YY=A(I,J,NF)* (CE+CW+CN+CS)
YY=YY-(CE*A(I+1,J,NF)+ CW+A(I-1,J,NF)+ CN*A(I,J+1,NF)
1 +CS*A(I,J=1,NF))
Z=A(I,J,NW)
FDEQCT
A(I,J,NW)=YY/R(J)
GO TO 222
110 CONTINUE
CALL SORCCT(I,J,K,SOURCE,ZQ)
CALL COFFCT(I,J,K,INDG)
IF(K.NE.NW) GO TO 2
ZW=0.
ZWMU=0.0
GG1=0.0
IF(I.EQ.IL.AND.(J.GE.JB.OR.(J.LE.JA1.AND.J.GE.JA))) GO TO 500
IF(I.EQ.IH.AND.J.GE.JC) GO TO 600
GO TO 700
500 CONTINUE
IF(J.EQ.JB.OR.J.EQ.JA1.OR.J.EQ.JA.AND.JRZ1.NE.1) GO TO 999
501 LJNS=J+1
WKNS=FLOAT(1-INDG)
IF(J.EQ.JNM) RADNS=RADN
IF(J.EQ.(JB-1)) RADNS=RADB
IF(J.EQ.(JA-1)) RADNS=RADA
IF(J.EQ.(JC-1)) RADNS=RADC
CQNS=CN
ZWHALF=ZW
GO TO 666
502 LJNS=J-1
WKNS=FLOAT(INDG-1)
RADNS=RADA1
CQNS-CS
GO TO 666
600 LIEW=I+1
CQEW=CE
GO TO 555
700 IF(J.EQ.JNM) GO TO 501
IF (J.EQ.(JA-1).AND.I.LE.IAB) GO TO 501
IF (J.EQ.(JC-1).AND.I.GE.IC) GO TO 501
IF (J.EQ.(JA1+1).AND.I.LE.IAB) GO TO 502
IF((JB-JA1). GT.2.AND.(J.EQ.(JB-1).AND.I.LE.IAB)) GO TO 50
GO TO 999
555 DX1=ABS(X1(LIEW)-X1(I))
IF(JTAUO.EQ,1) GG1=G2(LIEW,J)
IF(JTAUO.NE,1) GG1=0.0
IF(IWZP.EQ,1.AND.I.EQ.IH.AND.J.GT.JD) GO TO 1117
RAD=R(J)
A(LIEW,J,NW)=WC1V(A(I,J,NF)-A(LIEW,J,NF),RO(LIEW,J),
1 RO(I,J)-RD(LIEW,J),R(J))
ZW=CQEW*WCZV(RO(LIEW,J),RO(I,J)-RO(LIEW,J))

FDEQCT
ZWMU=ZWMU+ZW*ZMU(LIEW,J)
1117 CONTINUE
GO TO 800
666 DX2=ABS(X2(LJNS)-X2(J))
IF(JTAUO.EQ.1) GG1=G1(I,LJNS)
IF(JTAUO.NE.1) GG1=0.0
A(I,LJNS,NW)=WC1H(WKNS.A(I,J,NF)-A(I,LJNS,NF),RO(I,LJNS),
1 RO(I,J)-RO(I,LJNS),RADNS)
ZW=CQNS*WCZH(WKNS,RO(I,LJNS),RO(I,J)-RO(I,LJNS),RADNS)
ZWMU=ZWMU+ZW*ZMU(I,LJNS)
IF(LJNS.EQ.(J+1)) GO TO 999
GO TO 900
800 IF(I.EQ.(IAB+1).AND.J,EQ.JNM) GO TO 501
IF(I.EQ.(IC-1).AND.J,EQ.JNM) GO TO 501
GO TO 999
900 IF((JB-JA1),EQ.2.AND,(J.EQ,(JB-1).AND.I.LE,IAB)) GO TO 501
999 CONTINUE
2 IF(K.NE.NF.AND.K.NE.NP) GO TO 1
IF(JRZ1.EQ,1) GO TO 1
IF(K,EQ,NF,AND(J.EQ.JA .AND,I.EQ.(IAB+1))) GO TO 223
IF(K,EQ,NF,AND(J.EQ.JA1 .AND,I.EQ.(IAB+1))) GO TO 223
IF(K,EQ,NF,AND(J.EQ.JB .AND,I.EQ.(IAB+1))) GO TO 223
IF(K,EQ,NF,AND(J.EQ.JD .AND,I.EQ.(IC+1)).AND.JR23,EQ,1)
5 GO TO 2223
4565 GO TO 1
223 A(I,J,NF)=A(I-1,J,NF)
GO TO 10
2223 A(I,J,NF)=A(I+1,J,NF)
GO TO 10
1 CONTINUE
CALL CONVEC (I,J,K,AU,ZU)
IF(K.EQ.NW) AU=AU*R(J)*R(J)
IF(K.EQ.NW) ZU=ZU*R(J)*R(J)
ANUM=CE*C(2)+A(I+1,J,K)+CW*C(4)*A(I-1,J,K)
1 +CN*C(3)*A(I,J+1,K) + CS*C(5)*A(I,J-1,K)+ AU+ SOURCE
ADNM=(CE+CW+CN+CS)+C(1)+ZU+ZO
Z=A(I,J,K)
A(I,J,K)=ANUM/ADNM
222 CONTINUE
DIF=A(I,J,K)=Z
A(I,J,K)=Z+RP(K) *DIF
IF(A(I,J,K).EQ,0.0) GO TO 1100
SSRS=(0.25*(A(I+1,J,K)+A(I-1,J,K)+A(I,J+1,K)+A(I,J-1,K)))
IF(SSRS.EQ,0.0) GO TO 1100
RSM=DIF/SSRS
RS=DIF/A(I,J,K)
IF(ABS(RSM).LT.ABS(RS)) RS=RSM
FDEQCT

IF(ABS(RS).GT.ABS(RSDU(K))) RSDU(K)=RS
IF(ABS(DIF).GT.ABS(DIFMAX(K))) DIFMAX(K)=DIF
1100 CONTINUE
10 CONTINUE
CALL BOUNCT(K)
RETURN
END
FDEQCT

SUBROUTINE MVBCCT
COMMON/CFTAU/ TJN,TBJ,TVJN,TVBJ
COMMON/CIWZP1/ IWZP,JO,ID,IF,JB1,IB1
COMMON/CJ2/ JTM,JPL,JV0,JTAU0,IJB1,JRZ1,JRZ2,JRZ3
COMMON/CPARZP/ VIN2,VT2,VT
COMMON/CR1WZP/ ROS,ROP,RO2
COMMON/CVP/ A(21,21,9),G1(21,21),G2(21,21),RO(21,21),ZMU(21,21)
COMMON/CGRID/ IMIN(21),IMAX(21),IN,JN,INM,JNM,IE
COMMON/CPARAM/ VINP,VINS
COMMON/CREFER/ ROREF,PREF,ZMUREF,TREF,CPREF
COMMON/CNUMBR/ NW,NF,NZK,NHS,NVT,NZML,NZMM,NZMH,NP
COMMON/CYLI/ JA,JA1,JB,JC,IAB,UC,RADA,RADA1,RADB,RADC,RADN,
1 DA,DA1,DB,DCN,DN,X2AXIS
COMMON/CCONST/ PR(9),ZMW(3),CPJ(3),GCPM,GC,ZJC
COMMON/CX1X2/ X1(21),X2(21),DELX1(21), DELX2(21) /CVEXIT/ VEXIT
COMMON/CRAD/ R(21)
COMMON/CINDEX/ INDE(9),INDG,INDRO,INDZMU
YY(YN,YP)=1./(1-(YP/YN)**2)
DIMENSION AQ(21,21,14)
EQUIVALENCE(A(1,1,1,),AQ(1,1,1,))
RR=PREF/(GCPM*TREF)
ROP=ZMW(1)*RR
IF(JA.EQ.JB) ROP=RO(1,1)
ROS=ZMW(2)*RR
DO 1 J=1,JA
IF(JV0.EQ.0.AND.(J.EQ.1.OR.J.EQ.JA)) GO TO 1510
G1(1,J)=ROP*VINP
1510 CONTINUE
1 G2(1,J)=0.0
DO 2 J=JA1<JB
IF(JV0.EQ.AND.(J.EQ.JA1.OR.J.EQ.JB)) GO TO 1511
G1(1,J)=ROS*VINS
1511 CONTINUE
2 G2(1,J)=0.0
IF(IWZP.NE.1) GO TO 1190
DO 1020 J=JD,JN
IF (JV0.EQ.0.AND.(J.EQ.JD .OR.J.EQ.JN)) GO TO 1512
G1(IC,J)=-RO2*VIN2
1512 CONTINUE
G2(IC,J)=0.0
1020 CONTINUE
1190 CONTINUE
JJC=JC-1
DO 31 J=2,JJC
G1(IN,J)=ADF(IN,J,2,NF)/R(J)
G2(IN,J)0.0
31 CONTINUE
DX2=X2(2)-X2(1)
MVBCCT

BB=1./(1.-((X2(2)-X2(1))/(X2(3)-X2(1)))**2)
DO 10 I=2,IN
IF(X2AXIS.EQ.0.0) G1(I,1)=BB*G1(I,2)-(BB-1.)*G1(I,3)
IF(X2AXIS.GT.0.0) G1(I,1)=0.0
10 G2(I,1)=0.0
ROSUM=0.0
DO 20 J=1,JC
20 ROSUM=ROSUM+RO(IN,J)
ROMEAN=ROSUM/FLOAT(JC)
IF(INDG. EQ.1)VEXIT=-A(IN,1,NF)/(X2(JC)*ROMEAN)
IF(INDG. EQ.2)VEXIT=-2.*A(IN,1,NF)/((RADC*RADC-X2AXIS)
1 *ROMEAN)
G1(1,JA)=0.0
G1(1,JA1)=0.0
G1(1,JB)=0.0
RETURN
END

FUNCTION ADF(I,J,LX,KQ)
COMMON/CIWZP1/ IWZP,JD,ID,IF,JB1,IB1
COMMON/CJ2/ JTM,JPL,JV0,JTAU0,IJB1,JRZ1,JRZ2,JRZ3
COMMON/CVP/ A(21,21,9),G1(21,21),G2(21,21),RO(21,21),ZMU(21,21)
COMMON/CYLI/ JA,JA1,JB,JC,IAB,IC,RADA,RADA1,RADB,RADC,RADN,
1 DA,DA1,DB,DCN,DN,X2AXIS
COMMON/CX1X2/ X1(21),X2(21),DELX1(21),DELX2(21)
COMMON/CGRID/ IMIN(21),IMAX(21),IN,JN,INM,JNM,IE
DIMENSION AQ(21,21,14)
DIMENSION BENQ(20),BWSR(20),BP(20)
EQUIVALENCE (A(1,1,1), AQ(1,1,1))
DF(PN,BENQ,BWSR,BP,XENQ,XWSR)=((XENQ*XENQ-XWSR*XWSR)*BP+
1 XWSR*XWSR*BENQ-XENQ*XENQ*BWSR)/(XENQ*XWSR*(PN*XENQ+XWSR))
M=1
IF(I.LT.IAB.AND.I.GE.1) GO TO 300
IF(J.EQ.1) GO TO 92
IF(I.EQ.IAB) GO TO 93
IF(I.EQ.JN.OR.(J.EQ.JC.AND. I .LE.IN.AND.I.GE.IC)) GO TO 94
IF(I.EQ.IC.AND.J.LE.JN.AND .J.GE.JC).OR,
1 (I.EQ.IN.AND.J.LE.JC.AND .J.GE.1)) M=S
GO TO 99
300 IF(J.EQ.1.OR.J.EQ.JA1) GO TO 92
IF(IJB1.EQ.1.AND.J.EQ.JB1) GO TO 92
IF(I.EQ.1) GO TO 93
IF(J.EQ.JB.OR.J.EQ.JA) GO TO 94
GO TO 99
92 M=2
GO TO 99
93 M=3
GO TO 99
94 M=4
CONTINUE
GO TO (1,2,3,4,5),M
1 PN=1.
IF(LX.EQ.1) GO TO 12
13 IF(J.EQ.1.OR.J.EQ.JA1) GO TO 21
IF(J.EQ.JN.OR.(J.EQ.JC.AND.I.EQ.IN).OR.(I.LT.IAB.AND.
1 (J.EQ.JB.CR.J.EQ.JA))) GO TO 41
DO 10 L=1.14
BENQ(L)=AQ(I,J+1,L)
BWSR(L)=AQ(I,J-1,L)
10 BP(L)=AQ(I,J,L)
XENQ=X2(J+1)=X2(J)
XWSR=X2(J) =X2(J-1)
GO TO 100
12 IF(I.EQ.IAB.OR.I.EQ.1) GO TO 31
IF((I.EQ.IC.AND.J.EQ.JN).OR.I.EQ.IN) GO TO 51
DO 11 L=1.14

BENQ(L)=AQ(I+1,J,L)
BWSR(L)=AQ(I-1,J,L)
BP(L)=AQ(I,J,L)
XENQ=X1(I+1)-X1(I)
XWSE=X1(I)-X1(I-1)
GO TO 100
2 IF(LX.EQ.1) GO TO 1
21 PN=-1.
DO 20 L=1,14
BENQ(L)=AQ(I,J+1,L)
BWSE(L)=AQ(I,J+2,L)
20 BP(L)=AQ(I,J,L)
XENQ=X2(J+1)-X2(J)
XWSR=X2(J+2)-X2(J)
GO TO 100
3 IF(LX.NE.1) GO TO 1
31 PN=-1.
DO 30 L=1,14
BENQ(L)=AQ(I+1,J,L)
BWSE(L)=AQ(I+1,J,L)
30 BP(L)=AQ(I,J,L)
XENQ=X1(I+1)-X1(I)
XWSR=X2(I+2)-X2(I)
GO TO 100
4 IF(LX.EQ.1) GO TO 1
41 PN=-1.
DO 40 L-1,14
BENQ(L)=AQ(I,J-1,L)
BWSE(L)=AQ(I,J-2,L)
40 BP(L)=AQ(I,J,L)
XENQ=X2(J-1)-X2(J)
XWSR=X2(J-2)-X2(J)
GO TO 100
5 IF(LX.NE,1) GO TO 1
51 PN=-1.
DO 50 L=1,14
BENQ(L)=AQ(I-1,J,L)
BWSE(L)=AQ(I-2,J,L)
50 BP(L)=AQ(I,J,L)
XENQ=X1(J-1)-X1(J)
XWSR=X1(J-2)-X1(J)
100 CONTINUE
BENQ(15)=BENQ(10)/BENQ(12)
BENQ(16)=BENQ(11)/BENQ(12)
BWSR(15)=BWSR(10)/BWSR(12)
BWSR(16)=BWSR(11)/BWSR(12)
BP(15)=BP(10)/BP(12)
BP(16)=BP(11)/BP(12)
ADF=DF(PN,BENQ(KQ),BWSR(KQ),BP(KQ),XENQ,XWSR)
RETURN
END

SUBROUTINE COEFCT(I,J,K,INDG)
COMMON/CVP/ A(21,21,9), G1(21,21),G2(21,21),RO(21,21),ZMU(21,21)
COMMON/CXIX2/ X1(21),X2(21),DELX1(21),DELX2(21)
COMMON/CCOEF/ CE,CW,CN,CS,O(5)
COMMON/CCONST/ PR(9),ZMW(3),OPJ(3),GCPM,GC,ZJC
COMMON/CRAD/ R(21)
COMMON/CNUMBR/ NW,NF,NZK,NHS,NVT,NZML,NZMM,NZMN,NP
1 DA,DA1,DB,DON,DN,X2AXIS
DIMENSION B(5)
DO 16 L=1,5
16 C(L)=1.0
IF(K.NE.NW.AND.K.NE.NP) GO TO 2
B(1)=R(J)**3
B(2)=B(1)
B(3)=R(J+1)**3
B(4)=B(1)
B(5)=R(J-1)**3
C(1)=ZMU(I,J)
C(2)=ZMU(I+1,J)
C(3)=ZMU(I,J+1)
C(4)=ZMU(I-1,J)
C(5)=ZMU(I,J-1)
GO TO 20
2 IF(K.NE.NF) GO TO 3
B(1)=1./(R(J)*RO(I,J))
B(2)=1./(R(J)*RO(I+1,J))
B(3)=8.0/((R(J)+R(J+1))*(RO(I,J)+RO(I,J+1)))-B(1)
B(4)=1.0/(R(J)*RO(I-1,J))
B(5)=8.0/((R(J)+R(J-1))*(RO(I,J)+RO(I,J-1)))-B(1)
GO TO 20
3 IF (K.NE.NVT) GO TO 15
B(1)=ZMU(I,J)*R(J)**3
B(2)=ZMU(I+1,J)*R(J)**3
B(3)=ZMU(I,J+1)*R(J+1)**3
B(4)=ZMU(I-1,J)*R(J)**3
B(5)=ZMU(I,J-1)*R(J-1)**3
C(1)=1./(R(J)R(J))
C(2)=C(1)
C(3)=1./(R(J+1)*R(J+1))
C(4)=C(1)
IF(R(J-1).GT.0.0) C(5)=1./(R(J-1)*R(J-1))
IF(R(J-1).EQ.0.0) C(5)-0.0
GO TO 20
15 B(1)=ZMU(I,J)*R(J)/PR(K)
B(2)=ZMU(I+1,J)*R(J)/PR(K)
B(3)=ZMU(I,J+1)*R(J+1)/PR(K)
B(4)=ZMU(I-1,J)*R(J)/PR(K)
B(5)=ZMU(I,J-1)*R(J-1)/PR(K)
20 CONTINUE
CE=(B(2)+B(1))/((X1(I+1)-X1(I))*2.*DELX1(I))
CW=(B(1)+B(4))/((X1(I)-X1(I-1))*2.*DELX1(I))
CN=(B(3)+B(1))/((X2(J+1)-X2(J))*2.*DELX2(J))
CS=(B(1)+B(5))/((X1(J)-X2(J-1))*2.*DELX2(J))
RETURN
END

SUBROUTINE CONVEC(I,J,K,AU,ZU)
COMMON/CVP/ A(21,21,9),G1(21,21) G2(21,21),RO(21,21),ZMU(21,21)
COMMON/CX1X2/ X1(21),X2(21),DELX1(21),DELX2(21)
COMMON/CNUMBR/ NW,NF,NZK,NHS,NVT,NZML,NZMM,NZMN,NP
COMMON/CYLI/ JA,JA1,JB,JC,IAB,IC,RADA,RADA1,RADB,RADC,RADN,
1 DA,DA1,DB,DCN,DN,X2AXIS
COMMON/CRAD/ R(21)
COMMON/CZW/ ZW,ZWHALF,IL,IH
COMMON/CGRID/ IMIN(21),IMAX(21),IN,JN,INM,JNM,IE
AU=0.
ZU-0.
IF(K.EQ.NF) GO TO 12
GX12=1./(4.*DELX1(I)*DELX2(J))
GM1E= DX12*(A(I,J+1,NF)-A(I,J-1,NF)+A(I+1,J+1,NF)-A(I+1,J-1,NF))
GM1W= DX12*(A(I,J+1,NF)-A(I,J-1,NF)+A(I-1,J+1,NF)-A(I-1,J-1,NF))
GM2S=-DX12*(A(I+1,J,NF)-A(I-1,J,NF)+A(I+1,J-1,NF)-A(I-1,J-1,NF))
GM2N=-DX12*(A(I+1,J,NF)-A(I-1,J,NF)+A(I+1,J+1,NF)-A(I-1,J-1,NF))
IF(GM1W)1,2,3
1 ZU=-GM1W
GO TO 3
2 AU=GH1W*A(I-1,J,K)
IF(K.EQ.NW.AND.
1 (I.EQ.IL.AND.(J.GE.JB.OR.(J.LE.JA1.AND.J.GE.JA)))) GO TO 500
3 IF(GM2S)4,6,5
4 ZU=ZU-GM2S
GO TO 6
5 AU=AU+GM2S*A(I,J-1,K)
IF(K.EQ.NW.AND.J.EQ.(JA1+1).AND.I.LE.IAB) GO TO 501
6 IF(GM1E)7,9,8
7 AU=AU-GM1E*A(I+1,J,K)
IF(K.EQ.NW.AND.(I.EQ.IH.AND.J.GE.JC)) GO TO 502
GO TO 9
8 ZU=ZU+GM1E
9 IF(GM2N)10,11,11
10 AU=AU-GMZN*A(I,J+1,K)
IF(K.EQ.NW.AND.(J.EQ.JNM.OR.J.EQ.(JB-1).OR.J.E
1 J.EQ.(JC-1))) GO TO 503
RETURN
11 ZU=ZU+GM2N
12 RETURN
500 IF(J.NE.JNM) ZU=ZU-GM1W*ZW/CW
IF(J.EQ.JNM) ZU=ZU-GM1W*ZWHALF/CW
GO TO 3
501 IF((JB-JA1).GT.2) ZU=ZU-GM2S*ZW/CS
IF((JB-JA1).EQ.2) ZU=ZU-GM2S*ZWHALF/CS
GO TO 6
502 IF(J.NE.JNM ) ZU=ZU+GM1E*ZW /CE
IF(J.EQ.JNM ) ZU=ZU+GM1E*ZWHALF/CE
GO TO 9
503 IF(J.EQ.JNM ) GO TO 504
IF(J.EQ.(JC-1).OR.J.EQ. (JA-1)) GO TO 504
IF(J.EQ.(JB-1).AND.((JB-JA1).GE.2)) GO TO 504
RETURN
504 ZU=ZU+GM2N*ZW/CN
RETURN
END