module zpad1 {                       # Zero pad.  Surround data by zeros. 1-D
#%  _lop( data,  padd)    
integer                              p,  d
do d= 1, size(data) {                p = d + (size(padd)-size(data))/2
        if( adj)
                        data(d) = data(d) + padd(p)
	else
                        padd(p) = padd(p) + data(d)
        }
}