Page History

...

Go to the Src directory that instruction has provided
Code Block
language bash
$ cd /home/jonney/20100904_Darfield_m7p1_201705221509/Src

Create the config file that is needed for generating a Rupture Model

Code Block

language	bash

$config_src.py -h
usage: config_src.py [-h] {1,2,3,4} ...
positional arguments:
  {1,2,3,4}   Choose the type of srf you want to generate, usage of each type
              use: gen_src_params.py type -h
    1         Generate type 1 srf (Point source srf).
    2         Generate type 2 srf (Point source to Finite Fault srf).
    3         Generate type 3 srf (Finite Fault descriptor to Finite Fault
              srf).
    4         Generate type 4 srf (Multiple Segment Finite Fault srf).
optional arguments:
  -h, --help  show this help message and exit

The extra options may differ from type to type. Use: "config_src.py type -h" will show you all the options can be provided for that type of source.

example of showing options for type 2:

Code Block

language	bash

$ config_src.py 2 -h
usage: config_src.py 2 [-h] [--LAT [LAT]] [--LON [LON]] [--MAG [MAG]]
                       [--STK [STK]] [--RAK [RAK]] [--DIP [DIP]] [--DT [DT]]
                       [--DEPTH [DEPTH]] [--SEED [SEED]] [--RVFRAC [RVFRAC]]
                       [--ROUGH [ROUGH]] [--SLIP_COV [SLIP_COV]]
                       [--MWSR [MWSR]] [--STOCH [STOCH]] [--GENSLIP [GENSLIP]]
                       [--PREFIX [PREFIX]]
optional arguments:
  -h, --help            show this help message and exit
  --LAT [LAT]           latitude (float)
  --LON [LON]           longitude (float)
  --MAG [MAG]           magnitude (float)
  --STK [STK]           strike (int)
  --RAK [RAK]           rake (int)
  --DIP [DIP]           dip (int)
  --DT [DT]             rupture timestep
  --DEPTH [DEPTH]
  --SEED [SEED]
  --RVFRAC [RVFRAC]
  --ROUGH [ROUGH]
  --SLIP_COV [SLIP_COV]
  --MWSR [MWSR]         Mw Scaling Relation (string), one
                        of:HanksBakun2002,BerrymanEtAl2002,VillamorEtAl2001
  --STOCH [STOCH]       String to set stoch gener ation. Default vaule will be
                        used if not specified
  --GENSLIP [GENSLIP]   Version of genslip
  --PREFIX [PREFIX]     The prefix you want for SRF files. defualt= source.
                        Ending with _ will auto append extra info to file
                        name. Defualt value will be used if not specified.

And example of creating type 2 Source file:

Code Block

language	bash

$ config_src.py 2 --LAT -43.6 --LON 172.3 --MAG 7.1
TYPE =  2
LAT =  -43.6
LON =  172.3
MAG =  7.1
PREFIX =  Srf/source
SEED =  103245
STOCH =  Stoch
GENSLIP =  3.3
####################
the params below are not specified, please manually edit setSrfParams.py to specify the value.
####################
MAG =
STK =
RAK =
DIP =
DT =
RVFRAC =
ROUGH =
SLIP_COV =
DEPTH =
MWSR =
GENSLIP =
####################
Afte manually fill in all the required parameters, please cd to where setSrfParams.py is located at and run make_src.sh
cd====================================================================================================
Cannot find folder /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245
####################

Go to the folder provided by instrutions
e.g. $ cd

. 
Attemp to create...
Created Folder: /home/jonney/20100904_Darfield_m7p1_

201705181350

201705221509/Src/Model/20100904_Darfield_m7p1

/t2_103245_20170518_175214

Edit the setSrfParams.py
note: every variables with #!!!PLEASE SPECIFY!!!# must be provided
Create the Model using make_src.sh
e.g. $ make_src.sh
Follow the instruction given on the screen to provide a path to a admin.

4. How to create observation files

note: you can skip this part if you are running for a future event (not yet happened)

Go to the Obs directory that instruction has provided
e.g. $ cd /home/jonney/20100904_Darfield_m7p1_201705181350/GM/Obs
Create the config file that is needed
$ config_obs.py geoNet_ftp_url
e.g. $ config_obs.py ftp://ftp.geonet.org.nz/strong/processed/Proc/2010/09_Sep/2010-09-04_045556/Vol1/data/
Run make_obs.sh to generate all the files
e.g. $ make_obs.sh
Follow the instruction given on the screen to provide a path to a admin.

5. How to setup the simulation

After admin tells you that all files are uploaded to Fitzroy, Run setup_remote.sh
$ setup_remote.sh username
e.g $ setup_remote.sh jonney

