|
| 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 | lon_a = 135.2207105 |
| |
| float | lat_a = 34.6533937 |
| |
| float | lon_b = 139.7671248 |
| |
| float | lat_b = 35.6812362 |
| |
| int | dist0 = 20 |
| |
| list | xticks = [130,135,140] |
| |
| list | yticks = [30,35,40] |
| |
| | xticks_f = np.array([tmp for tmp in np.arange(np.floor(lon_a)-2,np.ceil(lon_b)+3,2)]) |
| |
| | yticks_f = np.array([tmp for tmp in np.arange(np.floor(lat_a)-2,np.ceil(lat_b)+3,2)]) |
| |
| | 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_ab = ll2dist(lat_a,lon_a,lat_b,lon_b) |
| |
| | dist_a = ll2dist(lat_a,lon_a,lat,lon) |
| |
| | dist_b = ll2dist(lat_b,lon_b,lat,lon) |
| |
| list | arg_lim_list = [] |
| |
| list | dist_along_ab_list = [] |
| |
| | fig = plt.figure(figsize=(14,6)) |
| |
| list | axes |
| |
| | color |
| |
| | levels1 = np.linspace(0,500,11) |
| |
| | xlocs |
| |
| | ylocs |
| |
| | crs |
| |
| | cax = ax.contour(lon,lat,dist_a,levels=levels1,cmap="Reds_r") |
| |
| | cbar = plt.colorbar(cax,ax=ax,fraction=0.02) |
| |
| | distsum = dist_a+dist_b |
| |
| tuple | arg_lim = (distsum.values<dist_ab+dist0)&(dist_a.values<dist0*(tmp+1))&(dist_a.values>dist0*tmp) |
| |
| | c = np.mean(uppvar_limmean,axis=1)*uppvarfact |
| |
| | s |
| |
| | cmap |
| |
| | vmin |
| |
| | vmax |
| |
| | marker |
| |
| | sfcvar = xrvar[sfcvarname] |
| |
| | uppvar = xrvar[uppvarname] |
| |
| list | sfcvar_limmean_list = [] |
| |
| list | uppvar_limmean_list = [] |
| |
| | sfcvar_limmean = np.array(sfcvar_limmean_list) |
| |
| | uppvar_limmean = np.array(uppvar_limmean_list) |
| |
| | ax = fig.add_subplot(111,projection=ccrs.PlateCarree()) |
| |
| | figsize |
| |
| int | y = np.mean(sfcvar_limmean,axis=1)*sfcvarfact |
| |
| | linestyle |
| |
水平指定線分近傍の地表面変数・高層大気変数
指定した水平2地点を繋ぐ線分の近くのデータを抽出し, その線分上の地表面変数・高層大気変数の変化を描画
- Author
- Tomoro Yanase, Team SCALE
- Note
