From b71f385f2118dd85aa33a22ec42f8e38f5e55f1e Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 10 Mar 2023 16:24:40 -0500 Subject: [PATCH] Got rid of quasimc stuff from master where it isn't implemented yet --- python/ctem/ctem/driver.py | 36 ------------------------------------ 1 file changed, 36 deletions(-) diff --git a/python/ctem/ctem/driver.py b/python/ctem/ctem/driver.py index 24f16c9c..bbc9a616 100644 --- a/python/ctem/ctem/driver.py +++ b/python/ctem/ctem/driver.py @@ -46,7 +46,6 @@ def __init__(self, param_file="ctem.in", isnew=True): 'ctemfile': param_file, 'impfile': None, 'sfdcompare': None, - 'quasimc': None, 'sfdfile' : None, 'realcraterlist': None, } @@ -135,41 +134,6 @@ def __init__(self, param_file="ctem.in", isnew=True): # Scale the production function to the simulation domain self.scale_production() - # Setup Quasi-MC run - - if (self.user['quasimc'] == 'T'): - - #Read list of real craters - print("quasi-MC mode is ON") - print("Generating the crater scaling data in CTEM") - rclist = util.read_formatted_ascii(self.user['realcraterlist'], skip_lines = 0) - tempfile = os.path.join(currentdir, 'temp.in') - - # Generate craterlist.dat - shutil.copy2(self.user['ctemfile'], tempfile ) - - #Write a temporary input file to generate the necessary quasimc files - util.write_temp_input(self.user, tempfile) - util.write_datfile(self.user, self.output_filenames['dat'], self.seedarr) - self.compute_one_interval(ctemin=tempfile) - os.remove(tempfile) - - #Interpolate craterscale.dat to get impactor sizes from crater sizes given - df = pandas.read_csv(self.output_filenames['craterscale'], sep='\s+') - df['log(Dc)'] = np.log(df['Dcrat(m)']) - df['log(Di)'] = np.log(df['#Dimp(m)']) - xnew = df['log(Dc)'].values - ynew = df['log(Di)'].values - interp = interp1d(xnew, ynew, fill_value='extrapolate') - rclist[:,0] = np.exp(interp(np.log(rclist[:,0]))) - - #Convert age in Ga to "interval time" - rclist[:,5] = (self.user['interval'] * self.user['numintervals']) - craterproduction.Tscale(rclist[:,5], 'NPF_Moon') - rclist = rclist[rclist[:,5].argsort()] - - #Export to dat file - util.write_realcraters(self.output_filenames['craterlist'], rclist) - util.write_datfile(self.user, self.output_filenames['dat'], self.seedarr) else: print('Continuing a previous run')