Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.

Topics - Anne Hulsey

Pages: [1] 2
Hello. I would like to use the uq module to create a truncated lognormal fragility curve. I currently use stats.lognorm(beta, scale=median).pdf(x) to get a lognormal fragility curve so I'm hoping for something that could give me a similar result but with a truncation. Any suggestions?

Regional Hazard Simulation (R2D, rWhale) / openquake selection widget
« on: February 07, 2022, 10:46:56 PM »
Many thanks to Steve for implementing the openquake selection widget in R2D and to Adam for recognizing how useful it would be to me and coordinating the effort. I appreciate the proactive attentiveness to both user needs and the underlying details.

Earthquake Engineering (EE-UQ) / record max base shear for AutoSDA
« on: October 07, 2021, 09:03:44 PM »

I want to check whether each AutoSDA gm causes yield in the building. (I am using it to cross-check the Taghavi-Miranda modal approximation, which is limited to linear responses.) The best way I can think of is to compare the max base shear to the pushover curve. Is it possible to include max base shear as an output?



I would like to be able to upload a downloaded PEER gm folder to EE-UQ without having to manually enter the scale factors. There is currently an option to include a .txt file ( but this would need me to run a script to create the file, every time I download the folder.

Could you make the upload search for a _SearchResults.csv file and include the scale factors automatically? I think this is better than automatically including the .txt file when downloading via EE-UQ, since that would mean that only folders downloaded via EE-UQ could be uploaded in this way.

Earthquake Engineering (EE-UQ) / number of records for AutoSDA
« on: October 06, 2021, 04:39:16 PM »

I've gotten some surprising results when changing the number of records for the AutoSDA analysis.

1) I tried to run 1 gm so that I could get the design/model files stored locally before running more gms on DesignSafe. This caused an error.

2) 2 gms ran locally but sampled from the same two records to populate 11 rows in the results.

3) I tried to run 16 gms on DesignSafe to see if it only sampled 11 of them. The job says it finished but when I retrieve the results I get a blank table in the Data Values tab. This might be because I accidentally selected SRSS instead of H1 but I'm having a hard time testing it because my jobs are getting stuck in a queue. I'm attaching the log file. The zipped tmp folder is too big and looks like it's just input files anyway.

Earthquake Engineering (EE-UQ) / retrieving tmp files from DesignSafe
« on: October 01, 2021, 12:42:32 AM »

When I run AutoSDA locally, I save the entire tmp.SimCenter folder so that I can retrieve the geometry, loading, periods, etc. Is it possible to retrieve this whole folder when running it on DesignSafe?

Also, I got an error when retrieving my results. Do you recognize the error: Error Running Dakota: -- Expected 14 function value(s) but found 5?


Earthquake Engineering (EE-UQ) / pushover results for AutoSDA
« on: September 28, 2021, 12:12:13 AM »

Is there a way to get pushover results from the AutoSDA analysis? It looks like they have been commented out of the main_program script. Speaking of which, I'm a little surprised to see that the Eigen Value Analysis has also been commented out, given that I am still getting eigen value results, such as the first few periods of the building.

If you could help me confirm how to get all the results, in particular the pushoever curve, I would appreciate it.


Regional Hazard Simulation (R2D, rWhale) / newest openquake implementation
« on: September 21, 2021, 08:45:19 PM »
Hi Kuanshi,

Now that I have confirmed that openquake is producing the same results, I'm switching to using the standard gmms that the New Zealand NSHM is going to use (not the custom GNS ones I recently sent you). In particular, the Kuehn, Parker, and Abrahamson&Gulerce subduction models from 2020, all available here:

They were not automatically included in the openquake version that was packaged in the R2D build so I copied them over. However, openquake recently reorganized somethings and it looks like the R2D build cannot accommodate them. (For example, they require the base.add_alias function from the same folder I linked above.)

Would it be possible to include the openquake version that is on github, instead of what we are currently using? (I recall at one point I asked you to go back to a previous version - maybe that is what is included now.) Unfortunately, replacing the openquake folder wasn't sufficient because FetchOpenQuake calls, which is no longer included.


