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)