|
| str | dir1 = "./sampledata/scale-5.4.5/scale-rm/test/tutorial/real/experiment/run/" |
| |
| str | domainlabel = "_d01" |
| |
| str | timelabel = "" |
| |
| int | PRC_NUM_X = 2 |
| |
| int | PRC_NUM_Y = 2 |
| |
| str | dir_out = "./fig/" |
| |
| bool | savefig = True |
| |
| str | sfcvarname = "MSLP" |
| |
| str | sfcvarunit = "(hPa)" |
| |
| int | sfcvarfact = 1e-2 |
| |
| str | uppvarname = "RH" |
| |
| str | uppvarunit = "(%)" |
| |
| int | uppvarfact = 1 |
| |
| | uppvarlevels = np.linspace(0,100,21) |
| |
| str | uppvarcmap = "BrBG" |
| |
| float | lon0 = 135.2207105 |
| |
| float | lat0 = 34.6533937 |
| |
| int | dist0 = 50 |
| |
| list | xticks = [130,135,140] |
| |
| list | yticks = [30,35,40] |
| |
| | xticks_f = np.array([tmp for tmp in np.linspace(np.floor(lon0)-2,np.floor(lon0)+2,5)]) |
| |
| | yticks_f = np.array([tmp for tmp in np.linspace(np.floor(lat0)-2,np.floor(lat0)+2,5)]) |
| |
| | fpathlist = get_fpathlist(dir1,"history",domainlabel,timelabel,PRC_NUM_X,PRC_NUM_Y) |
| |
| | xrvar = get_xrvar(fpathlist) |
| |
| | time = xrvar.coords["time"] |
| |
| | tsize = len(time) |
| |
| | lon = xrvar.coords["lon"] |
| |
| | lat = xrvar.coords["lat"] |
| |
| | z = xrvar.coords["z"] |
| |
| | dist = ll2dist(lat0,lon0,lat,lon) |
| |
| | lon_lim = lon.values[dist<dist0] |
| |
| | lat_lim = lat.values[dist<dist0] |
| |
| | argmin_dist = np.argmin(dist.values.reshape(-1)) |
| |
| | lon_nn = lon.values.reshape(-1)[argmin_dist] |
| |
| | lat_nn = lat.values.reshape(-1)[argmin_dist] |
| |
| | fig = plt.figure(figsize=(12,5)) |
| |
| list | axes = [fig.add_subplot(121,projection=ccrs.PlateCarree()), fig.add_subplot(122,projection=ccrs.PlateCarree())] |
| |
| | color |
| |
| | levels1 = np.linspace(0,1000,11) |
| |
| | xlocs |
| |
| | ylocs |
| |
| | crs |
| |
| | cax = ax.contourf(lon,lat,dist,levels=levels1) |
| |
| | cbar = plt.colorbar(cax,ax=ax,fraction=0.02) |
| |
| | marker |
| |
| | s |
| |
| | sfcvar = xrvar[sfcvarname] |
| |
| | sfcvar_nn = sfcvar.values.reshape(time.size,-1)[:,argmin_dist] |
| |
| | sfcvar_lim = sfcvar.values[:,dist<dist0] |
| |
| | uppvar = xrvar[uppvarname] |
| |
| | uppvar_nn = uppvar.values.reshape(time.size,z.size,-1)[:,:,argmin_dist] |
| |
| | uppvar_lim = uppvar.values[:,:,dist<dist0] |
| |
| | ax = axes[0] |
| |
| | label |
| |
| | linewidth |
| |
| | linestyle |
| |
| | alpha |
| |
| | facecolor |
| |
| | fontsize |
| |
| | figsize |
| |
水平指定地点近傍の地表面変数・高層大気変数
モデル出力からある指定した水平地点の近くのデータを抽出・描画
- Author
- Tomoro Yanase, Team SCALE
- Note
