diff --git a/HIV_model_2.m b/Legacy Code/HIV_model_2.m
similarity index 96%
rename from HIV_model_2.m
rename to Legacy Code/HIV_model_2.m
index 37ea783..88b7264 100644
--- a/HIV_model_2.m
+++ b/Legacy Code/HIV_model_2.m	
@@ -34,8 +34,8 @@
 %   -> NOTE: parameter pairs are for E and K, respectively
 
 %     P(01) = V i.c. [log(/mL)]
-%     P(02) = E i.c. [/µL]
-%     P(03) = K i.c. [/µL]
+%     P(02) = E i.c. [/ç„¡]
+%     P(03) = K i.c. [/ç„¡]
 
 %     P(04) = N i.c. [pmol/kg]
 %     P(05) = C absorption rate const. [/d]
@@ -44,10 +44,10 @@
 
 %     P(08) = death rate constant [/d]
 %     P(09)
-%     P(10) = maximum concentration [/µL]
+%     P(10) = maximum concentration [/ç„¡]
 %     P(11)
 
-%     P(12) = killing rate constant [µL/d]
+%     P(12) = killing rate constant [ç„¡/d]
 %     P(13)
 %     P(14) = Ce killing modifier []
 %     P(15)
@@ -56,7 +56,7 @@
 %     P(17)
 %     P(18) = constant IL15 factor [pM]
 %     P(19)
-%     P(20) = V-dep IL15 factor [pM*µL]
+%     P(20) = V-dep IL15 factor [pM*ç„¡]
 %     P(21)
 
 %     P(22) = Ct max level []
@@ -71,13 +71,13 @@
 
 %     P(30) = escape mutation rate []
 %     P(31) = initial escape proportion []
-%     P(32) = response activation threshold [/µL]
+%     P(32) = response activation threshold [/ç„¡]
 %     P(33) = response prolif. rate const. [/d]
 
-%     P(34) = T i.c. [/µL]
+%     P(34) = T i.c. [/ç„¡]
 %     P(35) = T death rate constant [/d]
-%     P(36) = T maximum concentration [/µL]
-%     P(37) = T infection rate constant [/µL-d]
+%     P(36) = T maximum concentration [/ç„¡]
+%     P(37) = T infection rate constant [/ç„¡-d]
 %     P(38) = infected T decay rate constant [*(P(35)]
 
 % EquatedPars = vector of indices in 'P' used for equating parameter values
@@ -91,9 +91,9 @@
 
 % Y_MAIN = model output at each time point in 'SolTimes'
 % Y_MAIN(:,1) = V [log(#/mL)]
-% Y_MAIN(:,2) = E [/µL]
-% Y_MAIN(:,3) = K [/µL]
-% Y_MAIN(:,4) = T [/µL]
+% Y_MAIN(:,2) = E [/ç„¡]
+% Y_MAIN(:,3) = K [/ç„¡]
+% Y_MAIN(:,4) = T [/ç„¡]
 
 % Y_EXTRA = cell array of extra outputs (see 'output' section)
 
@@ -131,7 +131,7 @@
 % -------------------------------------------------------------------------
 % Declare initial conditions & calculated parameters ----------------------
 
-V = 10^(P(01)-3) ;% converting V i.c. to [/µL]
+V = 10^(P(01)-3) ;% converting V i.c. to [/ç„¡]
 
 V        = V       *10^(3*mutanum) ;% scaling up V to avoid rounding errors
 P(20:21) = P(20:21)/10^(3*mutanum) ;% scaling parameters to match
@@ -280,7 +280,7 @@
 Y = abs(Y)      ;% removing negatives resulting from numerical error
 
 Y_MAIN(:,1) = log10(sum(Y(:,10:10+mutanum),2))+3 ;% V [log(#/mL)]
-Y_MAIN(:,2:4) = Y(:,7:9)                         ;% E:K:T [/µL]
+Y_MAIN(:,2:4) = Y(:,7:9)                         ;% E:K:T [/ç„¡]
 
 % -------------------------------------------------------------------------
 % Prepare extra outputs 'Y_EXTRA' -----------------------------------------
@@ -385,4 +385,4 @@
     Y_EXTRA{15}{3,1}(n) =  funOn(2) ;
 end
 
-end
\ No newline at end of file
+end