1 | from common import * |
---|
2 | |
---|
3 | def plot_T850(day, lon,lat,T850): # T850 at final time |
---|
4 | print 'Reading data ...' |
---|
5 | lon, lat, T850 = lon[:], lat[:], T850[day-1, :, :] |
---|
6 | print '... done.' |
---|
7 | plt.figure(figsize=(12,5)) |
---|
8 | plt.contourf(lon,lat,T850) |
---|
9 | plt.colorbar() |
---|
10 | plt.title('T850') |
---|
11 | axis_longitude() |
---|
12 | axis_latitude() |
---|
13 | plt.savefig('T850_day%d.png'%day) |
---|
14 | |
---|
15 | gridfile = 'netcdf/output_dcmip2016_regular.nc' |
---|
16 | nc = cdf.Dataset(gridfile, "r") |
---|
17 | llm, nlon, nlat, ntime = getdims(nc, 'lev','lon','lat','time_counter') |
---|
18 | lon, lat, T850, T, Phi = getvars(nc, 'lon','lat','T850', 'T', 'PHI') |
---|
19 | |
---|
20 | plot_T850(7,lon,lat,T850) |
---|
21 | plot_T850(10,lon,lat,T850) |
---|
22 | plot_T850(30,lon,lat,T850) |
---|