TALK=T;RUN( 1, 1)
** LOAD(x577) from the x Input Library
GROUP 1. Run title and other preliminaries
****Q1 created by V. Semin for WUA-95/96
3d,BFC,Steady.
TEXT( CAR body.
BOOLEAN(LVISC,LTURB,LGEM); LVISC=T; LTURB=T; LGEM=T
REAL(UIN,DTHYD,REYNO)
REAL(LCAR,HCAR,ZB,ZLL,RWL,DXWL,DZWL,XW1,XW2,YW1,YW2,YLL)
REAL(X0,X00,X1,X2,X3,X4,X5,X6,X7,X8,X9,X10,X11,X12,X13)
REAL(Z00,Z3T,Z4T,Z5T,Z6T,Z7T,Z8T,Z10T,Z11B,Z11T,Z12B,Z12T)
REAL(XCR,ZCRB,ZCRT,PWY,PWX,PWZ,PWZY,PWYZ)
INTEGER(NZ1,NZ2,NZ3,NZT,NX1,NX2,NX3,NX4,NX5,NX6,NX7,NX8)
INTEGER(IS,JS,IS2,JS2,NTT,NT2)
INTEGER(NY1,NY2,NY3,NY4,NYT,IIF,IIL,JJF,JJL,KKF,KKL)
REAL(TKEIN,EPSIN,FRIC,LMIX,KTDKP,CF,TU)
UIN = 50.0
REYNO = 2.64E6
IF(LTURB) THEN
+ CF= GRND2
ELSE
+ CF= 1.0
ENDIF
** Z-direction.
NZ1 = 6; NZ2 = 10; NZ3 = 16
NZ = NZ1+NZ2+NZ3
** X-direction.
** Extra regions.
INTEGER(NX31,NX32,NX33,NX51,NX52,NX53,NX54,NX61,NX62)
NX1 = 3; NX2 = 9
NX31= 1; NX32= 3; NX33= 1; NX3=NX31+NX32+NX33
NX4=3
NX51= 2; NX52= 3; NX53= 2; NX54= 4; NX5=NX51+NX52+NX53+NX54
NX61= 1; NX62= 2; NX6=NX61+NX62
NX7=2; NX8=20
NX = NX1+NX2+NX3+NX4+NX5+NX6+NX7+NX8
** Y-direction.
NY1=4; NY2=2; NY3=1; NY4=10
NY = NY1+NY2+NY3+NY4
** CAR geometry.
LCAR = 0.810; HCAR = 0.27; ZB = -0.0335; ZLL = ZB+4*HCAR
X0 = 0.002; YLL = 3*HCAR
** Wheel parameters.
RWL=0.130/2.0; DXWL=RWL/SQRT(2.0)
** Position in Z direction.
DZWL=RWL+ZB-0.004205; DZWL= -DXWL+DZWL
** Position in X direction.
XW1=X0+0.150; XW2 = XW1+0.540
** Position in Y direction.
YW1=0.097; YW2=0.097+0.026
** Nonuniformity parameters.
PWX= 1.2; PWY= 1.6; PWZ= 1.6; PWZY= 1.4; PWYZ=1.4
X00= 0.069999; Z00 = 0.09548
X1 = X0-0.44
X2 = X1-0.5*LCAR
X3 = 0.09301; Z3T = 0.107799
X4 = 0.129; Z4T = 0.124925
X5 = 0.175; Z5T = 0.14466
X6 = 0.21099; Z6T = 0.16021
X7 = 0.63301; Z7T = 0.22544
X8 = 0.665; Z8T = 0.22106
X10= 0.711; Z10T=0.021443
X11= 0.74699; Z11B=0.013841; Z11T =0.20772
X12= 0.81; Z12B=0.023221; Z12T =0.19222
X13= X0+3.0*LCAR
** Turbulence setting.
TU = 0.05
LMIX = 0.09*10*LCAR/SQRT(REYNO)
TKEIN= UIN**2*TU**2; EPSIN= 0.1643*TKEIN**1.5/LMIX
GROUP 2. Transience; time-step specification
GROUP 3. X-direction grid specification
GROUP 4. Y-direction grid specification
GROUP 5. Z-direction grid specification
GROUP 6. Body-fitted coordinates or grid distortion
BFC = T; NONORT= T
GSET(D,NX,NY,NZ)
IF(LGEM) THEN
NTT = NY1+NY2+NY3
** Section 1
IS = NX1+1
IS
GSET(P,P1,X1,0.0 ,ZB ); GSET(P,P2,X1,0.128,ZB )
GSET(P,P3,X1,0.128,0.0); GSET(P,P4,X1,0.0 ,0.0)
GSET(P,P6,X1,0.0 ,Z00); GSET(P,P5,X1,0.128,0.1)
GSET(P,P8,X1,0.0 ,ZLL); GSET(P,P7,X1,0.128,ZLL)
GSET(P,A1,X1,YLL,ZB ); GSET(P,A2,X1,YLL,0.0)
GSET(P,A4,X1,YLL,ZLL); GSET(P,A3,X1,YLL,0.1)
CALL STLN2
CALL STFM2
GSET(I)
** Section 2.
IS = NX1+NX2+1
IS
GSET(P,P1,X0,0.0 ,ZB )
GSET(P,P2,X0,0.5895498E-01,ZB )
GSET(P,P3,X0,0.5895498E-01,0.1325599E-01)
GSET(P,P4,X0,0.0 ,0.1273900E-01)
GSET(P,P6,X0,0.0 ,0.2516499E-01)
GSET(P,P5,X0,0.5891598E-01,0.2345898E-01)
GSET(P,P8,X0,0.0 ,ZLL )
GSET(P,P7,X0,0.5891598E-01,ZLL )
GSET(P,A1,X0,YLL,ZB ); GSET(P,A2,X0,YLL,0.1325599E-01)
GSET(P,A4,X0,YLL,ZLL); GSET(P,A3,X0,YLL,0.2345898E-01)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,X0,0.5421098E-01,0.1273900E-01)
+ GSET(V,X0,0.5713499E-01,0.1286799E-01)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NY1+NY2+NY3,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,X0,0.5706197E-01,0.2388699E-01)
+ GSET(V,X0,0.5452598E-01,0.2417799E-01)
+ GSET(V,X0,0.5216499E-01,0.2437498E-01)
+ GSET(V,X0,0.5002899E-01,0.2450999E-01)
+ GSET(V,X0,0.4809999E-01,0.2460599E-01)
+ GSET(V,X0,0.4640798E-01,0.2467499E-01)
+ GSET(V,X0,0.4493598E-01,0.2472600E-01)
+ GSET(V,X0,0.4360398E-01,0.2476599E-01)
+ GSET(V,X0,0.4227598E-01,0.2480198E-01)
+ GSET(V,X0,0.4076399E-01,0.2483799E-01)
+ GSET(V,X0,0.3882899E-01,0.2487799E-01)
+ GSET(V,X0,0.3617498E-01,0.2492599E-01)
+ GSET(V,X0,0.3245099E-01,0.2498098E-01)
+ GSET(V,X0,0.2765999E-01,0.2503698E-01)
+ GSET(V,X0,0.2221899E-01,0.2508598E-01)
+ GSET(V,X0,0.1670999E-01,0.2512098E-01)
+ GSET(V,X0,0.1180499E-01,0.2514399E-01)
+ GSET(V,X0,0.7887997E-02,0.2515599E-01)
+ GSET(V,X0,0.4789997E-02,0.2516199E-01)
+ GSET(V,X0,0.2252999E-02,0.2516398E-01)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NY1+NY2+NY3,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X0,0.5976399E-01,0.1389600E-01)
+ GSET(V,X0,0.6040998E-01,0.1477500E-01)
+ GSET(V,X0,0.6097898E-01,0.1585799E-01)
+ GSET(V,X0,0.6141898E-01,0.1708699E-01)
+ GSET(V,X0,0.6167999E-01,0.1839599E-01)
+ GSET(V,X0,0.6171398E-01,0.1970799E-01)
+ GSET(V,X0,0.6144798E-01,0.2093298E-01)
+ GSET(V,X0,0.6084998E-01,0.2199299E-01)
+ GSET(V,X0,0.5997598E-01,0.2283899E-01)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
** Section 3
IS = NX1+NX2+NX31+1
IS
XCR = 0.1799998E-01
GSET(P,P1,XCR,0.0 ,ZB )
GSET(P,P2,XCR,0.1168119E+00,ZB )
GSET(P,P3,XCR,0.1168119E+00,0.6082997E-02)
GSET(P,P4,XCR,0.0 ,0.3732998E-02)
GSET(P,P6,XCR,0.0 ,0.6019798E-01)
GSET(P,P5,XCR,0.1167349E+00,0.5244798E-01)
GSET(P,P8,XCR,0.0 ,ZLL )
GSET(P,P7,XCR,0.1167349E+00,ZLL )
GSET(P,A1,XCR,YLL,ZB ); GSET(P,A2,XCR,YLL,0.6082997E-02)
GSET(P,A4,XCR,YLL,ZLL); GSET(P,A3,XCR,YLL,0.5244798E-01)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,XCR,0.1074129E+00,0.3732998E-02)
+ GSET(V,XCR,0.1132069E+00,0.4321996E-02)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NY1+NY2+NY3,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,XCR,0.1130620E+00,0.5439299E-01)
+ GSET(V,XCR,0.1080359E+00,0.5571298E-01)
+ GSET(V,XCR,0.1033589E+00,0.5660798E-01)
+ GSET(V,XCR,0.9912694E-01,0.5722298E-01)
+ GSET(V,XCR,0.9530497E-01,0.5765798E-01)
+ GSET(V,XCR,0.9195095E-01,0.5797099E-01)
+ GSET(V,XCR,0.8903497E-01,0.5820199E-01)
+ GSET(V,XCR,0.8639693E-01,0.5838499E-01)
+ GSET(V,XCR,0.8376598E-01,0.5854698E-01)
+ GSET(V,XCR,0.8076894E-01,0.5871098E-01)
+ GSET(V,XCR,0.7693595E-01,0.5889498E-01)
+ GSET(V,XCR,0.7167596E-01,0.5911198E-01)
+ GSET(V,XCR,0.6429893E-01,0.5936299E-01)
+ GSET(V,XCR,0.5480598E-01,0.5961797E-01)
+ GSET(V,XCR,0.4402398E-01,0.5983797E-01)
+ GSET(V,XCR,0.3310899E-01,0.5999998E-01)
+ GSET(V,XCR,0.2339099E-01,0.6010098E-01)
+ GSET(V,XCR,0.1562900E-01,0.6015498E-01)
+ GSET(V,XCR,0.9490997E-02,0.6018199E-01)
+ GSET(V,XCR,0.4464999E-02,0.6019499E-01)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NY1+NY2+NY3,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,XCR,0.1184149E+00,0.8991998E-02)
+ GSET(V,XCR,0.1196949E+00,0.1298800E-01)
+ GSET(V,XCR,0.1208230E+00,0.1790498E-01)
+ GSET(V,XCR,0.1216940E+00,0.2349199E-01)
+ GSET(V,XCR,0.1222119E+00,0.2943799E-01)
+ GSET(V,XCR,0.1222789E+00,0.3539999E-01)
+ GSET(V,XCR,0.1217520E+00,0.4096899E-01)
+ GSET(V,XCR,0.1205670E+00,0.4578399E-01)
+ GSET(V,XCR,0.1188359E+00,0.4962998E-01)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
IS = NX1+NX2+NX31+NX32+1
IS
GSET(P,P1,X00,0.0 ,ZB )
GSET(P,P2,X00,0.1286420E+00,ZB )
GSET(P,P3,X00,0.1286420E+00,0.3971998E-02)
GSET(P,P4,X00,0.0 ,0.0 )
GSET(P,P6,X00,0.0 ,Z00 )
GSET(P,P5,X00,0.1285570E+00,0.8237797E-01)
GSET(P,P8,X00,0.0 ,ZLL )
GSET(P,P7,X00,0.1285570E+00,ZLL )
GSET(P,A1,X00,YLL,ZB ); GSET(P,A2,X00,YLL,0.3971998E-02)
GSET(P,A4,X00,YLL,ZLL); GSET(P,A3,X00,YLL,0.8237797E-01)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,X00,0.1182919E+00,0.0000000E+00)
+ GSET(V,X00,0.1246719E+00,0.9949997E-03)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NTT,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,X00,0.1245130E+00,0.8566594E-01)
+ GSET(V,X00,0.1189770E+00,0.8789897E-01)
+ GSET(V,X00,0.1138260E+00,0.8941197E-01)
+ GSET(V,X00,0.1091670E+00,0.9045196E-01)
+ GSET(V,X00,0.1049569E+00,0.9118795E-01)
+ GSET(V,X00,0.1012639E+00,0.9171695E-01)
+ GSET(V,X00,0.9805197E-01,0.9210795E-01)
+ GSET(V,X00,0.9514594E-01,0.9241796E-01)
+ GSET(V,X00,0.9224898E-01,0.9269196E-01)
+ GSET(V,X00,0.8894795E-01,0.9296894E-01)
+ GSET(V,X00,0.8472693E-01,0.9328097E-01)
+ GSET(V,X00,0.7893497E-01,0.9364694E-01)
+ GSET(V,X00,0.7081097E-01,0.9407097E-01)
+ GSET(V,X00,0.6035598E-01,0.9450293E-01)
+ GSET(V,X00,0.4848298E-01,0.9487396E-01)
+ GSET(V,X00,0.3646199E-01,0.9514797E-01)
+ GSET(V,X00,0.2575999E-01,0.9531897E-01)
+ GSET(V,X00,0.1721199E-01,0.9541094E-01)
+ GSET(V,X00,0.1045300E-01,0.9545696E-01)
+ GSET(V,X00,0.4916996E-02,0.9547794E-01)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X00,0.1304069E+00,0.8891996E-02)
+ GSET(V,X00,0.1318169E+00,0.1564899E-01)
+ GSET(V,X00,0.1330599E+00,0.2396499E-01)
+ GSET(V,X00,0.1340190E+00,0.3341099E-01)
+ GSET(V,X00,0.1345890E+00,0.4346699E-01)
+ GSET(V,X00,0.1346629E+00,0.5354899E-01)
+ GSET(V,X00,0.1340820E+00,0.6296593E-01)
+ GSET(V,X00,0.1327779E+00,0.7110995E-01)
+ GSET(V,X00,0.1308709E+00,0.7761294E-01)
+ GSET(V,X00,0.1285570E+00,0.8237797E-01)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
IS = NX1+NX2+NX31+NX32+NX33+1
IS
GSET(P,P1,X4,0.0 ,ZB )
GSET(P,P2,X4,0.1304629E+00,ZB )
GSET(P,P3,X3,0.1304629E+00,0.4483998E-02)
GSET(P,P4,X3,0.0 ,0.0 )
GSET(P,P6,X3,0.0 ,Z3T )
GSET(P,P5,X3,0.1303759E+00,0.9294093E-01)
GSET(P,P8,X3,0.0 ,ZLL )
GSET(P,P7,X3,0.1303759E+00,ZLL )
GSET(P,A1,X4,YLL,ZB ); GSET(P,A2,X3,YLL,0.4483998E-02)
GSET(P,A4,X3,YLL,ZLL); GSET(P,A3,X3,YLL,0.9294093E-01);
GSET(P,B1,X4,YW1,ZB ); GSET(P,B2,X4,YW2,ZB )
GSET(P,B4,X3,YW1,0.0); GSET(P,B3,X3,YW2,0.0)
CALL STLN3
GSET(P,PA,XW1-DXWL, YW1,DZWL)
GSET(P,PB,XW1-DXWL, YW2,DZWL)
GSET(P,PC,XW1-DXWL,0.1304629E+00,DZWL)
GSET(L,B14,B1,B4,NZ1,1,ARC,PA)
GSET(L,B23,B2,B3,NZ1,1,ARC,PB)
GSET(L,L23,P2,P3,NZ1,1,ARC,PC)
GSET(V,CV2,S,P5)
+ GSET(V,X3,0.1262749E+00,0.9665096E-01)
+ GSET(V,X3,0.1206610E+00,0.9916997E-01)
+ GSET(V,X3,0.1154370E+00,0.1008769E+00)
+ GSET(V,X3,0.1107119E+00,0.1020499E+00)
+ GSET(V,X3,0.1064420E+00,0.1028829E+00)
+ GSET(V,X3,0.1026970E+00,0.1034830E+00)
+ GSET(V,X3,0.9943897E-01,0.1039270E+00)
+ GSET(V,X3,0.9649295E-01,0.1042799E+00)
+ GSET(V,X3,0.9355396E-01,0.1045920E+00)
+ GSET(V,X3,0.9020692E-01,0.1049089E+00)
+ GSET(V,X3,0.8592695E-01,0.1052650E+00)
+ GSET(V,X3,0.8005297E-01,0.1056849E+00)
+ GSET(V,X3,0.7181293E-01,0.1061710E+00)
+ GSET(V,X3,0.6120998E-01,0.1066679E+00)
+ GSET(V,X3,0.4916899E-01,0.1070960E+00)
+ GSET(V,X3,0.3697799E-01,0.1074119E+00)
+ GSET(V,X3,0.2612399E-01,0.1076089E+00)
+ GSET(V,X3,0.1745599E-01,0.1077150E+00)
+ GSET(V,X3,0.1060000E-01,0.1077690E+00)
+ GSET(V,X3,0.4985996E-02,0.1077930E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X3,0.1322529E+00,0.1003800E-01)
+ GSET(V,X3,0.1336830E+00,0.1766398E-01)
+ GSET(V,X3,0.1349429E+00,0.2704699E-01)
+ GSET(V,X3,0.1359159E+00,0.3770599E-01)
+ GSET(V,X3,0.1364930E+00,0.4904899E-01)
+ GSET(V,X3,0.1365690E+00,0.6042198E-01)
+ GSET(V,X3,0.1359800E+00,0.7104397E-01)
+ GSET(V,X3,0.1346570E+00,0.8022994E-01)
+ GSET(V,X3,0.1327239E+00,0.8756596E-01)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM3
GSET(I)
** Section 4.
IS = NX1+NX2+NX3+NX4+1
IS
GSET(P,P1,X5,0.0 ,ZB )
GSET(P,P2,X5,0.1351299E+00,ZB )
GSET(P,P3,X6,0.1351299E+00,0.5889997E-02)
GSET(P,P4,X6,0.0 ,0.0 )
GSET(P,P6,X6,0.0 ,Z6T )
GSET(P,P5,X6,0.1351950E+00,0.1177869E+00)
GSET(P,P8,X6,0.0 ,ZLL )
GSET(P,P7,X6,0.1351950E+00,ZLL )
GSET(P,A1,X5,YLL,ZB ); GSET(P,A2,X6,YLL,0.5889997E-02)
GSET(P,A4,X6,YLL,ZLL); GSET(P,A3,X6,YLL,0.1177869E+00);
GSET(P,B1,X5,YW1,ZB ); GSET(P,B2,X5,YW2,ZB )
GSET(P,B4,X6,YW1,0.0); GSET(P,B3,X6,YW2,0.0)
CALL STLN3
** Wheels.
GSET(P,PA,XW1+DXWL,YW1, DZWL)
GSET(P,PB,XW1+DXWL,YW2, DZWL)
GSET(P,PC,XW1+DXWL,0.1351299E+00,DZWL)
GSET(L,B14,B1,B4,NZ1,1,ARC,PA)
GSET(L,B23,B2,B3,NZ1,1,ARC,PB)
GSET(L,L23,P2,P3,NZ1,1,ARC,PC)
GSET(V,CV2,S,P5)
+ GSET(V,X6,0.1313490E+00,0.1226130E+00)
+ GSET(V,X6,0.1261470E+00,0.1260729E+00)
+ GSET(V,X6,0.1210970E+00,0.1287400E+00)
+ GSET(V,X6,0.1163089E+00,0.1310909E+00)
+ GSET(V,X6,0.1120340E+00,0.1334940E+00)
+ GSET(V,X6,0.1082870E+00,0.1360260E+00)
+ GSET(V,X6,0.1048660E+00,0.1386729E+00)
+ GSET(V,X6,0.1016359E+00,0.1413890E+00)
+ GSET(V,X6,0.9837693E-01,0.1441090E+00)
+ GSET(V,X6,0.9477395E-01,0.1467630E+00)
+ GSET(V,X6,0.9039593E-01,0.1492869E+00)
+ GSET(V,X6,0.8469695E-01,0.1516320E+00)
+ GSET(V,X6,0.7702994E-01,0.1537600E+00)
+ GSET(V,X6,0.6731594E-01,0.1556109E+00)
+ GSET(V,X6,0.5624198E-01,0.1571239E+00)
+ GSET(V,X6,0.4474899E-01,0.1582769E+00)
+ GSET(V,X6,0.3379798E-01,0.1590899E+00)
+ GSET(V,X6,0.2403399E-01,0.1596299E+00)
+ GSET(V,X6,0.1536100E-01,0.1599699E+00)
+ GSET(V,X6,0.7474996E-02,0.1601579E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X6,0.1369870E+00,0.1313899E-01)
+ GSET(V,X6,0.1384739E+00,0.2301599E-01)
+ GSET(V,X6,0.1397809E+00,0.3504899E-01)
+ GSET(V,X6,0.1407849E+00,0.4857199E-01)
+ GSET(V,X6,0.1413779E+00,0.6282097E-01)
+ GSET(V,X6,0.1414599E+00,0.7699198E-01)
+ GSET(V,X6,0.1408550E+00,0.9019792E-01)
+ GSET(V,X6,0.1394979E+00,0.1016740E+00)
+ GSET(V,X6,0.1375329E+00,0.1109279E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM3
GSET(I)
** Section 5.
IS = NX1+NX2+NX3+NX4+NX51+1
IS
XCR = 0.2939997E+00
GSET(P,P1,XCR,0.0 ,ZB )
GSET(P,P2,XCR,0.1364899E+00,ZB )
GSET(P,P3,XCR,0.1364899E+00,0.5576998E-02)
GSET(P,P4,XCR,0.0 ,0.0 )
GSET(P,P6,XCR,0.0 ,0.1942750E+00)
GSET(P,P5,XCR,0.1369550E+00,0.1220250E+00)
GSET(P,P8,XCR,0.0 ,ZLL )
GSET(P,P7,XCR,0.1369550E+00,ZLL )
GSET(P,A1,XCR,YLL,ZB ); GSET(P,A2,XCR,YLL,0.5576998E-02)
GSET(P,A4,XCR,YLL,ZLL); GSET(P,A3,XCR,YLL,0.1220250E+00)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,XCR,0.1255149E+00,0.0000000E+00)
+ GSET(V,XCR,0.1322809E+00,0.1391999E-02)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NY1+NY2+NY3,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,XCR,0.1341159E+00,0.1275269E+00)
+ GSET(V,XCR,0.1304519E+00,0.1319060E+00)
+ GSET(V,XCR,0.1262890E+00,0.1360539E+00)
+ GSET(V,XCR,0.1217290E+00,0.1407869E+00)
+ GSET(V,XCR,0.1177930E+00,0.1464869E+00)
+ GSET(V,XCR,0.1143519E+00,0.1528400E+00)
+ GSET(V,XCR,0.1107720E+00,0.1594409E+00)
+ GSET(V,XCR,0.1070470E+00,0.1658959E+00)
+ GSET(V,XCR,0.1031989E+00,0.1718689E+00)
+ GSET(V,XCR,0.9920394E-01,0.1771159E+00)
+ GSET(V,XCR,0.9492695E-01,0.1815059E+00)
+ GSET(V,XCR,0.9016997E-01,0.1850100E+00)
+ GSET(V,XCR,0.8467793E-01,0.1876850E+00)
+ GSET(V,XCR,0.7815695E-01,0.1896809E+00)
+ GSET(V,XCR,0.7060295E-01,0.1911619E+00)
+ GSET(V,XCR,0.6193697E-01,0.1922539E+00)
+ GSET(V,XCR,0.5162898E-01,0.1930479E+00)
+ GSET(V,XCR,0.3973098E-01,0.1936129E+00)
+ GSET(V,XCR,0.2688598E-01,0.1939889E+00)
+ GSET(V,XCR,0.1354700E-01,0.1942049E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NY1+NY2+NY3,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,XCR,0.1383730E+00,0.1257500E-01)
+ GSET(V,XCR,0.1398900E+00,0.2233299E-01)
+ GSET(V,XCR,0.1412150E+00,0.3450299E-01)
+ GSET(V,XCR,0.1422200E+00,0.4845098E-01)
+ GSET(V,XCR,0.1428069E+00,0.6333894E-01)
+ GSET(V,XCR,0.1428980E+00,0.7822198E-01)
+ GSET(V,XCR,0.1423020E+00,0.9217393E-01)
+ GSET(V,XCR,0.1409649E+00,0.1044310E+00)
+ GSET(V,XCR,0.1390780E+00,0.1144680E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
IS = NX1+NX2+NX3+NX4+NX51+NX52+1
IS
XCR = 0.3939997E+00
GSET(P,P1,XCR,0.0 ,ZB )
GSET(P,P2,XCR,0.1370609E+00,ZB )
GSET(P,P3,XCR,0.1370609E+00,0.5563997E-02)
GSET(P,P4,XCR,0.0 ,0.0 )
GSET(P,P6,XCR,0.0 ,0.2279699E+00)
GSET(P,P5,XCR,0.1378459E+00,0.1242099E+00)
GSET(P,P8,XCR,0.0 ,ZLL )
GSET(P,P7,XCR,0.1378459E+00,ZLL )
GSET(P,A1,XCR,YLL,ZB ); GSET(P,A2,XCR,YLL,0.5563997E-02)
GSET(P,A4,XCR,YLL,ZLL); GSET(P,A3,XCR,YLL,0.1242099E+00)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,XCR,0.1260449E+00,0.0000000E+00)
+ GSET(V,XCR,0.1328369E+00,0.1385999E-02)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NY1+NY2+NY3,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,XCR,0.1358209E+00,0.1301799E+00)
+ GSET(V,XCR,0.1334019E+00,0.1352659E+00)
+ GSET(V,XCR,0.1299639E+00,0.1405999E+00)
+ GSET(V,XCR,0.1256040E+00,0.1472659E+00)
+ GSET(V,XCR,0.1219530E+00,0.1557499E+00)
+ GSET(V,XCR,0.1187710E+00,0.1654729E+00)
+ GSET(V,XCR,0.1150770E+00,0.1757069E+00)
+ GSET(V,XCR,0.1109709E+00,0.1857529E+00)
+ GSET(V,XCR,0.1066650E+00,0.1950169E+00)
+ GSET(V,XCR,0.1023729E+00,0.2030810E+00)
+ GSET(V,XCR,0.9819293E-01,0.2097249E+00)
+ GSET(V,XCR,0.9420693E-01,0.2149149E+00)
+ GSET(V,XCR,0.9047496E-01,0.2187679E+00)
+ GSET(V,XCR,0.8653396E-01,0.2215660E+00)
+ GSET(V,XCR,0.8182395E-01,0.2236099E+00)
+ GSET(V,XCR,0.7545197E-01,0.2251129E+00)
+ GSET(V,XCR,0.6569797E-01,0.2262189E+00)
+ GSET(V,XCR,0.5213699E-01,0.2270179E+00)
+ GSET(V,XCR,0.3600499E-01,0.2275569E+00)
+ GSET(V,XCR,0.1835299E-01,0.2278680E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NY1+NY2+NY3,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,XCR,0.1389579E+00,0.1258099E-01)
+ GSET(V,XCR,0.1404929E+00,0.2242599E-01)
+ GSET(V,XCR,0.1418279E+00,0.3475999E-01)
+ GSET(V,XCR,0.1428300E+00,0.4892999E-01)
+ GSET(V,XCR,0.1434100E+00,0.6404495E-01)
+ GSET(V,XCR,0.1435080E+00,0.7910597E-01)
+ GSET(V,XCR,0.1429200E+00,0.9325296E-01)
+ GSET(V,XCR,0.1416050E+00,0.1058009E+00)
+ GSET(V,XCR,0.1397870E+00,0.1162319E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
IS = NX1+NX2+NX3+NX4+NX51+NX52+NX53+1
IS
XCR = 0.4779997E+00
GSET(P,P1,XCR,0.0 ,ZB )
GSET(P,P2,XCR,0.1369129E+00,ZB )
GSET(P,P3,XCR,0.1369129E+00,0.5718995E-02)
GSET(P,P4,XCR,0.0 ,0.0 )
GSET(P,P6,XCR,0.0 ,0.2350439E+00)
GSET(P,P5,XCR,0.1376989E+00,0.1245729E+00)
GSET(P,P8,XCR,0.0 ,ZLL )
GSET(P,P7,XCR,0.1376989E+00,ZLL )
GSET(P,A1,XCR,YLL,ZB ); GSET(P,A2,XCR,YLL,0.5718995E-02)
GSET(P,A4,XCR,YLL,ZLL); GSET(P,A3,XCR,YLL,0.1245729E+00)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,XCR,0.1259090E+00,0.0000000E+00)
+ GSET(V,XCR,0.1326939E+00,0.1426999E-02)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NY1+NY2+NY3,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,XCR,0.1356800E+00,0.1305280E+00)
+ GSET(V,XCR,0.1332679E+00,0.1356089E+00)
+ GSET(V,XCR,0.1298379E+00,0.1409469E+00)
+ GSET(V,XCR,0.1254830E+00,0.1476330E+00)
+ GSET(V,XCR,0.1218380E+00,0.1561899E+00)
+ GSET(V,XCR,0.1186600E+00,0.1660849E+00)
+ GSET(V,XCR,0.1149700E+00,0.1766289E+00)
+ GSET(V,XCR,0.1108660E+00,0.1871469E+00)
+ GSET(V,XCR,0.1065630E+00,0.1970379E+00)
+ GSET(V,XCR,0.1022739E+00,0.2058449E+00)
+ GSET(V,XCR,0.9809995E-01,0.2132789E+00)
+ GSET(V,XCR,0.9412098E-01,0.2192329E+00)
+ GSET(V,XCR,0.9040093E-01,0.2237549E+00)
+ GSET(V,XCR,0.8647597E-01,0.2271059E+00)
+ GSET(V,XCR,0.8178294E-01,0.2295949E+00)
+ GSET(V,XCR,0.7542896E-01,0.2314500E+00)
+ GSET(V,XCR,0.6568795E-01,0.2328309E+00)
+ GSET(V,XCR,0.5213398E-01,0.2338359E+00)
+ GSET(V,XCR,0.3600499E-01,0.2345179E+00)
+ GSET(V,XCR,0.1835398E-01,0.2349139E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NY1+NY2+NY3,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,XCR,0.1388089E+00,0.1289700E-01)
+ GSET(V,XCR,0.1403419E+00,0.2290899E-01)
+ GSET(V,XCR,0.1416759E+00,0.3537298E-01)
+ GSET(V,XCR,0.1426759E+00,0.4960499E-01)
+ GSET(V,XCR,0.1432559E+00,0.6471097E-01)
+ GSET(V,XCR,0.1433539E+00,0.7971197E-01)
+ GSET(V,XCR,0.1427659E+00,0.9377795E-01)
+ GSET(V,XCR,0.1414530E+00,0.1062509E+00)
+ GSET(V,XCR,0.1396369E+00,0.1166269E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
IS = NX1+NX2+NX3+NX4+NX51+NX52+NX53+NX54+1
IS
GSET(P,P1,X8,0.0 ,ZB )
GSET(P,P2,X8,0.1351689E+00,ZB )
GSET(P,P3,X7,0.1351689E+00,0.5522996E-02)
GSET(P,P4,X7,0.0 ,0.0 )
GSET(P,P6,X7,0.0 ,Z7T )
GSET(P,P5,X7,0.1358309E+00,0.1203899E+00)
GSET(P,P8,X7,0.0 ,ZLL )
GSET(P,P7,X7,0.1358309E+00,ZLL )
GSET(P,A1,X8,YLL,ZB ); GSET(P,A2,X7,YLL,0.5522996E-02)
GSET(P,A4,X7,YLL,ZLL); GSET(P,A3,X7,YLL,0.1203899E+00);
GSET(P,B1,X8,YW1,ZB ); GSET(P,B2,X8,YW2,ZB )
GSET(P,B4,X7,YW1,0.0); GSET(P,B3,X7,YW2,0.0)
CALL STLN3
** Wheels.
GSET(P,PA,XW2-DXWL,YW1 ,DZWL)
GSET(P,PB,XW2-DXWL,YW2 ,DZWL)
GSET(P,PC,XW2-DXWL,0.1351689E+00,DZWL)
GSET(L,B14,B1,B4,NZ1,1,ARC,PA)
GSET(L,B23,B2,B3,NZ1,1,ARC,PB)
GSET(L,L23,P2,P3,NZ1,1,ARC,PC)
VIEW(I,IS)
GSET(V,CV2,S,P5)
+ GSET(V,X7,0.1332809E+00,0.1261519E+00)
+ GSET(V,X7,0.1294640E+00,0.1310669E+00)
+ GSET(V,X7,0.1247749E+00,0.1362290E+00)
+ GSET(V,X7,0.1191189E+00,0.1426899E+00)
+ GSET(V,X7,0.1147220E+00,0.1509479E+00)
+ GSET(V,X7,0.1114730E+00,0.1604729E+00)
+ GSET(V,X7,0.1078690E+00,0.1705889E+00)
+ GSET(V,X7,0.1039389E+00,0.1806359E+00)
+ GSET(V,X7,0.9982294E-01,0.1900359E+00)
+ GSET(V,X7,0.9565794E-01,0.1983550E+00)
+ GSET(V,X7,0.9155494E-01,0.2053350E+00)
+ GSET(V,X7,0.8762497E-01,0.2108889E+00)
+ GSET(V,X7,0.8396494E-01,0.2150819E+00)
+ GSET(V,X7,0.8022892E-01,0.2181759E+00)
+ GSET(V,X7,0.7596493E-01,0.2204639E+00)
+ GSET(V,X7,0.7026696E-01,0.2221630E+00)
+ GSET(V,X7,0.6129199E-01,0.2234249E+00)
+ GSET(V,X7,0.4855798E-01,0.2243419E+00)
+ GSET(V,X7,0.3335899E-01,0.2249639E+00)
+ GSET(V,X7,0.1686399E-01,0.2253239E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X7,0.1367589E+00,0.1242100E-01)
+ GSET(V,X7,0.1381409E+00,0.2205699E-01)
+ GSET(V,X7,0.1394029E+00,0.3407399E-01)
+ GSET(V,X7,0.1403559E+00,0.4781799E-01)
+ GSET(V,X7,0.1408539E+00,0.6242698E-01)
+ GSET(V,X7,0.1408179E+00,0.7694793E-01)
+ GSET(V,X7,0.1401510E+00,0.9056896E-01)
+ GSET(V,X7,0.1389239E+00,0.1026499E+00)
+ GSET(V,X7,0.1374009E+00,0.1126980E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM3
GSET(I)
** Section 6.
IS = IS+NX61
IS
XCR = 0.6699997E+00
GSET(P,P1,XCR,0.0 ,ZB )
GSET(P,P2,XCR,0.1335509E+00,ZB )
GSET(P,P3,XCR,0.1335509E+00,0.7587995E-02)
GSET(P,P4,XCR,0.0 ,0.2268999E-02)
GSET(P,P6,XCR,0.0 ,0.2205199E+00)
GSET(P,P5,XCR,0.1342170E+00,0.1179889E+00)
GSET(P,P8,XCR,0.0 ,ZLL )
GSET(P,P7,XCR,0.1342170E+00,ZLL )
GSET(P,A1,XCR,YLL,ZB ); GSET(P,A2,XCR,YLL,0.7587995E-02)
GSET(P,A4,XCR,YLL,ZLL); GSET(P,A3,XCR,YLL,0.1179889E+00)
GSET(P,B1,XCR,YW1,ZB ); GSET(P,B2,XCR,YW2,ZB )
GSET(P,B4,XCR,YW1,0.2268999E-02); GSET(P,B3,XCR,YW2,0.2268999E-02)
CALL STLN3
GSET(V,CV2,S,P5)
+ GSET(V,XCR,0.1311889E+00,0.1235160E+00)
+ GSET(V,XCR,0.1264479E+00,0.1282319E+00)
+ GSET(V,XCR,0.1208810E+00,0.1331859E+00)
+ GSET(V,XCR,0.1144350E+00,0.1393909E+00)
+ GSET(V,XCR,0.1096359E+00,0.1473320E+00)
+ GSET(V,XCR,0.1063560E+00,0.1565139E+00)
+ GSET(V,XCR,0.1028119E+00,0.1662979E+00)
+ GSET(V,XCR,0.9899694E-01,0.1760579E+00)
+ GSET(V,XCR,0.9500897E-01,0.1852369E+00)
+ GSET(V,XCR,0.9094197E-01,0.1934099E+00)
+ GSET(V,XCR,0.8690894E-01,0.2003099E+00)
+ GSET(V,XCR,0.8303994E-01,0.2058380E+00)
+ GSET(V,XCR,0.7944793E-01,0.2100359E+00)
+ GSET(V,XCR,0.7585794E-01,0.2131469E+00)
+ GSET(V,XCR,0.7188994E-01,0.2154589E+00)
+ GSET(V,XCR,0.6663996E-01,0.2171819E+00)
+ GSET(V,XCR,0.5822299E-01,0.2184640E+00)
+ GSET(V,XCR,0.4613198E-01,0.2193980E+00)
+ GSET(V,XCR,0.3165098E-01,0.2200319E+00)
+ GSET(V,XCR,0.1596100E-01,0.2203999E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,XCR,0.1351349E+00,0.1426399E-01)
+ GSET(V,XCR,0.1365299E+00,0.2357199E-01)
+ GSET(V,XCR,0.1378049E+00,0.3515599E-01)
+ GSET(V,XCR,0.1387509E+00,0.4837998E-01)
+ GSET(V,XCR,0.1392059E+00,0.6241198E-01)
+ GSET(V,XCR,0.1390979E+00,0.7634294E-01)
+ GSET(V,XCR,0.1383849E+00,0.8940297E-01)
+ GSET(V,XCR,0.1371779E+00,0.1009820E+00)
+ GSET(V,XCR,0.1357300E+00,0.1106139E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM3
GSET(I)
IS = IS+NX62
IS
GSET(P,P1,X10,0.0 ,ZB )
GSET(P,P2,X10,0.1265749E+00,ZB )
GSET(P,P3,X11,0.1265749E+00,0.1855699E-01)
GSET(P,P4,X11,0.0 ,Z11B )
GSET(P,P6,X11,0.0 ,Z11T )
GSET(P,P5,X11,0.1272919E+00,0.1130329E+00)
GSET(P,P8,X11,0.0 ,ZLL )
GSET(P,P7,X11,0.1272919E+00,ZLL )
GSET(P,A1,X10,YLL,ZB ); GSET(P,A2,X11,YLL,0.1855699E-01)
GSET(P,A4,X11,YLL,ZLL); GSET(P,A3,X11,YLL,0.1130329E+00)
GSET(P,B1,X10,YW1,ZB ); GSET(P,B2,X10,YW2,ZB )
GSET(P,B4,X11,YW1,Z11B); GSET(P,B3,X11,YW2,Z11B)
CALL STLN3
** Wheels.
GSET(P,PA,XW2+DXWL,YW1 ,DZWL)
GSET(P,PB,XW2+DXWL,YW2 ,DZWL)
GSET(P,PC,XW2+DXWL,0.1265749E+00,DZWL)
GSET(L,B14,B1,B4,NZ1,1,ARC,PA)
GSET(L,B23,B2,B3,NZ1,1,ARC,PB)
GSET(L,L23,P2,P3,NZ1,1,ARC,PC)
VIEW(I,IS)
GSET(V,CV2,S,P5)
+ GSET(V,X11,0.1230149E+00,0.1177300E+00)
+ GSET(V,X11,0.1155699E+00,0.1217449E+00)
+ GSET(V,X11,0.1074849E+00,0.1259720E+00)
+ GSET(V,X11,0.9951192E-01,0.1312789E+00)
+ GSET(V,X11,0.9414595E-01,0.1381209E+00)
+ GSET(V,X11,0.9093493E-01,0.1461269E+00)
+ GSET(V,X11,0.8764791E-01,0.1547980E+00)
+ GSET(V,X11,0.8421695E-01,0.1636249E+00)
+ GSET(V,X11,0.8065993E-01,0.1721280E+00)
+ GSET(V,X11,0.7698894E-01,0.1799020E+00)
+ GSET(V,X11,0.7331795E-01,0.1866469E+00)
+ GSET(V,X11,0.6979495E-01,0.1921920E+00)
+ GSET(V,X11,0.6655496E-01,0.1965029E+00)
+ GSET(V,X11,0.6346595E-01,0.1997609E+00)
+ GSET(V,X11,0.6029798E-01,0.2022189E+00)
+ GSET(V,X11,0.5622198E-01,0.2040739E+00)
+ GSET(V,X11,0.4940899E-01,0.2054679E+00)
+ GSET(V,X11,0.3933099E-01,0.2064899E+00)
+ GSET(V,X11,0.2709799E-01,0.2071880E+00)
+ GSET(V,X11,0.1372100E-01,0.2075949E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X11,0.1282299E+00,0.2443599E-01)
+ GSET(V,X11,0.1296809E+00,0.3257098E-01)
+ GSET(V,X11,0.1309810E+00,0.4260598E-01)
+ GSET(V,X11,0.1318970E+00,0.5396598E-01)
+ GSET(V,X11,0.1322650E+00,0.6593496E-01)
+ GSET(V,X11,0.1320279E+00,0.7776093E-01)
+ GSET(V,X11,0.1312469E+00,0.8881897E-01)
+ GSET(V,X11,0.1300809E+00,0.9861797E-01)
+ GSET(V,X11,0.1287529E+00,0.1067770E+00)
+ GSET(V,X11,0.1272919E+00,0.1130329E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM3
GSET(I)
** Section 7.
IS = NX1+NX2+NX3+NX4+NX5+NX6+NX7+1
IS
GSET(P,P1,X12,0.0 ,ZB )
GSET(P,P2,X12,0.1097680E+00,ZB )
GSET(P,P3,X12,0.1097680E+00,0.2730499E-01)
GSET(P,P4,X12,0.0 ,Z12B )
GSET(P,P6,X12,0.0 ,Z12T )
GSET(P,P5,X12,0.1104640E+00,0.1080009E+00)
GSET(P,P8,X12,0.0 ,ZLL )
GSET(P,P7,X12,0.1104640E+00,ZLL )
GSET(P,A1,X12,YLL,ZB ); GSET(P,A2,X12,YLL,0.2730499E-01)
GSET(P,A4,X12,YLL,ZLL); GSET(P,A3,X12,YLL,0.1080009E+00)
CALL STLN2
GSET(V,CV1,S,P4)
+ GSET(V,X12,0.1005369E+00,0.2322198E-01)
+ GSET(V,X12,0.1062729E+00,0.2424398E-01)
GSET(V,CV1,E,P3)
GSET(L,L43,P4,P3,NTT,1.0CRVCV1)
GSET(V,CV2,S,P5)
+ GSET(V,X12,0.1056160E+00,0.1120099E+00)
+ GSET(V,X12,0.9631395E-01,0.1154410E+00)
+ GSET(V,X12,0.8661193E-01,0.1190559E+00)
+ GSET(V,X12,0.7843494E-01,0.1236029E+00)
+ GSET(V,X12,0.7340598E-01,0.1294869E+00)
+ GSET(V,X12,0.7062298E-01,0.1364160E+00)
+ GSET(V,X12,0.6788194E-01,0.1439819E+00)
+ GSET(V,X12,0.6509393E-01,0.1517599E+00)
+ GSET(V,X12,0.6223498E-01,0.1593369E+00)
+ GSET(V,X12,0.5927899E-01,0.1663460E+00)
+ GSET(V,X12,0.5632098E-01,0.1724979E+00)
+ GSET(V,X12,0.5349598E-01,0.1776089E+00)
+ GSET(V,X12,0.5092099E-01,0.1816189E+00)
+ GSET(V,X12,0.4854798E-01,0.1846719E+00)
+ GSET(V,X12,0.4624798E-01,0.1869879E+00)
+ GSET(V,X12,0.4335699E-01,0.1887439E+00)
+ GSET(V,X12,0.3836899E-01,0.1900679E+00)
+ GSET(V,X12,0.3081199E-01,0.1910419E+00)
+ GSET(V,X12,0.2147898E-01,0.1917079E+00)
+ GSET(V,X12,0.1103500E-01,0.1920959E+00)
GSET(V,CV2,E,P6)
GSET(L,L65,P5,P6,NTT,:PWYZ:CRVCV2)
GSET(V,CV3,S,P3)
+ GSET(V,X12,0.1113449E+00,0.3238299E-01)
+ GSET(V,X12,0.1127070E+00,0.3938499E-01)
+ GSET(V,X12,0.1138999E+00,0.4799398E-01)
+ GSET(V,X12,0.1147100E+00,0.5770498E-01)
+ GSET(V,X12,0.1149949E+00,0.6790996E-01)
+ GSET(V,X12,0.1147189E+00,0.7797694E-01)
+ GSET(V,X12,0.1139799E+00,0.8738297E-01)
+ GSET(V,X12,0.1129400E+00,0.9571993E-01)
+ GSET(V,X12,0.1117859E+00,0.1026670E+00)
GSET(V,CV3,E,P5)
GSET(L,L35,P3,P5,NZ2,S:PWZY:CRVCV3)
CALL STFM2
GSET(I)
** Interpolation.
IIL=NX1+1
GSET(C,I1,F,I:IIL:,1,NY,1,NZ,+,X2-X1,0.0,0.0,INC,1.0)
IIF=IIL; IIL = IIF+NX2
GSET(T,I:IIL:,F,I:IIF:,0.6)
IIF=IIL; IIL = IIF+NX31
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX32
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX33
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX4
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX51
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX52
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX53
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX54
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX61
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX62
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX7
GSET(T,I:IIL:,F,I:IIF:)
IIF=IIL; IIL = IIF+NX8
GSET(C,I:IIL:,F,I:IIF:,1,NY,1,NZ,+,X13-X12,0.0,0.0,INC,:PWZ:)
DUMPC(AAA)
ELSE
READCO(AAA)
DOMAIN(1,NX+1,1,NY+1,NZ1+NZ2+1,NZ+1)
SLIDH=T
LIK =T; LJK=F; LIJ=F; MSWP =3; MAGIC(L)
DOMAIN(1,NX+1,NY1+NY2+NY3+1,NY+1,1,NZ+1)
SLIDH=F; SLIDN=T
LIJ =T; LIK=F; LJK=F; MSWP =3; MAGIC(L)
DOMAIN(1,NX1+NX2,1,NY+1,1,NZ+1)
SLIDN=T; SLIDS=T; SLIDH=T; SLIDL=T
LIJ =F; LIK=F; LJK=T; MSWP =3; MAGIC(L)
ENDIF
GROUP 7. Variables stored, solved & named
STORE(PRPS)
SOLVE(P1,U1,V1,W1)
SOLUTN(P1,Y,Y,Y,N,N,N)
IF(LTURB) THEN
+ TURMOD(KEMODL)
+ FIINIT(KE)= TKEIN; FIINIT(EP)= EPSIN
ENDIF
GROUP 8. Terms (in differential equations) & devices
GROUP 9. Properties of the medium (or media)
ENUL = UIN*LCAR/REYNO
RHO1 = 1.189
GROUP 10. Inter-phase-transfer processes and properties
GROUP 11. Initialization of variable or porosity fields
INIADD = F
** Define car body.
PATCH(CAR,INIVAL,NX1+NX2+1,NX-NX8,1,NY1+NY2+NY3,NZ1+1,NZ1+NZ2,1,1)
COVAL(CAR,PRPS,0.0,199.0)
IIF = NX1+NX2+NX3+1; IIL=IIF+NX4-1
PATCH(WL1,INIVAL,IIF,IIL,NY1+1,NY1+NY2,1,NZ1,1,1)
COVAL(WL1,PRPS,0.0,199.0)
IIF = NX1+NX2+NX3+NX4+NX5+1; IIL=IIF+NX6-1
PATCH(WL2,INIVAL,IIF,IIL,NY1+1,NY1+NY2,1,NZ1,1,1)
COVAL(WL2,PRPS,0.0,199.0)
IF(LVISC) THEN
*** Boundary conditions.
** Car body.
PATCH(WH,LWALL,NX1+NX2+1,NX-NX8,1,NY1+NY2+NY3,NZ1+NZ2,NZ1+NZ2,1,1)
PATCH(WL,HWALL,NX1+NX2+1,NX-NX8,1,NY1+NY2+NY3,NZ1 ,NZ1 ,1,1)
PATCH(WS,SWALL,NX1+NX2+1,NX-NX8,1,NY1+NY2+NY3,NZ1+1,NZ1+NZ2,1,1)
PATCH(WE,EWALL,NX1+NX2,NX1+NX2 ,1,NY1+NY2+NY3,NZ1+1,NZ1+NZ2,1,1)
PATCH(WW,WWALL,NX-NX8+1,NX-NX8+1,1,NY1+NY2+NY3,NZ1+1,NZ1+NZ2,1,1)
COVAL(WH,U1,:CF:,0.0);COVAL(WH,V1,:CF:,0.0)
COVAL(WL,U1,:CF:,0.0);COVAL(WL,V1,:CF:,0.0)
COVAL(WS,U1,:CF:,0.0); COVAL(WS,W1,:CF:,0.0)
+ COVAL(WE,V1,:CF:,0.0);COVAL(WE,W1,:CF:,0.0)
+ COVAL(WW,V1,:CF:,0.0);COVAL(WW,W1,:CF:,0.0)
** Wheels.
** Front.
IIF= NX1+NX2+NX3+1; IIL= IIF+NX4-1
PATCH(S1,SWALL,IIF,IIL,NY1+NY2+1,NY1+NY2+1,1,NZ1,1,1)
PATCH(N1,NWALL,IIF,IIL,NY1 ,NY1 ,1,NZ1,1,1)
IIF= IIF-1
PATCH(E1,EWALL,IIF ,IIF ,NY1+1,NY1+NY2,1,NZ1,1,1)
PATCH(W1,WWALL,IIL+1,IIL+1,NY1+1,NY1+NY2,1,NZ1,1,1)
** Rear.
IIF= NX1+NX2+NX3+NX4+NX5+1; IIL = IIF+NX6-1
PATCH(S2,SWALL,IIF,IIL,NY1+NY2+1,NY1+NY2+1,1,NZ1,1,1)
PATCH(N2,NWALL,IIF,IIL,NY1 ,NY1 ,1,NZ1,1,1)
IIF= IIF-1
PATCH(E2,EWALL,IIF ,IIF ,NY1+1,NY1+NY2,1,NZ1,1,1)
PATCH(W2,WWALL,IIL+1,IIL+1,NY1+1,NY1+NY2,1,NZ1,1,1)
+ COVAL(E1,V1,:CF:,0.0);COVAL(E1,W1,:CF:,0.0)
+ COVAL(W1,V1,:CF:,0.0);COVAL(W1,W1,:CF:,0.0)
+COVAL(S1,U1,:CF:,0.0); COVAL(S1,W1,:CF:,0.0)
+COVAL(N1,U1,:CF:,0.0); COVAL(N1,W1,:CF:,0.0)
+ COVAL(E2,V1,:CF:,0.0);COVAL(E2,W1,:CF:,0.0)
+ COVAL(W2,V1,:CF:,0.0);COVAL(W2,W1,:CF:,0.0)
+COVAL(S2,U1,:CF:,0.0); COVAL(S2,W1,:CF:,0.0)
+COVAL(N2,U1,:CF:,0.0); COVAL(N2,W1,:CF:,0.0)
** Floor.
PATCH(ZF,LWALL,1,NX,1,NY,1,1,1,1)
COVAL(ZF,U1,:CF:,0.0);COVAL(WH,V1,:CF:,0.0)
ENDIF
FIINIT(U1) = UIN; FIINIT(V1)=0.0; FIINIT(W1)=0.0
GROUP 12. Unused
GROUP 13. Boundary conditions and special sources
** Inlet.
INLET(INL,WEST,1,1,1,NY,1,NZ,1,1)
VALUE(INL,P1,RHO1*UIN)
IF(LTURB) THEN
+ VALUE(INL,KE,TKEIN); VALUE(INL,EP,EPSIN)
ENDIF
VALUE(INL,U1,UIN)
VALUE(INL,V1,0.0)
VALUE(INL,W1,0.0)
** Outlet.
PATCH(OUT1,EAST,NX,NX,1,NY,1,NZ,1,1)
COVAL(OUT1,P1,1.E3,0.0)
PATCH(OUT2,HIGH,1,NX,1,NY,NZ,NZ,1,1)
COVAL(OUT2,P1,1.E3,0.0)
PATCH(OUT3,NORTH,1,NX,NY,NY,1,NZ,1,1)
COVAL(OUT3,P1,1.E3,0.0)
GROUP 14. Downstream pressure for PARAB=.TRUE.
GROUP 15. Termination of sweeps
LSWEEP = 500; TSTSWP = -1
GROUP 16. Termination of iterations
SELREF = T; RESFAC = 1.E-3
GROUP 17. Under-relaxation devices
DIFCUT = 0.0
DTHYD = 10.*HCAR/((NZ1+NZ2)*UIN)
RELAX(P1,LINRLX,0.5)
RELAX(U1,FALSDT,DTHYD)
RELAX(V1,FALSDT,DTHYD)
RELAX(W1,FALSDT,DTHYD)
IF(LTURB) THEN
+ RELAX(KE,linrlx,0.25); RELAX(EP,linrlx,0.25)
+ KELIN = 3
ENDIF
GROUP 19. Data communicated by satellite to GROUND
GCV=T
GROUP 21. Print-out of variables
OUTPUT(P1 ,Y,N,N,Y,Y,Y)
OUTPUT(U1 ,Y,N,N,Y,Y,Y)
OUTPUT(V1 ,Y,N,N,Y,Y,Y)
OUTPUT(W1 ,Y,N,N,Y,Y,Y)
GROUP 22. Spot-value print-out
IXMON=2; IYMON=NY/2; IZMON=NZ/2
GROUP 23. Field print-out and plot control
TSTSWP = -1
IF(DISTIL) THEN
+ CALL TEST:IRUNN:
ENDIF
ENDMAIN !
SUBROUTINE STLN
+ GSET(L,L21,P2,P1,NTT,1.0); GSET(L,L43,P4,P3,NTT,1.0)
+ GSET(L,L65,P6,P5,NTT,1.0); GSET(L,L78,P8,P7,NTT,1.0)
+ GSET(L,L14,P1,P4,NZ1,1.0); GSET(L,L23,P2,P3,NZ1,1.0)
+ GSET(L,L46,P4,P6,NZ2,1.0); GSET(L,L35,P3,P5,NZ2,1.0)
+ GSET(L,L68,P6,P8,NZ3,1.0); GSET(L,L57,P5,P7,NZ3,1.0)
ENDSUB
SUBROUTINE STFM
+ GSET(F,F1,P1,-,P2,-,P3,-,P4,-); GSET(M,F1,+I+K,IS,JS, 1)
+ GSET(F,F2,P4,-,P3,-,P5,-,P6,-); GSET(M,F2,+I+K,IS,JS, NZ1+1)
+ GSET(F,F3,P6,-,P5,-,P7,-,P8,-); GSET(M,F3,+I+K,IS,JS,NZ1+NZ2+1)
ENDSUB
SUBROUTINE STLN2
+ NTT = NY1+NY2+NY3
+ GSET(L,L21,P2,P1,NTT, 1.0 ); GSET(L,L43,P4,P3,NTT, 1.0 )
+ GSET(L,L65,P6,P5,NTT, 1.0 ); GSET(L,L78,P8,P7,NTT, 1.0 )
+ GSET(L,L14,P1,P4,NZ1, S1.0 ); GSET(L,L23,P2,P3,NZ1, S1.0 )
+ GSET(L,L46,P4,P6,NZ2, 1.0 ); GSET(L,L35,P3,P5,NZ2,S:PWZY:)
+ GSET(L,L68,P6,P8,NZ3, :PWZ:); GSET(L,L57,P5,P7,NZ3, :PWZ:)
+ GSET(L,A21,P2,A1,NY4, :PWY:); GSET(L,L32,P3,A2,NY4, :PWY:)
+ GSET(L,A53,P5,A3,NY4, :PWY:); GSET(L,A74,P7,A4,NY4, :PWY:)
+ GSET(L,A12,A1,A2,NZ1, S1.0 ); GSET(L,A23,A2,A3,NZ2, 1.0 )
+ GSET(L,A34,A3,A4,NZ3, :PWZ:)
ENDSUB
SUBROUTINE STFM2
+ GSET(F,F1,P1,-,P2,-,P3,-,P4,-);GSET(M,F1,+J+K,IS,1, 1)
+ GSET(F,F2,P4,-,P3,-,P5,-,P6,-);GSET(M,F2,+J+K,IS,1, NZ1+1)
+ GSET(F,F3,P6,-,P5,-,P7,-,P8,-);GSET(M,F3,+J+K,IS,1,NZ1+NZ2+1)
+ JS2=NY1+NY2+NY3+1
+ GSET(F,F4,P2,-,A1,-,A2,-,P3,-);GSET(M,F4,+J+K,IS,JS2, 1)
+ GSET(F,F5,P3,-,A2,-,A3,-,P5,-);GSET(M,F5,+J+K,IS,JS2, NZ1+1)
+ GSET(F,F6,P5,-,A3,-,A4,-,P7,-);GSET(M,F6,+J+K,IS,JS2,NZ1+NZ2+1)
ENDSUB
SUBROUTINE STLN3
+ NTT = NY1+NY2+NY3
+ GSET(L,B11,P1,B1,NY1, 1.0 ); GSET(L,B44,P4,B4,NY1, 1.0 )
+ GSET(L,B43,B4,B3,NY2, 1.0 ); GSET(L,B12,B1,B2,NY2, 1.0 )
+ GSET(L,B33,B3,P3,NY3, 1.0 ); GSET(L,B22,B2,P2,NY3, 1.0 )
+ GSET(L,L65,P6,P5,NTT, 1.0 ); GSET(L,L78,P8,P7,NTT, 1.0 )
+ GSET(L,B14,B1,B4,NZ1, 1.0 ); GSET(L,B23,B2,B3,NZ1, 1.0 )
+ GSET(L,L14,P1,P4,NZ1, S1.0 ); GSET(L,L23,P2,P3,NZ1, S1.0 )
+ GSET(L,L46,P4,P6,NZ2, 1.0 ); GSET(L,L35,P3,P5,NZ2,S:PWZY:)
+ GSET(L,L68,P6,P8,NZ3, :PWZ:); GSET(L,L57,P5,P7,NZ3, :PWZ:)
+ GSET(L,A21,P2,A1,NY4, :PWY:); GSET(L,L32,P3,A2,NY4, :PWY:)
+ GSET(L,A53,P5,A3,NY4, :PWY:); GSET(L,A74,P7,A4,NY4, :PWY:)
+ GSET(L,A12,A1,A2,NZ1, S1.0 ); GSET(L,A23,A2,A3,NZ2, 1.0 )
+ GSET(L,A34,A3,A4,NZ3, :PWZ:)
ENDSUB
SUBROUTINE STFM3
+GSET(F,F11,P1,-,B1,-,B4,-,P4,-); GSET(M,F11,+J+K,IS, 1,1)
+GSET(F,F12,B1,-,B2,-,B3,-,B4,-); GSET(M,F12,+J+K,IS, NY1+1,1)
+GSET(F,F13,B2,-,P2,-,P3,-,B3,-); GSET(M,F13,+J+K,IS,NY1+NY2+1,1)
+GSET(F,F2,P4,B4.B3,P3,-,P5,-,P6,-); GSET(M,F2,+J+K,IS,1, NZ1+1)
+GSET(F,F3,P6,- ,P5,-,P7,-,P8,-); GSET(M,F3,+J+K,IS,1,NZ2+NZ1+1)
+JS2=NY1+NY2+NY3+1
+GSET(F,F4,P2,-,A1,-,A2,-,P3,-); GSET(M,F4,+J+K,IS,JS2, 1)
+GSET(F,F5,P3,-,A2,-,A3,-,P5,-); GSET(M,F5,+J+K,IS,JS2, NZ1+1)
+GSET(F,F6,P5,-,A3,-,A4,-,P7,-); GSET(M,F6,+J+K,IS,JS2,NZ1+NZ2+1)
ENDSUB
SUBROUTINE TEST1
ENDSUB
LIBREF = 577
STOP