File |
Line number |
Procedure |
Code |
./pkg/generic_advdiff/gad_osc_hat_r.F |
118 |
GAD_OSC_HAT_R |
integer ir |
./pkg/generic_advdiff/gad_osc_hat_r.F |
121 |
GAD_OSC_HAT_R |
do ir = +1-3, Nr+3 |
./pkg/generic_advdiff/gad_osc_hat_r.F |
123 |
GAD_OSC_HAT_R |
CALL GAD_OSC_LOC_R(ir,mask,fbar,ohat) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
13 |
GAD_OSC_LOC_R |
SUBROUTINE GAD_OSC_LOC_R(ir,mask,fbar,ohat) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
21 |
GAD_OSC_LOC_R |
integer ir |
./pkg/generic_advdiff/gad_osc_hat_r.F |
30 |
GAD_OSC_LOC_R |
if (ir.gt.+1-3 .and. |
./pkg/generic_advdiff/gad_osc_hat_r.F |
31 |
GAD_OSC_LOC_R |
& ir.lt.Nr+3) then |
./pkg/generic_advdiff/gad_osc_hat_r.F |
35 |
GAD_OSC_LOC_R |
floc(+0) = fbar(+0+ir) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
38 |
GAD_OSC_LOC_R |
& mask(ir-1)*(fbar(ir-1)-floc(+0)) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
40 |
GAD_OSC_LOC_R |
& mask(ir+1)*(fbar(ir+1)-floc(+0)) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
44 |
GAD_OSC_LOC_R |
ohat(+1,ir) = floc(+1)*0.25 _d 0 |
./pkg/generic_advdiff/gad_osc_hat_r.F |
47 |
GAD_OSC_LOC_R |
ohat(+2,ir) = floc(+1)*0.25 _d 0 |
./pkg/generic_advdiff/gad_osc_hat_r.F |
53 |
GAD_OSC_LOC_R |
if (ir.eq.+1-3) then |
./pkg/generic_advdiff/gad_osc_hat_r.F |
57 |
GAD_OSC_LOC_R |
floc(+0) = fbar(+0+ir) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
60 |
GAD_OSC_LOC_R |
& mask(ir+1)*(fbar(ir+1)-floc(+0)) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
64 |
GAD_OSC_LOC_R |
ohat(+1,ir) = floc(+1)*0.50 _d 0 |
./pkg/generic_advdiff/gad_osc_hat_r.F |
67 |
GAD_OSC_LOC_R |
ohat(+2,ir) = 0. _d 0 |
./pkg/generic_advdiff/gad_osc_hat_r.F |
71 |
GAD_OSC_LOC_R |
if (ir.eq.Nr+3) then |
./pkg/generic_advdiff/gad_osc_hat_r.F |
75 |
GAD_OSC_LOC_R |
floc(+0) = fbar(+0+ir) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
78 |
GAD_OSC_LOC_R |
& mask(ir-1)*(fbar(ir-1)-floc(+0)) |
./pkg/generic_advdiff/gad_osc_hat_r.F |
82 |
GAD_OSC_LOC_R |
ohat(+1,ir) = floc(+0)*0.50 _d 0 |
./pkg/generic_advdiff/gad_osc_hat_r.F |
85 |
GAD_OSC_LOC_R |
ohat(+2,ir) = 0. _d 0 |
./pkg/generic_advdiff/gad_osc_mul_r.F |
16 |
GAD_OSC_MUL_R |
integer ir,hh |
./pkg/generic_advdiff/gad_osc_mul_r.F |
36 |
GAD_OSC_MUL_R |
do ii = ir-hh, ir+hh |
./pkg/generic_advdiff/gad_osc_mul_r.F |
39 |
GAD_OSC_MUL_R |
dels = (ii - ir) * 2. _d 0 |
./pkg/generic_advdiff/gad_osc_mul_r.F |
6 |
GAD_OSC_MUL_R |
SUBROUTINE GAD_OSC_MUL_R(ir,hh,mask,ohat,scal) |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
101 |
GAD_PPM_ADV_R |
do ir = 2, Nr |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
104 |
GAD_PPM_ADV_R |
& + abs(velR(ix,iy,ir)) |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
109 |
GAD_PPM_ADV_R |
do ir = 1, Nr |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
111 |
GAD_PPM_ADV_R |
floc(ir) = fbar(ix,iy,ir) |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
112 |
GAD_PPM_ADV_R |
mloc(ir) = |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
113 |
GAD_PPM_ADV_R |
& maskC(ix,iy,ir,bi,bj) |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
151 |
GAD_PPM_ADV_R |
do ir = 2, Nr |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
153 |
GAD_PPM_ADV_R |
flux(ix,iy,ir) = 0.0 _d +0 |
./pkg/generic_advdiff/gad_ppm_adv_r.F |
61 |
GAD_PPM_ADV_R |
integer ix,iy,ir |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
100 |
GAD_PPM_FLX_R |
wCFL = wvel(ix,iy,ir) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
101 |
GAD_PPM_FLX_R |
& * delT(ir-0)*recip_drF(ir-0) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
113 |
GAD_PPM_FLX_R |
intF = ivec(1) * fhat(1,ir-0) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
114 |
GAD_PPM_FLX_R |
& + ivec(2) * fhat(2,ir-0) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
115 |
GAD_PPM_FLX_R |
& + ivec(3) * fhat(3,ir-0) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
122 |
GAD_PPM_FLX_R |
flux(ix,iy,ir) = + wfac(ix,iy,ir) * intF |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
53 |
GAD_PPM_FLX_R |
integer ir |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
67 |
GAD_PPM_FLX_R |
do ir = +2, Nr |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
69 |
GAD_PPM_FLX_R |
if (wvel(ix,iy,ir) .eq. 0. _d 0) then |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
71 |
GAD_PPM_FLX_R |
flux(ix,iy,ir) = 0. _d 0 |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
75 |
GAD_PPM_FLX_R |
if (wvel(ix,iy,ir) .lt. 0. _d 0) then |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
78 |
GAD_PPM_FLX_R |
wCFL = wvel(ix,iy,ir) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
79 |
GAD_PPM_FLX_R |
& * delT(ir-1)*recip_drF(ir-1) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
91 |
GAD_PPM_FLX_R |
intF = ivec(1) * fhat(1,ir-1) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
92 |
GAD_PPM_FLX_R |
& + ivec(2) * fhat(2,ir-1) |
./pkg/generic_advdiff/gad_ppm_flx_r.F |
93 |
GAD_PPM_FLX_R |
& + ivec(3) * fhat(3,ir-1) |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
113 |
GAD_PPM_HAT_R |
CALL GAD_OSC_MUL_R(ir,+2,mask, |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
131 |
GAD_PPM_HAT_R |
fhat(ii,ir) = lhat(ii) |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
56 |
GAD_PPM_HAT_R |
integer ii,ir |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
66 |
GAD_PPM_HAT_R |
do ir = +1, Nr |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
69 |
GAD_PPM_HAT_R |
ff00 = fbar(ir+0) |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
71 |
GAD_PPM_HAT_R |
& + mask(ir-1)*(fbar(ir-1)-ff00) |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
73 |
GAD_PPM_HAT_R |
& + mask(ir+1)*(fbar(ir+1)-ff00) |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
75 |
GAD_PPM_HAT_R |
fell = edge(ir-0) |
./pkg/generic_advdiff/gad_ppm_hat_r.F |
76 |
GAD_PPM_HAT_R |
ferr = edge(ir+1) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
28 |
GAD_PPM_P3E_R |
integer ir |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
33 |
GAD_PPM_P3E_R |
do ir = +1, Nr+1 |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
36 |
GAD_PPM_P3E_R |
mloc(-1) = mask(ir-1) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
37 |
GAD_PPM_P3E_R |
mloc(+0) = mask(ir+0) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
39 |
GAD_PPM_P3E_R |
floc(-1) = fbar(ir+0) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
40 |
GAD_PPM_P3E_R |
& + mloc(-1)*(fbar(ir-1)-fbar(ir+0)) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
41 |
GAD_PPM_P3E_R |
floc(+0) = fbar(ir-1) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
42 |
GAD_PPM_P3E_R |
& + mloc(+0)*(fbar(ir+0)-fbar(ir-1)) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
44 |
GAD_PPM_P3E_R |
mloc(-2) = mask(ir-2) * mloc(-1) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
48 |
GAD_PPM_P3E_R |
& + mloc(-2)*(fbar(ir-2)-ftmp) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
50 |
GAD_PPM_P3E_R |
mloc(+1) = mask(ir+1) * mloc(+0) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
54 |
GAD_PPM_P3E_R |
& + mloc(+1)*(fbar(ir+1)-ftmp) |
./pkg/generic_advdiff/gad_ppm_p3e_r.F |
57 |
GAD_PPM_P3E_R |
edge(ir) = |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
101 |
GAD_PQM_ADV_R |
do ir = +1, Nr |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
104 |
GAD_PQM_ADV_R |
& + abs(wfac(ix,iy,ir)) |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
106 |
GAD_PQM_ADV_R |
floc(ir) = fbar (ix,iy,ir) |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
107 |
GAD_PQM_ADV_R |
mloc(ir) = |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
108 |
GAD_PQM_ADV_R |
& maskC(ix,iy,ir,bi,bj) |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
148 |
GAD_PQM_ADV_R |
do ir = 2, Nr |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
150 |
GAD_PQM_ADV_R |
flux(ix,iy,ir) = 0. _d 0 |
./pkg/generic_advdiff/gad_pqm_adv_r.F |
61 |
GAD_PQM_ADV_R |
integer ix,iy,ir |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
106 |
GAD_PQM_FLX_R |
wCFL = wvel(ix,iy,ir) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
107 |
GAD_PQM_FLX_R |
& * delT(ir-0)*recip_drF(ir-0) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
123 |
GAD_PQM_FLX_R |
intF = ivec(1) * fhat(1,ir-0) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
124 |
GAD_PQM_FLX_R |
& + ivec(2) * fhat(2,ir-0) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
125 |
GAD_PQM_FLX_R |
& + ivec(3) * fhat(3,ir-0) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
126 |
GAD_PQM_FLX_R |
& + ivec(4) * fhat(4,ir-0) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
127 |
GAD_PQM_FLX_R |
& + ivec(5) * fhat(5,ir-0) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
134 |
GAD_PQM_FLX_R |
flux(ix,iy,ir) = + wfac(ix,iy,ir) * intF |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
53 |
GAD_PQM_FLX_R |
integer ir |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
67 |
GAD_PQM_FLX_R |
do ir = +2, Nr |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
69 |
GAD_PQM_FLX_R |
if (wfac(ix,iy,ir) .eq. 0. _d 0) then |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
71 |
GAD_PQM_FLX_R |
flux(ix,iy,ir) = 0. _d 0 |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
75 |
GAD_PQM_FLX_R |
if (wfac(ix,iy,ir) .lt. 0. _d 0) then |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
78 |
GAD_PQM_FLX_R |
wCFL = wvel(ix,iy,ir) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
79 |
GAD_PQM_FLX_R |
& * delT(ir-1)*recip_drF(ir-1) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
95 |
GAD_PQM_FLX_R |
intF = ivec(1) * fhat(1,ir-1) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
96 |
GAD_PQM_FLX_R |
& + ivec(2) * fhat(2,ir-1) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
97 |
GAD_PQM_FLX_R |
& + ivec(3) * fhat(3,ir-1) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
98 |
GAD_PQM_FLX_R |
& + ivec(4) * fhat(4,ir-1) |
./pkg/generic_advdiff/gad_pqm_flx_r.F |
99 |
GAD_PQM_FLX_R |
& + ivec(5) * fhat(5,ir-1) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
133 |
GAD_PQM_HAT_R |
CALL GAD_OSC_MUL_R(ir,+2,mask, |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
151 |
GAD_PQM_HAT_R |
fhat(ii,ir) = lhat(ii) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
62 |
GAD_PQM_HAT_R |
integer ii,ir |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
75 |
GAD_PQM_HAT_R |
do ir = +1, Nr |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
78 |
GAD_PQM_HAT_R |
rhat = drF(ir) * .5 _d 0 |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
81 |
GAD_PQM_HAT_R |
ff00 = fbar(ir+0) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
83 |
GAD_PQM_HAT_R |
& + mask(ir-1)*(fbar(ir-1)-ff00) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
85 |
GAD_PQM_HAT_R |
& + mask(ir+1)*(fbar(ir+1)-ff00) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
87 |
GAD_PQM_HAT_R |
fell = edge(+1,ir-0) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
88 |
GAD_PQM_HAT_R |
ferr = edge(+1,ir+1) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
90 |
GAD_PQM_HAT_R |
dell = edge(+2,ir-0) |
./pkg/generic_advdiff/gad_pqm_hat_r.F |
91 |
GAD_PQM_HAT_R |
derr = edge(+2,ir+1) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
29 |
GAD_PQM_P5E_R |
integer ir |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
34 |
GAD_PQM_P5E_R |
do ir = +1, Nr+1 |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
37 |
GAD_PQM_P5E_R |
mloc(-1) = mask(ir-1) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
38 |
GAD_PQM_P5E_R |
mloc(+0) = mask(ir+0) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
40 |
GAD_PQM_P5E_R |
floc(-1) = fbar(ir+0) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
41 |
GAD_PQM_P5E_R |
& + mloc(-1)*(fbar(ir-1)-fbar(ir+0)) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
42 |
GAD_PQM_P5E_R |
floc(+0) = fbar(ir-1) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
43 |
GAD_PQM_P5E_R |
& + mloc(+0)*(fbar(ir+0)-fbar(ir-1)) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
45 |
GAD_PQM_P5E_R |
mloc(-2) = mask(ir-2) * mloc(-1) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
46 |
GAD_PQM_P5E_R |
mloc(-3) = mask(ir-3) * mloc(-2) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
50 |
GAD_PQM_P5E_R |
& + mloc(-2)*(fbar(ir-2)-ftmp) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
53 |
GAD_PQM_P5E_R |
& + mloc(-3)*(fbar(ir-3)-ftmp) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
55 |
GAD_PQM_P5E_R |
mloc(+1) = mask(ir+1) * mloc(+0) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
56 |
GAD_PQM_P5E_R |
mloc(+2) = mask(ir+2) * mloc(+1) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
60 |
GAD_PQM_P5E_R |
& + mloc(+1)*(fbar(ir+1)-ftmp) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
63 |
GAD_PQM_P5E_R |
& + mloc(+2)*(fbar(ir+2)-ftmp) |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
66 |
GAD_PQM_P5E_R |
edge(1,ir) = |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
71 |
GAD_PQM_P5E_R |
edge(2,ir) = ( |
./pkg/generic_advdiff/gad_pqm_p5e_r.F |
75 |
GAD_PQM_P5E_R |
& ) * recip_drC(ir) |