6. Run the simulation on Low Frequency (LF)

Go to the Simulation Folder, located at /nesi/projects/nesi00213/RunFolder/user/run_name
e.g. $ cd /nesi/projects/nesi00213/RunFolder/ykh22/20100904_Darfield_m7p1_201705181350
Run submit_emod3d.sh
$ ./submit_emod3d.h
It will ask you to enter a wall_clock_limit and if it should auto submit the job.
note: if you prefer to change some detail configurations, press 2(no) for auto submit.*
use llsubmit to submit your job if you selected to not to auto submit.
Use llq to check if your job has finished. (job status: R= running, C= Complete, I = idle)
$ llq -u user
After Emod3d finished, run submit_post_emod3d.sh
note: if you prefer to change some detail configurations, press 2(no) for auto submit.*
$ ./submit_post_emod3d.sh

7. Run the simulation on High Frequency (HF)

note: this part can run before LF is finished

Run install_bb.sh, to setup the parameters
$ ./install_bb.sh
Choose from the options according to simulation requirement
Run submit_hf.sh
$ ./submit_hf.sh
note: if you prefer to change some detail configurations, press 2(no) for auto submit.*

*: use llsubmit to submit your job if you selected not to auto submit. $ llsubmit job.ll

8. Run the Simulation on Broadband (BB)

Note: This part MUST wait untill both LF and HF are finished

...


no existing folder with type 2, 
Using sub-folder name: t2_103245
####################
Please go to:
cd /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245
to manually edit all the parameters that is empty in setSrfParams.py and run make_src.sh
####################

If config_src.sh has been run more than once, it will ask you if you want to update/replace a current param or create a new one.

You can choose to create a new sub-folder with time stamp or replace the old files.

example of creating new folder:

Code Block

language	bash

$ config_src.py 2 --LAT -43.6 --LON 172.3 --MAG 7.1 --STK 150.0 --RAK 6.0
TYPE =  2
LAT =  -43.6
LON =  172.3
MAG =  7.1
STK =  150.0
RAK =  6.0
PREFIX =  Srf/source
SEED =  103245
STOCH =  Stoch
GENSLIP =  3.3
####################
the params below are not specified, please manually edit setSrfParams.py to specify the value.
####################
DIP =
DT =
RVFRAC =
ROUGH =
SLIP_COV =
DEPTH =
MWSR =
====================================================================================================
====================================================================================================
select the src folder you want to place the setSrfParams.py into
====================================================================================================
 1. Create new setSrfParams.py
 2. t2_103245
Enter the number you wish to select (1-2):1	
====================================================================================================
Do you want to create a sub-folder: t2_103245_20170524_171346?
====================================================================================================
 1. Yes
 2. No
Enter the number you wish to select (1-2):1
####################
Please go to:
cd /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/220100904_Darfield_m7p1/t2_103245_20170524_171346
to manually edit all the parameters that is empty in setSrfParams.py and run make_src.sh
####################

If you choose to replace existing params, it will show a comparison between two params(old on the left, new on the right).
Example of replacing a old param file:

Code Block

language	bash

$ config_src.py 2 --LAT -43.6 --LON 172.3 --MAG 7.1 --STK 1.0 --RAK 1.0 --DT 0.025
TYPE =  2
LAT =  -43.6
LON =  172.3
MAG =  7.1
STK =  1.0
RAK =  1.0
DT =  0.025
PREFIX =  Srf/source
SEED =  103245
STOCH =  Stoch
GENSLIP =  3.3
####################
the params below are not specified, please manually edit setSrfParams.py to specify the value.
####################
DIP =
RVFRAC =
ROUGH =
SLIP_COV =
DEPTH =
MWSR =
====================================================================================================
====================================================================================================
select the src folder you want to place the setSrfParams.py into
====================================================================================================
 1. Create new setSrfParams.py
 2. t2_103245
 3. t2_103245_20170524_171346
Enter the number you wish to select (1-3):2
diff -y /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245/setSrfParams.py /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p/params.tmp
## Sets Variables for SRF/Stoch Generation			## Sets Variables for SRF/Stoch Generation
####################						####################
#!!!Important!!							#!!!Important!!
####################						####################
#Please change the values accordingly if '#!!!PLEASE SPECIFY!	#Please change the values accordingly if '#!!!PLEASE SPECIFY!

