from pyhdf import VS
import pyhdf
import pyhdf.HDF as HDF
import pyhdf.SD
import numpy as np
from readcloudsat import get_geo
import matplotlib.pyplot as plt
if __name__=="__main__":
ecmwfname='2008291181813_13156_CS_ECMWF-AUX_GRANULE_P_R04_E02.hdf'
lat,lon,time_utc,prof_seconds,dem_elevation=get_geo(ecmwfname)
hdffile=HDF.HDF(ecmwfname,HDF.HC.READ)
vs=hdffile.vstart()
the_vs_data=vs.attach('EC_height')
nrecs=the_vs_data._nrecs
height=the_vs_data.read(nRec=nrecs)
height=np.array(height).squeeze()
the_vs_data.detach()
readonly=pyhdf.SD.SDC.READ
sdfile=pyhdf.SD.SD(ecmwfname,readonly)
temp_soundings=sdfile.select('Temperature')
temp_array=temp_soundings.get()
sdfile.end()
start=20000
stop=22000
temp_array=temp_array[start:stop,:]
temp_avg=temp_array.mean(axis=0)
fig1=plt.figure(1)
fig1.clf()
axis1=fig1.add_subplot(1,1,1)
axis1.plot(temp_avg,height/1.e3)
""" modified 2012/11/23 to eliminate need for matlab"""
import matplotlib
import datetime
import dateutil.tz as tz
import numpy as np
import matplotlib.dates as mpldates
import matplotlib.pyplot as plt
import h5py
import sys
from readcloudsat import get_geo,convert_field
if __name__=="__main__":
#http://www.cloudsat.cira.colostate.edu/dataSpecs.php?prodid=9
rainfilename='2010247105814_23156_CS_2C-RAIN-PROFILE_GRANULE_P_R04_E03.h5'
#rainfilename='2009197050922_17109_CS_2C-RAIN-PROFILE_GRANULE_P_R04_E02.h5'
lat,lon,time_utc,prof_seconds,dem_elevation=get_geo(rainfilename)
with h5py.File(rainfilename,'r') as f:
rain_rate=f['2C-RAIN-PROFILE']['Data Fields']['rain_rate'].value
rain_rate=[item[0] for item in rain_rate]
rain_rate=np.array(rain_rate)
#rain_rate[rain_rate < -9000]=np.nan
plt.close('all')
fig1=plt.figure(1)
fig1.clf()
axis1=fig1.add_subplot(1,1,1)
axis1.plot(rain_rate)
axis1.set_xlabel('index value')
axis1.set_ylabel('rain rate (mm/hour)')
axis1.set_title('rain rate for all times')
axis1.set_ylim([0,10])
