c----------------------------------------------------------------------
	subroutine GGALAM(nl,ngga)
c calculates and outputs (to eclam.d) the exact lambda
c decomposition of a given approximate GGA, denoted by ngga
c
c author: Kieron Burke, 8.14.99
c     in: densities
c	  nl = no of lambda intervals, uniformly spaced between 0 and 1
c	  ngga=which correlation gga to use
c	  1=LSD,2=PW91,3=PBE,4=BLYP
c    out: eclam = Ec[n_{1/lambda}]
c         tclam = Tc[n_{1/lambda}]
c         uclam = Uc[n_{1/lambda}]
c  calls: SCALEGGA
c  some variables used: egga(j,k)
c                       j : 1=>X, 2=>C,
c                       tcgga(j,k)
c                       j : 1=> Tc, 2=>Uc
c                       vdum(j,k,l)
c               	l : 1=>up, 2=>dn
	parameter(ml=100)
	include 'basics.inc'
	include 'dens.inc'
	include 'sys.inc'
	dimension eclam(0:ml),tclam(0:ml),uclam(0:ml)
	dimension der(3,4),egga(2,4),eggap(2,4),eggam(2,4)
	dimension dum(0:m,4,0:ml),vdum(2,4,2)