def plot_dg_hist_inset(evnum,sim_file,DG_MIN=-30.0,DG_MAX=30.0,NUM_BINS=100):
from pyvc import vcplots
from matplotlib import pyplot as plt
#******************************************
#evnum = 109382
#sim_file = 'ALLCAL2_1-7-11_no-creep_dyn-05_st-20.h5'
#*******************************************
out_file = 'local/dg_field_{}_hist_inset.png'.format(evnum)
# For this vvv to work you must go to vcplots and make this method return the figure
this_fig = vcplots.plot_event_field(sim_file,evnum,output_file=out_file,
save_file_prefix='animation_test_g/field_values/'+str(evnum)+'_',
field_type='gravity')
UNIT = pow(10,-8) #microgals
# Width and height are fixed
ph = 768.0
pw = 1024.0
width_frac = 200.0/pw
height_frac = 150.0/ph
left_frac = 130.0/pw
bottom_frac = 100.0/ph
event_field,mag = get_dg_field(evnum,sim_file)
dg = event_field.dG.flatten()/UNIT
dg_binned,counts = hist_1d(dg,bmin=DG_MIN,bmax=DG_MAX,num_bins=NUM_BINS,norm=True)
LABEL = r'$M={}$'.format(str(round(mag,1)))
hist_inset = this_fig.add_axes((left_frac,bottom_frac,width_frac,height_frac))
hist_inset.plot(dg_binned,counts,color='k')
hist_inset.tick_params(axis='y',bottom='off',top='off',left='off',right='off',
labelbottom='off',labeltop='off',labelleft='off',labelright='off')
hist_inset.tick_params(axis='x', labelsize='small')
#hist_inset.tick_params(axis='x',top='off',bottom='off')
# Make first and last ticks be <MIN and >MAX
#tick_labels = [item.get_text() for item in hist_inset.get_xticklabels()]
#tick_labels[0] = '<'+tick_labels[0]
#tick_labels[-1] = '>'+tick_labels[-1]
#hist_inset.set_xticklabels(tick_labels)
#hist_inset.legend(loc=2,frameon=False,numpoints=1,handlelength=0,handletextpad=0)
hist_inset.patch.set_alpha(0.0)
this_fig.savefig(out_file,dpi=200)
plt.clf()
#start_year = 0.0
#end_year = 100.0
#event_range={'type':'year','filter':(start_year,end_year)}
#vcplots.event_field_animation(sim_file, output_directory, event_range,
# field_type=FIELD, fringes=True, padding=0.08, cutoff=None,
# animation_target_length=LENGTH, animation_fps = FPS, fade_seconds = FADE,
# min_mag_marker = MIN_MAG_MARK, force_plot=PLOT)
#------------------------------------------------------------------------------
#------------------------------------------------------------------------------
#------------------------------------------------------------------------------
"""
#sim_file = '../VCModels/ALLCAL2_1-7-11_no-creep/ALLCAL2_1-7-11_no-creep_dyn-0-5_st-5.h5'
sim_file = '/home/kasey/Okubo-test/ALLCAL2_1-7-11_no-creep_dyn-05_st-20.h5'
evnum = 109382 # Mag = 8.0
#evnum = 1455 # Mag = 7.83
#evnum = 48467 # Mag = 7.99
#evnum = 49261 # Mag = 7.63
#evnum = 54526 # Mag = 7.73
#evnum = 54358 # Mag = 7.79
out_file = '/home/kasey/Okubo-test/dg_plots/dg_'+str(evnum)+'_match.png'
FRINGES = True
FIELD = 'gravity'
LAT_RANGE = (30.65,42.85)
LON_RANGE = (-125.52,-113.24)
#if FIELD == 'gravity':
