Updated ISARA to include hygroscopic growth impacts on IRI.

ISARA.py

@@ -243,13 +243,12 @@ def Retr_kappa(wvl,
     gf = np.power((1+Kappa[i1]*RH/(100-RH)),1/3)#D/Ddry = (1+kappa*RH/(100-RH))**(1/3)
     D_1 = np.squeeze(np.multiply(gf,dpg1))
     D_2 = np.squeeze(np.multiply(gf,dpg2))
-    #IRIf = (CRI[0]+((gf**3)-1)*IRIw)/(gf**3)
     dpg[1] = D_1
     RRI[1] = (CRI1[0]+((gf**3)-1)*RRIw)/(gf**3)
-    IRI[1] = CRI1[1]
+    IRI[1] = (CRI1[1]+((gf**3)-1)*IRIw)/(gf**3)#CRI1[1]
     dpg[2] = D_2
     RRI[2] = (CRI2[0]+((gf**3)-1)*RRIw)/(gf**3)
-    IRI[2] = CRI2[1]
+    IRI[2] = (CRI2[1]+((gf**3)-1)*IRIw)/(gf**3)
 
     if stop_indx == 0:
       results = mopsmap_wrapper.Model(wvl,size_equ,sd,dpg,RRI,IRI,nonabs_fraction,shape,rho,0,kappa,num_theta,optical_dataset,path_mopsmap_executable)

ISARA_ACTIVATE_Data_Retrieval.py

@@ -27,7 +27,7 @@ def RunISARA():
     sys.path.insert(0, os.path.abspath("../"))  
 
     # Number of cores you want to use
-    number_of_cores = 16    
+    number_of_cores = 32    
 
     # This should be at the start of the code to minimize the fork size
     pool = ProcessPool(ncpus=number_of_cores)   
@@ -82,7 +82,7 @@ def grab_keydata(key_starts_with):
             CRI[io, :] = [RRIp[i2], IRIp[i1]]
             io += 1 
 
-    kappa = np.arange(0.0, 1.40, 0.01).reshape(-1)  
+    kappa = np.arange(0.0, 1.40, 0.001).reshape(-1)  
 
     # set the non-absorbing fraction of the aerosol SD
     nonabs_fraction = 0