In [1]:

    

import os
import sys
import platform
import matplotlib

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.dates as dates
import matplotlib.ticker as tick
from matplotlib.backends.backend_pdf import PdfPages
from datetime import datetime,timedelta
from pylab import rcParams
rcParams['figure.figsize'] = 15, 10
%matplotlib inline


    

In [2]:

    

# custom functions for transducer data import
import Snake_Valley_Data_Import as svdi


    

In [ ]:

    

if not os.path.exists(checkFolder):
    print('Creating Check Folder')
    os.makedirs(checkFolder)
else:
    print('Check Folder Exists')


    

In [ ]:

    

def UDless(S,r,trans,time):
    return (S*r**2)/(4*Trans*time)


    

In [10]:

    

def WTheis(u):
    if u > 10:
        wth = 0
    elif u < 0.01:
        wth = -0.577216 - np.log(u)
    else:
        wth = -0.577216 - np.log(u)
    
        delta = 1
        I = 0
        Fact = 1
        while abs(delta) > 1E-99:
            I = I + 1
            Fact = Fact * I
            delta = (-1) ** (I + 1) * u ** I / I / Fact
            wth = wth + delta
    return wth


    

In [15]:

    

WTheis(3)


    

    
Out[15]:





-1.6758282886681097