<jupyter_start><jupyter_text>Regular Theis well<jupyter_code>ml = ModelMaq(kaq=25, z=[20, 0], Saq=S/20, tmin=0.01, tmax=100)
w = DischargeWell(ml, tsandQ=[(0, Q)], rw=1e-5)
ml.solve()
h1 = ml.head(r, 0, t)<jupyter_output>self.neq  0
No unknowns. Solution complete<jupyter_text>Hantush for several different $c$ values<jupyter_code>clist = [1e2, 1e3, 1e4]
hhantush = np.zeros((len(clist), len(t)))
for i, c in enumerate(clist):
    ml = ModelMaq(kaq=25, z=[21, 20, 0], c=c, Saq=S/20, topboundary='semi', tmin=0.01, tmax=100)
    w = DischargeWell(ml, tsandQ=[(0, Q)], rw=1e-5)
    ml.solve()
    hhantush[i] = ml.head(r, 0, t)[0]
plt.semilogx(t, h1[0], label='Theis')
for i in range(3):
    plt.semilogx(t, hhantush[i], label='c=' + str(int(clist[i])))
plt.legend(loc='best')
plt.title('head at r=20');<jupyter_output><empty_output><jupyter_text>Two-aquifer system with same $T$ in both layers and several different $c$ values<jupyter_code>clist = [1e2, 1e3, 1e4]
htwolayer = np.zeros((len(clist), len(t)))
for i, c in enumerate(clist):
    ml = ModelMaq(kaq=[25, 25], z=[20, 0, -5, -25], c=c, Saq=S/20, topboundary='conf', tmin=0.01, tmax=100)
    w = DischargeWell(ml, tsandQ=[(0, Q)], rw=1e-5)
    ml.solve()
    htwolayer[i] = ml.head(r, 0, t)[0]
plt.semilogx(t, h1[0], label='Theis')
for i in range(3):
    plt.semilogx(t, htwolayer[i], label='c=' + str(int(clist[i])))
plt.legend(loc='best')
plt.title('head at r=20');<jupyter_output><empty_output><jupyter_text>Model that starts at $t=10$ days<jupyter_code>ml = ModelMaq(kaq=25, z=[20, 0], Saq=S/20, tmin=0.01, tmax=100, tstart=10)
w = DischargeWell(ml, tsandQ=[(10, Q)], rw=1e-5)
ml.solve()
ht10 = ml.head(r, 0, t)
plt.semilogx(t, ht10[0], label='time shifted in model')
plt.semilogx(t + 10, h1[0], 'r--', label='time shifted by hand')
plt.legend(loc='best')
plt.title('note that heads are not computed at the same times');<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.20732984293193718, 'fraction_numerical': 0.060209424083769635, 'mean_word_length': 5.752650176678445, 'pattern_counts': {'":': 0, '<': 16, '<?xml version=': 0, '>': 16, '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': '4', 'dir': 'jupyter-structured-clean-dedup', 'id': '612487', 'license': 'MIT', 'n_tokens_mistral': 963, 'n_tokens_neox': 884, 'n_words': 182, 'path': 'notebooks/wells_in_different_systems.ipynb', 'repo_name': 'vcantarella/ttim', 'size': '1910'}