5 common /chb
/ pi
,sig
(0:N
)
8 c initialize some variables
21 c ******************************************************************
25 common /chb
/ pi
,sig
(0:N
)
27 dimension yzin1
(0:N
), yzin2
(0:N
)
29 c function subprograms
30 z
(i
,j
,k
) = 0.5*h
(i
,j
)*(sig
(k
)-1.)
32 c a four-way average of rhobar
35 & ( z
(2,2,k
) + z
(1,2,k
) +
36 & z
(2,1,k
) + z
(1,1,k
) )
39 c another four-way average of rhobar
46 & ( rtmp1
+ rtmp2
+ rtmp3
+ rtmp4
)
50 if (yzin1
(k
) .ne
. yzin2
(k
)) call abort
52 if (yzin1
(0) .ne
. -1371.) call abort
53 if (yzin1
(1) .ne
. -685.5) call abort
54 if (yzin1
(2) .ne
. 0.) call abort