Changeset Navigation

Changeset 104

Timestamp:

02/12/08 15:03:54

Author:

haines

Message:

Files:

r101	r104
1	1	#!/usr/bin/env python
2		# Last modified: Time-stamp: <2008-0~~1-03 09:50:53~~ haines>
	2	# Last modified: Time-stamp: <2008-02-12 14:39:18 haines>
3	3	"""Utilities to help data processing
4	4
…	…
7	7	TO DO:
8	8	check_configs()
	9	unit conversions (udunits?)
9	10	"""
10	11
…	…
277	278	return str
278	279
	280	# unit conversions
	281	def meters2feet(meters):
	282	"""Convert meters to feet: <feet> = <meters>*3.28084 """
	283	return meters*3.28084
	284
279	285	#

r101	r104
1	1	#!/usr/bin/env python
2		# Last modified: Time-stamp: <2008-0~~1-14 11:03:16~~ haines>
	2	# Last modified: Time-stamp: <2008-02-12 14:43:50 haines>
3	3	"""Process raw data to monthly netCDF data files
4	4
…	…
31	31
32	32	# define config file location to run under cron
33		defconfigs='/~~afs/isis.unc.edu/depts/marine/workspace/haines/nc-~~coos/raw2proc'
	33	defconfigs='/home/haines/nccoos/raw2proc'
34	34
35	35	import numpy
…	…
348	348	data = parse(pi, si, lines)
349	349	# determine which index of data is within the specified timeframe (usually the month)
350		data['in'] = data['dt']>si['proc_start_dt'] and data['dt']<=si['proc_end_dt']
	350	# data['in'] = data['dt']>si['proc_start_dt'] and data['dt']<=si['proc_end_dt']
	351	n = len(data['dt'])
	352	data['in'] = numpy.array([False for i in range(n)])
	353	for index, val in enumerate(data['dt']):
	354	if val>si['proc_start_dt'] and val<=si['proc_end_dt']:
	355	data['in'][index] = True
	356
351	357	# if any records are in the month then write to netcdf
352	358	if data['in'].any():