outRootGrp = Dataset(outNetCDF_File, 'w', format='NETCDF3_CLASSIC') if (inRootGrp.start_year != startDate.year): raise Exception("Specified simulation start year ({0}) does not match with input vapor pressure data file data year.".format(startDate.year)) exit() # DEBUG: print global attributes of the original file: ## for gAttribute in inRootGrp.ncattrs(): ## print 'Global attribute name:', gAttribute, '=', getattr(inRootGrp,gAttribute) # add global file level attributes to the new netcdf file outRootGrp.start_year = inRootGrp.start_year outRootGrp.data_variable_name = inRootGrp.data_variable_name outRootGrp.data_time_step = inTimeStep outRootGrp.orginal_data_source = 'Daymet Software Version 2.0' outRootGrp.conventions = 'CF-1.0' outRootGrp.modified_data_source = 'CI Water System' outRootGrp.spatial_reference = 'NAD83_UTM_Zone_12N' outRootGrp.datum = 'D_North_America_1983' # get dimension values from the original netcdf file inputTimeVar = inRootGrp.variables['time'] inputXvar = inRootGrp.variables['x'] inputYvar = inRootGrp.variables['y'] inputVpVar = inRootGrp.variables[inRootGrp.data_variable_name] print(inputTimeVar.shape[0]) print(inputXvar.shape[0]) print(inputYvar.shape[0]) print(inputVpVar.shape)
# open the netCDF file in readonly mode based on which we will be creating a new netcdf file inRootGrp = Dataset(inPrecpNetCDF_File, "r", format="NETCDF3_CLASSIC") # open a new blank netcdf file to which we will be writting data outRootGrp = Dataset(outWindNetCDF_File, "w", format="NETCDF3_CLASSIC") # DEBUG: print global attributes of the original file: ## for gAttribute in inRootGrp.ncattrs(): ## print 'Global attribute name:', gAttribute, '=', getattr(inRootGrp,gAttribute) # add global file level attributes to the new netcdf file outRootGrp.start_year = inRootGrp.start_year outRootGrp.data_variable_name = "V" outRootGrp.data_time_step = inRootGrp.data_time_step outRootGrp.orginal_data_source = "CI Water Generated Data Source" outRootGrp.conventions = "CF-1.0" outRootGrp.modified_data_source = "CI Water System" outRootGrp.spatial_reference = "NAD83_UTM_Zone_12N" outRootGrp.datum = "D_North_America_1983" # get dimension values from the original netcdf file inputTimeVar = inRootGrp.variables["time"] inputXvar = inRootGrp.variables["x"] inputYvar = inRootGrp.variables["y"] inputPrecVar = inRootGrp.variables["Prec"] # DEBUG: print dimensions of variables in input netcdf file print(inputTimeVar.shape[0]) print(inputXvar.shape[0]) print(inputYvar.shape[0])