# TYPE:								# TYPE:
# 1: point source to point source srf				# 1: point source to point source srf
# 2: point source to finite fault srf				# 2: point source to finite fault srf
# 3: finite fault to finite fault srf				# 3: finite fault to finite fault srf
# 4: multi-segment finite fault srf				# 4: multi-segment finite fault srf
TYPE = 2							TYPE = 2
# specify basename for gsf, srf and stoch file created		# specify basename for gsf, srf and stoch file created
# PREFIX for gsf/srf/stoch files				# PREFIX for gsf/srf/stoch files
# if prefix ends with '_', automatic naming follows		# if prefix ends with '_', automatic naming follows
PREFIX = 'Srf/source'						PREFIX = 'Srf/source'
# directory for Stoch file(s)					# directory for Stoch file(s)
# set to None to not produce the stoch file			# set to None to not produce the stoch file
STOCH = 'Stoch'							STOCH = 'Stoch'
# rupture timestep (float), e.g. DT = 0.025			# rupture timestep (float), e.g. DT = 0.025
DT = #!!!PLEASE SPECIFY!!!#				      |	DT = 0.025
# latitude (float), e.g. LAT = -43.5095				# latitude (float), e.g. LAT = -43.5095
LAT = -43.6							LAT = -43.6
# longitude (float), e.g. LON = 172.1811			# longitude (float), e.g. LON = 172.1811
LON = 172.3							LON = 172.3
...
...
SEED = 103245							SEED = 103245
# Mw Scaling Relation (string), one of:				# Mw Scaling Relation (string), one of:
# HanksBakun2002, BerrymanEtAl2002, VillamorEtAl2001		# HanksBakun2002, BerrymanEtAl2002, VillamorEtAl2001
MWSR = ''#!!!PLEASE SPECIFY!!!#					MWSR = ''#!!!PLEASE SPECIFY!!!#

====================================================================================================
do you want to replace this param file?
====================================================================================================
 1. Yes
 2. No
Enter the number you wish to select (1-2):1
####################
Please go to:
cd /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245
to manually edit all the parameters that is empty in setSrfParams.py and run make_src.sh
####################

Go to the folder provided by instrutions

Code Block

language	bash

$ cd /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245

Edit the setSrfParams.py
note: every variables with #!!!PLEASE SPECIFY!!!# must be provided

Code Block

language	bash

$ vim setSrfParams.py

1 ## Sets Variables for SRF/Stoch Generation
 2 ####################
 3 #!!!Important!!
 4 ####################
 5 #Please change the values accordingly if '#!!!PLEASE SPECIFY!!!#'
 6 
 7 # TYPE:
 8 # 1: point source to point source srf
 9 # 2: point source to finite fault srf
10 # 3: finite fault to finite fault srf
11 # 4: multi-segment finite fault srf
12 TYPE = 2
13 # specify basename for gsf, srf and stoch file created
14 # PREFIX for gsf/srf/stoch files
15 # if prefix ends with '_', automatic naming follows
16 PREFIX = 'Srf/source'
17 # directory for Stoch file(s)
18 # set to None to not produce the stoch file
19 STOCH = 'Stoch'
20 # rupture timestep (float), e.g. DT = 0.025
21 DT = 0.025
22 # latitude (float), e.g. LAT = -43.5095
23 LAT = -43.6
24 # longitude (float), e.g. LON = 172.1811
25 LON = 172.3
26 # depth (float), ex. DEPTH = 2.0
27 DEPTH = 2.0#!!!PLEASE SPECIFY!!!#
28 # magnitude (float), e.g. MAG = 5.5
...
...
40 ## roughness of fault only for genslip 5.0+
41 ## 0.1 is a good value to use - Rob Graves
42 ## 0.0050119 = (10^(-2.3)) - Shi & Day 2014. Used as default in 5.2.3a
43 #ROUGH = #!!!PLEASE SPECIFY!!!#
44 ## (float), e.g. RVFRAC = 0.8
45 #RVFRAC = #!!!PLEASE SPECIFY!!!#
46 ## (float), e.g. SLIP_COV = 0.85
47 #SLIP_COV = #!!!PLEASE SPECIFY!!!#
48 %s# seed for stoch, e.g. SEED = 103245
49 SEED = 103245
50 # Mw Scaling Relation (string), one of:
51 # HanksBakun2002, BerrymanEtAl2002, VillamorEtAl2001
52 MWSR = ''#!!!PLEASE SPECIFY!!!#

Create the Model using make_src.sh

Code Block

language	bash

$ make_src.sh

Ececuting createSRF.py
Executing /nesi/projects/nesi00213/tools/genslip-v3.3read_erf=0write_srf=1read_gsf=1write_gsf=0infile=Srf/source.gsfmag=5.500000nx=58ny=58ns=1nh=1seed=103245velfile=/home/jonney/code/SrfGen//lp_generic1d-gp01_v1.vmodshypo=0.000000dhypo=2.877200dt=0.025000plane_header=1srf_version=2.0rvfrac=Noneshal_vrup=0.6deep_vrup=0.6fdrup_time=1alpha_rough=Noneslip_sigma=None
seed=103245
mag= 5.50 median mag= 5.52 nslip= 1 nhypo= 1 random_seed =103245 
nx= 58 ny= 58 dx=     0.0992 dy=     0.0992
  0:   172.30055   -43.59977 58 58     5.7543     5.7543 150.0 54.0     0.0000
