Piinshin

import pandas as pd
import os
import matplotlib.pyplot as plt
import numpy as np
#%% Read csv file
files = os.listdir('.')
for file_name in files:
    if file_name.startswith('dat') and file_name.endswith('.csv'):
        dat = pd.read_csv(file_name)
dat_xr = dat.set_index(['speed', 'seed', 'backoff', 'downstream'])

Piinshin

import pandas as pd
import os
import matplotlib.pyplot as plt
import numpy as np
#%% Read csv file
files = os.listdir('.')
for file_name in files:
    if file_name.startswith('dat') and file_name.endswith('.csv'):
        dat = pd.read_csv(file_name)
dat_xr = dat.set_index(['speed', 'seed', 'backoff', 'downstream'])

Piinshin

import pandas as pd
import csv
import os
import matplotlib.pyplot as plt
import numpy as np
#%% Read csv file
files = os.listdir('.')
for file_name in files:
    if file_name.startswith('dat') and file_name.endswith('.csv'):
        dat = pd.read_csv(file_name)

Piinshin

Piinshin

import netCDF4
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import datetime
import calendar
from scipy.interpolate import griddata
from mpl_toolkits.basemap import Basemap
from scipy.interpolate import RegularGridInterpolator
plt.rcParams['figure.dpi'] = 600

Piinshin

Piinshin

Piinshin

import numpy as np
import matplotlib.pyplot as plt

def k(wavelength):
    return (2*np.pi/wavelength)
def c(g,k,h):
    return np.sqrt(g*np.tanh(k*h)/k)
def G(k,h):
    return 2*k*h/np.sinh(2*k*h)
def cg(g,k,h):

Piinshin

Piinshin