forked from autoplot/scripts
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathreadNoaaPoes.jyds
17 lines (11 loc) · 1.21 KB
/
readNoaaPoes.jyds
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
# correct poes files from CDAWeb, which have faulty timetags. Here we
# read the timetags from Year,Day,msec instead of time.
resourceURI= getParam( 'resourceURI', 'ftp://satdat.ngdc.noaa.gov/sem/poes/data/raw/ngdc/2013/noaa19/poes_n19_20130409_raw.nc', 'example file to load' )
parm= getParam( 'parm', 'mep_pro_tel0_cps_p3', 'parameter to load', [ 'alt', 'lat', 'lon', 'mep_pro_tel0_cps_p1', 'mep_pro_tel0_cps_p2', 'mep_pro_tel0_cps_p3', 'mep_pro_tel0_cps_p4', 'mep_pro_tel0_cps_p5', 'mep_pro_tel0_cps_p6', 'mep_pro_tel90_cps_p1', 'mep_pro_tel90_cps_p2', 'mep_pro_tel90_cps_p3', 'mep_pro_tel90_cps_p4', 'mep_pro_tel90_cps_p5', 'mep_pro_tel90_cps_p6', 'mep_ele_tel0_cps_e1', 'mep_ele_tel0_cps_e2', 'mep_ele_tel0_cps_e3', 'mep_ele_tel90_cps_e1', 'mep_ele_tel90_cps_e2', 'mep_ele_tel90_cps_e3', 'mep_omni_cps_p6', 'mep_omni_cps_p7', 'mep_omni_cps_p8', 'mep_omni_cps_p9' ] )
data= getDataSet( '%s?%s' % ( resourceURI, parm ), monitor )
data.putProperty( QDataSet.LABEL, parm )
tyr= getDataSet( '%s?%s' % ( resourceURI, 'year' ) )
tday= getDataSet( '%s?%s' % ( resourceURI, 'day' ) )
tmsec= getDataSet( '%s?%s' % ( resourceURI, 'msec' ) )
ttag= toTimeDataSet(tyr,None,tday,None,None,None,tmsec*1e6)
result= link( ttag, data )