forked from rwharvey/cql3d
-
Notifications
You must be signed in to change notification settings - Fork 2
/
coeffpad.f
39 lines (34 loc) · 1.17 KB
/
coeffpad.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
c
c
subroutine coeffpad(k)
implicit integer (i-n), real*8 (a-h,o-z)
c real*8,dimension(iy):: prnt1,prnt2,prnt3
c real*8,dimension(iy):: prnt4,prnt5,prnt6
c..................................................................
c This routine adds in the collisional contribution to the
c coefficients employed in time advancement..
c scatfrac=0. disables pitch angle scattering (along with
c mx=0, see cqlinput_help. scatfrac=1. by default).
c..................................................................
include 'param.h'
include 'comm.h'
do 10 j=1,jx
do 11 i=1,iy
da(i,j)=da(i,j)+cal(i,j,k,l_)
db(i,j)=db(i,j)+cbl(i,j,k,l_)
dc(i,j)=dc(i,j)+scatfrac*ccl(i,j,k,l_)
dd(i,j)=dd(i,j)+scatfrac*cdl(i,j,k,l_) !YuP[2019-11-18] Added scatfrac
de(i,j)=de(i,j)+scatfrac*cel(i,j,k,l_)
df(i,j)=df(i,j)+scatfrac*cfl(i,j,k,l_)
11 continue
10 continue
c do i=1,iy
c prnt1(i)=da(i,2)
c prnt2(i)=db(i,2)
c prnt3(i)=dc(i,2)
c prnt4(i)=dd(i,2)
c prnt5(i)=de(i,2)
c prnt6(i)=df(i,2)
c enddo
return
end