<jupyter_start><jupyter_text>Deterministic point jet<jupyter_code>import matplotlib.pyplot as plt
import pandas as pd

import numpy as np
import matplotlib.pylab as pl<jupyter_output><empty_output><jupyter_text>\begin{equation}    \partial_t \zeta = \frac{\zeta_{jet}}{\tau} - \mu \zeta + \nu_\alpha \nabla^{2\alpha} - \beta \partial_x \psi - J(\psi, \zeta) \zeta\end{equation} Here $\zeta_{jet}$ is the profile of a prescribed jet and $\tau = 1/\mu$ is the relaxation time. Also $\beta = 2\Omega cos \theta$ and $\theta = 0$ corresponds to the equator. Parameters> $\Omega = 2\pi$> $\mu = 0.05$> $\nu = 0$> $\nu_4 = 0$> $\Xi = 1.0$> $\Delta \theta = 0.1$> $\tau$ = 20 days<jupyter_code>dn = "pointjet/8x8/"
# Parameters: μ = 0.05, τ = 20.0, Ξ = 1.0*Ω and Δθ = 0.1
M = 8
N = 8

colors = pl.cm.nipy_spectral(np.linspace(0,1,M))<jupyter_output><empty_output><jupyter_text>NL v GQL v GCE2<jupyter_code>nl = np.load(dn+"nl.npz",allow_pickle=True) 
gql = np.load(dn+"gql.npz",allow_pickle=True) 
gce2 = np.load(dn+"gce2.npz",allow_pickle=True) 
# fig,ax = plt.subplots(1,2,figsize=(14,5))
    
# # Energy
# ax[0].plot(nl['t'],nl['Etav'],label=r'$\langle NL \rangle$')
# ax[0].plot(gql['t'],gql['Etav'],label=r'$\langle GQL(M) \rangle$')
# ax[0].plot(gce2['t'],gce2['Etav'],label=r'$\langle GCE2(M) \rangle$')

# ax[0].set_xlabel(r'$t$',fontsize=14)
# ax[0].set_ylabel(r'$E$',fontsize=14)
# # ax[0].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

# # Enstrophy
# ax[1].plot(nl['t'],nl['Ztav'],label=r'$\langle NL \rangle$')
# ax[1].plot(gql['t'],gql['Ztav'],label=r'$\langle GQL(M) \rangle$')
# ax[1].plot(gce2['t'],gce2['Ztav'],label=r'$\langle GCE2(M) \rangle$')


# ax[1].set_xlabel(r'$t$',fontsize=14)
# ax[1].set_ylabel(r'$Z$',fontsize=14)
# ax[1].legend(loc=4,fontsize=14)

# plt.show()
fig,ax = plt.subplots(1,3,figsize=(14,5))

ax[0].set_title(f'NL')
for i,x in enumerate(nl['Emt'].T):    
    ax[0].plot(nl['t'],x,label=i,c=colors[i])

ax[1].set_title(f'GQL(M)')
for i,x in enumerate(gql['Emt'].T):
    ax[1].plot(gql['t'],x,label=i,c=colors[i])

ax[2].set_title(f'GCE2(M)')
for i,x in enumerate(gce2['Emt'].T):
    ax[2].plot(gce2['t'],x,label=i,c=colors[i])

for a in ax:
    
    a.set_xlabel(r'$t$',fontsize=14)
    a.set_yscale('log')
    a.set_ylim(1e-12,1e2)

ax[0].set_ylabel(r'$E(m)$',fontsize=14)
ax[2].legend(bbox_to_anchor=(1.01,0.85),ncol=1)

plt.show()
# fig,ax = plt.subplots(1,3,figsize=(16,4))

# im = ax[0].imshow((nl['Vxy'][:,:,0]),interpolation="bicubic",cmap="RdBu_r",origin="lower")
# fig.colorbar(im, ax=ax[0])
# ax[0].set_title(r'NL: $\zeta(x,y,t = 0 )$',fontsize=14)

# im = ax[1].imshow((gql['Vxy'][:,:,0]),interpolation="bicubic",cmap="RdBu_r",origin="lower")
# fig.colorbar(im, ax=ax[1])
# ax[1].set_title(r'GQL: $\zeta(x,y,t = 0 )$',fontsize=14)

