forked from lgcrego/Dynemol
-
Notifications
You must be signed in to change notification settings - Fork 0
/
abcpes-Coul.F
635 lines (635 loc) · 38.2 KB
/
abcpes-Coul.F
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
BLOCK DATA abpeGP5
implicit real*8 (a-h,o-z)
common / abcpesGP5 / abc(5), pes(5)
DATA (abc(i),i=1,5) / -9.06179845938664D-1,-5.384693101056831D-1,0
&.D0,5.384693101056831D-1,9.06179845938664D-1 /
DATA (pes(i),i=1,5) / 2.369268850561891D-1,4.786286704993665D-1,5.
&688888888888889D-1,4.786286704993665D-1,2.369268850561891D-1 /
end
SUBROUTINE abcpes5(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP5 / abc(5), pes(5)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,5
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end
!
!
!
BLOCK DATA abpeGP11
implicit real*8 (a-h,o-z)
common / abcpesGP11 / abc(11), pes(11)
DATA (abc(i),i=1,11) / -9.06179845938664D-1,-5.384693101056831D-1,
&0.D0,5.384693101056831D-1,9.06179845938664D-1,-9.840853600948425D-
&1,-7.541667265708492D-1,-2.796304131617832D-1,2.796304131617832D-1
&,7.541667265708492D-1,9.840853600948425D-1 /
DATA (pes(i),i=1,11) / 1.152333166224734D-1,2.410403392286476D-1,2
&.829874178574912D-1,2.410403392286476D-1,1.152333166224734D-1,4.25
&8203675108183D-2,1.868007965564927D-1,2.728498019125589D-1,2.72849
&8019125589D-1,1.868007965564927D-1,4.258203675108183D-2 /
end
SUBROUTINE abcpes11(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP11 / abc(11), pes(11)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,11
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end
!
!
!
BLOCK DATA abpeGP23
implicit real*8 (a-h,o-z)
common / abcpesGP23 / abc(23), pes(23)
DATA (abc(i),i=1,23) / -9.06179845938664D-1,-5.384693101056831D-1,
&0.D0,5.384693101056831D-1,9.06179845938664D-1,-9.840853600948425D-
&1,-7.541667265708492D-1,-2.796304131617832D-1,2.796304131617832D-1
&,7.541667265708492D-1,9.840853600948425D-1,-9.975454975868673D-1,-
&9.54373144738509D-1,-8.39137030805526D-1,-6.531199781049166D-1,-4.
&130159172460904D-1,-1.410830454016985D-1,1.410830454016985D-1,4.13
&0159172460904D-1,6.531199781049166D-1,8.39137030805526D-1,9.543731
&44738509D-1,9.975454975868673D-1 /
DATA (pes(i),i=1,23) / 5.767783624685351D-2,1.205322764830333D-1,1
&.415024762815157D-1,1.205322764830333D-1,5.767783624685351D-2,2.10
&1171242196147D-2,9.337876279217436D-2,1.36415530721228D-1,1.364155
&30721228D-1,9.337876279217436D-2,2.101171242196147D-2,6.7888161722
&1286D-3,3.876209361566127D-2,7.624831129039682D-2,1.08297714633535
&1D-1,1.298931290615022D-1,1.402425784206833D-1,1.402425784206833D-
&1,1.298931290615022D-1,1.082977146335351D-1,7.624831129039682D-2,3
&.876209361566127D-2,6.78881617221286D-3 /
end
SUBROUTINE abcpes23(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP23 / abc(23), pes(23)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,23
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end
!
!
!
BLOCK DATA abpeGP47
implicit real*8 (a-h,o-z)
common / abcpesGP47 / abc(47), pes(47)
DATA (abc(i),i=1,47) / -9.06179845938664D-1,-5.384693101056831D-1,
&0.D0,5.384693101056831D-1,9.06179845938664D-1,-9.840853600948425D-
&1,-7.541667265708492D-1,-2.796304131617832D-1,2.796304131617832D-1
&,7.541667265708492D-1,9.840853600948425D-1,-9.975454975868673D-1,-
&9.54373144738509D-1,-8.39137030805526D-1,-6.531199781049166D-1,-4.
&130159172460904D-1,-1.410830454016985D-1,1.410830454016985D-1,4.13
&0159172460904D-1,6.531199781049166D-1,8.39137030805526D-1,9.543731
&44738509D-1,9.975454975868673D-1,-9.996428940629623D-1,-9.92609474
&2582342D-1,-9.714532078050934D-1,-9.326448966646406D-1,-8.74982562
&3646846D-1,-7.987951932551814D-1,-7.055094496670671D-1,-5.97324442
&3227357D-1,-4.769125924857435D-1,-3.471380549365452D-1,-2.10834758
&239592D-1,-7.069888485151317D-2,7.069888485151317D-2,2.10834758239
&592D-1,3.471380549365452D-1,4.769125924857435D-1,5.973244423227357
&D-1,7.055094496670671D-1,7.987951932551814D-1,8.749825623646846D-1
&,9.326448966646406D-1,9.714532078050934D-1,9.926094742582342D-1,9.
&996428940629623D-1 /
DATA (pes(i),i=1,47) / 2.883895682580973D-2,6.026613906811125D-2,7
&.075123853871084D-2,6.026613906811125D-2,2.883895682580973D-2,1.05
&0800141129907D-2,4.668938456191255D-2,6.820776582310313D-2,6.82077
&6582310313D-2,4.668938456191255D-2,1.050800141129907D-2,3.36900409
&1470427D-3,1.938081126665124D-2,3.812414630481552D-2,5.41488558820
&8636D-2,6.494656395501522D-2,7.012128879802731D-2,7.01212887980273
&1D-2,6.494656395501522D-2,5.414885588208636D-2,3.812414630481552D-
&2,1.938081126665124D-2,3.369004091470427D-3,1.009319084223881D-3,6
&.626149655510141D-3,1.481224370075563D-2,2.409015843488085D-2,3.35
&4142415014526D-2,4.252442368225165D-2,5.057557853415348D-2,5.73843
&2328883686D-2,6.278691566576171D-2,6.675055504840345D-2,6.93282442
&2931844D-2,7.059412726810144D-2,7.059412726810144D-2,6.93282442293
&1844D-2,6.675055504840345D-2,6.278691566576171D-2,5.73843232888368
&6D-2,5.057557853415348D-2,4.252442368225165D-2,3.354142415014526D-
&2,2.409015843488085D-2,1.481224370075563D-2,6.626149655510141D-3,1
&.009319084223881D-3 /
end
SUBROUTINE abcpes47(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP47 / abc(47), pes(47)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,47
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end
!
!
!
BLOCK DATA abpeGP95
implicit real*8 (a-h,o-z)
common / abcpesGP95 / abc(95), pes(95)
DATA (abc(i),i=1,95) / -9.06179845938664D-1,-5.384693101056831D-1,
&0.D0,5.384693101056831D-1,9.06179845938664D-1,-9.840853600948425D-
&1,-7.541667265708492D-1,-2.796304131617832D-1,2.796304131617832D-1
&,7.541667265708492D-1,9.840853600948425D-1,-9.975454975868673D-1,-
&9.54373144738509D-1,-8.39137030805526D-1,-6.531199781049166D-1,-4.
&130159172460904D-1,-1.410830454016985D-1,1.410830454016985D-1,4.13
&0159172460904D-1,6.531199781049166D-1,8.39137030805526D-1,9.543731
&44738509D-1,9.975454975868673D-1,-9.996428940629623D-1,-9.92609474
&2582342D-1,-9.714532078050934D-1,-9.326448966646406D-1,-8.74982562
&3646846D-1,-7.987951932551814D-1,-7.055094496670671D-1,-5.97324442
&3227357D-1,-4.769125924857435D-1,-3.471380549365452D-1,-2.10834758
&239592D-1,-7.069888485151317D-2,7.069888485151317D-2,2.10834758239
&592D-1,3.471380549365452D-1,4.769125924857435D-1,5.973244423227357
&D-1,7.055094496670671D-1,7.987951932551814D-1,8.749825623646846D-1
&,9.326448966646406D-1,9.714532078050934D-1,9.926094742582342D-1,9.
&996428940629623D-1,-9.999498338042262D-1,-9.988910202139731D-1,-9.
&954852422297365D-1,-9.888331807398711D-1,-9.78307721609879D-1,-9.6
&34845653255284D-1,-9.440979512415834D-1,-9.200061936727053D-1,-8.9
&11692033084103D-1,-8.576327998942287D-1,-8.195162842820891D-1,-7.7
&70016931497028D-1,-7.303239512945866D-1,-6.797614383331393D-1,-6.2
&56266847329246D-1,-5.682571221004128D-1,-5.080060481051527D-1,-4.4
&52341960318211D-1,-3.803024612130273D-1,-3.1356635867478D-1,-2.453
&726209119229D-1,-1.760580133874071D-1,-1.059500582206001D-1,-3.536
&907741210191D-2,3.536907741210191D-2,1.059500582206001D-1,1.760580
&133874071D-1,2.453726209119229D-1,3.1356635867478D-1,3.80302461213
&0273D-1,4.452341960318211D-1,5.080060481051527D-1,5.68257122100412
&8D-1,6.256266847329246D-1,6.797614383331393D-1,7.303239512945866D-
&1,7.770016931497028D-1,8.195162842820891D-1,8.576327998942287D-1,8
&.911692033084103D-1,9.200061936727053D-1,9.440979512415834D-1,9.63
&4845653255284D-1,9.78307721609879D-1,9.888331807398711D-1,9.954852
&422297365D-1,9.988910202139731D-1,9.999498338042262D-1 /
DATA (pes(i),i=1,95) / 1.441947841293669D-2,3.013306953405563D-2,3
&.537561926935542D-2,3.013306953405563D-2,1.441947841293669D-2,5.25
&4000981537219D-3,2.334469228095641D-2,3.410388291155157D-2,3.41038
&8291155157D-2,2.334469228095641D-2,5.254000981537219D-3,1.68462244
&2460738D-3,9.690405635101188D-3,1.906207315240915D-2,2.70744279410
&432D-2,3.247328197750761D-2,3.506064439901366D-2,3.506064439901366
&D-2,3.247328197750761D-2,2.70744279410432D-2,1.906207315240915D-2,
&9.690405635101188D-3,1.684622442460738D-3,5.019096272361336D-4,3.3
&13069672920139D-3,7.406121831040587D-3,1.204507921722963D-2,1.6770
&71207506663D-2,2.126221184112543D-2,2.528778926707669D-2,2.8692161
&64441842D-2,3.139345783288085D-2,3.337527752420172D-2,3.4664122114
&65922D-2,3.529706363405072D-2,3.529706363405072D-2,3.4664122114659
&22D-2,3.337527752420172D-2,3.139345783288085D-2,2.869216164441842D
&-2,2.528778926707669D-2,2.126221184112543D-2,1.677071207506663D-2,
&1.204507921722963D-2,7.406121831040587D-3,3.313069672920139D-3,5.0
&19096272361336D-4,1.435403334315947D-4,1.026466183499721D-3,2.4526
&12875446336D-3,4.25132994730987D-3,6.309094660181704D-3,8.53592235
&3880551D-3,1.086231940055549D-2,1.323265211877946D-2,1.56004179943
&4175D-2,1.792595499520407D-2,2.017530939295059D-2,2.23196264066824
&1D-2,2.43348397098274D-2,2.620155177412264D-2,2.790500675618432D-2
&,2.943504346402941D-2,3.07859080875392D-2,3.195582266511717D-2,3.2
&94625744378213D-2,3.376094124020843D-2,3.440473951757565D-2,3.4882
&59232918837D-2,3.519869263654232D-2,3.535599212846248D-2,3.5355992
&12846248D-2,3.519869263654232D-2,3.488259232918837D-2,3.4404739517
&57565D-2,3.376094124020843D-2,3.294625744378213D-2,3.1955822665117
&17D-2,3.07859080875392D-2,2.943504346402941D-2,2.790500675618432D-
&2,2.620155177412264D-2,2.43348397098274D-2,2.231962640668241D-2,2.
&017530939295059D-2,1.792595499520407D-2,1.560041799434175D-2,1.323
&265211877946D-2,1.086231940055549D-2,8.535922353880551D-3,6.309094
&660181704D-3,4.25132994730987D-3,2.452612875446336D-3,1.0264661834
&99721D-3,1.435403334315947D-4 /
end
SUBROUTINE abcpes95(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP95 / abc(95), pes(95)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,95
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end
!
!
!
BLOCK DATA abpeGP191
implicit real*8 (a-h,o-z)
common / abcpesGP191 / abc(191), pes(191)
DATA (abc(i),i=1,150) / -9.06179845938664D-1,-5.384693101056831D-1
&,0.D0,5.384693101056831D-1,9.06179845938664D-1,-9.840853600948425D
&-1,-7.541667265708492D-1,-2.796304131617832D-1,2.796304131617832D-
&1,7.541667265708492D-1,9.840853600948425D-1,-9.975454975868673D-1,
&-9.54373144738509D-1,-8.39137030805526D-1,-6.531199781049166D-1,-4
&.130159172460904D-1,-1.410830454016985D-1,1.410830454016985D-1,4.1
&30159172460904D-1,6.531199781049166D-1,8.39137030805526D-1,9.54373
&144738509D-1,9.975454975868673D-1,-9.996428940629623D-1,-9.9260947
&42582342D-1,-9.714532078050934D-1,-9.326448966646406D-1,-8.7498256
&23646846D-1,-7.987951932551814D-1,-7.055094496670671D-1,-5.9732444
&23227357D-1,-4.769125924857435D-1,-3.471380549365452D-1,-2.1083475
&8239592D-1,-7.069888485151317D-2,7.069888485151317D-2,2.1083475823
&9592D-1,3.471380549365452D-1,4.769125924857435D-1,5.97324442322735
&7D-1,7.055094496670671D-1,7.987951932551814D-1,8.749825623646846D-
&1,9.326448966646406D-1,9.714532078050934D-1,9.926094742582342D-1,9
&.996428940629623D-1,-9.999498338042262D-1,-9.988910202139731D-1,-9
&.954852422297365D-1,-9.888331807398711D-1,-9.78307721609879D-1,-9.
&634845653255284D-1,-9.440979512415834D-1,-9.200061936727053D-1,-8.
&911692033084103D-1,-8.576327998942287D-1,-8.195162842820891D-1,-7.
&770016931497028D-1,-7.303239512945866D-1,-6.797614383331393D-1,-6.
&256266847329246D-1,-5.682571221004128D-1,-5.080060481051527D-1,-4.
&452341960318211D-1,-3.803024612130273D-1,-3.1356635867478D-1,-2.45
&3726209119229D-1,-1.760580133874071D-1,-1.059500582206001D-1,-3.53
&6907741210191D-2,3.536907741210191D-2,1.059500582206001D-1,1.76058
&0133874071D-1,2.453726209119229D-1,3.1356635867478D-1,3.8030246121
&30273D-1,4.452341960318211D-1,5.080060481051527D-1,5.6825712210041
&28D-1,6.256266847329246D-1,6.797614383331393D-1,7.303239512945866D
&-1,7.770016931497028D-1,8.195162842820891D-1,8.576327998942287D-1,
&8.911692033084103D-1,9.200061936727053D-1,9.440979512415834D-1,9.6
&34845653255284D-1,9.78307721609879D-1,9.888331807398711D-1,9.95485
&2422297365D-1,9.988910202139731D-1,9.999498338042262D-1,-9.9999308
&00957606D-1,-9.998413733868334D-1,-9.993326036260333D-1,-9.9830059
&04912557D-1,-9.966114125082292D-1,-9.941549532502683D-1,-9.9083864
&5014144D-1,-9.865846355418413D-1,-9.81328454820258D-1,-9.750176168
&805138D-1,-9.67610132633725D-1,-9.590731842806808D-1,-9.4938205432
&55846D-1,-9.385192846195623D-1,-9.264740062449535D-1,-9.1324138646
&75595D-1,-8.988221544560521D-1,-8.832221810638305D-1,-8.6645209709
&84496D-1,-8.485269399006613D-1,-8.294658210610491D-1,-8.0929160978
&71924D-1,-7.880306274725631D-1,-7.657123497864587D-1,-7.4236911329
&73484D-1,-7.180358243555749D-1,-6.927496687383048D-1,-6.6654982141
&33578D-1,-6.39477156702661D-1,-6.115739600970961D-1,-5.82883643952
&0516D-1,-5.534504702200102D-1,-5.233192841812306D-1,-4.92535263734
&6438D-1,-4.611436891258013D-1,-4.291897379448173D-1,-3.96718309776
&3532D-1,-3.637738840148685D-1,-3.304004131090205D-1,-2.96641251960
&2923D-1,-2.625391225146616D-1,-2.281361109327633D-1,-1.93473693299
&7303D-1,-1.585927848221899D-1,-1.23533806989227D-1,-8.833676729901
&243D-2,-5.304134682592874D-2,-1.768699197607406D-2,1.7686991976074
&06D-2,5.304134682592874D-2,8.833676729901243D-2,1.23533806989227D-
&1,1.585927848221899D-1,1.934736932997303D-1,2.281361109327633D-1 /
&
DATA (abc(i),i=151,191) / 2.625391225146616D-1,2.966412519602923D-
&1,3.304004131090205D-1,3.637738840148685D-1,3.967183097763532D-1,4
&.291897379448173D-1,4.611436891258013D-1,4.925352637346438D-1,5.23
&3192841812306D-1,5.534504702200102D-1,5.828836439520516D-1,6.11573
&9600970961D-1,6.39477156702661D-1,6.665498214133578D-1,6.927496687
&383048D-1,7.180358243555749D-1,7.423691132973484D-1,7.657123497864
&587D-1,7.880306274725631D-1,8.092916097871924D-1,8.294658210610491
&D-1,8.485269399006613D-1,8.664520970984496D-1,8.832221810638305D-1
&,8.988221544560521D-1,9.132413864675595D-1,9.264740062449535D-1,9.
&385192846195623D-1,9.493820543255846D-1,9.590731842806808D-1,9.676
&10132633725D-1,9.750176168805138D-1,9.81328454820258D-1,9.86584635
&5418413D-1,9.90838645014144D-1,9.941549532502683D-1,9.966114125082
&292D-1,9.983005904912557D-1,9.993326036260333D-1,9.998413733868334
&D-1,9.999930800957606D-1 /
DATA (pes(i),i=1,150) / 7.209739206468344D-3,1.506653476702782D-2,
&1.768780963467771D-2,1.506653476702782D-2,7.209739206468344D-3,2.6
&27000490768626D-3,1.16723461404782D-2,1.705194145577578D-2,1.70519
&4145577578D-2,1.16723461404782D-2,2.627000490768626D-3,8.423112289
&532929D-4,4.845202817550594D-3,9.531036576204576D-3,1.353721397052
&16D-2,1.623664098875381D-2,1.753032219950683D-2,1.753032219950683D
&-2,1.623664098875381D-2,1.35372139705216D-2,9.531036576204576D-3,4
&.845202817550594D-3,8.423112289532929D-4,2.50965635591339D-4,1.656
&534836468448D-3,3.703060915520294D-3,6.022539608614813D-3,8.385356
&037533313D-3,1.063110592056272D-2,1.264389463353834D-2,1.434608082
&220921D-2,1.569672891644043D-2,1.668763876210086D-2,1.733206105732
&961D-2,1.764853181702536D-2,1.764853181702536D-2,1.733206105732961
&D-2,1.668763876210086D-2,1.569672891644043D-2,1.434608082220921D-2
&,1.264389463353834D-2,1.063110592056272D-2,8.385356037533313D-3,6.
&022539608614813D-3,3.703060915520294D-3,1.656534836468448D-3,2.509
&65635591339D-4,7.143137772267442D-5,5.132328103896912D-4,1.2263064
&37485539D-3,2.125664973654593D-3,3.154547330090851D-3,4.2679611769
&40275D-3,5.431159700277743D-3,6.616326059389729D-3,7.8002089971708
&73D-3,8.962977497602036D-3,1.008765469647529D-2,1.115981320334121D
&-2,1.21674198549137D-2,1.310077588706132D-2,1.395250337809216D-2,1
&.47175217320147D-2,1.53929540437696D-2,1.597791133255859D-2,1.6473
&12872189106D-2,1.688047062010422D-2,1.720236975878782D-2,1.7441296
&16459418D-2,1.759934631827116D-2,1.767799606423124D-2,1.7677996064
&23124D-2,1.759934631827116D-2,1.744129616459418D-2,1.7202369758787
&82D-2,1.688047062010422D-2,1.647312872189106D-2,1.597791133255859D
&-2,1.53929540437696D-2,1.47175217320147D-2,1.395250337809216D-2,1.
&310077588706132D-2,1.21674198549137D-2,1.115981320334121D-2,1.0087
&65469647529D-2,8.962977497602036D-3,7.800208997170873D-3,6.6163260
&59389729D-3,5.431159700277743D-3,4.267961176940275D-3,3.1545473300
&90851D-3,2.125664973654593D-3,1.226306437485539D-3,5.1323281038969
&12D-4,7.143137772267442D-5,1.991706885703549D-5,1.496107785399853D
&-4,3.728569938879705D-4,6.702583391738669D-4,1.028038174940154D-3,
&1.436115964703238D-3,1.886670513946234D-3,2.372699521551761D-3,2.8
&87842289651325D-3,3.426483604750896D-3,3.983726580652995D-3,4.5552
&74633238511D-3,5.137304648633606D-3,5.726366236714079D-3,6.3193109
&4183439D-3,6.913244587476468D-3,7.505495016816082D-3,8.09358953789
&2668D-3,8.675238553145434D-3,9.248323378736119D-3,9.81088715268888
&2D-3,1.036112819269042D-2,1.089739537407941D-2,1.141818517380656D-
&2,1.192214003393026D-2,1.24080476741203D-2,1.28748409457053D-2,1.3
&32159778154144D-2,1.374754076526337D-2,1.415203582870152D-2,1.4534
&58959590346D-2,1.489484493432945D-2,1.523257435476985D-2,1.5547671
&02436904D-2,1.584013731978726D-2,1.611007104151424D-2,1.6357649620
&221D-2,1.658311285001439D-2,1.678674485546577D-2,1.6968856112407D-
&2,1.712976637354502D-2,1.726978928459762D-2,1.738921931409695D-2,1
&.748832137611585D-2,1.756732323277551D-2,1.762641046871624D-2,1.76
&6572358475625D-2,1.768535660983149D-2,1.768535660983149D-2,1.76657
&2358475625D-2,1.762641046871624D-2,1.756732323277551D-2,1.74883213
&7611585D-2,1.738921931409695D-2,1.726978928459762D-2 /
DATA (pes(i),i=151,191) / 1.712976637354502D-2,1.6968856112407D-2,
&1.678674485546577D-2,1.658311285001439D-2,1.6357649620221D-2,1.611
&007104151424D-2,1.584013731978726D-2,1.554767102436904D-2,1.523257
&435476985D-2,1.489484493432945D-2,1.453458959590346D-2,1.415203582
&870152D-2,1.374754076526337D-2,1.332159778154144D-2,1.287484094570
&53D-2,1.24080476741203D-2,1.192214003393026D-2,1.141818517380656D-
&2,1.089739537407941D-2,1.036112819269042D-2,9.810887152688882D-3,9
&.248323378736119D-3,8.675238553145434D-3,8.093589537892668D-3,7.50
&5495016816082D-3,6.913244587476468D-3,6.31931094183439D-3,5.726366
&236714079D-3,5.137304648633606D-3,4.555274633238511D-3,3.983726580
&652995D-3,3.426483604750896D-3,2.887842289651325D-3,2.372699521551
&761D-3,1.886670513946234D-3,1.436115964703238D-3,1.028038174940154
&D-3,6.702583391738669D-4,3.728569938879705D-4,1.496107785399853D-4
&,1.991706885703549D-5 /
end
SUBROUTINE abcpes191(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP191 / abc(191), pes(191)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,191
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end
!
!
!
BLOCK DATA abpeGP383
implicit real*8 (a-h,o-z)
common / abcpesGP383 / abc(383), pes(383)
DATA (abc(i),i=1,150) / -9.06179845938664D-1,-5.384693101056831D-1
&,0.D0,5.384693101056831D-1,9.06179845938664D-1,-9.840853600948425D
&-1,-7.541667265708492D-1,-2.796304131617832D-1,2.796304131617832D-
&1,7.541667265708492D-1,9.840853600948425D-1,-9.975454975868673D-1,
&-9.54373144738509D-1,-8.39137030805526D-1,-6.531199781049166D-1,-4
&.130159172460904D-1,-1.410830454016985D-1,1.410830454016985D-1,4.1
&30159172460904D-1,6.531199781049166D-1,8.39137030805526D-1,9.54373
&144738509D-1,9.975454975868673D-1,-9.996428940629623D-1,-9.9260947
&42582342D-1,-9.714532078050934D-1,-9.326448966646406D-1,-8.7498256
&23646846D-1,-7.987951932551814D-1,-7.055094496670671D-1,-5.9732444
&23227357D-1,-4.769125924857435D-1,-3.471380549365452D-1,-2.1083475
&8239592D-1,-7.069888485151317D-2,7.069888485151317D-2,2.1083475823
&9592D-1,3.471380549365452D-1,4.769125924857435D-1,5.97324442322735
&7D-1,7.055094496670671D-1,7.987951932551814D-1,8.749825623646846D-
&1,9.326448966646406D-1,9.714532078050934D-1,9.926094742582342D-1,9
&.996428940629623D-1,-9.999498338042262D-1,-9.988910202139731D-1,-9
&.954852422297365D-1,-9.888331807398711D-1,-9.78307721609879D-1,-9.
&634845653255284D-1,-9.440979512415834D-1,-9.200061936727053D-1,-8.
&911692033084103D-1,-8.576327998942287D-1,-8.195162842820891D-1,-7.
&770016931497028D-1,-7.303239512945866D-1,-6.797614383331393D-1,-6.
&256266847329246D-1,-5.682571221004128D-1,-5.080060481051527D-1,-4.
&452341960318211D-1,-3.803024612130273D-1,-3.1356635867478D-1,-2.45
&3726209119229D-1,-1.760580133874071D-1,-1.059500582206001D-1,-3.53
&6907741210191D-2,3.536907741210191D-2,1.059500582206001D-1,1.76058
&0133874071D-1,2.453726209119229D-1,3.1356635867478D-1,3.8030246121
&30273D-1,4.452341960318211D-1,5.080060481051527D-1,5.6825712210041
&28D-1,6.256266847329246D-1,6.797614383331393D-1,7.303239512945866D
&-1,7.770016931497028D-1,8.195162842820891D-1,8.576327998942287D-1,
&8.911692033084103D-1,9.200061936727053D-1,9.440979512415834D-1,9.6
&34845653255284D-1,9.78307721609879D-1,9.888331807398711D-1,9.95485
&2422297365D-1,9.988910202139731D-1,9.999498338042262D-1,-9.9999308
&00957606D-1,-9.998413733868334D-1,-9.993326036260333D-1,-9.9830059
&04912557D-1,-9.966114125082292D-1,-9.941549532502683D-1,-9.9083864
&5014144D-1,-9.865846355418413D-1,-9.81328454820258D-1,-9.750176168
&805138D-1,-9.67610132633725D-1,-9.590731842806808D-1,-9.4938205432
&55846D-1,-9.385192846195623D-1,-9.264740062449535D-1,-9.1324138646
&75595D-1,-8.988221544560521D-1,-8.832221810638305D-1,-8.6645209709
&84496D-1,-8.485269399006613D-1,-8.294658210610491D-1,-8.0929160978
&71924D-1,-7.880306274725631D-1,-7.657123497864587D-1,-7.4236911329
&73484D-1,-7.180358243555749D-1,-6.927496687383048D-1,-6.6654982141
&33578D-1,-6.39477156702661D-1,-6.115739600970961D-1,-5.82883643952
&0516D-1,-5.534504702200102D-1,-5.233192841812306D-1,-4.92535263734
&6438D-1,-4.611436891258013D-1,-4.291897379448173D-1,-3.96718309776
&3532D-1,-3.637738840148685D-1,-3.304004131090205D-1,-2.96641251960
&2923D-1,-2.625391225146616D-1,-2.281361109327633D-1,-1.93473693299
&7303D-1,-1.585927848221899D-1,-1.23533806989227D-1,-8.833676729901
&243D-2,-5.304134682592874D-2,-1.768699197607406D-2,1.7686991976074
&06D-2,5.304134682592874D-2,8.833676729901243D-2,1.23533806989227D-
&1,1.585927848221899D-1,1.934736932997303D-1,2.281361109327633D-1 /
&
DATA (abc(i),i=151,300) / 2.625391225146616D-1,2.966412519602923D-
&1,3.304004131090205D-1,3.637738840148685D-1,3.967183097763532D-1,4
&.291897379448173D-1,4.611436891258013D-1,4.925352637346438D-1,5.23
&3192841812306D-1,5.534504702200102D-1,5.828836439520516D-1,6.11573
&9600970961D-1,6.39477156702661D-1,6.665498214133578D-1,6.927496687
&383048D-1,7.180358243555749D-1,7.423691132973484D-1,7.657123497864
&587D-1,7.880306274725631D-1,8.092916097871924D-1,8.294658210610491
&D-1,8.485269399006613D-1,8.664520970984496D-1,8.832221810638305D-1
&,8.988221544560521D-1,9.132413864675595D-1,9.264740062449535D-1,9.
&385192846195623D-1,9.493820543255846D-1,9.590731842806808D-1,9.676
&10132633725D-1,9.750176168805138D-1,9.81328454820258D-1,9.86584635
&5418413D-1,9.90838645014144D-1,9.941549532502683D-1,9.966114125082
&292D-1,9.983005904912557D-1,9.993326036260333D-1,9.998413733868334
&D-1,9.999930800957606D-1,-9.999990540780864D-1,-9.999779201180271D
&-1,-9.999053843277264D-1,-9.997548099292319D-1,-9.995029904670225D
&-1,-9.991293637004486D-1,-9.986154377889328D-1,-9.979445541563222D
&-1,-9.97101673838173D-1,-9.960731134876332D-1,-9.94846326350319D-1
&,-9.934097684367996D-1,-9.917528288320241D-1,-9.898657893182227D-1
&,-9.877397894749672D-1,-9.853667873289969D-1,-9.827395144462354D-1
&,-9.798514279898772D-1,-9.766966628152694D-1,-9.732699859447113D-1
&,-9.695667547791078D-1,-9.655828795794257D-1,-9.613147901930057D-1
&,-9.567594066806098D-1,-9.519141133552408D-1,-9.467767357107768D-1
&,-9.413455197489509D-1,-9.356191132738967D-1,-9.295965487935348D-1
&,-9.232772277347039D-1,-9.166609057383601D-1,-9.097476788503845D-1
&,-9.025379704626385D-1,-8.950325188890622D-1,-8.872323654843452D-1
&,-8.791388432296039D-1,-8.707535657219533D-1,-8.620784165140526D-1
&,-8.531155387565449D-1,-8.43867325101528D-1,-8.343364078293007D-1,
&-8.245256491640314D-1,-8.144381317469599D-1,-8.040771492384722D-1,
&-7.934461970230241D-1,-7.825489629935193D-1,-7.713893183944502D-1,
&-7.599713087059173D-1,-7.482991445535909D-1,-7.363771926327818D-1,
&-7.242099666380492D-1,-7.118021181931931D-1,-6.99158427780045D-1,-
&6.862837956681647D-1,-6.731832328513438D-1,-6.598618520006799D-1,-
&6.463248584478648D-1,-6.325775412161859D-1,-6.186252641204988D-1,-
&6.044734569610379D-1,-5.901276068392971D-1,-5.755932496272772D-1,-
&5.608759616240512D-1,-5.459813514357845D-1,-5.309150521169581D-1,-
&5.156827136115164D-1,-5.002899955329195D-1,-4.847425603215671D-1,-
&4.690460668167286D-1,-4.532061642779477D-1,-4.372284868878711D-1,-
&4.211186487646215D-1,-4.048822395072182D-1,-3.88524820292248D-1,-3
&.720519205340819D-1,-3.5546903511457D-1,-3.387816221814662D-1,-3.2
&19951015080254D-1,-3.05114853399454D-1,-2.881462181253911D-1,-2.71
&0944958515495D-1,-2.539649470382427D-1,-2.36762793268957D-1,-2.194
&932184685372D-1,-2.021613704680775D-1,-1.847723628723129D-1,-1.673
&312771852347D-1,-1.49843165150782D-1,-1.323130512677196D-1,-1.1474
&59354410743D-1,-9.714679573659039D-2,-7.952059120936545D-2,-6.1872
&2647828821D-2,-4.420674615979685D-2,-2.652895475079536D-2,-8.84380
&2612311925D-3,8.843802612311925D-3,2.652895475079536D-2,4.42067461
&5979685D-2,6.18722647828821D-2,7.952059120936545D-2,9.714679573659
&039D-2,1.147459354410743D-1,1.323130512677196D-1,1.49843165150782D
&-1,1.673312771852347D-1,1.847723628723129D-1,2.021613704680775D-1,
&2.194932184685372D-1 /
DATA (abc(i),i=301,383) / 2.36762793268957D-1,2.539649470382427D-1
&,2.710944958515495D-1,2.881462181253911D-1,3.05114853399454D-1,3.2
&19951015080254D-1,3.387816221814662D-1,3.5546903511457D-1,3.720519
&205340819D-1,3.88524820292248D-1,4.048822395072182D-1,4.2111864876
&46215D-1,4.372284868878711D-1,4.532061642779477D-1,4.6904606681672
&86D-1,4.847425603215671D-1,5.002899955329195D-1,5.156827136115164D
&-1,5.309150521169581D-1,5.459813514357845D-1,5.608759616240512D-1,
&5.755932496272772D-1,5.901276068392971D-1,6.044734569610379D-1,6.1
&86252641204988D-1,6.325775412161859D-1,6.463248584478648D-1,6.5986
&18520006799D-1,6.731832328513438D-1,6.862837956681647D-1,6.9915842
&7780045D-1,7.118021181931931D-1,7.242099666380492D-1,7.36377192632
&7818D-1,7.482991445535909D-1,7.599713087059173D-1,7.71389318394450
&2D-1,7.825489629935193D-1,7.934461970230241D-1,8.040771492384722D-
&1,8.144381317469599D-1,8.245256491640314D-1,8.343364078293007D-1,8
&.43867325101528D-1,8.531155387565449D-1,8.620784165140526D-1,8.707
&535657219533D-1,8.791388432296039D-1,8.872323654843452D-1,8.950325
&188890622D-1,9.025379704626385D-1,9.097476788503845D-1,9.166609057
&383601D-1,9.232772277347039D-1,9.295965487935348D-1,9.356191132738
&967D-1,9.413455197489509D-1,9.467767357107768D-1,9.519141133552408
&D-1,9.567594066806098D-1,9.613147901930057D-1,9.655828795794257D-1
&,9.695667547791078D-1,9.732699859447113D-1,9.766966628152694D-1,9.
&798514279898772D-1,9.827395144462354D-1,9.853667873289969D-1,9.877
&397894749672D-1,9.898657893182227D-1,9.917528288320241D-1,9.934097
&684367996D-1,9.94846326350319D-1,9.960731134876332D-1,9.9710167383
&8173D-1,9.979445541563222D-1,9.986154377889328D-1,9.99129363700448
&6D-1,9.995029904670225D-1,9.997548099292319D-1,9.999053843277264D-
&1,9.999779201180271D-1,9.999990540780864D-1 /
DATA (pes(i),i=1,150) / 3.604869603234172D-3,7.533267383513908D-3,
&8.843904817338855D-3,7.533267383513908D-3,3.604869603234172D-3,1.3
&13500245384313D-3,5.836173070239102D-3,8.525970727887892D-3,8.5259
&70727887892D-3,5.836173070239102D-3,1.313500245384313D-3,4.2115561
&44766466D-4,2.422601408775297D-3,4.765518288102288D-3,6.7686069852
&608D-3,8.118320494376903D-3,8.765161099753414D-3,8.765161099753414
&D-3,8.118320494376903D-3,6.7686069852608D-3,4.765518288102288D-3,2
&.422601408775297D-3,4.211556144766466D-4,1.254828182226315D-4,8.28
&2674182342238D-4,1.851530457760147D-3,3.011269804307406D-3,4.19267
&8018766656D-3,5.315552960281358D-3,6.321947316769172D-3,7.17304041
&1104604D-3,7.848364458220213D-3,8.34381938105043D-3,8.666030528664
&805D-3,8.824265908512679D-3,8.824265908512679D-3,8.666030528664805
&D-3,8.34381938105043D-3,7.848364458220213D-3,7.173040411104604D-3,
&6.321947316769172D-3,5.315552960281358D-3,4.192678018766656D-3,3.0
&11269804307406D-3,1.851530457760147D-3,8.282674182342238D-4,1.2548
&28182226315D-4,3.571684231728753D-5,2.566164051950194D-4,6.1315321
&87427696D-4,1.062832486827296D-3,1.577273665045425D-3,2.1339805884
&70138D-3,2.715579850138872D-3,3.308163029694864D-3,3.9001044985854
&37D-3,4.481488748801018D-3,5.043827348237647D-3,5.579906601670604D
&-3,6.08370992745685D-3,6.55038794353066D-3,6.976251689046081D-3,7.
&358760866007352D-3,7.6964770218848D-3,7.988955666279293D-3,8.23656
&4360945532D-3,8.440235310052109D-3,8.601184879393912D-3,8.72064808
&2297092D-3,8.799673159135581D-3,8.83899803211562D-3,8.838998032115
&62D-3,8.799673159135581D-3,8.720648082297092D-3,8.601184879393912D
&-3,8.440235310052109D-3,8.236564360945532D-3,7.988955666279293D-3,
&7.6964770218848D-3,7.358760866007352D-3,6.976251689046081D-3,6.550
&38794353066D-3,6.08370992745685D-3,5.579906601670604D-3,5.04382734
&8237647D-3,4.481488748801018D-3,3.900104498585437D-3,3.30816302969
&4864D-3,2.715579850138872D-3,2.133980588470138D-3,1.57727366504542
&5D-3,1.062832486827296D-3,6.131532187427696D-4,2.566164051950194D-
&4,3.571684231728753D-5,9.91523227811318D-6,7.480536669060081D-5,1.
&864284969356289D-4,3.351291695869296D-4,5.140190874700771D-4,7.180
&579823516192D-4,9.433352569731171D-4,1.18634976077588D-3,1.4439211
&44825663D-3,1.713241802375448D-3,1.991863290326498D-3,2.2776373166
&19255D-3,2.568652324316803D-3,2.86318311835704D-3,3.15965547091719
&5D-3,3.456622293738234D-3,3.752747508408041D-3,4.046794768946334D-
&3,4.337619276572717D-3,4.624161689368059D-3,4.905443576344441D-3,5
&.180564096345209D-3,5.448697687039704D-3,5.70909258690328D-3,5.961
&07001696513D-3,6.20402383706015D-3,6.437420472852648D-3,6.66079889
&0770722D-3,6.873770382631686D-3,7.076017914350761D-3,7.26729479795
&1729D-3,7.447422467164724D-3,7.616287177384924D-3,7.77383551218452
&D-3,7.920068659893632D-3,8.05503552075712D-3,8.1788248101105D-3,8.
&291556425007196D-3,8.393372427732884D-3,8.484428056203499D-3,8.564
&883186772511D-3,8.634894642298812D-3,8.694609657048475D-3,8.744160
&688057923D-3,8.783661616387756D-3,8.813205234358121D-3,8.832861792
&378126D-3,8.842678304915744D-3,8.842678304915744D-3,8.832861792378
&126D-3,8.813205234358121D-3,8.783661616387756D-3,8.744160688057923
&D-3,8.694609657048475D-3,8.634894642298812D-3 /
DATA (pes(i),i=151,300) / 8.564883186772511D-3,8.484428056203499D-
&3,8.393372427732884D-3,8.291556425007196D-3,8.1788248101105D-3,8.0
&5503552075712D-3,7.920068659893632D-3,7.77383551218452D-3,7.616287
&177384924D-3,7.447422467164724D-3,7.267294797951729D-3,7.076017914
&350761D-3,6.873770382631686D-3,6.660798890770722D-3,6.437420472852
&648D-3,6.20402383706015D-3,5.96107001696513D-3,5.70909258690328D-3
&,5.448697687039704D-3,5.180564096345209D-3,4.905443576344441D-3,4.
&624161689368059D-3,4.337619276572717D-3,4.046794768946334D-3,3.752
&747508408041D-3,3.456622293738234D-3,3.159655470917195D-3,2.863183
&11835704D-3,2.568652324316803D-3,2.277637316619255D-3,1.9918632903
&26498D-3,1.713241802375448D-3,1.443921144825663D-3,1.1863497607758
&8D-3,9.433352569731171D-4,7.180579823516192D-4,5.140190874700771D-
&4,3.351291695869296D-4,1.864284969356289D-4,7.480536669060081D-5,9
&.91523227811318D-6,2.730900453802486D-6,2.10271345352118D-5,5.3714
&81098362814D-5,9.878770118479197D-5,1.547400687862095D-4,2.2042637
&0643783D-4,2.948859700437433D-4,3.772485523447D-4,4.66770178996769
&2D-4,5.62834655593204D-4,6.649139459856811D-4,7.725279023570785D-4
&,8.852212736095659D-4,1.0025562808784D-3,1.124113348989334D-3,1.24
&949409642859D-3,1.378323620704722D-3,1.510251216758654D-3,1.644949
&690510391D-3,1.782113750729078D-3,1.921457949623798D-3,2.062714518
&609035D-3,2.205631313704426D-3,2.34996998093712D-3,2.4955043815145
&53D-3,2.642019274021699D-3,2.789309228651758D-3,2.937177739477908D
&-3,3.085436499658327D-3,3.233904807520907D-3,3.382409076310557D-3,
&3.530782425586379D-3,3.678864337075955D-3,3.826500361898042D-3,3.9
&73541869369421D-3,4.119845830164231D-3,4.265274628501346D-3,4.4096
&95899419508D-3,4.552982388177057D-3,4.695011829482839D-3,4.8356668
&44707901D-3,4.974834855506346D-3,5.112408012435225D-3,5.2482831372
&42267D-3,5.382361677511354D-3,5.514549672337077D-3,5.6447577276540
&98D-3,5.772900999783749D-3,5.898899185686309D-3,6.022676518328035D
&-3,6.144161765492019D-3,6.263288230285296D-3,6.379993751525392D-3,
&6.494220702131653D-3,6.605915983604318D-3,6.715031014651607D-3,6.8
&21521712026309D-3,6.925348461662744D-3,7.026476078266477D-3,7.1248
&73751606503D-3,7.220514977895769D-3,7.313377474823167D-3,7.4034430
&79019696D-3,7.490697625003239D-3,7.575130804948921D-3,7.6567360089
&71964D-3,7.735510145982581D-3,7.811453445570876D-3,7.8845692417955
&49D-3,7.954863740173162D-3,8.022345769583236D-3,8.087026521205611D
&-3,8.1489192769768D-3,8.208039130377487D-3,8.264402702630289D-3,8.
&318027857582548D-3,8.368933418661979D-3,8.417138891314384D-3,8.462
&664194256029D-3,8.505529402695916D-3,8.545754506406038D-3,8.583359
&185146405D-3,8.618362603495841D-3,8.650783226613983D-3,8.680638657
&882861D-3,8.707945498770452D-3,8.732719230648456D-3,8.754974117709
&213D-3,8.774723129589149D-3,8.791977881844435D-3,8.80674859206207D
&-3,8.819044049145431D-3,8.82887159320115D-3,8.836237103480898D-3,8
&.841144991996892D-3,8.843598200725715D-3,8.843598200725715D-3,8.84
&1144991996892D-3,8.836237103480898D-3,8.82887159320115D-3,8.819044
&049145431D-3,8.80674859206207D-3,8.791977881844435D-3,8.7747231295
&89149D-3,8.754974117709213D-3,8.732719230648456D-3,8.7079454987704
&52D-3,8.680638657882861D-3,8.650783226613983D-3 /
DATA (pes(i),i=301,383) / 8.618362603495841D-3,8.583359185146405D-
&3,8.545754506406038D-3,8.505529402695916D-3,8.462664194256029D-3,8
&.417138891314384D-3,8.368933418661979D-3,8.318027857582548D-3,8.26
&4402702630289D-3,8.208039130377487D-3,8.1489192769768D-3,8.0870265
&21205611D-3,8.022345769583236D-3,7.954863740173162D-3,7.8845692417
&95549D-3,7.811453445570876D-3,7.735510145982581D-3,7.6567360089719
&64D-3,7.575130804948921D-3,7.490697625003239D-3,7.403443079019696D
&-3,7.313377474823167D-3,7.220514977895769D-3,7.124873751606503D-3,
&7.026476078266477D-3,6.925348461662744D-3,6.821521712026309D-3,6.7
&15031014651607D-3,6.605915983604318D-3,6.494220702131653D-3,6.3799
&93751525392D-3,6.263288230285296D-3,6.144161765492019D-3,6.0226765
&18328035D-3,5.898899185686309D-3,5.772900999783749D-3,5.6447577276
&54098D-3,5.514549672337077D-3,5.382361677511354D-3,5.2482831372422
&67D-3,5.112408012435225D-3,4.974834855506346D-3,4.835666844707901D
&-3,4.695011829482839D-3,4.552982388177057D-3,4.409695899419508D-3,
&4.265274628501346D-3,4.119845830164231D-3,3.973541869369421D-3,3.8
&26500361898042D-3,3.678864337075955D-3,3.530782425586379D-3,3.3824
&09076310557D-3,3.233904807520907D-3,3.085436499658327D-3,2.9371777
&39477908D-3,2.789309228651758D-3,2.642019274021699D-3,2.4955043815
&14553D-3,2.34996998093712D-3,2.205631313704426D-3,2.06271451860903
&5D-3,1.921457949623798D-3,1.782113750729078D-3,1.644949690510391D-
&3,1.510251216758654D-3,1.378323620704722D-3,1.24949409642859D-3,1.
&124113348989334D-3,1.0025562808784D-3,8.852212736095659D-4,7.72527
&9023570785D-4,6.649139459856811D-4,5.62834655593204D-4,4.667701789
&967692D-4,3.772485523447D-4,2.948859700437433D-4,2.20426370643783D
&-4,1.547400687862095D-4,9.878770118479197D-5,5.371481098362814D-5,
&2.10271345352118D-5,2.730900453802486D-6 /
end
SUBROUTINE abcpes383(u,p,epsi)
implicit real*8 (a-h,o-z)
dimension u(383),p(383)
common / abcpesGP383 / abc(383), pes(383)
epsim=epsi*0.5D0
uepsi2=1.D0-epsi**2
Do i=1,383
v=1.D0/(abc(i)*epsi+1.D0)
u(i)=0.5D0+(0.5D0*abc(i)+epsim)*v
p(i)=0.5D0*pes(i)*uepsi2*v*v
end do
return
end