]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
7ef50f50 | 5 | * Revision 1.4 2000/09/18 10:02:35 hristov |
6 | * Makefile added to PDF8 | |
7 | * | |
21886bb6 | 8 | * Revision 1.1.1.2 1996/10/30 08:32:14 cernlib |
9 | * Version 7.04 | |
10 | * | |
11 | * Revision 1.2 1996/10/29 09:30:56 pdflib | |
12 | * more LINNUX corrections in DATA statements | |
13 | * | |
14 | * Revision 1.1.1.1 1996/04/12 15:29:07 plothow | |
15 | * Version 7.01 | |
16 | * | |
17 | * | |
18 | #include "pdf/pilot.h" | |
19 | FUNCTION Ctq1Pd (ISET, IPRTN, XX, QQ, IRT) | |
20 | ||
21 | C Returns xf(x,Q) -- the momentum fraction distribution !! | |
22 | C Returns valence and sea rather than combined flavor distr. | |
23 | ||
24 | C Iset : PDF set label | |
25 | ||
26 | C Iprtn : Parton label: 2, 1 = d_ and u_ valence | |
27 | C 0 = gluon | |
28 | C -1, ... -6 = u, d, s, c, b, t sea quarks | |
29 | ||
30 | C XX : Bjorken-x | |
31 | C QQ : scale parameter "Q" | |
32 | C Irt : Return code | |
33 | ||
34 | C IMPLICIT DOUBLE PRECISION (A-H, O-Z) | |
35 | C+SEQ, IMPDP. | |
36 | PARAMETER (D0=0D0, D1=1D0, D2=2D0, D3=3D0, D4=4D0, D10=1D1) | |
37 | PARAMETER (NEX = 5, MXFL = 6, NPN = 3, NST = 30, NEXPT=20) | |
38 | PARAMETER (NST4 = NST*4) | |
39 | ||
7ef50f50 | 40 | #if defined(CERNLIB_DOUBLE) |
41 | DOUBLE PRECISION | |
42 | #endif | |
43 | #if defined(CERNLIB_SINGLE) | |
44 | REAL | |
45 | #endif | |
21886bb6 | 46 | + P012, P34, P5, PDFD |
47 | REAL PDF | |
48 | ||
49 | DIMENSION | |
50 | > IORD(NST), ISCH(NST), NQRK(NST),ALM(NST) | |
51 | > , VLM(4:6,NST), QMS(4:6, NST) | |
52 | > , XMN(NST), QMN(NST), QMX(NST), NEXP(NEXPT) | |
53 | > , MEX(NST), MPN(NST), EXPN(NEXPT, NST), EXPNOR(NEXPT) | |
54 | ||
55 | C CTEQ1M | |
56 | DATA | |
57 | > ISCH(1), IORD(1), NQRK(1), ALM(1) / 1, 2, 6, .152 / | |
58 | > ,(VLM(I,1), I=4,6) / .231, .152, .059 / | |
59 | > ,(QMS(I,1), I=4,6) / 1.50, 5.00, 180.0 / | |
60 | > ,XMN(1), QMN(1), QMX(1) / 1.E-5, 2.00, 1.E3 / | |
61 | > ,MEX(1), MPN(1), NEXP(1) / 5, 3, 8 / | |
62 | > ,(EXPN(I, 1), I=1,8) | |
63 | > / 0.989, 1.00, 1.02, 0.978, 1.10, 0.972, 0.987, 0.846 / | |
64 | C CTEQ1MS | |
65 | DATA | |
66 | > ISCH(2), IORD(2), NQRK(2), ALM(2) / 1, 2, 6, .152 / | |
67 | > ,(VLM(I,2), I=4,6) / .231, .152, .059 / | |
68 | > ,(QMS(I,2), I=4,6) / 1.50, 5.00, 180.0 / | |
69 | > ,XMN(2), QMN(2), QMX(2) / 1.E-5, 2.00, 1.E3 / | |
70 | > ,MEX(2), MPN(2), NEXP(2) / 5, 3, 8 / | |
71 | > ,(EXPN(I, 2), I=1,8 ) | |
72 | > / 0.989, 1.00, 1.02, 0.984, 1.05, 0.891, 0.923, 0.824 / | |
73 | C CTEQ1ML | |
74 | DATA | |
75 | > ISCH(3), IORD(3), NQRK(3), ALM(3) / 1, 2, 6, .220 / | |
76 | > ,(VLM(I,3), I=4,6) / .322, .220, .088 / | |
77 | > ,(QMS(I,3), I=4,6) / 1.50, 5.00, 180.0 / | |
78 | > ,XMN(3), QMN(3), QMX(3) / 1.E-5, 2.00, 1.E3 / | |
79 | > ,MEX(3), MPN(3), NEXP(3) / 5, 3, 8 / | |
80 | > ,(EXPN(I, 3), I=1,8 ) | |
81 | > / 0.985, 1.00, 1.01, 0.977, 1.07, 1.31, 1.19, 1.09 / | |
82 | ||
83 | C CTEQ1D | |
84 | DATA | |
85 | > ISCH(4), IORD(4), NQRK(4), ALM(4) / 2, 2, 6, .164 / | |
86 | > ,(VLM(I,4), I=4,6) / .247, .164, .064 / | |
87 | > ,(QMS(I,4), I=4,6) / 1.50, 5.00, 180.0 / | |
88 | > ,XMN(4), QMN(4), QMX(4) / 1.E-5, 2.00, 1.E3 / | |
89 | > ,MEX(4), MPN(4), NEXP(4) / 5, 3, 8 / | |
90 | > ,(EXPN(I, 4), I=1,8 ) | |
91 | > / 0.983, 1.00, 1.01, 0.975, 0.964, 1.23, 1.00, 1.12 / | |
92 | C CTEQ1L | |
93 | DATA | |
94 | > ISCH(5), IORD(5), NQRK(5), ALM(5) / 1, 1, 6, .125 / | |
95 | > ,(VLM(I,5), I=4,6) / .168, .125, .063 / | |
96 | > ,(QMS(I,5), I=4,6) / 1.50, 5.00, 180.0 / | |
97 | > ,XMN(5), QMN(5), QMX(5) / 1.E-5, 2.00, 1.E3 / | |
98 | > ,MEX(5), MPN(5), NEXP(5) / 5, 3, 8 / | |
99 | > ,(EXPN(I, 5), I=1,8 ) | |
100 | > / 0.982, 1.01, 1.00, 0.972, 0.840, 0.959, 0.930, 0.861 / | |
101 | ||
102 | DATA IST, LP, QSTO, ALN2 / 0, -10, 1.2345, 0.6931 / | |
103 | ||
104 | X = XX | |
105 | QQ = MAX(QQ,QMN(ISET)) | |
106 | IF(ISET.EQ.IST.AND.IPRTN.EQ.LP.AND.QSTO.EQ.QQ) GOTO 100 | |
107 | ||
108 | IRT = 0 | |
109 | IP = ABS(IPRTN) | |
110 | IF (IP.GE.5.AND.QQ.LE.QMS(IP, ISET)) THEN | |
111 | Ctq1Pd = 0.0 | |
112 | RETURN | |
113 | ENDIF | |
114 | ||
115 | C if heavy parton, different logarithmic scale | |
116 | ||
117 | IF(IP.GE.5) THEN | |
118 | QI = QMS(IP,ISET) | |
119 | ELSE | |
120 | QI = QMN (ISET) | |
121 | ENDIF | |
122 | ||
123 | ALAM = ALM (ISET) | |
124 | STA = LOG(QQ/ALAM) | |
125 | STB = LOG(QI/ALAM) | |
126 | ||
127 | C SBL = LOG(QQ/Alam) / LOG(Qi/Alam) | |
128 | SBL = STA/STB | |
129 | SB = LOG (SBL) | |
130 | SB2 = SB*SB | |
131 | SB3 = SB2*SB | |
132 | ||
133 | IFLV = 3 - IPRTN | |
134 | ||
135 | GOTO (1, 2, 3, 4, 5), ISET | |
136 | ||
137 | 1 GOTO(11,12,13,14,15,16,17,18,19)IFLV | |
138 | C ifl = 2 | |
139 | 11 A0=0.3636E+01*(1.0 + 0.3122E+00*SB+0.1396E+00*SB2+0.4251E+00*SB3) | |
140 | A1=0.6930E+00-.2574E-01*SB+0.1047E+00*SB2-.2794E-01*SB3 | |
141 | A2=0.3195E+01+0.4045E+00*SB-.3737E+00*SB2-.1677E+00*SB3 | |
142 | A3=0.1009E+00*(1.0 -.1784E+01*SB+0.6263E+00*SB2+0.7337E-01*SB3) | |
143 | $ -1.0 | |
144 | A4=0.2910E+00-.2793E+00*SB+0.6155E-01*SB2+0.5150E-02*SB3 | |
145 | A5=0.0000E+00+0.3185E+00*SB+0.1953E+00*SB2+0.4184E-01*SB3 | |
146 | GOTO 100 | |
147 | C ifl = 1 | |
148 | 12 A0=0.2851E+00*(1.0 + 0.3617E+00*SB-.4526E+00*SB2+0.5787E-01*SB3) | |
149 | A1=0.2690E+00+0.1104E-01*SB+0.1888E-01*SB2-.1031E-01*SB3 | |
150 | A2=0.3766E+01+0.7850E+00*SB-.3053E+00*SB2+0.1822E+00*SB3 | |
151 | A3=0.2865E+02*(1.0 -.9774E+00*SB+0.5958E+00*SB2-.1234E+00*SB3) | |
152 | $ -1.0 | |
153 | A4=0.8230E+00-.3612E+00*SB+0.5520E-01*SB2+0.1571E-01*SB3 | |
154 | A5=0.0000E+00+0.2145E-01*SB+0.2289E+00*SB2-.4947E-01*SB3 | |
155 | GOTO 100 | |
156 | C ifl = 0 | |
157 | 13 A0=0.2716E+01*(1.0 -.2092E+01*SB+0.1500E+01*SB2-.3703E+00*SB3) | |
158 | A1=-.3100E-01-.7963E+00*SB+0.1129E+01*SB2-.4191E+00*SB3 | |
159 | A2=0.8015E+01+0.1168E+01*SB-.1625E+01*SB2-.1130E+01*SB3 | |
160 | A3=0.4813E+02*(1.0 -.4951E+00*SB-.8715E+00*SB2+0.5893E+00*SB3) | |
161 | $ -1.0 | |
162 | A4=0.2773E+01-.6329E+00*SB-.1048E+01*SB2+0.1418E+00*SB3 | |
163 | A5=0.0000E+00+0.5048E+00*SB+0.2390E+01*SB2-.4159E+00*SB3 | |
164 | GOTO 100 | |
165 | C ifl = -1 | |
166 | 14 A0=0.3085E+00*(1.0 + 0.9422E+00*SB-.2606E+01*SB2+0.1364E+01*SB3) | |
167 | A1=0.5000E-02-.6433E+00*SB+0.4980E+00*SB2-.1780E+00*SB3 | |
168 | A2=0.7490E+01+0.9112E+00*SB-.2047E+01*SB2+0.1456E+01*SB3 | |
169 | A3=0.1145E-01*(1.0 + 0.4610E+01*SB+0.1699E+01*SB2+0.1296E+00*SB3) | |
170 | $ -1.0 | |
171 | A4=0.6030E+00-.8081E+00*SB+0.9410E+00*SB2-.4458E+00*SB3 | |
172 | A5=0.0000E+00-.1736E+01*SB+0.2863E+01*SB2-.1268E+01*SB3 | |
173 | GOTO 100 | |
174 | C ifl = -2 | |
175 | 15 A0=0.1324E+00*(1.0 -.1050E+01*SB+0.4844E+00*SB2-.1043E+00*SB3) | |
176 | A1=-.1580E+00+0.1672E+00*SB-.4100E+00*SB2+0.1793E+00*SB3 | |
177 | A2=0.8559E+01-.7351E-01*SB+0.5898E+00*SB2-.2655E+00*SB3 | |
178 | A3=0.2378E+02*(1.0 -.1108E+00*SB-.1646E-01*SB2+0.1129E-01*SB3) | |
179 | $ -1.0 | |
180 | A4=0.1477E+01+0.3312E-01*SB-.2191E+00*SB2+0.9588E-01*SB3 | |
181 | A5=0.0000E+00+0.1850E+01*SB-.1481E+01*SB2+0.6222E+00*SB3 | |
182 | GOTO 100 | |
183 | C ifl = -3 | |
184 | 16 A0=0.3208E+00*(1.0 -.4755E+00*SB-.4003E+00*SB2+0.2300E+00*SB3) | |
185 | A1=-.3200E-01-.3357E+00*SB+0.3222E-01*SB2+0.5011E-01*SB3 | |
186 | A2=0.1164E+02+0.1048E+01*SB-.1097E+01*SB2-.4431E+00*SB3 | |
187 | A3=0.5065E+02*(1.0 + 0.2484E+00*SB-.9235E+00*SB2+0.1935E+00*SB3) | |
188 | $ -1.0 | |
189 | A4=0.3300E+01-.6785E+00*SB+0.5337E+00*SB2-.4035E+00*SB3 | |
190 | A5=0.0000E+00-.2496E+00*SB+0.3903E+00*SB2+0.1392E+00*SB3 | |
191 | GOTO 100 | |
192 | C ifl = -4 | |
193 | 17 A0=0.7967E-06*(1.0 + 0.1587E+01*SB+0.1812E+02*SB2-.1333E+02*SB3) | |
194 | $ *SQRT(STA - STB) | |
195 | A1=0.1096E+01-.1236E+01*SB+0.1014E+02*SB2+0.1940E+01*SB3 | |
196 | A2=0.4366E+00+0.1197E+02*SB-.5471E+00*SB2-.5427E+01*SB3 | |
197 | A3=0.4650E+03*(1.0 + 0.1310E+02*SB-.1918E+02*SB2+0.6791E+01*SB3) | |
198 | $ -1.0 | |
199 | A4=-.8486E+00+0.7457E+00*SB-.1083E+02*SB2-.1210E+01*SB3 | |
200 | A5=0.3494E+01-.3511E+01*SB-.1766E+01*SB2+0.3442E+01*SB3 | |
201 | GOTO 100 | |
202 | C ifl = -5 | |
203 | 18 A0=0.1713E-03*(1.0 + 0.2562E+02*SB-.2988E+02*SB2+0.4798E+01*SB3) | |
204 | $ *SQRT(STA - STB) | |
205 | A1=-.5276E-01+0.4105E+00*SB-.1079E+01*SB2+0.6278E+00*SB3 | |
206 | A2=0.4515E+01+0.8369E+01*SB-.1192E+02*SB2+0.3403E+01*SB3 | |
207 | A3=0.1756E+01*(1.0 + 0.1325E+02*SB-.2997E+02*SB2+0.1758E+02*SB3) | |
208 | $ -1.0 | |
209 | A4=0.3557E-01+0.4159E+01*SB-.6947E+01*SB2+0.2982E+01*SB3 | |
210 | A5=0.2551E+01+0.2168E+01*SB-.5119E+01*SB2+0.3739E+01*SB3 | |
211 | GOTO 100 | |
212 | C ifl = -6 | |
213 | 19 A0=0.7510E-04*(1.0 + 0.2836E+02*SB-.3000E+02*SB2-.2979E+02*SB3) | |
214 | $ *SQRT(STA - STB) | |
215 | A1=-.1855E+00+0.4543E+00*SB-.1448E+01*SB2+0.2009E-01*SB3 | |
216 | A2=0.6775E+01-.4210E+01*SB-.1221E+01*SB2+0.1199E+02*SB3 | |
217 | A3=0.1070E+01*(1.0 + 0.8356E+01*SB-.2992E+02*SB2+0.2433E+02*SB3) | |
218 | $ -1.0 | |
219 | A4=-.4601E-01+0.4248E+01*SB-.1736E+01*SB2+0.1187E+02*SB3 | |
220 | A5=0.2771E+01+0.1382E+01*SB-.4797E+01*SB2+0.1273E+01*SB3 | |
221 | GOTO 100 | |
222 | ||
223 | ||
224 | 2 GOTO(21,22,23,24,25,26,27,28,29)IFLV | |
225 | C CTEQ1MS | |
226 | C ifl = 2 | |
227 | 21 A0=0.1828E+01*(1.0 -.8698E+00*SB+0.2906E+00*SB2-.2003E-01*SB3) | |
228 | A1=0.6060E+00+0.8595E-01*SB-.4934E-01*SB2+0.2221E-01*SB3 | |
229 | A2=0.3454E+01-.3115E+00*SB+0.1321E+01*SB2-.3490E+00*SB3 | |
230 | A3=0.2616E+00*(1.0 -.1670E+01*SB+0.2333E+01*SB2+0.7730E-01*SB3) | |
231 | $ -1.0 | |
232 | A4=0.8920E+00-.8500E-02*SB+0.4960E+00*SB2-.4045E-01*SB3 | |
233 | A5=0.0000E+00+0.1091E+01*SB-.1613E+00*SB2+0.3773E-01*SB3 | |
234 | GOTO 100 | |
235 | C ifl = 1 | |
236 | 22 A0=0.2885E+00*(1.0 + 0.3388E+00*SB-.4550E+00*SB2+0.6005E-01*SB3) | |
237 | A1=0.2730E+00+0.1198E-01*SB+0.1880E-01*SB2-.1077E-01*SB3 | |
238 | A2=0.3736E+01+0.7687E+00*SB-.2731E+00*SB2+0.1638E+00*SB3 | |
239 | A3=0.2741E+02*(1.0 -.9585E+00*SB+0.5925E+00*SB2-.1239E+00*SB3) | |
240 | $ -1.0 | |
241 | A4=0.8040E+00-.3546E+00*SB+0.6123E-01*SB2+0.1086E-01*SB3 | |
242 | A5=0.0000E+00+0.4277E-01*SB+0.2187E+00*SB2-.4646E-01*SB3 | |
243 | GOTO 100 | |
244 | C ifl = 0 | |
245 | 23 A0=0.8416E-01*(1.0 -.1996E+01*SB+0.1903E+01*SB2-.6722E+00*SB3) | |
246 | A1=-.4790E+00-.5459E+00*SB+0.1638E+01*SB2-.4342E+00*SB3 | |
247 | A2=0.5071E+01+0.1470E+01*SB-.2401E+01*SB2+0.1273E+01*SB3 | |
248 | A3=0.2847E+02*(1.0 + 0.1124E+00*SB-.1338E+01*SB2+0.7115E+00*SB3) | |
249 | $ -1.0 | |
250 | A4=0.4990E+00-.7208E+00*SB+0.3333E-03*SB2-.2354E+00*SB3 | |
251 | A5=0.0000E+00-.4480E+00*SB+0.3720E+01*SB2-.1838E+01*SB3 | |
252 | GOTO 100 | |
253 | C ifl = -1 | |
254 | 24 A0=0.4378E+00*(1.0 -.1244E+01*SB+0.3278E+01*SB2-.2098E+01*SB3) | |
255 | A1=0.3500E-01-.1298E+01*SB+0.1229E+01*SB2-.3665E+00*SB3 | |
256 | A2=0.6781E+01+0.4078E+01*SB-.9711E+00*SB2-.1536E+01*SB3 | |
257 | A3=0.1527E-03*(1.0 + 0.1430E+02*SB+0.3000E+02*SB2+0.2771E+02*SB3) | |
258 | $ -1.0 | |
259 | A4=0.3060E+00+0.1011E+01*SB-.2045E+01*SB2+0.9422E+00*SB3 | |
260 | A5=0.0000E+00-.3205E+01*SB+0.2683E+01*SB2-.1746E+00*SB3 | |
261 | GOTO 100 | |
262 | C ifl = -2 | |
263 | 25 A0=0.7413E-01*(1.0 + 0.1291E+01*SB-.2667E+01*SB2+0.1076E+01*SB3) | |
264 | A1=-.2730E+00-.1206E+00*SB+0.1828E+00*SB2-.1001E+00*SB3 | |
265 | A2=0.7719E+01+0.1537E+01*SB-.6410E+00*SB2-.3920E-01*SB3 | |
266 | A3=0.1799E+02*(1.0 -.1334E+01*SB+0.1916E+01*SB2-.8878E+00*SB3) | |
267 | $ -1.0 | |
268 | A4=0.1167E+01-.9176E-01*SB+0.5132E+00*SB2-.3460E+00*SB3 | |
269 | A5=0.0000E+00-.5023E+00*SB+0.1951E+01*SB2-.8427E+00*SB3 | |
270 | GOTO 100 | |
271 | C ifl = -3 | |
272 | 26 A0=0.6551E+00*(1.0 -.5968E-01*SB+0.5621E-02*SB2-.2074E+00*SB3) | |
273 | A1=0.2800E-01-.1138E+01*SB+0.1178E+01*SB2-.4425E+00*SB3 | |
274 | A2=0.7553E+01+0.3996E+01*SB-.4448E+01*SB2+0.1673E+01*SB3 | |
275 | A3=0.9264E-01*(1.0 -.1760E+01*SB+0.1634E+01*SB2-.4067E+00*SB3) | |
276 | $ -1.0 | |
277 | A4=0.1970E+00+0.5256E+00*SB-.9775E+00*SB2+0.4488E+00*SB3 | |
278 | A5=0.0000E+00-.3668E+01*SB+0.4757E+01*SB2-.1717E+01*SB3 | |
279 | GOTO 100 | |
280 | C ifl = -4 | |
281 | 27 A0=0.1486E-03*(1.0 + 0.2107E+01*SB-.1056E+02*SB2+0.1403E+02*SB3) | |
282 | $ * SQRT(STA - STB) | |
283 | A1=0.2115E+00-.1702E+01*SB+0.2571E+01*SB2-.1177E+01*SB3 | |
284 | A2=0.3533E+01+0.1367E+01*SB-.3397E+01*SB2+0.6260E+01*SB3 | |
285 | A3=0.1096E+02*(1.0 + 0.9213E+01*SB-.2020E+02*SB2+0.1084E+02*SB3) | |
286 | $ -1.0 | |
287 | A4=0.7041E+00-.7236E+00*SB+0.2766E-01*SB2+0.7352E+00*SB3 | |
288 | A5=0.3904E+01-.4398E+01*SB+0.7056E+01*SB2-.3722E+01*SB3 | |
289 | GOTO 100 | |
290 | C ifl = -5 | |
291 | 28 A0=0.1201E-03*(1.0 + 0.5408E+01*SB-.1489E+02*SB2+0.1667E+02*SB3) | |
292 | $ * SQRT(STA - STB) | |
293 | A1=0.1420E-01-.1525E+01*SB+0.2408E+01*SB2-.1154E+01*SB3 | |
294 | A2=0.4254E+01+0.2836E+01*SB-.6018E+00*SB2+0.4133E+00*SB3 | |
295 | A3=0.5696E+01*(1.0 + 0.9451E+01*SB-.2029E+02*SB2+0.1033E+02*SB3) | |
296 | $ -1.0 | |
297 | A4=0.4775E+00-.6695E+00*SB+0.2747E+00*SB2-.1051E+00*SB3 | |
298 | A5=0.3330E+01-.5133E+01*SB+0.6921E+01*SB2-.3283E+01*SB3 | |
299 | GOTO 100 | |
300 | C ifl = -6 | |
301 | 29 A0=0.7697E-04*(1.0 + 0.2801E+02*SB-.1901E+02*SB2-.2880E+02*SB3) | |
302 | $ *SQRT(STA - STB) | |
303 | A1=-.2249E+00+0.4432E+00*SB-.1454E+01*SB2+0.3509E-01*SB3 | |
304 | A2=0.6642E+01-.2702E+01*SB+0.8229E+01*SB2+0.8243E+01*SB3 | |
305 | A3=0.1146E+01*(1.0 + 0.8104E+01*SB-.2998E+02*SB2+0.2812E+02*SB3) | |
306 | $ -1.0 | |
307 | A4=-.6421E-01+0.4246E+01*SB-.2908E+01*SB2+0.9686E-02*SB3 | |
308 | A5=0.2606E+01+0.1261E+01*SB-.4933E+01*SB2+0.3476E+00*SB3 | |
309 | GOTO 100 | |
310 | ||
311 | ||
312 | 3 GOTO(31,32,33,34,35,36,37,38,39)IFLV | |
313 | C CTEQ1ML | |
314 | C ifl = 2 | |
315 | 31 A0=0.3777E+01*(1.0 + 0.6986E+00*SB-.20655E+01*SB2+.10334E+01*SB3) | |
316 | A1=0.7100E+00+.2880E-01*SB-.7930E-01*SB2+0.5600E-01*SB3 | |
317 | A2=0.3259E+01+0.1508E+01*SB-.3932E+01*SB2+0.20613E+01*SB3 | |
318 | A3=0.1304E+00*(1.0 -.2016E+00*SB-.30015E+01*SB2+0.19118E+01*SB3) | |
319 | $ -1.0 | |
320 | A4=0.2890E+00-0.4311E+00*SB+0.7387E+00*SB2-.3697E+00*SB3 | |
321 | A5=0.0000E+00+0.4320E+00*SB+0.2449E+00*SB2-0.6670E-01*SB3 | |
322 | GOTO 100 | |
323 | C ifl = 1 | |
324 | 32 A0=0.2780E+00*(1.0 + 0.4355E+00*SB-0.4584E+00*SB2+0.4390E-01*SB3) | |
325 | A1=0.2760E+00+0.1420E-01*SB+0.1480E-01*SB2-.9800E-02*SB3 | |
326 | A2=0.3710E+01+0.8250E+00*SB-.3581E+00*SB2+0.1978E+00*SB3 | |
327 | A3=0.2928E+02*(1.0 -.10154E+01*SB+0.6037E+00*SB2-.1175E+00*SB3) | |
328 | $ -1.0 | |
329 | A4=0.8070E+00-.3575E+00*SB+0.4920E-01*SB2+0.1584E-01*SB3 | |
330 | A5=0.0000E+00+0.1860E-01*SB+0.2080E+00*SB2-.450E-01*SB3 | |
331 | GOTO 100 | |
332 | C ifl = 0 | |
333 | 33 A0=0.2924E+01*(1.0 -.18916E+01*SB+0.1191E+01*SB2-.2492E+00*SB3) | |
334 | A1=0.0000E+00-.9167E+00*SB+0.11147E+01*SB2-.3329E+00*SB3 | |
335 | A2=0.8529E+01+0.7080E+00*SB-.11345E+01*SB2-.10563E+01*SB3 | |
336 | A3=0.1420E+03*(1.0 -.15346E+01*SB+0.7261E+00*SB2-.5730E-01*SB3) | |
337 | $ -1.0 | |
338 | A4=0.3396E+01-.11541E+01*SB-.8834E+00*SB2+0.2430E+00*SB3 | |
339 | A5=0.0000E+00+0.1645E+00*SB+0.19041E+01*SB2+0.1474E+00*SB3 | |
340 | GOTO 100 | |
341 | C ifl = -1 | |
342 | 34 A0=0.3471E+00*(1.0- 0.1753E+00*SB-.9189E+00*SB2+0.6211E+00*SB3) | |
343 | A1=0.1900E-01-.4579E+00*SB+0.2112E+00*SB2-.6180E-01*SB3 | |
344 | A2=0.7301E+01-.17308E+01*SB+.13666E+01*SB2-.6400E-02*SB3 | |
345 | A3=0.1853E-04*(1.0 -.18260E+02*SB-.2872E+02*SB2-.23456E+02*SB3) | |
346 | $ -1.0 | |
347 | A4=0.4400E+00-.4672E+00*SB+0.6532E+00*SB2-.3222E+00*SB3 | |
348 | A5=0.0000E+00-.4679E+00*SB+0.10741E+01*SB2-.5663E+00*SB3 | |
349 | GOTO 100 | |
350 | C ifl = -2 | |
351 | 35 A0=0.1702E+00*(1.0 -.1041E+01*SB+0.4064E+00*SB2-.5888E-01*SB3) | |
352 | A1=-.9300E-01-.4742E-01*SB-.1959E+00*SB2+0.1039E+00*SB3 | |
353 | A2=0.9119E+01-.7331E-01*SB+0.3506E+00*SB2-.2081E+00*SB3 | |
354 | A3=0.2981E+02*(1.0 -.1912E+00*SB-.8947E-02*SB2+0.8805E-02*SB3) | |
355 | $ -1.0 | |
356 | A4=0.1668E+01-.6678E-02*SB-.2894E+00*SB2+0.1221E+00*SB3 | |
357 | A5=0.0000E+00+0.1245E+01*SB-.7843E+00*SB2+0.3724E+00*SB3 | |
358 | GOTO 100 | |
359 | C ifl = -3 | |
360 | 36 A0=0.3910E+00*(1.0 -.1103E+01*SB+0.5383E+00*SB2-.1083E+00*SB3) | |
361 | A1=-.1400E-01-.2471E+00*SB-.8042E-01*SB2+0.7193E-01*SB3 | |
362 | A2=0.9812E+01-.4860E+01*SB+0.5958E+01*SB2-.2342E+01*SB3 | |
363 | A3=0.3749E+00*(1.0 -.3569E+01*SB+0.5456E+01*SB2-.2344E+01*SB3) | |
364 | $ -1.0 | |
365 | A4=0.4940E+00+0.2772E+00*SB-.2732E+00*SB2+0.6466E-01*SB3 | |
366 | A5=0.0000E+00+0.3927E+00*SB-.3216E+00*SB2+0.2164E+00*SB3 | |
367 | GOTO 100 | |
368 | C ifl = -4 | |
369 | 37 A0=0.3815E-02*(1.0 + 0.2039E+02*SB-.2834E+02*SB2+0.1070E+02*SB3) | |
370 | $ * SQRT(STA - STB) | |
371 | A1=-.2789E-01-.7345E-03*SB-.3251E+00*SB2+0.1946E+00*SB3 | |
372 | A2=0.3223E+01-.4268E+00*SB+0.4387E+01*SB2-.2401E+01*SB3 | |
373 | A3=0.3338E-01*(1.0 -.1163E+02*SB+0.2995E+02*SB2-.1471E+02*SB3) | |
374 | $ -1.0 | |
375 | A4=0.3646E+00-.5767E+00*SB+0.6088E+00*SB2-.2514E+00*SB3 | |
376 | A5=0.1200E+01+0.2178E+00*SB-.4230E+00*SB2+0.4739E+00*SB3 | |
377 | GOTO 100 | |
378 | C ifl = -5 | |
379 | 38 A0=0.1666E-02*(1.0 + 0.9518E+01*SB-.4715E+01*SB2-.1060E+01*SB3) | |
380 | $ * SQRT(STA - STB) | |
381 | A1=-.1231E+00+0.1656E+00*SB-.5219E+00*SB2+0.2750E+00*SB3 | |
382 | A2=0.3693E+01+0.4922E+01*SB-.1200E+02*SB2+0.7929E+01*SB3 | |
383 | A3=0.1778E+00*(1.0 + 0.3036E+01*SB-.1184E+02*SB2+0.7940E+01*SB3) | |
384 | $ -1.0 | |
385 | A4=0.5353E+00-.1401E+01*SB+0.1970E+01*SB2-.9405E+00*SB3 | |
386 | A5=0.1590E+01+0.1025E+01*SB-.2318E+01*SB2+0.1380E+01*SB3 | |
387 | GOTO 100 | |
388 | C ifl = -6 | |
389 | 39 A0=0.4319E-03*(1.0 + 0.1100E+02*SB-.9520E+00*SB2+0.1434E+02*SB3) | |
390 | $ * SQRT(STA - STB) | |
391 | A1=-.2512E+00+0.3554E+00*SB-.4120E+00*SB2+0.1328E+00*SB3 | |
392 | A2=0.4764E+01-.3513E+00*SB+0.1199E+02*SB2-.8290E+01*SB3 | |
393 | A3=0.8458E-01*(1.0 + 0.2618E+01*SB+0.4407E+01*SB2+0.2991E+02*SB3) | |
394 | $ -1.0 | |
395 | A4=0.3991E+00-.1363E+01*SB+0.1526E+01*SB2-.3179E+01*SB3 | |
396 | A5=0.1981E+01+0.1496E+01*SB-.1501E+01*SB2+0.3880E+01*SB3 | |
397 | GOTO 100 | |
398 | ||
399 | 4 GOTO(41,42,43,44,45,46,47,48,49)IFLV | |
400 | C CTEQ1D | |
401 | C ifl = 2 | |
402 | 41 A0=0.1634E+01*(1.0 -.8336E+00*SB+0.1640E+00*SB2+0.1530E+00*SB3) | |
403 | A1=0.5790E+00+0.8587E-01*SB-.6087E-01*SB2+0.1361E-01*SB3 | |
404 | A2=0.2839E+01+0.3720E+00*SB+0.5264E+00*SB2+0.3538E-01*SB3 | |
405 | A3=0.1095E+00*(1.0 -.4830E+00*SB+0.3708E+01*SB2-.6165E+00*SB3) | |
406 | $ -1.0 | |
407 | A4=0.8010E+00-.1432E+00*SB+0.1442E+01*SB2-.1286E+01*SB3 | |
408 | A5=0.0000E+00+0.1035E+01*SB-.5910E-01*SB2-.1982E+00*SB3 | |
409 | GOTO 100 | |
410 | C ifl = 1 | |
411 | 42 A0=0.3535E+00*(1.0 + 0.4352E+00*SB-.2095E+00*SB2-.8455E-02*SB3) | |
412 | A1=0.2660E+00-.4096E-03*SB+0.1502E-01*SB2-.1163E-01*SB3 | |
413 | A2=0.3514E+01+0.8219E+00*SB-.2330E+00*SB2+0.1055E+00*SB3 | |
414 | A3=0.2200E+02*(1.0 -.9716E+00*SB+0.4552E+00*SB2-.8202E-01*SB3) | |
415 | $ -1.0 | |
416 | A4=0.9000E+00-.3207E+00*SB-.4808E-01*SB2+0.3492E-01*SB3 | |
417 | A5=0.0000E+00-.6273E-01*SB+0.1497E+00*SB2-.5683E-01*SB3 | |
418 | GOTO 100 | |
419 | C ifl = 0 | |
420 | 43 A0=0.2743E+01*(1.0 -.2027E+01*SB+0.1517E+01*SB2-.4145E+00*SB3) | |
421 | A1=0.7000E-02-.9431E+00*SB+0.1231E+01*SB2-.4834E+00*SB3 | |
422 | A2=0.8200E+01+0.1827E+01*SB-.3453E+01*SB2+0.6763E+00*SB3 | |
423 | A3=0.4975E+02*(1.0 -.2329E+00*SB-.1245E+01*SB2+0.7194E+00*SB3) | |
424 | $ -1.0 | |
425 | A4=0.2387E+01-.4077E+00*SB-.5542E+00*SB2-.9677E-02*SB3 | |
426 | A5=0.0000E+00+0.2702E+00*SB+0.2389E+01*SB2-.8274E+00*SB3 | |
427 | GOTO 100 | |
428 | C ifl = -1 | |
429 | 44 A0=0.2015E+00*(1.0 -.2133E+00*SB-.6770E+00*SB2+0.5011E+00*SB3) | |
430 | A1=-.7700E-01-.7104E-01*SB-.3720E+00*SB2+0.2159E+00*SB3 | |
431 | A2=0.8008E+01-.2049E+01*SB+0.1800E+01*SB2-.4660E+00*SB3 | |
432 | A3=0.2923E-05*(1.0 + 0.2327E+02*SB+0.1500E+02*SB2+0.2633E+02*SB3) | |
433 | $ -1.0 | |
434 | A4=0.9020E+00-.9191E+00*SB+0.1104E+01*SB2-.5863E+00*SB3 | |
435 | A5=0.0000E+00+0.5840E+00*SB-.8720E+00*SB2+0.4234E+00*SB3 | |
436 | GOTO 100 | |
437 | C ifl = -2 | |
438 | 45 A0=0.9117E-01*(1.0 -.4089E+00*SB-.4361E+00*SB2+0.2512E+00*SB3) | |
439 | A1=-.2370E+00+0.2492E+00*SB-.3267E+00*SB2+0.1055E+00*SB3 | |
440 | A2=0.8447E+01+0.6009E+00*SB+0.1003E+01*SB2-.1287E+01*SB3 | |
441 | A3=0.3106E+02*(1.0 -.3901E-01*SB+0.1443E+00*SB2-.3433E+00*SB3) | |
442 | $ -1.0 | |
443 | A4=0.1629E+01+0.7855E-01*SB-.1573E+00*SB2-.8595E-01*SB3 | |
444 | A5=0.0000E+00+0.1558E+01*SB-.6295E+00*SB2+0.1847E+00*SB3 | |
445 | GOTO 100 | |
446 | C ifl = -3 | |
447 | 46 A0=0.3997E+00*(1.0 -.1046E+01*SB+0.6194E+00*SB2-.1342E+00*SB3) | |
448 | A1=0.2000E-02-.2544E+00*SB-.1958E+00*SB2+0.1458E+00*SB3 | |
449 | A2=0.9613E+01-.3919E+01*SB+0.9573E+01*SB2-.5623E+01*SB3 | |
450 | A3=0.3620E+00*(1.0 -.1858E+01*SB+0.8312E+01*SB2-.5900E+01*SB3) | |
451 | $ -1.0 | |
452 | A4=0.3840E+00+0.3572E+00*SB-.1191E+01*SB2+0.7310E+00*SB3 | |
453 | A5=0.0000E+00+0.3351E+00*SB-.7709E+00*SB2+0.4296E+00*SB3 | |
454 | GOTO 100 | |
455 | C ifl = -4 | |
456 | 47 A0=0.2156E-03*(1.0 + 0.2879E+02*SB-.2310E+02*SB2+0.9812E+01*SB3) | |
457 | $ * SQRT(STA - STB) | |
458 | A1=0.9086E-01-.1250E+00*SB-.7373E-01*SB2-.2201E-01*SB3 | |
459 | A2=0.3588E+01+0.4518E+01*SB-.8930E-01*SB2+0.9163E-02*SB3 | |
460 | A3=0.5216E+01*(1.0 + 0.5912E+00*SB-.4111E+00*SB2+0.7330E+00*SB3) | |
461 | $ -1.0 | |
462 | A4=0.3145E+00+0.1233E+01*SB-.7478E+00*SB2+0.4657E+00*SB3 | |
463 | A5=0.2723E+01-.4110E+00*SB+0.4868E-01*SB2-.3075E+00*SB3 | |
464 | GOTO 100 | |
465 | C ifl = -5 | |
466 | 48 A0=0.7476E-03*(1.0 + 0.1454E+02*SB-.2509E+02*SB2+0.1184E+02*SB3) | |
467 | $ * SQRT(STA - STB) | |
468 | A1=-.1955E-01-.1712E+00*SB-.1686E+00*SB2+0.2339E+00*SB3 | |
469 | A2=0.4616E+01-.6859E+00*SB-.3959E+01*SB2+0.5530E+01*SB3 | |
470 | A3=0.9881E+01*(1.0 -.1239E+02*SB+0.2721E+02*SB2-.1850E+02*SB3) | |
471 | $ -1.0 | |
472 | A4=0.1200E+02-.1133E+02*SB+0.8138E+01*SB2+0.1199E+02*SB3 | |
473 | A5=0.2226E+01-.5738E+00*SB+0.5239E+00*SB2+0.3825E+00*SB3 | |
474 | GOTO 100 | |
475 | C ifl = -6 | |
476 | 49 A0=0.8392E-06*(1.0 + 0.1844E+02*SB-.1110E+02*SB2-.2504E+02*SB3) | |
477 | $ * SQRT(STA - STB) | |
478 | A1=0.2127E+00-.5602E+00*SB+0.4777E+01*SB2-.1014E+02*SB3 | |
479 | A2=0.1229E+01+0.7495E+01*SB-.5024E+01*SB2-.1200E+02*SB3 | |
480 | A3=0.2868E+02*(1.0 + 0.7634E+01*SB-.2916E+02*SB2+0.2953E+02*SB3) | |
481 | $ -1.0 | |
482 | A4=0.5970E+00+0.1138E+01*SB-.1439E+01*SB2-.1966E+01*SB3 | |
483 | A5=0.6429E+01-.6673E+00*SB+0.7008E+01*SB2-.1157E+02*SB3 | |
484 | GOTO 100 | |
485 | ||
486 | 5 GOTO(51,52,53,54,55,56,57,58,59)IFLV | |
487 | C CTEQ1L | |
488 | C ifl = 2 | |
489 | 51 A0= 1.791*(1.0 -0.449*SB-0.445*SB2+ 0.401*SB3) | |
490 | A1= 0.608+ 0.069*SB+ 0.005*SB2-0.037*SB3 | |
491 | A2= 3.470-0.375*SB+ 2.267*SB2-1.261*SB3 | |
492 | A3= 0.315*(1.0 -2.628*SB+ 6.481*SB2-3.834*SB3)-1.0 | |
493 | A4= 1.007-0.732*SB+ 1.490*SB2-0.966*SB3 | |
494 | A5= 0.000+ 0.741*SB+ 0.563*SB2-0.525*SB3 | |
495 | GOTO 100 | |
496 | C ifl = 1 | |
497 | 52 A0= 0.513*(1.0 + 0.032*SB-0.120*SB2+ 0.013*SB3) | |
498 | A1= 0.276+ 0.052*SB+ 0.000*SB2-0.006*SB3 | |
499 | A2= 3.579+ 0.763*SB-0.135*SB2+ 0.083*SB3 | |
500 | A3= 17.993*(1.0 -0.725*SB+ 0.241*SB2-0.020*SB3)-1.0 | |
501 | A4= 1.120-0.357*SB+ 0.008*SB2+ 0.028*SB3 | |
502 | A5= 0.000+ 0.311*SB+ 0.029*SB2-0.010*SB3 | |
503 | GOTO 100 | |
504 | C ifl = 0 | |
505 | 53 A0= 2.710*(1.0 -1.773*SB+ 0.970*SB2-0.149*SB3) | |
506 | A1= -0.010-1.636*SB+ 2.087*SB2-0.637*SB3 | |
507 | A2= 7.174+ 2.102*SB-2.209*SB2-0.420*SB3 | |
508 | A3= 29.904*(1.0 -0.756*SB-0.506*SB2+ 0.605*SB3)-1.0 | |
509 | A4= 2.572-0.437*SB-0.968*SB2+ 0.243*SB3 | |
510 | A5= 0.000-1.776*SB+ 4.266*SB2-0.335*SB3 | |
511 | GOTO 100 | |
512 | C ifl = -1 | |
513 | 54 A0= 0.278*(1.0 - 1.022*SB+ 0.6457*SB2-0.1824*SB3) | |
514 | A1= 0.0862*SB-0.8657*SB2+ 0.4185*SB3 | |
515 | A2= 11.000-1.2809*SB+ 1.2516*SB2+0.061*SB3 | |
516 | A3= 37.338*(1.0 - 0.9404*SB+ 0.2517*SB2+0.1364*SB3)-1.0 | |
517 | A4= 1.960- 0.3385*SB-0.3422*SB2+0.3653*SB3 | |
518 | A5= 0.000+1.424*SB-2.7503*SB2+ 1.2226*SB3 | |
519 | GOTO 100 | |
520 | C ifl = -2 | |
521 | 55 A0= 0.154*(1.0 -0.659*SB+ 0.005*SB2+ 0.061*SB3) | |
522 | A1= -0.128+ 0.279*SB-0.786*SB2+ 0.363*SB3 | |
523 | A2= 8.649+ 0.071*SB+ 0.351*SB2-0.051*SB3 | |
524 | A3= 43.685*(1.0 -0.603*SB+ 0.037*SB2+ 0.134*SB3)-1.0 | |
525 | A4= 2.238-0.338*SB-0.199*SB2+ 0.157*SB3 | |
526 | A5= 0.000+ 1.681*SB-2.068*SB2+ 0.975*SB3 | |
527 | GOTO 100 | |
528 | C ifl = -3 | |
529 | 56 A0= 0.372*(1.0 -1.939*SB+ 1.504*SB2-0.440*SB3) | |
530 | A1= 0.009+ 0.610*SB-1.387*SB2+ 0.579*SB3 | |
531 | A2= 10.273-4.833*SB+ 6.583*SB2-2.633*SB3 | |
532 | A3= 0.160*(1.0 + 10.325*SB-2.027*SB2+ 1.571*SB3)-1.0 | |
533 | A4= 0.819-1.660*SB+ 1.845*SB2-0.829*SB3 | |
534 | A5= 0.000+ 3.558*SB-3.940*SB2+ 1.302*SB3 | |
535 | GOTO 100 | |
536 | C ifl = -4 | |
537 | 57 A0= (7.5242E-5)*(1.0+22.0905*SB+7.1209*SB2-8.303*SB3)* | |
538 | $ SQRT(STA - STB) | |
539 | A1= 0.125-0.3027*SB+0.1564*SB2-0.091*SB3 | |
540 | A2= 2.0388+1.2161*SB+11.5296*SB2-8.0659*SB3 | |
541 | A3= 14.849*(1.0 -2.556*SB+3.5268*SB2-1.6353*SB3)-1.0 | |
542 | A4= 0.3061-0.0901*SB+0.953*SB2-0.4871*SB3 | |
543 | A5= 2.7352+0.1811*SB-0.5167*SB2+0.0543*SB3 | |
544 | GOTO 100 | |
545 | C ifl = -5 | |
546 | 58 A0= (3.751E-4)*(1.0 + 21.5993*SB+3.1379*SB2-18.8328*SB3)* | |
547 | $ SQRT(STA - STB) | |
548 | A1= -0.0256-0.7717*SB+ 1.1499*SB2-0.5037*SB3 | |
549 | A2= 4.9241+4.0107*SB-4.7012*SB2+0.1097*SB3 | |
550 | A3= 2.842*(1.0 -2.2184*SB+ 2.0293*SB2-0.6907*SB3)-1.0 | |
551 | A4= -0.1352+ 0.8753*SB-1.2626*SB2+ 0.667*SB3 | |
552 | A5= 1.5627-0.4917*SB+ 1.5927*SB2-0.351*SB3 | |
553 | GOTO 100 | |
554 | C ifl = -6 | |
555 | 59 A0=(2.725E-4)*(1.0 + 18.8497*SB-26.5797*SB2-29.0774*SB3)* | |
556 | $ SQRT(STA - STB) | |
557 | A1= -0.2204-1.0048*SB+0.9415*SB2-0.4274*SB3 | |
558 | A2= 11.034-9.8362*SB-11.1034*SB2-9.1977*SB3 | |
559 | A3= 2.084*(1.0 -2.881*SB+1.2778*SB2-2.9328*SB3)-1.0 | |
560 | A4= -0.0872+ 0.200*SB-1.6187*SB2-1.6058*SB3 | |
561 | A5= 0.8684+4.7047*SB-1.4614*SB2-5.2309*SB3 | |
562 | GOTO 100 | |
563 | ||
564 | 100 CONTINUE | |
565 | C 100 CtqPdf = A0*(x**A1)*((1.-x)**A2)*(1.+A3*(x**A4)) | |
566 | C $ *(log(1.+1./x))**A5 | |
567 | P012 = A0*(X**A1)*((1.-X)**A2) | |
568 | P34 = (1.+A3*(X**A4)) | |
569 | P5 = (LOG(1.+1./X))**A5 | |
570 | C | |
571 | PDFD = P012 * P34 * P5 | |
572 | PDF = PDFD | |
573 | C | |
574 | Ctq1Pd = PDF | |
575 | IF(Ctq1Pd.LT.0.0) Ctq1Pd = 0.0 | |
576 | ||
577 | IST = ISET | |
578 | ||
579 | LP = IPRTN | |
580 | QSTO = QQ | |
581 | ||
582 | RETURN | |
583 | C ----------------------- | |
584 | ENTRY Wlamd1 (ISET, IORDER) | |
585 | ||
586 | IORDER = IORD (ISET) | |
587 | Wlamd1 = ALM (ISET) | |
588 | ||
589 | Ctq1Pd = 0.0 | |
590 | RETURN | |
591 | C ----------------------- | |
592 | ENTRY PrCtq1 | |
593 | > (ISET, IORDR, ISCHEME, MXFLV, | |
594 | > ALAM4, ALAM5, ALAM6, AMAS4, AMAS5, AMAS6, | |
595 | > XMIN, QINI, QMAX, EXPNOR) | |
596 | ||
597 | C Return QCD parameters and Fitting parameters | |
598 | C associated with parton distribution set Iset. | |
599 | C Iord : Order Of Fit | |
600 | C Ischeme : (0, 1, 2) for (LO, MS-bar-NLO, DIS-NLO) resp. | |
601 | C MxFlv : Maximum number of flavors included | |
602 | C Alam_i : i = 4,5,6 Effective lambda for i-flavors | |
603 | ||
604 | C Amas_i : i = 4,5,6 Mass parameter for flavor i | |
605 | C Xmin, Qini, Qmax : self explanary | |
606 | C ExpNor(I) : Normalization factor for the experimental data set used in | |
607 | C obtaining the best global fit for parton distributions Iset: | |
608 | C I = 1, 2, 3, 4, 5, 6, 7, 8 | |
609 | C BCDMS NMC90 NMC280 CCFR E605 WA70 E706 UA6 | |
610 | ||
611 | IORDR = IORD (ISET) | |
612 | ISCHEME= ISCH (ISET) | |
613 | MXFLV = NQRK (ISET) | |
614 | ||
615 | ALAM4 = VLM(4,ISET) | |
616 | ALAM5 = VLM(5,ISET) | |
617 | ALAM6 = VLM(6,ISET) | |
618 | ||
619 | AMAS4 = QMS(4,ISET) | |
620 | AMAS5 = QMS(5,ISET) | |
621 | AMAS6 = QMS(6,ISET) | |
622 | ||
623 | XMIN = XMN (ISET) | |
624 | QINI = QMN (ISET) | |
625 | QMAX = QMX (ISET) | |
626 | ||
627 | DO 101 IEXP = 1, NEXP(ISET) | |
628 | EXPNOR(IEXP) = EXPN(IEXP, ISET) | |
629 | 101 CONTINUE | |
630 | ||
631 | Ctq1Pd = 0.0 | |
632 | RETURN | |
633 | C ************************* | |
634 | END |