;+
; NAME:
;  rt_szss96_zgrid
; PURPOSE: (one line)
;  Create an altitude grid
; DESCRIPTION:
;  Create an altitude grid
; CATEGORY:
;  RT
; CALLING SEQUENCE:
;  z = rt_szss96_zgrd(nz, delz0_cm, fac1)
; INPUTS:
;  nz = number of altitude layers
;  delz0_cm = delta z of the lowest layer (cm)
;  fac1 = fractional increase in subsequent delta z (uniteless)
; OPTIONAL INPUT PARAMETERS:
; KEYWORD INPUT PARAMETERS:
; KEYWORD OUTPUT PARAMETERS:
;  none
; OUTPUTS:
;    z   altitudes [nz]  (cm)
; COMMON BLOCKS:
; SIDE EFFECTS:
; RESTRICTIONS:
;  Modifies atm
; PROCEDURE:
; 
; EXAMPLE:
;
;  nz = rt_szss96_zgrid(106, 1d5, 1.03)
; 
; MODIFICATION HISTORY:
;  Written 2007 Aug 31, Leslie Young SwRI
;  based on plutot.f, written by Xun Zhu
;-

function rt_szss96_zgrid, nz, delz0_cm, fac1

    z = dblarr(nz)
    z[0] = 0.0d
    z[1] = z[0] + delz0_cm
    for k = 2, nz-1 do begin
        z[k] = z[k-1] + (z[k-1] - z[k-2]) * fac1
    end

    return, z

end