Skip to content

Commit

Permalink
Reverted rrup perturbation
Browse files Browse the repository at this point in the history
  • Loading branch information
sungeunbae committed Jul 19, 2024
1 parent f161fd7 commit df8275c
Showing 1 changed file with 0 additions and 30 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -148,33 +148,6 @@ def generate_source_params(
else:
print("Didn't get vs30")

#if realisation["vs30"]["vs30"] <= 0:
# realisation["vs30"]["vs30"] = 0
if 'sites_info' in realisation['params']:
sites_info_csv = Path(realisation['params']['sites_info'])
sites_info = pd.read_csv(realisation['params']['sites_info'], index_col=0)

rrup_data = sites_info.rrup
rrup_data_perturbed = distributions.truncated_log_normal(rrup_data , 0.1, 4)
sites_info["rrup"] = rrup_data_perturbed
if "vs30" in realisation["vs30"]:
sites_info["vs30"] = realisation["vs30"]["vs30"] # update vs30 values if vs30 values were perturbed.
new_sites_info_csv = sites_info_csv.parent / f"{sites_info_csv.stem}_perturbed.csv"
sites_info.to_csv(new_sites_info_csv, index=False)
else:
print("Didn't get sites_info CSV")

if sites_info is not None:
rrup_data = sites_info.rrup
rrup_data_perturbed = distributions.truncated_log_normal(rrup_data , 0.1, 4)
sites_info["rrup"] = rrup_data_perturbed
new_sites_info_csv = sites_info_csv.parent / f"{sites_info_csv.stem}_rrup_perturbed.csv"
sites_info.to_csv(new_sites_info_csv, index=False)
else:
print("Didn't get sites_info CSV")



### End of custom code area

# verify_realisation_params(realisation["params"]) #modified
Expand Down Expand Up @@ -210,9 +183,6 @@ def generate_from_gcmt(
#rdist_sigma = 2
#depth_sigma = additional_source_parameters["depth_sigma"]

sites_info_csv = Path(additional_source_parameters["sites_info"])
assert sites_info_csv.exists()

### the following parameters feed into srfgen

mag = distributions.truncated_normal(sources_line.mag, 0.2, 4) # magnitude
Expand Down

0 comments on commit df8275c

Please sign in to comment.