;+
; NAME:
;  rt_atmread2
; PURPOSE: (one line)
;  read a model atmosphere
; DESCRIPTION:
;  read a model atmosphere, return atmospheric structure
; CATEGORY:
;  RT
; CALLING SEQUENCE:
;  atm = rt_atmread2(fname, pmax)
; INPUTS:
;  fname = file name of atmospheric structure
;  pmax = maximum pressure (ubar) 
; OPTIONAL INPUT PARAMETERS:
; KEYWORD INPUT PARAMETERS:
;  none
; KEYWORD OUTPUT PARAMETERS:
;  none
; OUTPUTS:
;    returns an atmosphere structure, with the following fields
;    nz  number of atmospheric levels
;    ns  number of species
;    z   altitudes [nz]  (cm)
;    p   pressures [nz] (cm)
;    t   temperatures [nz] (K)
;    n   number density [nz,ns] (cm^-3)
;    ncol   number column density [nz,ns] (cm^-2)
; COMMON BLOCKS:
; SIDE EFFECTS:
; RESTRICTIONS:
;  Modifies atm
; PROCEDURE:
; 
; EXAMPLE:
;
;  atm = rt_atmread2('titan.atm',1.d+2)
; 
; MODIFICATION HISTORY:
;  Written 2004 May 8, Leslie Young SwRI
;  based on atmread2.f, written by Roger Yelle, 1990
;-

function rt_atmread2, fname, pmax
    
    atm1 = rt_atmread1(fname)
    nz = rt_locate(atm1.p, pmax)+2
    ns = atm1.ns

    z = atm1.z[0:nz-1] 
    p = atm1.p[0:nz-1] 
    t = atm1.t[0:nz-1] 
    n = atm1.n[0:nz-1,*] 
    ncol = atm1.ncol[0:nz-1,*] 
    
    atm = {  nz:nz, ns:ns,  z:z, p:p, t:t, n:n, ncol:ncol }

    return, atm

end