dogukankaratas / scalepy Goto Github PK

#%%
import pandas as pd , os , numpy as np
from aad_TBDY2018_ResponseSpectrumCreator_class import tbdy2018_spectra
import datetime
import json

#%%

fault_list = [ "Strike - Slip" ,  "Normal" , "Reverse" , 'Reverse - Oblique' , 'Normal - Oblique']
range_list = ["0,50","0,100","0,200","0,300",",50,100","50,200","50,300"]
mangitude_list = ["4,9"]
vs_list= ["0,180","180,360","360,760","760,1500","1500,10000"]
T1_list = [.25 , .5 , .75 , 1.0 , 1.25 , 1.5]
intensity_list = ["DD1","DD2","DD3","DD4"]

sonuc = {}

print( datetime.datetime.now())
for magnitude_range in mangitude_list[:] : 
    for vs_range in vs_list[:]:
        for rjb_range in range_list[:]:
            for fault_mechanism in fault_list[:]:
                for intensity in intensity_list[:]:
                    for T1 in T1_list[:]:
                        #pass
                        try: 
                            vs30 = 0.5*( float(vs_range.split(',')[0]) + float(vs_range.split(',')[1]) )
                            """
                            """
                            spectrum_object = tbdy2018_spectra()
                            spectrum_object.getSpectraValue( 39.39062046395253, 32.307622830661735 , intensity )
                            #spectrum_object.show_spectral_values()
                            spectrum_object.get_spectral_ordinates( vs30 )
                            spectrum_object.select_records_aad( magnitude_range, vs_range, rjb_range, fault_mechanism, T1)
                            spectrum_object.BothComponentScaling(T1 , plot_save = True )
                            
                            print( f"{datetime.datetime.now() }-{intensity}-{magnitude_range}- { vs_range}- { rjb_range}- { fault_mechanism}- { T1} ")
                            
                            if len( spectrum_object.ScaleFactors ) == 11 : 
                                sonuc[f"{intensity}-{magnitude_range}- { vs_range}- { rjb_range}- { fault_mechanism}- { T1}"] = spectrum_object.ScaleFactors
                            else : 

                                sonuc[f"{intensity}-{magnitude_range}- { vs_range}- { rjb_range}- { fault_mechanism}- { T1}"] = "FAIL"


                        except Exception as err:

                            print( err)
# %%
json.dump( sonuc, open( "6ICEES_Sonuc.json", 'w' ) )