;####################################################################
function dotp,a,b

    ;print,'IN dotp'
    sza = size(a)
    szb = size(b)

    if sza[1] eq szb[1] then begin
        dot = 0.0d
        for i=0,sza[1]-1 do begin
            dot = dot + a[i]*b[i]
        endfor;i
    endif
    return,dot

end;dotp
