update

pyatmos/StandardAtmosphere/StandardAtmosphere.py

@@ -100,4 +100,4 @@ def isa(altitude,altitude_type='geometric'):
             h0 = geopotential_alt[i+1]
 
     density = pressure / (R * temperature)
-    return {'temperature':temperature,'pressure':pressure,'density':density}
+    return {'temperature[K]':temperature,'pressure[Pa]':pressure,'density[kg/m^3]':density}

pyatmos/__init__.py

@@ -3,3 +3,7 @@ pyatmos package
 
 This package is an archive of scientific routines that implements the estimation of atmospheric properties for various atmosphere models.  
 '''    
+
+from .StandardAtmosphere.StandardAtmosphere import isa
+from .msise.spaceweather import download_sw,read_sw
+from .msise.nrlmsise00 import nrlmsise00

pyatmos/msise/nrlmsise00.py

@@ -49,6 +49,7 @@ nrlmsise00
 '''
 
 import numpy as np
+from astropy.time import Time
 from scipy.interpolate import CubicSpline
 from pyshtools.legendre import PLegendreA,PlmIndex
 import pkg_resources
@@ -1021,7 +1022,8 @@ def gts7(inputp,flags,gsurf,re):
 
 # ============================ nrlmsise00 =========================== #
 
-def nrlmsise00(t,lat,lon,alt,SW_OBS_PRE,omode,aphmode):
+def nrlmsise00(t,lat,lon,alt,SW_OBS_PRE,omode='Oxygen',aphmode='NoAph'):
+    t = Time(t)
     lon_wrap = wraplon(lon)
     t_yday = t.yday.split(':')
     t_ymd = t.iso.split()[0].split('-')