Examples added.

gstinoco · Sep 18, 2023 · c44aaf5 · c44aaf5
1 parent 36f29fe
commit c44aaf5
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/Example_1.py b/Example_1.py
@@ -15,7 +15,7 @@ def fADV(x, y, t, c):
     fun = 0.2*np.exp((-(x-.5-c[0]*t)**2-(y-.5-c[1]*t)**2)/.001)
     return fun
 
-mat = loadmat('Regions/Clouds/CAB_2.mat')
+mat = loadmat('Regions/Clouds/CAB_3.mat')
 p   = mat['p']
 tt  = mat['tt']
 if tt.min() == 1:

diff --git a/GFD.py b/GFD.py
@@ -27,12 +27,13 @@ def __init__(self, points, triangulation, time_interval, time_steps):
         self.tt   = triangulation
         self.t    = time_steps
         self.m    = len(points[:,0])
+        self.L    = np.vstack([[0], [0], [0], [0], [0]])
+        self.K    = np.zeros([self.m,self.m])
         self.T    = np.linspace(time_interval[0],time_interval[1],time_steps)
         self.dt   = self.T[1] - self.T[0]
-        self.L    = np.vstack([[0], [0], [0], [0], [0]])
         self.u_ap = np.zeros([self.m,time_steps])
         self.u_ex = np.zeros([self.m,time_steps])
-        self.K    = np.zeros([self.m,self.m])
+
 
     def Boundary(self, f, coefficient):
         for k in np.arange(self.t):