source: PSPA/parmelaPSPA/trunk/blockdata.f @ 418

Last change on this file since 418 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 2.0 KB
Line 
1      block data
2c----------------------------------------------------------------------
3c
4      include 'ucom.h'
5      include 'upawcom.h'
6c
7      common/splwall/rr,s,sig,sx,sy,sxy
8c     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/
19c---set up spline functions needed for computing
20c---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
55c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.