1 subroutine r8util_bcherm2(fherm,idimx1,idimx2,
2 > jbcx1a,jbcx1b, jbcx2a,jbcx2b,
3 > zbcx1a,zbcx1b, zbcx2a,zbcx2b,
10 INTEGER,
PARAMETER :: R8=selected_real_kind(12,100)
11 integer :: idimx1,idimx2
12 real*8 :: fherm(0:3,idimx1,idimx2)
14 integer :: jbcx1a,jbcx1b
15 real*8 :: zbcx1a(idimx2)
16 real*8 :: zbcx1b(idimx2)
18 integer :: jbcx2a,jbcx2b
19 real*8 :: zbcx2a(idimx1)
20 real*8 :: zbcx2b(idimx1)
22 real*8 :: x1(idimx1),x2(idimx2)
29 if((jbcx1a.eq.1).or.(jbcx1b.eq.1))
then
32 fherm(1,1,1:idimx2)=zbcx1a(1:idimx2)
36 fherm(1,1,ix)=(fherm(0,2,ix)-fherm(0,1,ix))/zdx
41 fherm(1,idimx1,1:idimx2)=zbcx1b(1:idimx2)
43 zdx = x1(idimx1)-x1(idimx1-1)
46 > (fherm(0,idimx1,ix)-fherm(0,idimx1-1,ix))/zdx
52 if((jbcx2a.eq.1).or.(jbcx2b.eq.1))
then
55 fherm(2,1:idimx1,1)=zbcx2a(1:idimx1)
59 fherm(2,ix,1)=(fherm(0,ix,2)-fherm(0,ix,1))/zdx
64 fherm(2,1:idimx1,idimx2)=zbcx2b(1:idimx1)
66 zdx=x2(idimx2)-x2(idimx2-1)
69 > (fherm(0,ix,idimx2)-fherm(0,ix,idimx2-1))/zdx
74 end subroutine r8util_bcherm2