pro opscaleheight,nz,ndensop,tp,zp,tn,rairv,sch,sch0

pop = ndensop * 1.e6 * 1.38e-23 * tp

sch = fltarr(nz)

for k=1,nz-2 do begin
   dz = (zp(k+1)-zp(k-1))*rairv(k)*tn(k)/9.8
   sch(k) = dz/(alog(pop(k+1))-alog(pop(k-1)))
endfor

sch0 = 1.38e-23 * 2.*tp * 6.023e26 / 9.8 / 16.

return
end