RTModel 2.4__tar.gz → 3.1__tar.gz

{rtmodel-2.4 → rtmodel-3.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: RTModel
-Version: 2.4
+Version: 3.1
 Summary: RTModel is a tool for microlensing event interpretation.
 Keywords: Microlensing analysis and fitting
 Author-Email: Valerio Bozza <valboz@sa.infn.it>
@@ -16,11 +16,11 @@ Requires-Dist: pytest
 Requires-Dist: matplotlib
 Requires-Dist: tqdm
 Requires-Dist: tabulate
-Requires-Dist: VBMicrolensing>=4.1.0
+Requires-Dist: VBMicrolensing>=5.0
 Description-Content-Type: text/markdown
 
 # RTModel
-`RTModel` is a package for modeling and interpretation of microlensing events. It uses photometric time series collected from ground and/or space telescopes to propose one or more possible models among the following:
+`RTModel` is a package for modeling and interpretation of microlensing events. It uses photometric and/or astrometric time series collected from ground and/or space telescopes to propose one or more possible models among the following:
 - Single-lens-single-source microlensing (i.e. Paczynski)
 - Single-lens-binary-source microlensing (with or without xallarap)
 - Binary-lens-single-source microlensing (including planetary microlensing, parallax and orbital motion)
@@ -37,11 +37,11 @@ A second subpackage **`RTModel.templates`** helps the user in the visualization
 
 ## Attribution
 
-`RTModel` has been created by Valerio Bozza (University of Salerno) as a product of many years of direct experience on microlensing modeling (see [RTModel webpage](http://www.fisica.unisa.it/GravitationAstrophysics/RTModel.htm)). 
+`RTModel` has been created by Valerio Bozza (University of Salerno) as a product of many years of direct experience on microlensing modeling (see [RTModel webpage](https://projects.phys.unisa.it/GravitationAstrophysics/RTModel.htm)). 
 
 Any scientific use of `RTModel` should be acknowledged by citing the paper [V.Bozza, A&A 688 (2024) 83](https://ui.adsabs.harvard.edu/abs/2024A%26A...688A..83B/abstract), describing all the algorithms behind the code.
 
-We are grateful to Greg Olmschenk, who revised the package installation in order to make it as cross-platform as possible. We also thank all the users who are providing suggestions, reporting bugs or failures: Etienne Bachelet, David Bennett, Jonathan Brashear, Stela Ishitani Silva, Yiannis Tsapras, Keto Zhang.
+We are grateful to Greg Olmschenk, who revised the package installation in order to make it as cross-platform as possible. We also thank all the users who are providing suggestions, reporting bugs or failures: Etienne Bachelet, David Bennett, Jonathan Brashear, Laura Salmeri, Stela Ishitani Silva, Yiannis Tsapras, Sigfried Vanaverbeke, Keto Zhang.
 
 ## Installation
 
@@ -61,9 +61,9 @@ In alternative, you may directly install it from PyPI without cloning this repos
 pip install RTModel
 ```
 
-Currently, `RTModel` works on Linux, Windows and MacOS, requiring Python >= 3.6. 
+Currently, `RTModel` works on Linux, Windows and MacOS, requiring Python >= 3.7. 
 A C++ compiler compatible with C++17 standard is needed for installation.
-`RTModel` also incorporates version 4.1 of [`VBMicrolensing`](https://github.com/valboz/VBMicrolensing). You are encouraged to cite the relevant papers listed in that repository as well.
+`RTModel` uses [`VBMicrolensing`](https://github.com/valboz/VBMicrolensing) for all calculations. You are encouraged to cite the relevant papers listed in that repository as well.
 
 ## Documentation
 Full [documentation for the use of RTModel](docs/README.md) is available.

{rtmodel-2.4 → rtmodel-3.1}/README.md

@@ -1,5 +1,5 @@
 # RTModel
-`RTModel` is a package for modeling and interpretation of microlensing events. It uses photometric time series collected from ground and/or space telescopes to propose one or more possible models among the following:
+`RTModel` is a package for modeling and interpretation of microlensing events. It uses photometric and/or astrometric time series collected from ground and/or space telescopes to propose one or more possible models among the following:
 - Single-lens-single-source microlensing (i.e. Paczynski)
 - Single-lens-binary-source microlensing (with or without xallarap)
 - Binary-lens-single-source microlensing (including planetary microlensing, parallax and orbital motion)
@@ -16,11 +16,11 @@ A second subpackage **`RTModel.templates`** helps the user in the visualization
 
 ## Attribution
 
-`RTModel` has been created by Valerio Bozza (University of Salerno) as a product of many years of direct experience on microlensing modeling (see [RTModel webpage](http://www.fisica.unisa.it/GravitationAstrophysics/RTModel.htm)). 
+`RTModel` has been created by Valerio Bozza (University of Salerno) as a product of many years of direct experience on microlensing modeling (see [RTModel webpage](https://projects.phys.unisa.it/GravitationAstrophysics/RTModel.htm)). 
 
 Any scientific use of `RTModel` should be acknowledged by citing the paper [V.Bozza, A&A 688 (2024) 83](https://ui.adsabs.harvard.edu/abs/2024A%26A...688A..83B/abstract), describing all the algorithms behind the code.
 
-We are grateful to Greg Olmschenk, who revised the package installation in order to make it as cross-platform as possible. We also thank all the users who are providing suggestions, reporting bugs or failures: Etienne Bachelet, David Bennett, Jonathan Brashear, Stela Ishitani Silva, Yiannis Tsapras, Keto Zhang.
+We are grateful to Greg Olmschenk, who revised the package installation in order to make it as cross-platform as possible. We also thank all the users who are providing suggestions, reporting bugs or failures: Etienne Bachelet, David Bennett, Jonathan Brashear, Laura Salmeri, Stela Ishitani Silva, Yiannis Tsapras, Sigfried Vanaverbeke, Keto Zhang.
 
 ## Installation
 
@@ -40,9 +40,9 @@ In alternative, you may directly install it from PyPI without cloning this repos
 pip install RTModel
 ```
 
-Currently, `RTModel` works on Linux, Windows and MacOS, requiring Python >= 3.6. 
+Currently, `RTModel` works on Linux, Windows and MacOS, requiring Python >= 3.7. 
 A C++ compiler compatible with C++17 standard is needed for installation.
-`RTModel` also incorporates version 4.1 of [`VBMicrolensing`](https://github.com/valboz/VBMicrolensing). You are encouraged to cite the relevant papers listed in that repository as well.
+`RTModel` uses [`VBMicrolensing`](https://github.com/valboz/VBMicrolensing) for all calculations. You are encouraged to cite the relevant papers listed in that repository as well.
 
 ## Documentation
 Full [documentation for the use of RTModel](docs/README.md) is available.

{rtmodel-2.4 → rtmodel-3.1}/RTModel/RTModel.py

@@ -32,7 +32,10 @@ class RTModel:
             self.eventname = os.getcwd()
         else:
             self.eventname = os.path.realpath(event)
-            print("Event name: " + self.eventname)
+            if(os.path.exists(self.eventname)):
+                print("Event name: " + self.eventname)
+            else:
+                print("! Invalid path for event: " + self.eventname)
         self.inidir = "ini"
         self.modelcodes = ['PS', 'PX', 'BS', 'BO', 'LS', 'LX', 'LO', 'LK', 'TS', 'TX']
         self.endphase =  len(self.modelcodes)*2+3
@@ -44,15 +47,40 @@ class RTModel:
         self.config_LevMar()
         self.config_ModelSelector()
         self.satellitedir = '.'
-
+        self.astrometric = False
+        self.parameters_ranges = {'PS': [[-11.,1.0, 1.0],[-4.6, 7.6, 1.0],[-300,300,5.0],[-11.5,2.3,2.3]],
+                                  'PX': [[-3.0,3.0, 0.5],[-4.6, 7.6, 1.0],[-300,300,5.0],[-11.5,2.3,2.3],[-3.0,3.0,0.1],[-3.0,3.0,0.1]],
+                                  'BS': [[-4.6,7.6,1.0],[-11.5,0.0,0.5],[0,3.0,0.5],[0,3.0,0.5],[-300,300,1.0],[-300,300,1.0],[-11.5,2.3,2.3]],
+                                  'BO': [[-4.6,7.6,1.0],[-11.5,0.0,0.5],[0,3.0,0.5],[0,3.0,0.5],[-300,300,1.0],[-300,300,1.0],[-11.5,2.3,2.3],
+                                         [-3.0,3.0,0.03],[-3.0,3.0,0.03],[-1.0,1.0,0.01],[-1.0,1.0,0.01],[1.e-7,1.0,0.01]],
+                                  'LS': [[-4.0,3.0,.1],[-16.1,16.1,0.5],[-3.0,3.0,0.1], [-12.56,12.56,0.1],[-11.5,-2.5,0.3],[-4.6,7.6,0.6],
+                                         [-300,300,5.0]],
+                                  'LX': [[-4.0,3.0,.1],[-16.1,16.1,0.5],[-3.0,3.0,0.1], [-12.56,12.56,0.1],[-11.5,-2.5,0.3],[-4.6,7.6,0.6],
+                                         [-300,300,5.0],[-3.0,3.0,0.03],[-3.0,3.0,0.03]],
+                                  'LO': [[-4.0,3.0,.1],[-16.1,16.1,0.5],[-3.0,3.0,0.1], [-12.56,12.56,0.1],[-11.5,-2.5,0.3],[-4.6,7.6,0.6],
+                                         [-300,300,5.0],[-3.0,3.0,0.03],[-3.0,3.0,0.03],[-1.0,1.0,0.01],[-1.0,1.0,0.01],[1.e-7,1.0,0.01]],
+                                  'LK': [[-4.0,3.0,.1],[-16.1,16.1,0.5],[-3.0,3.0,0.1], [-12.56,12.56,0.1],[-11.5,-2.5,0.3],[-4.6,7.6,0.6],
+                                         [-300,300,5.0],[-3.0,3.0,0.03],[-3.0,3.0,0.03],[-1.0,1.0,0.01],[-1.0,1.0,0.01],[1.e-7,1.0,0.01],
+                                        [-10,10,0.1],[0.5001,10,0.1]],
+                                  'TS': [[-4.0,3.0,.1],[-16.1,16.1,0.5],[-3.0,3.0,0.1], [-12.56,12.56,0.1],[-11.5,-2.5,0.3],[-4.6,7.6,0.6],
+                                         [-300,300,5.0],[-4.0,3.0,0.3],[-11.5,11.5,0.5], [-12.56,12.56,0.3]],
+                                  'TX': [[-4.0,3.0,.1],[-16.1,16.1,0.5],[-3.0,3.0,0.1], [-12.56,12.56,0.1],[-11.5,-2.5,0.3],[-4.6,7.6,0.6],
+                                         [-300,300,5.0],[-4.0,3.0,0.3],[-11.5,11.5,0.5], [-12.56,12.56,0.3],[-3.0,3.0,0.03],[-3.0,3.0,0.03]],
+                                  'astrometry': [[-30.0,30.0,1.0],[-30.0,30.0,1.0],[0.05,1.0,0.1],[0.001,30.0,0.2]]
+                                 }		                                  
+	
     def set_processors(self, nprocessors):
         self.nprocessors = nprocessors
 
     def set_event(self, event):
         self.eventname = os.path.realpath(event)
+        if(not os.path.exists(self.eventname)):
+            print("! Invalid path for event: " + self.eventname)
 
     def set_satellite_dir(self, satellitedir):
-        self.satellitedir = satellitedir
+        self.satellitedir = os.path.realpath(satellitedir)
+        if(not os.path.exists(self.satellitedir)):
+            print("! Invalid path for satellite directory: " + self.satellitedir)
 
     def set_constraints(self, constraints = None):
         self.constraints = constraints
@@ -62,7 +90,21 @@ class RTModel:
             for cons in constraints:
                 f.write(cons[0] + ' = '+ str(cons[1]) + ' '+ str(cons[2]) + ' '+ str(cons[3]) + ' ' + '\n')
 
-    def config_Reader(self, tau = 0.1, binning = 4000, otherseasons = 1, renormalize = 1, thresholdoutliers = 10):
+    def set_parameter_ranges(self):
+        if(not os.path.exists(self.eventname + '/' + self.inidir)):
+            os.makedirs(self.eventname + '/' + self.inidir)
+        with open(self.eventname + '/' + self.inidir + '/Parameters_Ranges.ini','w') as f:
+            for modelcode in self.modelcodes:
+                f.write(modelcode + '\n')
+                for par in self.parameters_ranges[modelcode]:
+                    f.write(str(par[0]) + ' ' + str(par[1]) + ' ' + str(par[2]) + '\n')
+            modelcode = 'astrometry'    
+            f.write(modelcode + '\n')
+            for par in self.parameters_ranges[modelcode]:
+                f.write(str(par[0]) + ' ' + str(par[1]) + ' ' + str(par[2]) + '\n')
+        
+
+    def config_Reader(self, tau = 0.1, binning = 4000, otherseasons = 100, renormalize = 1, thresholdoutliers = 10):
         self.Reader_tau= tau # conventional correlation time for consecutive points
         self.Reader_binning = binning # maximum number of points left after re-binning
         self.Reader_otherseasons = otherseasons # How to use other seasons (0 = Yes, 1 = decrease significance, 2 = remove)
@@ -82,6 +124,15 @@ class RTModel:
         print('  Pre-processing data...')
         try:
             completedprocess=subprocess.run([self.bindir+self.readerexe,self.eventname], cwd = self.bindir, stdout=subprocess.PIPE, stderr=subprocess.PIPE, check=True, text = True)
+            with open(self.eventname + '/LCToFit.txt') as f:
+                lines = f.readlines()
+                del(lines[0])
+                self.astrometric = False
+                for line in lines:
+                    if(float(line.split()[6])>0):
+                        self.astrometric = True
+                        print('  Astrometric data found: static models will be skipped')
+                        break                    
             print('  OK')
         except subprocess.CalledProcessError as e:
             print('\033[30;41m! Error in pre-processing. Please check your data!\033[m')
@@ -130,9 +181,9 @@ class RTModel:
         try:
             completedprocess=subprocess.run([self.bindir+self.initcondexe,self.eventname], cwd = self.bindir, stdout=subprocess.PIPE, stderr=subprocess.PIPE, check=True, text = True)
             peaksearch = True
-            i=0
+            imod=0
             while(peaksearch):
-                initfils=glob.glob(self.eventname + '/InitCond/InitCond'+ self.modelcodes[i] + '*')
+                initfils=glob.glob(self.eventname + '/InitCond/InitCond'+ self.modelcodes[imod] + '*')
                 if(len(initfils)!=0):
                     peaksearch = False
                     with open(initfils[0], 'r') as f:
@@ -140,6 +191,7 @@ class RTModel:
                         print('Peaks:  ',end ='')
                         for i in range(0,npeaks):
                             print(f'{float(f.readline().split()[0]):.4f}',end = '  ')
+                imod+=1
             print('\n  OK')
         except subprocess.CalledProcessError as e:
             print('\033[30;41m! Error in setting initial conditions!\033[m')
@@ -148,12 +200,20 @@ class RTModel:
             print('\033[30;41m! Program stopped here!\033[m')
             self.done = True
 
-    def config_LevMar(self, nfits = 6, offsetdegeneracy = 3, timelimit = 600.0, maxsteps = 50, bumperpower = 2.0):
+    def config_LevMar(self, nfits = 6, offsetdegeneracy = 3, timelimit = 600.0, maxsteps = 50, bumperpower = 2.0, \
+                      mass_luminosity_exponent = None, mass_radius_exponent = None, lens_mass_luminosity_exponent = None, \
+                     turn_off_secondary_source = False, turn_off_secondary_lens = False):
         self.LevMar_nfits = nfits # Number of models to be calculated from the same initial condition using the bumper method
         self.LevMar_offsetdegeneracy = offsetdegeneracy # Number of models to be fit after applying offset degeneracy to best model found so far
         self.LevMar_maxsteps = maxsteps # Maximum number of steps in each fit
         self.LevMar_timelimit = timelimit # Maximum time in seconds for total execution
         self.LevMar_bumperpower = bumperpower # Repulsion factor of bumpers
+        self.LevMar_mass_luminosity_exponent = mass_luminosity_exponent # mass-luminosity exponent for binary sources
+        self.LevMar_mass_radius_exponent = mass_radius_exponent # mass-radius exponent for binary sources
+        self.LevMar_lens_mass_luminosity_exponent = lens_mass_luminosity_exponent # mass-luminosity exponent for binary lenses
+        self.LevMar_turn_off_secondary_lens = turn_off_secondary_lens # Option for dark secondary lenses
+        self.LevMar_turn_off_secondary_source = turn_off_secondary_source # Option for dark secondary sources
+        self.LevMar_stepchainsave = False # If True, step chains are saved
     
     def LevMar(self,strmodel, parameters_file = None, parameters = None):
         if(not os.path.exists(self.eventname + '/' + self.inidir)):
@@ -165,14 +225,27 @@ class RTModel:
                 for fl in parameters:
                     line = line + str(fl) + ' '
                 f.write(line)
+        self.set_parameter_ranges()
         with open(self.eventname + '/' + self.inidir + '/LevMar.ini','w') as f:
             f.write('nfits = ' + str(self.LevMar_nfits) + '\n')
             f.write('offsetdegeneracy = ' + str(self.LevMar_offsetdegeneracy) + '\n')
             f.write('maxsteps = ' + str(self.LevMar_maxsteps) + '\n')
             f.write('timelimit = ' + str(self.LevMar_timelimit) + '\n')
             f.write('bumperpower = ' + str(self.LevMar_bumperpower) + '\n')
+            if(self.LevMar_mass_luminosity_exponent != None):
+                f.write('mass_luminosity_exponent = ' + str(self.LevMar_mass_luminosity_exponent) + '\n')
+            if(self.LevMar_mass_radius_exponent != None):
+                f.write('mass_radius_exponent = ' + str(self.LevMar_mass_radius_exponent) + '\n')
+            if(self.LevMar_mass_luminosity_exponent != None):
+                f.write('lens_mass_luminosity_exponent = ' + str(self.LevMar_lens_mass_luminosity_exponent) + '\n')
+            if(self.LevMar_turn_off_secondary_lens):
+                f.write('turn_off_secondary_lens = True\n')
+            if(self.LevMar_turn_off_secondary_source):
+                f.write('turn_off_secondary_source = True\n')
             if(parameters_file != None):
                 f.write('parametersfile = ' + parameters_file)
+            if(self.LevMar_stepchainsave):
+                f.write('stepchainsave = True\n')
         print('- Launching: LevMar')
         print('  Fitting ' + strmodel + ' ...')
         try:
@@ -194,6 +267,18 @@ class RTModel:
             f.write('maxsteps = ' + str(self.LevMar_maxsteps) + '\n')
             f.write('timelimit = ' + str(self.LevMar_timelimit) + '\n')
             f.write('bumperpower = ' + str(self.LevMar_bumperpower) + '\n')
+            if(self.LevMar_mass_luminosity_exponent != None):
+                f.write('mass_luminosity_exponent = ' + str(self.LevMar_mass_luminosity_exponent) + '\n')
+            if(self.LevMar_mass_radius_exponent != None):
+                f.write('mass_radius_exponent = ' + str(self.LevMar_mass_radius_exponent) + '\n')
+            if(self.LevMar_lens_mass_luminosity_exponent != None):
+                f.write('lens_mass_luminosity_exponent = ' + str(self.LevMar_lens_mass_luminosity_exponent) + '\n')
+            if(self.LevMar_turn_off_secondary_lens):
+                f.write('turn_off_secondary_lens = True\n')
+            if(self.LevMar_turn_off_secondary_source):
+                f.write('turn_off_secondary_source = True\n')
+            if(self.LevMar_stepchainsave):
+                f.write('stepchainsave = True\n')
         stringfits = {'PS' : '- Single-lens-Single-source fits',
                       'PX' : '- Single-lens-Single-source fits with parallax',
                       'BS' : '- Single-lens-Binary-source fits',
@@ -204,6 +289,7 @@ class RTModel:
                       'LK' : '- Binary-lens-Single-source fits with eccentric orbital motion',
                       'TS' : '- Triple-lens-Single-source fits',
                       'TX' : '- Triple-lens-Single-source fits with parallax'}       
+        self.set_parameter_ranges()
         print(stringfits[modelcode])
         initcondfile = self.eventname + '/InitCond/' + 'InitCond'+ modelcode + '.txt'
         if(os.path.exists(initcondfile)):
@@ -227,13 +313,22 @@ class RTModel:
                         processes[i].kill()
                         timeouts += 1
                         crashes -= 1
-                        premodfiles = glob.glob(self.eventname +'/PreModels/*.txt')
-                        strmodel =  modelcode + '{:0>4}'.format(str(procnumbers[i]))
-                        with open(self.eventname +'/PreModels/' + strmodel + '/t' + strmodel + '.dat','w') as f:
-                            f.write(f'{len(premodfiles)} {self.LevMar_nfits}')
+                        #premodfiles = glob.glob(self.eventname +'/PreModels/*.txt')
+                        #strmodel =  modelcode + '{:0>4}'.format(str(procnumbers[i]))
+                        #with open(self.eventname +'/PreModels/' + strmodel + '/t' + strmodel + '.dat','w') as f:
+                        #    f.write(f'{len(premodfiles)} {self.LevMar_nfits}')
                     if(processes[i].poll() != None):
                         if(processes[i].returncode!=0):
                             crashes +=1
+                        # Here we have to append results to main model file
+                        if(not self.LevMar_stepchainsave):
+                            strmodel = modelcode + '{:0>4}'.format(str(procnumbers[i])) + ".txt"
+                            if(os.path.exists(self.eventname +'/PreModels/' + strmodel)):
+                                with open(self.eventname +'/PreModels/' + strmodel) as f:
+                                    content = f.read()
+                                with open(self.eventname +'/PreModels/'+ modelcode + ".txt","a") as f:
+                                    f.write(content)
+                                os.remove(self.eventname +'/PreModels/' + strmodel)
                         processes.pop(i)
                         procnumbers.pop(i)
                         procepochs.pop(i)
@@ -242,12 +337,12 @@ class RTModel:
                         i += 1
                 while(iinitcond < ninitconds and len(processes) < self.nprocessors):
                     strmodel =  modelcode + '{:0>4}'.format(str(iinitcond))
-                    if(glob.glob(self.eventname +'/PreModels/' + strmodel + '/t' + strmodel + '.dat')==[]):
-                        processes.append(subprocess.Popen([self.bindir+self.levmarexe,self.eventname, strmodel,self.satellitedir], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL))
-                        procnumbers.append(iinitcond)
-                        procepochs.append(time.time())
-                    else:
-                        finitcond += 1
+                    #if(glob.glob(self.eventname +'/PreModels/' + strmodel + '/t' + strmodel + '.dat')==[]):
+                    processes.append(subprocess.Popen([self.bindir+self.levmarexe,self.eventname, strmodel,self.satellitedir], cwd = self.bindir, shell = False, stdout=subprocess.DEVNULL))
+                    procnumbers.append(iinitcond)
+                    procepochs.append(time.time())
+                    #else:
+                    #    finitcond += 1
                     iinitcond += 1
                 if(finitcond != finitcondold):
                     #print('  Fits launched: {}; completed: {}/{}'.format(iinitcond, finitcond, ninitconds))
@@ -312,7 +407,9 @@ class RTModel:
             print('\033[30;41m! Program stopped here!\033[m')
             self.done = True
 
-    def run(self, event = None, cleanup = False):
+    def run(self, event = None, cleanup = True):
+        if(not cleanup):
+            self.LevMar_stepchainsave = True
         phase =0
         if(event!= None):
             self.eventname = os.path.realpath(event)
@@ -472,6 +569,16 @@ class RTModel:
                         self.LevMar_timelimit = float(chunks[2])
                     elif(chunks[0]=='bumperpower'):
                         self.LevMar_bumperpower = float(chunks[2])
+                    elif(chunks[0] == 'turn_off_secondary_source' and chunks[2] == 'True'):
+                        self.LevMar_turn_off_secondary_source = True
+                    elif(chunks[0] == 'turn_off_secondary_lens' and chunks[2] == 'True'):
+                        self.LevMar_turn_off_secondary_lens = True
+                    elif(chunks[0] == 'mass_luminosity_exponent'):
+                        self.LevMar_mass_luminosity_exponent = float(chunks[2])
+                    elif(chunks[0] == 'mass_radius_exponent'):
+                        self.LevMar_mass_radius_exponent = float(chunks[2])
+                    elif(chunks[0] == 'lens_mass_luminosity_exponent'):
+                        self.LevMar_lens_mass_luminosity_exponent = float(chunks[2])
         if(os.path.exists(pathname + '/' + self.inidir + '/ModelSelector.ini')):
             with open(pathname + '/' + self.inidir + '/ModelSelector.ini','r') as f:
                 lines = f.read().splitlines()

{rtmodel-2.4 → rtmodel-3.1}/RTModel/__init__.py

@@ -1,4 +1,4 @@
-__version__ = "2.4"
+__version__ = "3.0"
 __author__ = 'Valerio Bozza'
 __credits__ = 'University of Salerno, Italy'