ETS  \$Id: Doxyfile 2162 2020-02-26 14:16:09Z g2dpc $
 All Classes Files Functions Variables Pages
_bnd.f
Go to the documentation of this file.
1  subroutine f_bndmat(rk,zk,nk,rlop,zlop,nlop,rprob,zprob,nprob)
2 
3  include 'double.inc'
4  include 'parevo.inc'
5  parameter(nkp=njlim)
6  include 'dim.inc'
7  include 'compol.inc'
8  include 'compol_add.inc'
9  real*8 rk(*),zk(*),rlop(*),zlop(*),rprob(*),zprob(*)
10 
11 
12  do j=2,nt1
13  rr=r(nr,j)
14  zz=z(nr,j)
15  do jb=2,nt1
16 
17  r0=r(nr,jb)
18  z0=z(nr,jb)
19 
20  r1=r(nr,jb+1)
21  z1=z(nr,jb+1)
22 
23  call bint(rr,zz,r0,z0,r1,z1,fint,1)
24 
25  binadg(jb,j)=fint
26 
27  enddo
28  enddo
29 
30  return
31 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
32 
33  ik_out=0
34 
35  do ik=1,nk
36 
37  rr=rk(ik)
38  zz=zk(ik)
39 
40  call numcel(rr,zz,ic,jc)
41 
42  if(ic.lt.nr) then !!!
43 
44  iprcon(ik)=1
45 
46  else !!!!!!!!!!!!!!!!!!!!!
47 
48  ik_out=ik_out+1
49  iprcon(ik)=0
50 
51  endif !!!!!!!!!!!!!!!!!!!!
52 
53  do 110 j=2,nt1
54 
55  r0=r(nr,j)
56  z0=z(nr,j)
57 
58  r1=r(nr,j+1)
59  z1=z(nr,j+1)
60 
61  call bint(rr,zz,r0,z0,r1,z1,fint,1)
62 
63  pinadg(ik,j)=fint
64  !pinadg(ik_out,j)=Fint
65 
66  !psitok(ik)=psitok(ik)-Fint*(dgdn(j)+dgdn(j+1))*0.5d0
67 
68  110 continue
69 
70  ! endif !!!!!!!!!!!!!!!!!!!!
71 
72  enddo
73 
74  nk_out=ik_out
75 
76 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
77 
78  ik_out=0
79 
80  do ik=1,nlop
81 
82  rr=rlop(ik)
83  zz=zlop(ik)
84 
85  call numcel(rr,zz,ic,jc)
86 
87  if(ic.lt.nr) then !!!
88 
89  iprlop(ik)=1
90 
91  else !!!!!!!!!!!!!!!!!!!!!
92 
93  ik_out=ik_out+1
94  iprlop(ik)=0
95 
96  do 120 j=2,nt1
97 
98  r0=r(nr,j)
99  z0=z(nr,j)
100 
101  r1=r(nr,j+1)
102  z1=z(nr,j+1)
103 
104  call bint(rr,zz,r0,z0,r1,z1,fint,1)
105 
106  adginl(ik_out,j)=fint
107 
108  !psitok(ik)=psitok(ik)-Fint*(dgdn(j)+dgdn(j+1))*0.5d0
109 
110  120 continue
111 
112  endif !!!!!!!!!!!!!!!!!!!!
113 
114  enddo
115 
116  nlop_out=ik_out
117 
118 
119 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
120 
121  ik_out=0
122 
123  do ik=1,nprob
124 
125  rr=rprob(ik)
126  zz=zprob(ik)
127 
128  call numcel(rr,zz,ic,jc)
129 
130  if(ic.lt.nr) then !!!
131 
132  iprprob(ik)=1
133 
134  else !!!!!!!!!!!!!!!!!!!!!
135 
136  ik_out=ik_out+1
137  iprprob(ik)=0
138 
139  do 130 j=2,nt1
140 
141  r0=r(nr,j)
142  z0=z(nr,j)
143 
144  r1=r(nr,j+1)
145  z1=z(nr,j+1)
146 
147  call bint_d(rr,zz,r0,z0,r1,z1,fint_r,fint_z,1)
148 
149  adginr(ik_out,j)=fint_r
150  adginz(ik_out,j)=fint_z
151 
152  !psitok(ik)=psitok(ik)-Fint*(dgdn(j)+dgdn(j+1))*0.5d0
153 
154  130 continue
155 
156  endif !!!!!!!!!!!!!!!!!!!!
157 
158  enddo
159 
160  nprob_out=ik_out
161 
162  return
163  end
164 
165 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
166 
167 c********************************************************************
168 c********************************************************************
169 c********************************************************************
170 c********************************************************************
171 
bint_d
subroutine bint_d(X, Y, R0, Z0, r1, z1, Fr, Fz, i)
Definition: scu.f:637
numcel
subroutine numcel(rrk, zzk, icell, jcell)
Definition: _fluxt.f:259
bint
subroutine bint(X, Y, R0, Z0, r1, z1, F, I)
Definition: scu.f:557
f_bndmat
subroutine f_bndmat(rk, zk, nk, rlop, zlop, nlop, rprob, zprob, nprob)
Definition: _bnd.f:1