Tsunami-wave-heights-from-fault-slip-distribution-with-linear-combinations.md

@@ -85,12 +85,21 @@ Tsunami initial conditions for this fault structure are created by running:
 ./run_ps_init_cond.sh < list_slip_hellenic_arc.txt
 ```
 
-This bash script\* create and launch (and then also remove) a job script on the cluster for each slip distribution scenario in the list. Tsunami initial conditions will be saved as `netcdf` files and stored in the `ps_initial_conditions/` folder, converting the folder path structure of the slip distribution into the name of the tsunami initial condition file (for example: `HeA/hellenic_Hazard_slip_HeA_025_variable_mu_6_5000_Strasser_Slip4cells_02635_000.nc`).
+This bash script create and launch (and then also remove) a job script on the cluster for each slip distribution scenario in the list. Tsunami initial conditions will be saved as `netcdf` files and stored in the `ps_initial_conditions/` folder, converting the folder path structure of the slip distribution into the name of the tsunami initial condition file (for example: `HeA/hellenic_Hazard_slip_HeA_025_variable_mu_6_5000_Strasser_Slip4cells_02635_000.nc`).
 
-_\* At the moment, the script `run_ps_init_cond` is set for the auriga cluster only, because of mercalli cluster has different library versions._
+Create ps_links folder and set the following variables in link_initial_conditions.py: 
 
-### 3\. Run linear combinations
+```
+dir_links = "/scratch/projects/cat/stress_drop_experiment/ps_links" 
+dir_ini_con = "/scratch/projects/cat/stress_drop_experiment/ps_initial_conditions"
+```
 
-Before running linear combinations, the following operations are needed to set up the procedure.
+Then run the script:
 
+```
+python link_initial_conditions.py CaA > init_cond_CaA.txt
+```
+
+### 3\. Run linear combinations
 
+Before running linear combinations, the following operations are needed to set up the procedure.
\ No newline at end of file