#Marissa E Cameron Tips and tricks for how the command line Satstress process works. This uses GMT4 and old SatStress command line. If you are running these codes on GMT5, you will need to update to the appropriate GMT syntax. Alternatively, use gmtswitch to go back and forth between GMT4 and GMT5 as needed. __________________________________________ Coordinate info and file format, TERMINAL ------------------------------------------ 1) Create/format longitude and latitude files. Currently SatStress wants EAST POSITIVE LONGITUDE and CO-LATITUDE. INPUT FILE: FILENAME_elon_lat.txt OUTPUT FILE: FILENAME_elon_plat.txt ________________________________ Calculate beta angles, MATLAB -------------------------------- calculate_beta.m 2) This will append a column of beta angles to your elon_plat file. INPUT FILE: FILENAME_elon_plat.txt OUTPUT FILE: FILENAME_elon_plat_beta.txt ________________________________ SatStress + GMT, TERMINAL -------------------------------- grdtrack_stress_Ganymede.com run_grdtrack_stress_Ganymede_NSR.com 3) File can now be run through SatStress. NOTES: -o flag is for Principal stresses, -d flag is for TAU. INPUT FILE: FILENAME_elon_plat_beta.txt OUTPUT FILE WHEN RUN WITH run_grdtrack_stress.com : m0_midptlonlatbetaTTPTPP.dat : m360_midptlonlatbetaTTPTPP.dat ___________________________________ Shear and Normal stresses, TERMINAL ----------------------------------- compute_shear_normal.com run_compute_shear_normal.com 4) Conmpute shear and normal stresses. INPUT FILES WHEN RUN WITH run_grdtrack_stress.com for varying m: m0_midptlonlatbetaTTPTPP.dat : m360_midptlonlatbetaTTPTPP.dat OUTPUT FILES: m0_midptlonlatbetaTTPTPPTsTn.dat : m360_midptlonlatbetaTTPTPPTsTn.dat _________________________ Coulomb stresses, MATLAB ------------------------- Complete_get_coulomb_sense_of_shear_Ganymede.m 5) Run in MATLAB. There is A LOT to experiment with here. Code currently contains a few nestled loops to allow me to quickly loop through ALL orbital postions (0:360, in increments of 10 degrees) and also compute at depths of z=0:2000 in increments of 100 meters. User can change coefficient of friction as needed, H/d if necessary. INPUT FILES: m0_midptlonlatbetaTTPTPPTsTn.dat : m360_midptlonlatbetaTTPTPPTsTn.dat OUTPUT FILES: failure_m0_midptlonlatbetaTTPTPP_z0.out : failure_m360_midptlonlatbetaTTPTPP_z0.out WITH z also going from 0:2000. * Remember, this will result in a lot of files since you're not only looping over the orbital positions but also the depths. _________________________ Plot, GMT and TERMINAL ------------------------- plot_NSR_compfail.com run_plot_NSR_compfail.com 6) Plot in GMT4.