subroutine tridi1	(	integer	l,
		integer	n,
		real(kind=kind_phys), dimension(l,2:n)	cl,
		real(kind=kind_phys), dimension(l,n)	cm,
		real(kind=kind_phys), dimension(l,n-1)	cu,
		real(kind=kind_phys), dimension(l,n)	r1,
		real(kind=kind_phys), dimension(l,n-1)	au,
		real(kind=kind_phys), dimension(l,n)	a1
	)