ISOSIMpy/Calculate.py at master · jvollhueter/ISOSIMpy · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
# -*- coding: utf-8 -*-
"""
Created on Tue Mar 29 15:23:22 2022

@author: Vollhüter
"""

import numpy as np
import Multis


class Tracer:
    """Calculate Tracer."""

    def __init__(self, par, Cin, rain):

        Multis.Multis()

        par.Thalf = par.Thalf_1

        res = Multis.Multis.runMultis(Cin, rain[1], par)

        self.result = res[:rain.shape[0]]


class TracerTracer:
    """Calculate Tracer Tracer."""

    def __init__(self, par, Cin, Cin_2, rain, rain_2, TTs):

        result_tt = np.zeros((2, len(TTs), rain.shape[0]))

        k = 0
        for n in TTs:

            Multis.Multis()

            par.TT = n
            par.Thalf = par.Thalf_1
            res = Multis.Multis.runMultis(Cin, rain[1], par)
            result_tt[0, k] = res[:rain.shape[0]]

            par.Thalf = par.Thalf_2
            res = Multis.Multis.runMultis(Cin_2, rain_2[1], par)
            result_tt[1, k] = res[:rain_2.shape[0]]

            k += 1

        self.result_tt = result_tt


class TriHe:
    """Calculate tritium helium."""

    def __init__(self, par, Cin, rain, TTs):

        result_tt = np.zeros((2, len(TTs), rain.shape[0]))

        k = 0
        for n in TTs:

            Multis.Multis()

            par.TT = n
            par.Thalf = par.Thalf_1
            res = Multis.Multis.runMultis(Cin, rain[1], par)
            result_tt[0, k] = res[:rain.shape[0]]

            par.Thalf = par.Thalf_2
            res = Multis.Multis.runMultis(Cin, rain[1], par)
            result_tt[1, k] = res[:rain.shape[0]]

            k += 1

        self.result_tt = result_tt