{"id":16455748,"url":"https://github.com/natgeo-wong/sumatra-modeq","last_synced_at":"2025-06-20T03:32:20.764Z","repository":{"id":103609979,"uuid":"196195985","full_name":"natgeo-wong/Sumatra-ModEQ","owner":"natgeo-wong","description":"Code used in modelling moderate Sumatran earthquakes from SuGAr data.","archived":false,"fork":false,"pushed_at":"2019-08-30T22:33:51.000Z","size":6861,"stargazers_count":0,"open_issues_count":0,"forks_count":0,"subscribers_count":1,"default_branch":"master","last_synced_at":"2025-01-09T23:41:36.470Z","etag":null,"topics":["earthquakes","gps","modelling","sumatra"],"latest_commit_sha":null,"homepage":"","language":"MATLAB","has_issues":true,"has_wiki":null,"has_pages":null,"mirror_url":null,"source_name":null,"license":null,"status":null,"scm":"git","pull_requests_enabled":true,"icon_url":"https://github.com/natgeo-wong.png","metadata":{"files":{"readme":"README.md","changelog":null,"contributing":null,"funding":null,"license":null,"code_of_conduct":null,"threat_model":null,"audit":null,"citation":null,"codeowners":null,"security":null,"support":null,"governance":null,"roadmap":null,"authors":null,"dei":null,"publiccode":null,"codemeta":null}},"created_at":"2019-07-10T11:52:08.000Z","updated_at":"2019-08-30T22:33:52.000Z","dependencies_parsed_at":null,"dependency_job_id":"b80a6c08-d811-4759-8f09-15bf6d6883e8","html_url":"https://github.com/natgeo-wong/Sumatra-ModEQ","commit_stats":{"total_commits":15,"total_committers":1,"mean_commits":15.0,"dds":0.0,"last_synced_commit":"afeeaee57bc6d39a325dab4dbe7a8b6f206a77fd"},"previous_names":[],"tags_count":0,"template":false,"template_full_name":null,"repository_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/natgeo-wong%2FSumatra-ModEQ","tags_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/natgeo-wong%2FSumatra-ModEQ/tags","releases_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/natgeo-wong%2FSumatra-ModEQ/releases","manifests_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories/natgeo-wong%2FSumatra-ModEQ/manifests","owner_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners/natgeo-wong","download_url":"https://codeload.github.com/natgeo-wong/Sumatra-ModEQ/tar.gz/refs/heads/master","host":{"name":"GitHub","url":"https://github.com","kind":"github","repositories_count":240959115,"owners_count":19884911,"icon_url":"https://github.com/github.png","version":null,"created_at":"2022-05-30T11:31:42.601Z","updated_at":"2022-07-04T15:15:14.044Z","host_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub","repositories_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repositories","repository_names_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/repository_names","owners_url":"https://repos.ecosyste.ms/api/v1/hosts/GitHub/owners"}},"keywords":["earthquakes","gps","modelling","sumatra"],"created_at":"2024-10-11T10:23:32.103Z","updated_at":"2025-02-27T01:18:04.115Z","avatar_url":"https://github.com/natgeo-wong.png","language":"MATLAB","funding_links":[],"categories":[],"sub_categories":[],"readme":"# Sumatran-ModEQ\n\n`EQ` is an earthquake forward-modelling framework based on the `GTdef` (Georgia Tech \nDeformation) code.  It can be configured to model earthquakes along a given slab geometry \n(in our case Slab1.0), or it may otherwise be used to create a fault plane based on given \ninput of (longitude,latitude,depth,strike,dip), or lastly along a horizontal plane \ncentered around a point (longitude,latitude,depth).  This file is designed to help the \nuser set up the code for their own use, with some tips on how to customise model input \nparameters.\n\nTests indicate that despite the limits of forward modelling data, we are able to relocate \nevents with sparse GPS data closer to regional seismic estimates better than their initial \nteleseismic inputs.  This therefore allows those interested to quickly obtain improved \nconstrained estimates of event location without requiring detailed seismic data.\n\nAny questions related to the main `EQ` code may be directed to me.  Questions relating to \nthe `GTdef` code dependency should be directed to Dr. Lujia Feng, currently at NTU.\n\nAuthors:\n* Nathanael Zhixin Wong (Main `EQ` code): nathanaelwong@fas.harvard.edu\n* Dr. Lujia Feng (`GTdef` dependency): lfeng@ntu.edu.sg\n\n\n## Download and Requirements\n\nThe first step is always to clone/fork the repository\n```{bash}\ngit clone https://github.com/natgeo-wong/Sumatra-ModEQ.git\n```\n\nRunning `EQ` requires MATLAB 2016b and above, at least.  It also requires the following \ntoolboxes:\n* Parallel toolbox\n* Optimization toolbox\n\n\n## Setup of `EQ` input files via `EQ_create.m`.\n\nThe basis of our forward models is that we must already have a set of predefined inputs.  \nWe generally take this from gCMT teleseismic estimates.  Based on these estimates, we \ncreate a `.eq` input file by running the `EQ_create.m` script.  This script creates a \nfolder in the current working directory that contains the relevant `.eq` script.  Both the \nworking directory and the `.eq` input script are named based on the date, magnitude, \nregression used and the rigidity.\n\nIt is advised that the user create a separate directory that contains all these input \nfiles and directories, because `EQ` will run the modelling and analysis inside these \ndirectories.\n\nAfter this, it is necessary to copy the contents of `EQ_sub` into the working directory.  \n`EQ_sub` contains both the master script `EQ.m` and submission scripts (e.g. `*.pbs`) that \nare necessary to run these scripts on a server.  The submission scripts will vary \ndepending on the different servers and user style of inputs.  It is also necessary to \nmodify the `startup.m` file, to add the paths of the main `EQ` codes, slab files and GPS \nvelocity input files.\n\nFor example, in the server, my codes are found in `/home/nwong002/Codes/EQ/`.  Therefore, \nin my startup folder, I add the following paths:\n\n```{bash}\n/home/nwong002/Codes/EQ/         # Folder containing main EQ codes\n/home/nwong002/Codes/EQ/Slab/    # Folder containing slab/fault input information\n/home/nwong002/Codes/EQ/GPS_vel/ # Folder containing GPS velocity information\n/home/nwong002/Codes/EQ/GTdef/   # GTdef code directory, requires a separate download\n```\n\n\n## Running the `EQ` model and analysing the results\n\nIn the directory, you should have at least the following files:\n```{bash}\n*.eq\nstartup.m\n*submission script\n```\n\nOnce the script has been run, you should see at least the following additional files and \nfolders:\n```{bash}\n./GTdef_L_Input/   # Directory containing input information for L loop\n./GTdef_L_Output/  # Directory containing output information for L loop\n./GTdef_L_Misfit/  # Directory containing misfit information for L loop\n./GTdef_SV_Input/  # Directory containing input information for SV loop\n./GTdef_SV_Output/ # Directory containing output information for SV loop\n./GTdef_SV_Misfit/ # Directory containing misfit information for SV loop\n./GTdef_SW_Input/  # Directory containing input information for SW loop\n./GTdef_SW_Output/ # Directory containing output information for SW loop\n./GTdef_SW_Misfit/ # Directory containing misfit information for SW loop\n*bfit*.txt         # Textfile containing bestfit information\n```\n\nBased on this information, you can then run the `analysis.m` script.  Once the script has \nbeen run, you should see the following additional folders:\n```{bash}\n./Analysis_ve1/ # Contains information required to plot patchfiles in GMT4 (from SV loop)\n./Analysis_ve2/ # Contains information required to plot patchfiles in GMT4 (from SW loop)\n```\n\nGenerally, I choose to use the SW (GPS-weighted variance explained) information in my plotting and analysis, as it tends to put more weightage on the horizontal displacement observed, and measurements of  horizontal displacements are more accurate than those of vertical displacements for GPS.\n\n\n## Downloading the GTdef code.\n\nThe GTdef code was created by my mentor Dr. Lujia Feng and her advisor Prof. Andrew Newman at Georgia Technological University.  Though the `EQ` code requires the `GTdef` code as a dependency, the `GTdef` code is not opensource and therefore I do not have the rights to distribute it in public.  However, those who are interested in the GTdef code may email Dr. Feng at lfeng@ntu.edu.sg to attain a copy of the code for your use.  If `GTdef` is used in your work, please acknowledge Dr. Feng and her advisor Prof. Andrew Newman in your work.  Any questions relating to the `GTdef` code may be directed towards her.","project_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fnatgeo-wong%2Fsumatra-modeq","html_url":"https://awesome.ecosyste.ms/projects/github.com%2Fnatgeo-wong%2Fsumatra-modeq","lists_url":"https://awesome.ecosyste.ms/api/v1/projects/github.com%2Fnatgeo-wong%2Fsumatra-modeq/lists"}