if dim_name is None: dim_name='dim_%d' % dim_length try: ncfile.createDimension(dim_name,dim_length) except: pass return dim_name def add_var(ncfile,var_name,var_val): the_dim='dim_%d' % var_val.shape[0] new_var=ncfile.createVariable(var_name,var_val.dtype,(the_dim,'var_cols')) new_var[:,:]=var_val[...] return new_var dimnames=[add_dim(ncfile,item[0]) for item in keep_arrays] ncfile.createDimension('var_cols',6) varnames=[add_var(ncfile,item[2],item[1]) for item in keep_arrays] col_names=array_list[0][2].dtype.names col_units=[array_list[0][2].dtype.fields[name][2] for name in col_names] col_units=[item.split('_')[1] for item in col_units] ncfile.units=','.join(col_units) ncfile.col_names=','.join(col_names) ## array_lengths.sort() ## array_lengths=set(array_lengths) ncfile.close()
keep_arrays.append((np_array.shape[0],np_array,date_name)) array_lengths.append(np_array.shape[0]) dimnames=[add_dim(ncfile,item[0]) for item in keep_arrays] ncfile.createDimension('var_cols',8) varnames=[add_var(ncfile,item[2],item[1]) for item in keep_arrays] col_names=array_list[0][2].dtype.names a=[(item[0],item[1]) for item in array_list[0][2].dtype.fields.items()] col_names=[item[0] for item in a] col_dtypes=[item[1] for item in a] col_units=[item[2] for item in col_dtypes] print "col_dtypes: ",col_dtypes print "col_units: ",col_units col_names=('press','height','temp','dewpt','relh','mixr','drct',\ 'sknt') titles=('hPa','height_m','temp_C','dewpt_C','relh_percent','mixr_gperkg',\ 'deg','knot') ncfile.units=','.join(titles) ncfile.col_names=','.join(col_names) ## array_lengths.sort() ## array_lengths=set(array_lengths) ncfile.close()