Plot of radar reflectivity from AUS2200 with a black and white terrain underneath.ΒΆ

Plots terrain height in black and white by reading in the AUS2200 terrain data and then overlays a plot of radar reflectivity (DBz) for a single time.

[9]:

import xarray as xr
import numpy as np
import datetime
import matplotlib.pyplot as plt
import matplotlib.cm as cm
import cmocean.cm as cmo
import os, datetime
from sys import exit, argv
import matplotlib.pyplot as plt
%matplotlib inline

[10]:

from dask.distributed import Client
from glob import glob
import xarray as xr

[11]:

flist=sorted(glob('/g/data/v45/cc6171/experiments/ECL-2016_smallshift/netcdf/u-cs142/2016060*T*Z/aus2200/d0198/RA3/um/umnsa_spec_*.nc'))

[ ]:

client=Client()

[13]:

client

[13]:

Client

Client-77e94147-636a-11ef-8840-000007d7fe80

Connection method: Cluster object Cluster type: distributed.LocalCluster
Dashboard: /proxy/34795/status

Cluster Info

 
[14]:

ds = xr.open_mfdataset(flist,parallel=True,preprocess=lambda x: x['fld_s04i111'])
dbz=ds['fld_s04i111']
#try opening a second dataset containing landmask for plotting coasts
datafile2 = '/g/data/v45/cc6171/experiments/ECL-2016_smallshift/netcdf/u-cs142/20160603T0000Z/aus2200/d0198/RA3/um/umnsaa_pa000.nc'
#bring in landmask from second file
landmask = ds2.fld_s00i030
#bring in topography from second file
topo = ds2.fld_s00i033
#
x = ds.lat
y = ds.lon
# TERRAIN PLOT UNDERNEATH
vmn = -500
vmx = 1500
interval = 100
tcks = np.arange(vmn,vmx+interval,interval)
lvl =  np.arange(vmn,vmx+interval,interval)
colormap = cm.gray
fig,ax=plt.subplots(1,1)
cp = ax.contourf(y,x,topo[1,:,:],cmap=colormap,levels=lvl,vmin=vmn,vmax=vmx,extend='max')
#END TERRAIN PLOT
vmn = 0
vmx = 70
interval = 5
tcks = np.arange(vmn,vmx+interval,interval)
lvl =  np.arange(vmn,vmx+interval,interval)
colormap = cm.hsv_r
cp = ax.contourf(y,x,dbz[100,:,:],cmap=colormap,levels=lvl,vmin=vmn,vmax=vmx,extend='max')
plt.rcParams['lines.linewidth'] = 0.1
cp2 = ax.contour(y,x,landmask[1,:,:],colors='k')
cbar_ax = fig.add_axes([0.25, 0.01, 0.5, 0.04])
cbar = fig.colorbar(cp, cax=cbar_ax, orientation='horizontal', label= 'Dbz', ticks=tcks, shrink=0.8)
plt.show()
../_images/Recipes_AUS2200_radar_terrain_6_0.png 
[ ]: