From 67a45b0d617469c3899935627cf5160e9d38284f Mon Sep 17 00:00:00 2001
From: alexjamesgarza <alexjamesgarza@gmail.com>
Date: Sun, 3 Mar 2024 16:20:32 -0600
Subject: [PATCH] implemented skeleton logic for solver

---
 ...c_box_species_concentration.cpython-39.pyc | Bin 0 -> 991 bytes
 src/box_model.py                              |  42 +++++++++++++++++-
 2 files changed, 41 insertions(+), 1 deletion(-)
 create mode 100644 src/__pycache__/music_box_species_concentration.cpython-39.pyc

diff --git a/src/__pycache__/music_box_species_concentration.cpython-39.pyc b/src/__pycache__/music_box_species_concentration.cpython-39.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..277677efe28bfdadb58502706e8e329a2f381cc5
GIT binary patch
literal 991
zcmbVL!D`$v5S45vYnCpQ(nBu=K`!y3YbZUGP)ftLlwS7GO)tI(Mbah$ueFtAlA5KL
z=4bi~{S;q%$S35MQ%AOTyM|oaff;G!dGj<6nN22R1a|c6xBMO<bnBCSiy?UjZl6H-
zhy)V&2K+}<AQ3H47IQcs7<B0o&tM}_fdwiefs0r~BECUIA`-wIv`-gn$(6J-U2$1C
z!<^Dpb2lbnpXJ&}3#XPaV*5J6_Y29Da+icTXViJ^q$RfIkXzy|B`j8R#hes7FssC{
z!F1d)o$flEbh&5W5*~*eW+&wA!LHLCWO%pZWvQ79u;C}Y>~0t+`Uc#Kyfwe~63)Qw
zV+b3(K@@Ewh*)5dE^68x@25&R#Y)wH0*qAhim1vu28xov6lNHxOn3<r1qJ-Xr{==`
zo&P^I@UY>q|EPm+FwdfP6x3-agLYZmCRUb9<K;14r=Q3iz;eIHzgTH(&PsWG!A!&Q
z_nOx$|9<k>V)o(H>zDbewu-+v*Vnl<yw50PG;q-Y)pY%>O(|8NE~Wklu?Yl}9^&Tt
zLxN7XwR3pb0Ib+bQreCwU1?F5e*cKluQe;fml35xb4tw!+IB<ig+qUaA$aeYQ=FvT
Z$vnf|x>xNLFdovTZ+f=39(M<h{s6;961@Nb

literal 0
HcmV?d00001

diff --git a/src/box_model.py b/src/box_model.py
index 6e189b1a..f403143d 100644
--- a/src/box_model.py
+++ b/src/box_model.py
@@ -204,7 +204,47 @@ def solve(self):
         """
         # TODO: Implement the logic to solve the box model simulation.
         # Update the internal state of the BoxModel instance to reflect the simulation results.
-        pass
+        
+
+        # TODO: Implement micm config setup
+        #solver = MICM(config)
+
+        #simulation time in minutes
+        sim_length_min = self.box_model_options.simulation_length * 60
+
+        #sets up initial conditions to be current conditions
+        curr_conditions = self.initial_conditions
+
+
+        #sets up next condition if evolving conditions is not empty
+        next_conditions = None
+        next_conditions_time = 0
+        next_conditions_index = 0
+        if(len(self.evolving_conditions.conditions) != 0):
+            next_conditions_index = 0
+            next_conditions = self.evolving_conditions.conditions[0]
+            next_conditions_time = self.evolving_conditions.times[0]
+            
+
+        #runs the simulation at each timestep
+        for curr_time in range(0, sim_length_min + 1, self.box_model_options.chem_step_time):
+
+            #iterates evolvings conditons if enough time has elapsed
+            if(next_conditions != None and next_conditions_time <= curr_time):
+                curr_conditions = next_conditions
+
+                #iterates next_conditions if there are remaining evolving conditions
+                if(len(self.evolving_conditions.conditions) > next_conditions_index + 1):
+                    next_conditions_index += 1
+                    next_conditions = self.evolving_conditions.conditions[next_conditions_index]
+                    next_conditions_time = self.evolving_conditions.times[next_conditions_index]
+                else:
+                    next_conditions = None
+                    
+
+            # TODO: Implement micm solver solve
+            #solver.solve(curr_time, curr_conditions)
+
 
     def readFromJson(self, path_to_json):
         """