# im = ax[2].imshow((gce2['Vxy'][:,:,0]),interpolation="bicubic",cmap="RdBu_r",origin="lower")
# fig.colorbar(im, ax=ax[2])
# ax[2].set_title(r'GCE2:$\zeta(x,y,t = 0 )$',fontsize=14)

# for a in ax:
#     a.set_xticks([0,M-1,2*M-2])
#     a.set_xticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)
#     a.set_yticks([0,M-1,2*M-2])
#     a.set_yticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)


# plt.show()
fig,ax = plt.subplots(1,3,figsize=(16,4))

im = ax[0].imshow((nl['Vxy'][:,:,-1]),interpolation="bicubic",cmap="RdBu_r",origin="lower")
fig.colorbar(im, ax=ax[0])
ax[0].set_title(r'NL: $\zeta(x,y,t = 0 )$',fontsize=14)

im = ax[1].imshow((gql['Vxy'][:,:,-1]),interpolation="bicubic",cmap="RdBu_r",origin="lower")
fig.colorbar(im, ax=ax[1])
ax[1].set_title(r'GQL: $\zeta(x,y,t = 0 )$',fontsize=14)

im = ax[2].imshow((gce2['Vxy'][:,:,-1]),interpolation="bicubic",cmap="RdBu_r",origin="lower")
fig.colorbar(im, ax=ax[2])
ax[2].set_title(r'GCE2:$\zeta(x,y,t = 0 )$',fontsize=14)

for a in ax:
    
    a.set_xticks([0,M-1,2*M-2])
    a.set_xticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)
    a.set_yticks([0,M-1,2*M-2])
    a.set_yticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)


plt.show()<jupyter_output><empty_output><jupyter_text>QL v CE2 v GCE2(0)<jupyter_code>ql = np.load(dn+"ql.npz",allow_pickle=True) 
ce2 = np.load(dn+"ce2.npz",allow_pickle=True)
gce2 = np.load(dn+"gce2_0.npz",allow_pickle=True)
fig,ax = plt.subplots(1,3,figsize=(18,5))
    
# Energy
ax[0].set_title(f'QL',fontsize=14)
for i,x in enumerate(ql['Emtav'].T):
    ax[0].plot(ql['t'],x,label=i,c=colors[i])

ax[1].set_title(f'CE2',fontsize=14)
for i,x in enumerate(ce2['Emtav'].T):
    ax[1].plot(ce2['t'],x,label=i,c=colors[i])

ax[2].set_title(f'GCE2(0)',fontsize=14)
for i,x in enumerate(gce2['Emtav'].T):
    ax[2].plot(gce2['t'],x,label=i,c=colors[i])

ax[2].legend(bbox_to_anchor=(1.01,0.5),ncol=1)


for a in ax:
    a.set_xlabel(r'$t$',fontsize=14)
    a.set_ylabel(r'$E$',fontsize=14)
    a.set_yscale('log')
    a.set_ylim(1e-12,1e2)
    

# plt.show()
# plt.savefig(dn+"ze_tau20_qlce2gce2_0.png",bbox_inches='tight',)<jupyter_output><empty_output><jupyter_text>GQL(1) v GCE2(1)<jupyter_code>gql = np.load(dn+"gql_1.npz",allow_pickle=True) 
gce2 = np.load(dn+"gce2_1.npz",allow_pickle=True) 
fig,ax = plt.subplots(1,2,figsize=(14,5))
    
# Energy
ax[0].set_title(f'GQL(1)',fontsize=14)
for i,x in enumerate(gql['Emtav'].T):
    ax[0].plot(gql['t'],x,label=i,c=colors[i])

ax[1].set_title(f'GCE2(1)',fontsize=14)
for i,x in enumerate(gce2['Emtav'].T):
    ax[1].plot(gce2['t'],x,label=i,c=colors[i])

for a in ax:
    a.set_xlabel(r'$t$',fontsize=14)
    a.set_ylabel(r'$E$',fontsize=14)
    a.set_yscale('log')
#     a.set_ylim(1e-1,1e0)
ax[1].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

plt.show()
fig,ax = plt.subplots(1,2,figsize=(15,6))

