| (4) |
| (5) |
# Zero pad. Surround data by zeros. 1-D
#
subroutine zpad1( adj,add, data,nd, padd,np )
integer adj,add, d, nd, p, np
real data(nd), padd(np)
call adjnull( adj,add, data,nd, padd,np)
do d= 1, nd { p = d + (np-nd)/2
if( adj == 0 )
padd(p) = padd(p) + data(d)
else
data(d) = data(d) + padd(p)
}
return; end