Taco Tuesday,gigigigigigi......happy all day!!!
CMIP5-rcp85可视化:
注:部分CMIP6可视化教程请参考以下往期推文
1 《Python基础地图构建(十)》
2 《Python基础地图构建(11)》
3 《Python基础地图构建(27)》
4 《Python基础地图构建(28)》
5 《Python基础地图构建(30)》
import numpy as npimport matplotlib.pyplot as pltimport matplotlib.ticker as mtickerimport xarray as xrimport maskout31import netCDF4 as ncimport cartopy.crs as ccrsimport cartopy.feature as cfeaturefrom cartopy.io.shapereader import Readerfrom cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatterfrom matplotlib import rcParamsconfig = {
"font.family":'Times New Roman', "font.size": 12, "mathtext.fontset":'stix'}rcParams.update(config)region=[70, 140, 15, 55]df= nc.Dataset(r'F:/Rpython/lp28/data/ts_Amon_GISS-E2-H_rcp85_r1i1p1_200601-205012.nc')print(df)lat0=df.variables['lat'][:]lon0=df.variables['lon'][:]z=(df.variables['ts'][:]-273.15)[12]#scale='10m'fig=plt.figure(dpi = 600, figsize = (12, 6))proj=ccrs.PlateCarree()ax = plt.axes(projection = proj)ax.set_extent(region, crs = proj)#ax.coastlines(scale, linewidth = 0.5)ax.set_xticks(np.arange(region[0], region[1] + 1, 10), crs &