4 common /chb
/ pi
,sig
(0:N
)
7 c initialize some variables
20 c ******************************************************************
24 common /chb
/ pi
,sig
(0:N
)
26 dimension yzin1
(0:N
), yzin2
(0:N
)
28 c function subprograms
29 z
(i
,j
,k
) = 0.5*h
(i
,j
)*(sig
(k
)-1.)
31 c a four-way average of rhobar
34 & ( z
(2,2,k
) + z
(1,2,k
) +
35 & z
(2,1,k
) + z
(1,1,k
) )
38 c another four-way average of rhobar
45 & ( rtmp1
+ rtmp2
+ rtmp3
+ rtmp4
)
49 if (yzin1
(k
) .ne
. yzin2
(k
)) call abort
51 if (yzin1
(0) .ne
. -1371.) call abort
52 if (yzin1
(1) .ne
. -685.5) call abort
53 if (yzin1
(2) .ne
. 0.) call abort