timeRefUnix = timeRefUnix / 10.**9
try:
rootgrpOut = writeData.createNetCdf(outPutFilePath, prefixL1)
time_ref = writeData.createTimeDimension(rootgrpOut, timeRefUnix, prefixL1)
range_ref = writeData.createRangeDimension(rootgrpOut, rangeRef, prefixL1)
timeDeviation = np.array(timeDeviation.astype(np.float32))
time_dev = writeData.createDeviation(rootgrpOut, timeDeviation,
'delta_time', radar, prefixL1)
rangeDeviation = np.array(rangeDeviation.astype(np.float32))
range_dev = writeData.createDeviation(rootgrpOut, rangeDeviation,
'delta_altitude', radar, prefixL1)
if varNameOutput == 'Ze_' + radar:
resampledTimeRange = resampledTimeRange + zeOffset
resampledTimeRange = np.array(resampledTimeRange.T.astype(np.float32))
var_resampled = writeData.createVariable(rootgrpOut, resampledTimeRange,
varFinalName, varNameOutput,
radar, prefixL1, np.float32)
except:
#except IOError:
print 'No Data To Write :('
rootgrpOut.close()
print 'Done %s' % (varNameOutput)
varListName = varNameOut.split('_')
if len(varListName) > 1:
varFinalName = '_'.join(varListName[:-1])
radar = varListName[-1]
else:
varFinalName = varListName[0]
radar = ''
dataDF = variableOutPut[varNameOut]['data']
if varFinalName == 'quality_flag_offset':
dataToWrite = np.array(dataDF.astype(np.uint16))
var_Written = writeData.createVariable(rootgrpOut, dataToWrite,
varFinalName, varNameOut, radar,
prefixL2, np.uint16)
else:
dataToWrite = np.array(dataDF.astype(np.float32))
var_Written = writeData.createVariable(rootgrpOut, dataToWrite,
varFinalName, varNameOut, radar,
prefixL2, np.float32)
#It writes the data from L1 in L2 file
for varNameOut in sorted(variableToCopy.keys()):
#it removes the noise from v_Ka
if varNameOut == 'rv_ka':
#indexV_Ka = variableToCopy.keys().index('rv_ka')
resampledTime = resampledTime.T
#resample in range
resampledTime['ranges'] = ranges
rangeIndexList = trLib.getIndexList(resampledTime, rangeRef, rangeTolerance)
emptyDataFrame = pd.DataFrame(index=rangeRef, columns=resampledTime.columns)
resampledTimeRange = trLib.getResampledDataPd(emptyDataFrame, resampledTime,
rangeIndexList)
rangeDeviation = trLib.getDeviationPd(rangeRef, resampledTimeRange)
del resampledTimeRange['ranges']
#To write the data
rootgrpOut = writeData.createNetCdf(outPutFilePath)
time_ref = writeData.createTimeDimension(rootgrpOut, timeRefUnix)
range_ref = writeData.createRangeDimension(rootgrpOut, rangeRef)
timeDeviation = np.array(timeDeviation.astype(np.float32))
time_dev = writeData.createDeviation(rootgrpOut, timeDeviation, 'delta_time',
radar)
rangeDeviation = np.array(rangeDeviation.astype(np.float32))
range_dev = writeData.createDeviation(rootgrpOut, rangeDeviation,
'delta_altitude', radar)
resampledTimeRange = np.array(resampledTimeRange.T.astype(np.float32))
var_resampled = writeData.createVariable(rootgrpOut, resampledTimeRange,
varFinalName, varNameOutput, radar)
rootgrpOut.close()
print 'Done %s' % (varNameOutput)