Regional Hazard Simulation (R2D, rWhale) / openquake poe values
« on: September 21, 2021, 12:24:58 AM »
Hi Kuanshi,

I just discovered that the poe terms in the oq .ini files is for the investigation time, not annual. Could you update line 238 in to be (something like):

1 - np.exp(-n_years * (1/RP))


Hello Kuanshi,

I've been crosschecking the new oq results against previous results. It turns out that the backarc term that automatically gets included in the OpenQuakeSiteFile.csv is messing up the Abrahamson GMM for subduction interface ruptures. In future builds, could you remove line 217 from


Earthquake Engineering (EE-UQ) / AutoSDA results for 10+ stories
« on: September 18, 2021, 12:22:48 AM »

It looks like the AutoSDA edps above the 10th story are getting collected incorrectly. The attached results have identical edp for the 11-15th and the 1-5th stories. See this (unrestricted) folder for the backend files:

Also, would it be possible to collect the peak drifts, in addition to the displacements and accelerations?


Thanks, Kuanshi, for your ongoing work on extracting the realizations.csv file from openquake. In case it is helpful, I am attaching the Jupyter Notebook that I currently use to extract it.

Regional Hazard Simulation (R2D, rWhale) / Running classical PSHA in R2D
« on: September 16, 2021, 10:54:06 PM »
Thanks to Kuanshi for fixing the Z1.0 and Z2.5 conversion from the SiteFile.csv to the OpenQuakeSiteModel.csv!

The new file looks good but OpenQuake crashes when it tries to run it. Here's a link to the working directory, accessible by Kuanshi and Adam:

I reran it all in PowerShell (see the uploaded .txt file in the linked folder), using:

python "C:/Users/ahul697/Desktop/R2D/build_date\\applications\\performRegionalEventSimulation\\regionalGroundMotion\\"   --hazard_config   "C:/Users/ahul697/Documents/R2D/LocalWorkDir\\HazardSimulation\\GroundMotions\\Input\\EQHazardConfiguration.json"

The 09_15 build gets to the same spot and crashes at :

FetchOpenQuake: OpenQuake configured.
ERROR 1: PROJ: proj_create_from_database: Cannot find proj.db
PROJ: proj_create_from_database: Cannot find proj.db
INFO:root:Upgrade completed in 0.000997304916381836 seconds
WARNING:root:DB server started with C:\Users\ahul697\Anaconda3\python.exe on tcp://, pid 23532
[2021-09-15 16:12:08 #4 WARNING] The parameter 'trucation_level' is unknown, ignoring
FetchOpenQuake: Classical PSHA failed.

I then reran the files in the LocalWorkDir using the 08/23 build (also in the .txt file), hoping that it would retrieve the correct OpenQuakeSiteModel.csv file. It was able to run and produce results, however, it recreated the .csv file, using the old scripts that didnít correctly include Z1.0 and Z2.5.

Note that both the 09/15 and 08/23 builds include the line:
[2021-09-15 16:12:08 #4 WARNING] The parameter 'trucation_level' is unknown, ignoring.
Iím curious about what that means. But, regardless, the unique issue for the 09/15 build seems to be the missing proj.db file.



I would like to use the R2D Hazard tool to simulate IM maps, rather than select ground motions. Additionally, I would like to be able to upload openquake files to identify the rupture scenario and include a gmpe logic tree.

I think the most important feature is being able to extract the predicted median and std dev at each location (perhaps a vector of them, representing the realizations of the logic tree). Is it possible to export this in the same format as EQHazard? Actually, even better than that, to export it as the .h5 files that are used in (See the extract_eqhazard_data function in the file.)

Beyond that, it could also be nice to run the whole simulation in R2D, then extract the IM maps.



The EVT page of EE-UQ includes a pulldown for the acceleration component of interest. As this is part of the record selection, it would be more intuitive if it were placed above the "Select Records" button. I'm attaching a suggestion.

Pages: [1] 2