# Import readsmart package 
import readsmart as rs

# Specify file locations
atm_file = 'icrccm_62.atm'
rad_file = 'earth_standard_hitran08_50_100000cm-1_clearsky_toa.rad'

# Get atm file data
pressure, temperature, abundance_profiles, gas_codes = rs.atm(atm_file, skiprows=11, plot=False, getdata=True)

# Plot atm
rs.atm(atm_file, skiprows=11)

# Get rad file data
wl, wno, solar_spec, TOA_flux, rad_streams = rs.rad(rad_file, plot=False, getdata=True)

# Plot flux
rs.rad(rad_file, plot=True, getdata=False)

# Plot radiance streams
rs.rad(rad_file, plot=True, getdata=False, ptype='rad', xran=[1.0,15])

# Specify jacobian file extension tag
jext = '.j_O3'

# Get Radiance Jacobians
wl, wno, pressure, radiance, jacobian, xprofile, xdelta = rs.jacobians(rad_file, jext, plot=False, getdata=True)

# Plot Radiance Jacobians
rs.jacobians(rad_file, jext, plot=True, getdata=False, xran=[3.0,5.0], pran=[100.,10000.])