<SYSTEM_TASK:>
Given the following text description, write Python code to implement the functionality described below step by step
<END_TASK>
<USER_TASK:>
Description:
 Compare Two NR Simulations from Different Groups
Step1:  <center>Find Simulations and Load data
Step2:  <center>Recompose the Waveforms
Step3:  <center>Plot the amplitudes to verify correct scaling between GT and SXS waveforms
<ASSISTANT_TASK:>
Python Code:
# Setup ipython environment
%load_ext autoreload
%autoreload 2
%matplotlib inline
# Setup plotting backend
import matplotlib as mpl
mpl.rcParams['lines.linewidth'] = 0.8
mpl.rcParams['font.family'] = 'serif'
mpl.rcParams['font.size'] = 12
mpl.rcParams['axes.labelsize'] = 20
from matplotlib.pyplot import *
# Import useful things
from nrutils import scsearch,gwylm
from numpy import array
Explanation: Compare Two NR Simulations from Different Groups: SXS, GaTech, q=2 Nonspinning
End of explanation
# Search for a simulation. Hopefully the results will be from different codes. NOTE that this could be done more manually so that we don't "hope" but know.
A = scsearch( q=[1,4], nonspinning=True, verbose=True )
# Select which of the search results we wish to keep
U,V = A[77],A[131]
# Load the modes 
u = gwylm(U,lmax=2,verbose=True)
v = gwylm(V,lmax=2,verbose=True)
# Plot the waveforms 
u.plot(); v.plot()
Explanation: <center>Find Simulations and Load data
End of explanation
#
theta,phi = 0,0
a,b = u.recompose(theta,phi,kind='strain'),v.recompose(theta,phi,kind='strain')
Explanation: <center>Recompose the Waveforms
End of explanation
figure( figsize=2*array([5,3]) ) 
plot( a.t - a.intrp_t_amp_max, a.amp )
plot( b.t - b.intrp_t_amp_max, b.amp )
gca().set_yscale("log", nonposy='clip')
ylim([1e-5,1e-1])
xlim([-400,100])
title('the amplitudes should be approx. the same')
a.plot();b.plot();
Explanation: <center>Plot the amplitudes to verify correct scaling between GT and SXS waveforms
End of explanation
<END_TASK>