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]
