1 | block data |
---|
2 | c---------------------------------------------------------------------- |
---|
3 | c |
---|
4 | include 'ucom.h' |
---|
5 | include 'upawcom.h' |
---|
6 | c |
---|
7 | common/splwall/rr,s,sig,sx,sy,sxy |
---|
8 | c units of input/output |
---|
9 | data nbz,nfld,nin,nsav,nschef /11,12,13,14,15/ |
---|
10 | data nav,ndes1,ndes2,ndiag,ndpout,nemit /16,17,18,19,20,21/ |
---|
11 | data nfield,nhist,nnout,nout1,nout2 /22,23,24,25,26/ |
---|
12 | data nsemit,nstat,ntrw /27,28,29/ |
---|
13 | data nuch/30/ ! for ip=3 |
---|
14 | data ninput/31/ |
---|
15 | data lun1,lun2,lun3,lun4 /31,32,33,34/ |
---|
16 | data nsnap,nback,nback1/35,36,37/ |
---|
17 | data nimp,npois/38,39/ ! nimp for ip=999 |
---|
18 | data lolun1,lolun2,lolun3,lolun4 /0,0,0,0/ |
---|
19 | c---set up spline functions needed for computing |
---|
20 | c---effects of image charge on wall. |
---|
21 | real rr(10),s(8),sig(8,10),sx(8,10), |
---|
22 | 1 sy(8,10),sxy(8,10) |
---|
23 | data rr(1),(sig(i,1),i=1,8) |
---|
24 | 1 / 0.0, .6619, .6566, .6154, .4989, .2477, |
---|
25 | 2 .0668, .0062, 0./ |
---|
26 | data rr(2),(sig(i,2),i=1,8) |
---|
27 | 1 / .1, .6684, .6624, .6200, .4998, .2458, |
---|
28 | 2 .0660, .0061, 0./ |
---|
29 | data rr(3),(sig(i,3),i=1,8) |
---|
30 | 1 / .2, .6888, .6804, .6341, .5022, .2402, |
---|
31 | 2 .0633, .0059, 0./ |
---|
32 | data rr(4),(sig(i,4),i=1,8) |
---|
33 | 1 / .3, .7257, .7132, .65845, .5052, .2302, |
---|
34 | 2 .0590, .0054, 0./ |
---|
35 | data rr(5),(sig(i,5),i=1,8) |
---|
36 | 1 / .4, .7840, .7657, .6945, .5069, .2151, |
---|
37 | 2 .0530, .0048, 0./ |
---|
38 | data rr(6),(sig(i,6),i=1,8) |
---|
39 | 1 / .5, .8735, .84745, .7433, .5033, .1940, |
---|
40 | 2 .0457, .0041, 0./ |
---|
41 | data rr(7),(sig(i,7),i=1,8) |
---|
42 | 1 / .6,1.0147, .9765, .8039, .4871, .1658, |
---|
43 | 2 .0373, .0034, 0./ |
---|
44 | data rr(8),(sig(i,8),i=1,8) |
---|
45 | 1 / .7,1.2604, 1.1867, .8662, .4447, .1301, |
---|
46 | 2 .0282, .0025, 0./ |
---|
47 | data rr(9),(sig(i,9),i=1,8) |
---|
48 | 1 / .8,1.78570, 1.5332, .8865, .3533, .0885, |
---|
49 | 2 .0187, .0016, 0./ |
---|
50 | data rr(10),(sig(i,10),i=1,8) |
---|
51 | 1 / .9,3.3579, 2.0374, .6850, .1953, .0443, |
---|
52 | 2 .0091, .0008, 0./ |
---|
53 | data (s(i),i=1,8)/0.,.08,.2,.4,.8,1.4,2.4,4./ |
---|
54 | end |
---|
55 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|