|
Revision 326
(checked in by haines, 14 years ago)
|
import first verison proc2plot
|
| Line | |
|---|
| 1 |
#!/usr/bin/env python |
|---|
| 2 |
# Last modified: Time-stamp: <2008-10-16 13:29:04 haines> |
|---|
| 3 |
"""test_concat_data""" |
|---|
| 4 |
|
|---|
| 5 |
import os, sys, glob |
|---|
| 6 |
import datetime, time, dateutil, dateutil.tz |
|---|
| 7 |
import pycdf |
|---|
| 8 |
import numpy |
|---|
| 9 |
|
|---|
| 10 |
sys.path.append('/home/haines/nccoos/raw2proc') |
|---|
| 11 |
del(sys) |
|---|
| 12 |
|
|---|
| 13 |
import procutil |
|---|
| 14 |
|
|---|
| 15 |
# test with jpier adcp |
|---|
| 16 |
print 'jpier_adcp ...' |
|---|
| 17 |
proc_dir ='/seacoos/data/nccoos/level1/jpier/adcp/' |
|---|
| 18 |
fns = glob.glob((os.path.join(proc_dir, '*.nc'))) |
|---|
| 19 |
fns.sort() |
|---|
| 20 |
nc = pycdf.CDFMF(fns) |
|---|
| 21 |
|
|---|
| 22 |
ncvars = nc.variables() |
|---|
| 23 |
# print ncvars |
|---|
| 24 |
es = nc.var('time')[:] |
|---|
| 25 |
units = nc.var('time').units |
|---|
| 26 |
dt = [procutil.es2dt(e) for e in es] |
|---|
| 27 |
# set timezone info to UTC (since data from level1 should be in UTC!!) |
|---|
| 28 |
dt = [e.replace(tzinfo=dateutil.tz.tzutc()) for e in dt] |
|---|
| 29 |
# return new datetime based on computer local |
|---|
| 30 |
dt_local = [e.astimezone(dateutil.tz.tzlocal()) for e in dt] |
|---|
| 31 |
dn = date2num(dt) |
|---|
| 32 |
z = nc.var('z')[:] |
|---|
| 33 |
wd = nc.var('water_depth')[:] |
|---|
| 34 |
u = nc.var('u')[:] |
|---|
| 35 |
v = nc.var('v')[:] |
|---|
| 36 |
nc.close() |
|---|
| 37 |
|
|---|
| 38 |
# averaged water depth over whole recorded deployment!! |
|---|
| 39 |
print wd.mean() |
|---|
