[502] | 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) |
---|