ax[0].set_title(f'GQL(1)',fontsize=14)
im = ax[0].imshow((gql['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
fig.colorbar(im,ax=ax[0])

ax[1].set_title(f'GCE2(1)',fontsize=14)
im = ax[1].imshow((gce2['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
fig.colorbar(im,ax=ax[1])

for a in ax:
    a.set_xticks([0,M-1,2*M-2])
    a.set_xticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)
    a.set_yticks([0,M-1,2*M-2])
    a.set_yticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)

plt.show()
fig,ax = plt.subplots(1,2,figsize=(12,6))

ax[0].set_title('GQL(1)')
im = ax[0].imshow((gql['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im)

ax[1].set_title('GCE2(1)')
im = ax[1].imshow((gce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im)

for a in ax:
    a.set_xticks([0,M-1,2*M-2])
    a.set_xticklabels([r'$-M$',r'$0$',r'$M$'],fontsize=14)
    a.set_yticks([0,M-1,2*M-2])
    a.set_yticklabels([r'$-M$',r'$0$',r'$M$'],fontsize=14)

plt.show()
gql = np.load(dn+"gql_3.npz",allow_pickle=True) 
gce2 = np.load(dn+"gce2_3.npz",allow_pickle=True) 
fig,ax = plt.subplots(1,2,figsize=(14,5))
    
# Energy
ax[0].set_title(f'GQL(3)',fontsize=14)
for i,x in enumerate(gql['Emtav'].T):
    ax[0].plot(gql['t'],x,label=i,c=colors[i])

ax[1].set_title(f'GCE2(3)',fontsize=14)
for i,x in enumerate(gce2['Emtav'].T):
    ax[1].plot(gce2['t'],x,label=i,c=colors[i])

for a in ax:
    a.set_xlabel(r'$t$',fontsize=14)
    a.set_ylabel(r'$E$',fontsize=14)
    a.set_yscale('log')
#     a.set_ylim(1e-1,1e0)

ax[1].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

plt.show()
fig,ax = plt.subplots(1,2,figsize=(15,6))

ax[0].set_title(f'GQL(3)',fontsize=14)
im = ax[0].imshow((gql['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
fig.colorbar(im,ax=ax[0])

ax[1].set_title(f'GCE2(3)',fontsize=14)
im = ax[1].imshow((gce2['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
fig.colorbar(im,ax=ax[1])

for a in ax:
    a.set_xticks([0,M-1,2*M-2])
    a.set_xticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)
    a.set_yticks([0,M-1,2*M-2])
    a.set_yticklabels([r'$0$',r'$\pi$',r'$2\pi$'],fontsize=14)

plt.show()
fig,ax = plt.subplots(1,2,figsize=(12,6))

ax[0].set_title('GQL(3)')
im = ax[0].imshow((gql['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im)

ax[1].set_title('GCE2(3)')
im = ax[1].imshow((gce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im)

for a in ax:
    a.set_xticks([0,M-1,2*M-2])
    a.set_xticklabels([r'$-M$',r'$0$',r'$M$'],fontsize=14)
    a.set_yticks([0,M-1,2*M-2])
    a.set_yticklabels([r'$-M$',r'$0$',r'$M$'],fontsize=14)

plt.show()
# gql = np.load(dn+"gql_5.npz",allow_pickle=True) 
# gce2 = np.load(dn+"gce2_5.npz",allow_pickle=True) 
# fig,ax = plt.subplots(1,2,figsize=(14,5))
    
# # Energy
# ax[0].set_title(f'GQL(5)',fontsize=14)
# for i,x in enumerate(gql['Emtav'].T):
#     ax[0].plot(gql['t'],x,label=i,c=colors[i])

# ax[1].set_title(f'GCE2(5)',fontsize=14)
# for i,x in enumerate(gce2['Emtav'].T):
#     ax[1].plot(gce2['t'],x,c=colors[i])

# for a in ax:
#     a.set_xlabel(r'$t$',fontsize=14)
#     a.set_ylabel(r'$E$',fontsize=14)
#     a.set_yscale('log')
# #     a.set_ylim(1e-1,1e0)
# # ax[1].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,2,figsize=(12,6))

# ax[0].set_title('GQL(5)')
# im = ax[0].imshow((gql['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# # fig.colorbar(im)

# ax[1].set_title('GCE2(5)')
# im = ax[1].imshow((gce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# # fig.colorbar(im)

# for a in ax:
#     a.set_xticks([0,M-1,2*M-2])
#     a.set_xticklabels([r'$-M$',r'$0$',r'$M$'],fontsize=14)
#     a.set_yticks([0,M-1,2*M-2])
#     a.set_yticklabels([r'$-M$',r'$0$',r'$M$'],fontsize=14)

# plt.show()
# dn = "pointjet/12x12/"
# Parameters: μ = 0.05, τ = 20.0, Ξ = 1.0*Ω and Δθ = 0.1
Nx = 12
Ny = 12
nl = np.load(dn+"nl.npz",allow_pickle=True) 
gql = np.load(dn+"gql.npz",allow_pickle=True) 
gce2 = np.load(dn+"gce2.npz",allow_pickle=True) 
# fig,ax = plt.subplots(1,2,figsize=(14,5))
    
# # Energy
# ax[0].plot(nl['t'],nl['Etav'],label=r'$\langle NL \rangle$')
# ax[0].plot(gql['t'],gql['Etav'],label=r'$\langle GQL(M) \rangle$')
# ax[0].plot(gce2['t'],gce2['Etav'],label=r'$\langle GCE2(M) \rangle$')

# ax[0].set_xlabel(r'$t$',fontsize=14)
# ax[0].set_ylabel(r'$E$',fontsize=14)
# # ax[0].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

# # Enstrophy
# ax[1].plot(nl['t'],nl['Ztav'],label=r'$\langle NL \rangle$')
# ax[1].plot(gql['t'],gql['Ztav'],label=r'$\langle GQL(M) \rangle$')
# ax[1].plot(gce2['t'],gce2['Ztav'],label=r'$\langle GCE2(M) \rangle$')


# ax[1].set_xlabel(r'$t$',fontsize=14)
# ax[1].set_ylabel(r'$Z$',fontsize=14)
# ax[1].legend(loc=4,fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,3,figsize=(14,5))

# ax[0].set_title(f'NL')
# for i,x in enumerate(nl['Emt'].T):    
#     ax[0].plot(nl['t'],x,label=i,c=colors[i])

# ax[1].set_title(f'GQL(M)')
# for i,x in enumerate(gql['Emt'].T):
#     ax[1].plot(gql['t'],x,label=i,c=colors[i])

# ax[2].set_title(f'GCE2(M)')
# for i,x in enumerate(gce2['Emt'].T):
#     ax[2].plot(gce2['t'],x,label=i,c=colors[i])

# for a in ax:
    
#     a.set_xlabel(r'$t$',fontsize=14)
#     a.set_yscale('log')
# #     a.set_ylim(1e-10,1e1)

# ax[0].set_ylabel(r'$E(m)$',fontsize=14)
# ax[2].legend(bbox_to_anchor=(1.01,0.85),ncol=1)

# plt.show()
# fig,ax = plt.subplots(1,3,figsize=(24,6))

# ax[0].set_title(f'NL',fontsize=14)
# im = ax[0].imshow((nl['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'GQL(M)',fontsize=14)
# im = ax[1].imshow((gql['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# ax[2].set_title(f'GCE2(M)',fontsize=14)
# im = ax[2].imshow((gce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[2])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,3,figsize=(24,6))

# ax[0].set_title(f'NL',fontsize=14)
# im = ax[0].imshow((nl['Vxy'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'GQL(M)',fontsize=14)
# im = ax[1].imshow((gql['Vxy'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# ax[2].set_title(f'GCE2(M)',fontsize=14)
# im = ax[2].imshow((gce2['Vxy'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[2])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()
ql = np.load(dn+"ql.npz",allow_pickle=True) 
ce2 = np.load(dn+"ce2.npz",allow_pickle=True)
gce2 = np.load(dn+"gce2_0.npz",allow_pickle=True)
# fig,ax = plt.subplots(1,3,figsize=(14,5))
    
# # Energy
# ax[0].set_title(f'QL',fontsize=14)
# for i,x in enumerate(ql['Emtav'].T):
#     ax[0].plot(ql['t'],x,label=i,c=colors[i])

# ax[1].set_title(f'CE2',fontsize=14)
# for i,x in enumerate(ce2['Emtav'].T):
#     ax[1].plot(ce2['t'],x,c=colors[i])

# ax[2].set_title(f'GCE2(0)',fontsize=14)
# for i,x in enumerate(gce2['Emtav'].T):
#     ax[2].plot(gce2['t'],x,c=colors[i])


# for a in ax:
#     a.set_xlabel(r'$t$',fontsize=14)
#     a.set_ylabel(r'$E$',fontsize=14)
#     a.set_yscale('log')
# #     a.set_ylim(1e-1,1e0)
# # ax[1].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,3,figsize=(24,6))

# ax[0].set_title(f'QL',fontsize=14)
# im = ax[0].imshow((ql['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'CE2',fontsize=14)
# im = ax[1].imshow((ce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# ax[2].set_title(f'GCE2(0)',fontsize=14)
# im = ax[2].imshow((gce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[2])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,3,figsize=(24,6))

# ax[0].set_title(f'QL',fontsize=14)
# im = ax[0].imshow((ql['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'CE2',fontsize=14)
# im = ax[1].imshow((ce2['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# ax[2].set_title(f'GCE2(0)',fontsize=14)
# im = ax[2].imshow((gce2['Vxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[2])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,3,figsize=(24,6))

# ax[0].set_title(f'QL',fontsize=14)
# im = ax[0].imshow((ql['Uxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'CE2',fontsize=14)
# im = ax[1].imshow((ce2['Uxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# ax[2].set_title(f'GCE2(0)',fontsize=14)
# im = ax[2].imshow((gce2['Uxy'][:,:,-1]),cmap="RdBu_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[2])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()
gql = np.load(dn+"gql_3.npz",allow_pickle=True) 
gce2 = np.load(dn+"gce2_3.npz",allow_pickle=True) 
# fig,ax = plt.subplots(1,2,figsize=(14,5))
    
# # Energy
# ax[0].set_title(f'GQL',fontsize=14)
# for i,x in enumerate(gql['Emtav'].T):
#     ax[0].plot(gql['t'],x,label=i,c=colors[i])

# ax[1].set_title(f'GCE2',fontsize=14)
# for i,x in enumerate(gce2['Emtav'].T):
#     ax[1].plot(gce2['t'],x,c=colors[i])

# for a in ax:
#     a.set_xlabel(r'$t$',fontsize=14)
#     a.set_ylabel(r'$E$',fontsize=14)
#     a.set_yscale('log')
# #     a.set_ylim(1e-1,1e0)
# # ax[1].legend(bbox_to_anchor=(1.01,0.85),fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,2,figsize=(15,6))

# ax[0].set_title(f'GQL(1)',fontsize=14)
# im = ax[0].imshow((gql['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'GCE2(1)',fontsize=14)
# im = ax[1].imshow((gce2['Emn'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()
# fig,ax = plt.subplots(1,2,figsize=(15,6))

# ax[0].set_title(f'GQL(1)',fontsize=14)
# im = ax[0].imshow((gql['Vxy'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[0])

# ax[1].set_title(f'GCE2(1)',fontsize=14)
# im = ax[1].imshow((gce2['Vxy'][:,:,-1]),cmap="nipy_spectral_r",origin="lower",interpolation="bicubic")
# fig.colorbar(im,ax=ax[1])

# for a in ax:
#     a.set_xticks([0,Nx-1,2*Nx-2])
#     a.set_xticklabels([r'$-N_x$',r'$0$',r'$N_x$'],fontsize=14)
#     a.set_yticks([0,Ny-1,2*Ny-2])
#     a.set_yticklabels([r'$-N_y$',r'$0$',r'$N_y$'],fontsize=14)

# plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.283661373294565, 'fraction_numerical': 0.054965332140460746, 'mean_word_length': 6.527356902356902, 'pattern_counts': {'":': 0, '<': 21, '<?xml version=': 0, '>': 28, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '5', 'dir': 'jupyter-structured-clean-dedup', 'id': '10707', 'license': 'MIT', 'n_tokens_mistral': 9396, 'n_tokens_neox': 9167, 'n_words': 753, 'path': 'validate/pointjet/pointjet.ipynb', 'repo_name': 'gvn22/ZonalFlow.jl', 'size': '17884'}