SPECFEM / specfem3d

SPECFEM3D_Cartesian simulates acoustic (fluid), elastic (solid), coupled acoustic/elastic, poroelastic or seismic wave propagation in any type of conforming mesh of hexahedra (structured or not).
GNU General Public License v3.0
403 stars 225 forks source link

How to set dynamic source in CMTSOLUTION file? #1469

Open xuhanqing98 opened 2 years ago

xuhanqing98 commented 2 years ago

Dear all I tried to use CUBIT to build a model and simulate the propagation of seismic wave, but after I changed the values of latorUTM and longorUTM in the CMTSOLUTION file, the sema file presented would not change. I think I set the dynamic source wrong., Here is my CMTSOLUTION FILE:

PDE 1999 01 01 00 00 00.00 0 0 -120 4.2 4.2 test1 event name: test1 time shift: 0.0000 half duration: 0.04 latorUTM: -50000.0 longorUTM: -50000.0 depth: -100000.0 Mrr: 5.690000e+27 Mtt: 1.040000e+26 Mpp: -5.790000e+27 Mrt: -1.990000e+27 Mrp: 5.610000e+27 Mtp: -3.480000e+27 PDE 1999 01 01 00 00 00.00 0 0 -120 4.2 4.2 test2 event name: test2 time shift: 0.0000 half duration: 0.37 latorUTM: -49000.0 longorUTM: -49000.0 depth: -110000.0 Mrr: -5.600000e+24 Mtt: 1.390000e+25 Mpp: -8.300000e+24 Mrt: -3.510000e+24 Mrp: 2.050000e+24 Mtp: 1.120000e+25 PDE 1999 01 01 00 00 00.00 0 0 -120 4.2 4.2 test3 event name: test3 time shift: 0.0000 half duration: 0.625 latorUTM: -49500.0 longorUTM: -49500.0 depth: -150000.0 Mrr: 2.400000e+24 Mtt: 3.780000e+23 Mpp: -2.780000e+24 Mrt: -6.860000e+23 Mrp: 1.430000e+24 Mtp: -1.410000e+24 PDE 1999 01 01 00 00 00.00 0 0 -120 4.2 4.2 test4 event name: test4 time shift: 0.0000 half duration: 1.00 latorUTM: -49300.0 longorUTM: -49300.0 depth: -120000.0 Mrr: -8.090000e+22 Mtt: 6.990001e+23 Mpp: -6.180000e+23 Mrt: 1.380000e+23 Mrp: 4.770000e+23 Mtp: 4.420000e+23 PDE 1999 01 01 00 00 00.00 0 0 -120 4.2 4.2 test5 event name: test5 time shift: 0.0000 half duration: 1.25 latorUTM: -49200.0 longorUTM: -49200.0 depth: -140000.0 Mrr: -4.000000e+23 Mtt: 4.490000e+23 Mpp: -4.970000e+22 Mrt: 1.430000e+23 Mrp: 2.510000e+22 Mtp: 2.280000e+23

jpampuero commented 2 years ago

To set a dynamic source (which means solving the spontaneous rupture process on a fault overned by a friction law) follow the guidelines in section 6 of the manual. The CMTSOLUTION file is not used in dynamic rupture simulations.

xuhanqing98 commented 2 years ago

@jpampuero So change FAULT_STATIONS and Par_file_faults?