

	call prem

	vrho=vrho*1000.		
	vmu=vmu*1000.
	vlam=vlam*1000.

	vmu=vrho*vmu**2
	vlam=vrho*vlam**2-2*vmu

        forall(j=1:nz)rho(:,j)=vrho(j)
	forall(j=1:nz)mu(:,j)=vmu(j)
	forall(j=1:nz)lam(:,j)=vlam(j)