****    16.5000	   16.5000
ratio (negative slip)/(positive slip)= 0.002216
mom=   1.99526e+24 avgslip= 48 maxslip= 98
orig_sigma= 0.403256 ... new_sigma= 0.396138
rt_scalefac= 7.388430
ravg=   1.00000e+00 rmed=   9.86539e-01 rmin=   3.09826e-01 rmax=   1.93412e+00
seg= 0
nstk= 58 nx= 3 nxdiv= 3 nxsum= 58
ndip= 58 ny= 3 nydiv= 3 nysum= 58
Plotting SRF as square plot...
...
/home/jonney/20100904_Darfield_m7p1_201705221509/gmsim.cfg
making a symbolic link from source_map.png to /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Figures/source_map.png
making a symbolic link from source_square.png to /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Figures/source_square.png
attempt_upload.sh /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245 /nesi/projects/nesi00213/RupModel/20100904_Darfield_m7p1 baes seb56
!!!!!  Ask Admin to upload: Give him the path: /home/jonney/20100904_Darfield_m7p1_201705221509/Src !!!!!

Under the working directory Src, a subdirectory is created

Code Block

language	bash

$ ls /home/jonney/20100904_Darfield_m7p1_201705221509/Src/Model/20100904_Darfield_m7p1/t2_103245
cnrs.txt  Info  params_plot.py  params_plot.pyc  setSrfParams.py  setSrfParams.pyc  Srf  Stoch

Give the admin (Sung) the path to the Src working directory and request to upload the model. The admin will need to execute sendme.sh in the Src directory.

4. How to create observation files

note: you can skip this part if you are running for a future event (not yet happened)

Go to the Obs directory that instruction has provided
Code Block
language bash
cd /home/jonney/20100904_Darfield_m7p1_201705221509/GM/Obs
Create the config file that is needed
$ config_obs.py geoNet_ftp_url
e.g. $ config_obs.py ftp://ftp.geonet.org.nz/strong/processed/Proc/2010/09_Sep/2010-09-04_045556/Vol1/data/
Run make_obs.sh to generate all the files
e.g. $ make_obs.sh
Follow the instruction given on the screen to provide a path to a admin.

5. How to setup the simulation

After admin tells you that all files are uploaded to Fitzroy, Run setup_remote.sh
$ setup_remote.sh username
e.g $ setup_remote.sh jonney

6. Run the simulation on Low Frequency (LF)

Go to the Simulation Folder, located at /nesi/projects/nesi00213/RunFolder/user/run_name
e.g. $ cd /nesi/projects/nesi00213/RunFolder/ykh22/20100904_Darfield_m7p1_201705181350
Run submit_emod3d.sh
$ ./submit_emod3d.h
It will ask you to enter a wall_clock_limit and if it should auto submit the job.
note: if you prefer to change some detail configurations, press 2(no) for auto submit.*
use llsubmit to submit your job if you selected to not to auto submit.
Use llq to check if your job has finished. (job status: R= running, C= Complete, I = idle)
$ llq -u user
After Emod3d finished, run submit_post_emod3d.sh
note: if you prefer to change some detail configurations, press 2(no) for auto submit.*
$ ./submit_post_emod3d.sh

7. Run the simulation on High Frequency (HF)

note: this part can run before LF is finished

Run install_bb.sh, to setup the parameters
$ ./install_bb.sh
Choose from the options according to simulation requirement
Run submit_hf.sh
$ ./submit_hf.sh
note: if you prefer to change some detail configurations, press 2(no) for auto submit.*

*: use llsubmit to submit your job if you selected not to auto submit. $ llsubmit job.ll

8. Run the Simulation on Broadband (BB)

Note: This part MUST wait untill both LF and HF are finished

Run submit_bb.sh
$ ./submit_bb.sh

Info

Child pages

Versions Compared

Old Version 26

New Version 27

Key

4. How to create observation files

5. How to setup the simulation

6. Run the simulation on Low Frequency (LF)

7. Run the simulation on High Frequency (HF)

8. Run the Simulation on Broadband (BB)

Related articles

4. How to create observation files

5. How to setup the simulation

6. Run the simulation on Low Frequency (LF)

7. Run the simulation on High Frequency (HF)

8. Run the Simulation on Broadband (BB)

Related articles