multiplicator=1, timestep=None, zlevel=0, dpi=300, drawfunction='contourf', ) (options, args) = parser.parse_args() if len(args)<2: parser.print_help() raise SystemExit bufferzone=options.bufferzone fh = ncf(args[0],'r') if options.gridfilename is not None: gridfh = ncf(options.gridfilename,'r') else: gridfh = fh if bufferzone>0: mbufferzone=-bufferzone else: mbufferzone=None if len(gridfh.variables[options.longitude_name].shape)==3: lon = gridfh.variables[options.longitude_name][0][bufferzone:mbufferzone,bufferzone:mbufferzone] lat = gridfh.variables[options.latitude_name][0][bufferzone:mbufferzone,bufferzone:mbufferzone] elif len(gridfh.variables[options.longitude_name].shape)==1: lon,lat = np.meshgrid(gridfh.variables[options.longitude_name][:], gridfh.variables[options.latitude_name][:])
#!/usr/bin/env python import sys try: from netCDF4 import Dataset as ncf except: from Scientific.IO.NetCDF import NetCDFFile as ncf infile = ncf(sys.argv[1],'r') print len(infile.dimensions[sys.argv[2]])
#!/usr/bin/env python import sys from netCDF4 import Dataset as ncf ff = ncf(sys.argv[1],'a') ff.variables['m2'][:] = 2 ff.sync() ff.close()
year2=m.group(2) if len(year1)>4: year1=year1[:4] if len(year2)>4: year2=year2[:4] else: raise Exception except: year1=int(sys.argv[2]) year2=int(sys.argv[3]) print year1, year2 ff = ncf(inputfilename,'a') ff.createDimension('bnds',2) tb = ff.createVariable('time_bnds','d',('time','bnds')) t = ff.variables['time'] tb.units=t.units tb.calendar=t.calendar t.time_bnds="time_bnds" t.axis="T" t.long_name="time" tbdata = [] tdata = [] for y in range(year1,year2+1): for m in range(1,13): daysinmonth=monthrange(y,m)[1]