From ff796a786de2ebe489e15c485d54fcdf51188cae Mon Sep 17 00:00:00 2001 From: buildbot Date: Sat, 1 May 2021 14:12:39 -0400 Subject: [PATCH] Buildbot ran the code beautifier. --- .../REF_SEIS/output_generate_databases.txt | 54 +- .../REF_SEIS/output_meshfem3D.txt | 80 +- .../sep_bathymetry/REF_SEIS/output_solver.txt | 696 +++++++++--------- .../REF_SEIS/output_generate_databases.txt | 58 +- .../REF_SEIS/output_meshfem3D.txt | 68 +- .../REF_SEIS/output_solver.txt | 270 +++---- src/shared/utm_geo.f90 | 4 +- .../write_output_ASCII_or_binary.f90 | 2 +- src/specfem3D/write_seismograms.f90 | 8 +- utils/convert_utm2geo.py | 2 +- 10 files changed, 621 insertions(+), 621 deletions(-) diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_generate_databases.txt b/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_generate_databases.txt index 2221060af..4253917d3 100644 --- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_generate_databases.txt +++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_generate_databases.txt @@ -46,7 +46,7 @@ defined materials : 2 undefined materials : 0 total number of spectral elements: 9600 - absorbing boundaries: + absorbing boundaries: xmin,xmax : 200 200 ymin,ymax : 2400 2400 bottom,top: 192 192 @@ -60,14 +60,14 @@ create regions: - ...allocating arrays + ...allocating arrays NGLLX = 5 NGLLY = 5 NGLLZ = 5 NGNOD = 8 NGNOD2D = 4 - main process setup: + main process setup: nspec = 3200 separating regular/irregular element shapes @@ -84,16 +84,16 @@ File DATA/Par_file_faults not found: assuming that there are no faults - ...setting up jacobian + ...setting up jacobian ...indexing global points - creating ibool indexing : x min/max = 0.00000000 / 10000.0000 + creating ibool indexing : x min/max = 0.00000000 / 10000.0000 creating indirect addressing: nglob = 222105 creating unique point locations - ...preparing MPI interfaces + ...preparing MPI interfaces number of interfaces : 1 - creating MPI indexing : x min/max = 0.00000000 / 10000.0000 + creating MPI indexing : x min/max = 0.00000000 / 10000.0000 tolerance = 9.9999999999999995E-007 total MPI interface points: 13668 @@ -125,7 +125,7 @@ 90 % time remaining: 5.2642145371919441E-008 s 100 % time remaining: 0.0000000000000000 s - ...detecting acoustic-elastic-poroelastic surfaces + ...detecting acoustic-elastic-poroelastic surfaces total acoustic elements : 384 total elastic elements : 9216 total poroelastic elements: 0 @@ -135,18 +135,18 @@ elastic - poroelastic coupling : total number of faces = 0 - ...element inner/outer separation + ...element inner/outer separation for overlapping of communications with calculations: percentage of edge elements 6.25000000 % percentage of volume elements 93.7500000 % - ...element mesh coloring + ...element mesh coloring use coloring = F - ...external binary models + ...external binary models using SEP model from directory: ./sep_model - ...creating mass matrix + ...creating mass matrix ...saving databases using binary file format @@ -179,29 +179,29 @@ ******** - Model: P velocity min,max = 1499.99988 4500.00000 - Model: S velocity min,max = 0.00000000 2590.00000 + Model: P velocity min,max = 1499.99988 4500.00000 + Model: S velocity min,max = 0.00000000 2590.00000 - Model: Poisson's ratio min,max = 0.250554323 0.500000000 + Model: Poisson's ratio min,max = 0.250554323 0.500000000 ******** ********************************************* *** Verification of simulation parameters *** ********************************************* - *** Xmin and Xmax of the model = 0.00000000 10000.0000 - *** Ymin and Ymax of the model = 1600.00000 2400.00000 - *** Zmin and Zmax of the model = -5000.00000 0.00000000 + *** Xmin and Xmax of the model = 0.00000000 10000.0000 + *** Ymin and Ymax of the model = 1600.00000 2400.00000 + *** Zmin and Zmax of the model = -5000.00000 0.00000000 - *** Max GLL point distance = 81.8317108 - *** Min GLL point distance = 16.1879883 - *** Max/min ratio = 5.05508852 + *** Max GLL point distance = 81.8317108 + *** Min GLL point distance = 16.1879883 + *** Max/min ratio = 5.05508852 - *** Max element size = 250.000000 - *** Min element size = 93.7500000 - *** Max/min ratio = 2.66666675 + *** Max element size = 250.000000 + *** Min element size = 93.7500000 + *** Max/min ratio = 2.66666675 - *** Minimum period resolved = 0.226449981 + *** Minimum period resolved = 0.226449981 *** Maximum suggested time step = 1.79866538E-03 Elapsed time for checking mesh resolution in seconds = 3.9490000000004244E-003 @@ -210,7 +210,7 @@ mesh regions done - min and max of elevation (i.e. height of the upper surface of the mesh) included in mesh in m is 0.0000000000000000 0.0000000000000000 + min and max of elevation (i.e. height of the upper surface of the mesh) included in mesh in m is 0.0000000000000000 0.0000000000000000 done mesh setup @@ -235,7 +235,7 @@ smallest and largest possible floating-point numbers are: 1.17549435E-38 3.40282347E+38 - Elapsed time for mesh generation and buffer creation in seconds = 4.5278070000000001 + Elapsed time for mesh generation and buffer creation in seconds = 4.5278070000000001 End of mesh generation done diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_meshfem3D.txt b/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_meshfem3D.txt index 533852295..88a922b77 100644 --- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_meshfem3D.txt +++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_meshfem3D.txt @@ -41,7 +41,7 @@ Creating global slice addressing Spatial distribution of slice numbers: - 0 1 2 + 0 1 2 This is process 0 There are 3 MPI processes Processes are numbered from 0 to 2 @@ -63,19 +63,19 @@ region selected: - latitude min = 1600.0000000000000 - latitude max = 2400.0000000000000 + latitude min = 1600.0000000000000 + latitude max = 2400.0000000000000 - longitude min = 0.0000000000000000 - longitude max = 10000.000000000000 + longitude min = 0.0000000000000000 + longitude max = 10000.000000000000 this is given directly as UTM - UTM X min = 0.0000000000000000 - UTM X max = 10000.000000000000 + UTM X min = 0.0000000000000000 + UTM X max = 10000.000000000000 - UTM Y min = 1600.0000000000000 - UTM Y max = 2400.0000000000000 + UTM Y min = 1600.0000000000000 + UTM Y max = 2400.0000000000000 UTM size of model along X is 10.000000000000000 km UTM size of model along Y is 0.80000000000000004 km @@ -95,46 +95,46 @@ number of interfaces: 2 mesh: - origin UTM minimum x/y (m) = 0.00000000 1600.00000 - origin UTM maximum x/y (m) = 10000.0000 2400.00000 + origin UTM minimum x/y (m) = 0.00000000 1600.00000 + origin UTM maximum x/y (m) = 10000.0000 2400.00000 reading interface 1 interface file : bathy.dat number of points x/y = 501 21 - origin x/y (m) = 0.00000000 1600.00000 - spacing x/y (m) = 20.0000000 40.0000000 + origin x/y (m) = 0.00000000 1600.00000 + spacing x/y (m) = 20.0000000 40.0000000 - dimension x-direction (m) = 0.00000000 / 10000.0000 - dimension y-direction (m) = 1600.00000 / 2400.00000 + dimension x-direction (m) = 0.00000000 / 10000.0000 + dimension y-direction (m) = 1600.00000 / 2400.00000 total number of file points = 10521 should be 10521 this point total is okay - original elevation min/max = -500.000000 -500.000000 - interpolated mesh elevation min/max = -500.000000 -500.000000 + original elevation min/max = -500.000000 -500.000000 + interpolated mesh elevation min/max = -500.000000 -500.000000 - interpolated mesh UTM minimum x/y (m) = 0.00000000 1600.00000 - interpolated mesh UTM maximum x/y (m) = 10000.0000 2400.00000 + interpolated mesh UTM minimum x/y (m) = 0.00000000 1600.00000 + interpolated mesh UTM maximum x/y (m) = 10000.0000 2400.00000 reading interface 2 interface file : topo.dat number of points x/y = 2 2 - origin x/y (m) = 0.00000000 1600.00000 - spacing x/y (m) = 10000.0000 800.000000 + origin x/y (m) = 0.00000000 1600.00000 + spacing x/y (m) = 10000.0000 800.000000 - dimension x-direction (m) = 0.00000000 / 10000.0000 - dimension y-direction (m) = 1600.00000 / 2400.00000 + dimension x-direction (m) = 0.00000000 / 10000.0000 + dimension y-direction (m) = 1600.00000 / 2400.00000 total number of file points = 4 should be 4 this point total is okay - original elevation min/max = 0.00000000 0.00000000 - interpolated mesh elevation min/max = 0.00000000 0.00000000 + original elevation min/max = 0.00000000 0.00000000 + interpolated mesh elevation min/max = 0.00000000 0.00000000 - interpolated mesh UTM minimum x/y (m) = 0.00000000 1600.00000 - interpolated mesh UTM maximum x/y (m) = 10000.0000 2400.00000 + interpolated mesh UTM minimum x/y (m) = 0.00000000 1600.00000 + interpolated mesh UTM maximum x/y (m) = 10000.0000 2400.00000 ************************** @@ -159,16 +159,16 @@ creating mesh region 1 (regular mesh) mesh dimensions: - Xmin and Xmax of the model = 0.00000000 10000.0000 - Ymin and Ymax of the model = 1600.00000 2400.00000 - Zmin and Zmax of the model = -5000.00000 0.00000000 + Xmin and Xmax of the model = 0.00000000 10000.0000 + Ymin and Ymax of the model = 1600.00000 2400.00000 + Zmin and Zmax of the model = -5000.00000 0.00000000 exact area = 8000000.00 (m^2) = 8.00000000 (km^2) Max element size = 250.000000 (m) Min element size = 93.7500000 (m) - Max/min ratio = 2.66666675 + Max/min ratio = 2.66666675 creating indirect addressing for unstructured mesh @@ -186,7 +186,7 @@ ************************** start computing the minimum and maximum edge size - done processing + done processing ------------ mesh quality parameter definitions: @@ -197,9 +197,9 @@ diagonal aspect ratio: 1. perfect, above 1. gives stretching factor ------------ - minimum length of an edge in the whole mesh (m) = 93.749999999999091 + minimum length of an edge in the whole mesh (m) = 93.749999999999091 - maximum length of an edge in the whole mesh (m) = 250.00000000000000 + maximum length of an edge in the whole mesh (m) = 250.00000000000000 *** *** max equiangle skewness = 5.5129614478700954E-015 in element 8 of slice 0 @@ -207,20 +207,20 @@ max deviation angle from a right angle (90 degrees) is therefore = 4.9616653030830860E-013 - worst angle in the mesh is therefore 89.999999999999503 + worst angle in the mesh is therefore 89.999999999999503 or 90.000000000000497 degrees - max edge aspect ratio = 2.2222222222222503 + max edge aspect ratio = 2.2222222222222503 - max diagonal aspect ratio = 1.0000000000000036 + max diagonal aspect ratio = 1.0000000000000036 *** *** Maximum suggested time step for simulation = 0.00169220 *** - *** Max CFL stability condition of the time scheme (must be below about 0.55 or so) = 0.47039999999999998 - *** computed using the maximum P wave velocity = 4500.0000000000000 + *** Max CFL stability condition of the time scheme (must be below about 0.55 or so) = 0.47039999999999998 + *** computed using the maximum P wave velocity = 4500.0000000000000 *** - that value is below the upper CFL limit of 0.55000000000000004 + that value is below the upper CFL limit of 0.55000000000000004 therefore the run should be stable creating histogram of mesh quality diff --git a/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_solver.txt b/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_solver.txt index 1ee85b8f3..68520d79b 100644 --- a/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_solver.txt +++ b/EXAMPLES/meshfem3D_examples/sep_bathymetry/REF_SEIS/output_solver.txt @@ -60,33 +60,33 @@ ******** - Model: P velocity min,max = 1499.99988 4500.00000 - Model: S velocity min,max = 0.00000000 2590.00000 + Model: P velocity min,max = 1499.99988 4500.00000 + Model: S velocity min,max = 0.00000000 2590.00000 - Model: Poisson's ratio min,max = 0.250554323 0.500000000 + Model: Poisson's ratio min,max = 0.250554323 0.500000000 ******** ********************************************* *** Verification of simulation parameters *** ********************************************* - *** Xmin and Xmax of the model = 0.00000000 10000.0000 - *** Ymin and Ymax of the model = 1600.00000 2400.00000 - *** Zmin and Zmax of the model = -5000.00000 0.00000000 + *** Xmin and Xmax of the model = 0.00000000 10000.0000 + *** Ymin and Ymax of the model = 1600.00000 2400.00000 + *** Zmin and Zmax of the model = -5000.00000 0.00000000 - *** Max GLL point distance = 81.8317108 - *** Min GLL point distance = 16.1879883 - *** Max/min ratio = 5.05508852 + *** Max GLL point distance = 81.8317108 + *** Min GLL point distance = 16.1879883 + *** Max/min ratio = 5.05508852 - *** Max element size = 250.000000 - *** Min element size = 93.7500000 - *** Max/min ratio = 2.66666675 + *** Max element size = 250.000000 + *** Min element size = 93.7500000 + *** Max/min ratio = 2.66666675 - *** Minimum period resolved = 0.226449981 + *** Minimum period resolved = 0.226449981 *** Maximum suggested time step = 1.79866538E-03 *** for DT : 1.8000000000000000E-003 - *** Max stability for wave velocities = 0.500370979 + *** Max stability for wave velocities = 0.500370979 Elapsed time for checking mesh resolution in seconds = 3.5490000000000001E-003 saving VTK files for Courant number and minimum period @@ -123,48 +123,48 @@ in element 3165 in acoustic domain - using force point source: - xi coordinate of source in that element: -1.0000000000000000 - eta coordinate of source in that element: 1.0000000000000000 - gamma coordinate of source in that element: 0.92000000000000004 + using force point source: + xi coordinate of source in that element: -1.0000000000000000 + eta coordinate of source in that element: 1.0000000000000000 + gamma coordinate of source in that element: 0.92000000000000004 - component of direction vector in East direction: 0.00000000 - component of direction vector in North direction: 0.00000000 - component of direction vector in Vertical direction: -1.00000000 + component of direction vector in East direction: 0.00000000 + component of direction vector in North direction: 0.00000000 + component of direction vector in Vertical direction: -1.00000000 - nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 - nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 - nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 + nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 + nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 + nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 - at (x,y,z) coordinates = 2500.0000000000000 2000.0000000000000 -9.9999999999999947 + at (x,y,z) coordinates = 2500.0000000000000 2000.0000000000000 -9.9999999999999947 source time function: using Ricker source time function - using a source of dominant frequency 4.0000000000000000 - t0_ricker = 0.29999999999999999 + using a source of dominant frequency 4.0000000000000000 + t0_ricker = 0.29999999999999999 Ricker frequency: 4.0000000000000000 Hz time shift: 0.0000000000000000 seconds magnitude of the source: - factor = 10000000000.000000 + factor = 10000000000.000000 original (requested) position of the source: - latitude: 2000.0000000000000 - longitude: 2500.0000000000000 + latitude: 2000.0000000000000 + longitude: 2500.0000000000000 - x: 2500.0000000000000 - y: 2000.0000000000000 + x: 2500.0000000000000 + y: 2000.0000000000000 depth: 1.0000000000000000E-002 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 2500.0000000000000 - y: 2000.0000000000000 + x: 2500.0000000000000 + y: 2000.0000000000000 depth: 9.9999999999999950E-003 km - z: -9.9999999999999947 + z: -9.9999999999999947 error in location of the source: 5.32907052E-15 m @@ -175,48 +175,48 @@ in element 3167 in acoustic domain - using force point source: - xi coordinate of source in that element: -0.20000046875018310 - eta coordinate of source in that element: 1.0000000000000000 - gamma coordinate of source in that element: 0.92000000000000004 + using force point source: + xi coordinate of source in that element: -0.20000046875018310 + eta coordinate of source in that element: 1.0000000000000000 + gamma coordinate of source in that element: 0.92000000000000004 - component of direction vector in East direction: 0.00000000 - component of direction vector in North direction: 0.00000000 - component of direction vector in Vertical direction: -1.00000000 + component of direction vector in East direction: 0.00000000 + component of direction vector in North direction: 0.00000000 + component of direction vector in Vertical direction: -1.00000000 - nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 - nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 - nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 + nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 + nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 + nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 - at (x,y,z) coordinates = 3000.0000000000000 2000.0000000000000 -9.9999999999999947 + at (x,y,z) coordinates = 3000.0000000000000 2000.0000000000000 -9.9999999999999947 source time function: using Ricker source time function - using a source of dominant frequency 4.0000000000000000 - t0_ricker = 0.29999999999999999 + using a source of dominant frequency 4.0000000000000000 + t0_ricker = 0.29999999999999999 Ricker frequency: 4.0000000000000000 Hz time shift: 0.20000000000000001 seconds magnitude of the source: - factor = 10000000000.000000 + factor = 10000000000.000000 original (requested) position of the source: - latitude: 2000.0000000000000 - longitude: 3000.0000000000000 + latitude: 2000.0000000000000 + longitude: 3000.0000000000000 - x: 3000.0000000000000 - y: 2000.0000000000000 + x: 3000.0000000000000 + y: 2000.0000000000000 depth: 1.0000000000000000E-002 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 3000.0000000000000 - y: 2000.0000000000000 + x: 3000.0000000000000 + y: 2000.0000000000000 depth: 9.9999999999999950E-003 km - z: -9.9999999999999947 + z: -9.9999999999999947 error in location of the source: 5.32907052E-15 m @@ -227,48 +227,48 @@ in element 3153 in acoustic domain - using force point source: - xi coordinate of source in that element: 0.59999953125036853 - eta coordinate of source in that element: 1.0000000000000000 - gamma coordinate of source in that element: 0.92000000000000004 + using force point source: + xi coordinate of source in that element: 0.59999953125036853 + eta coordinate of source in that element: 1.0000000000000000 + gamma coordinate of source in that element: 0.92000000000000004 - component of direction vector in East direction: 0.00000000 - component of direction vector in North direction: 0.00000000 - component of direction vector in Vertical direction: -1.00000000 + component of direction vector in East direction: 0.00000000 + component of direction vector in North direction: 0.00000000 + component of direction vector in Vertical direction: -1.00000000 - nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 - nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 - nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 + nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 + nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 + nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 - at (x,y,z) coordinates = 3500.0000000000000 2000.0000000000000 -9.9999999999999964 + at (x,y,z) coordinates = 3500.0000000000000 2000.0000000000000 -9.9999999999999964 source time function: using Ricker source time function - using a source of dominant frequency 4.0000000000000000 - t0_ricker = 0.29999999999999999 + using a source of dominant frequency 4.0000000000000000 + t0_ricker = 0.29999999999999999 Ricker frequency: 4.0000000000000000 Hz time shift: 0.40000000000000002 seconds magnitude of the source: - factor = 10000000000.000000 + factor = 10000000000.000000 original (requested) position of the source: - latitude: 2000.0000000000000 - longitude: 3500.0000000000000 + latitude: 2000.0000000000000 + longitude: 3500.0000000000000 - x: 3500.0000000000000 - y: 2000.0000000000000 + x: 3500.0000000000000 + y: 2000.0000000000000 depth: 1.0000000000000000E-002 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 3500.0000000000000 - y: 2000.0000000000000 + x: 3500.0000000000000 + y: 2000.0000000000000 depth: 9.9999999999999967E-003 km - z: -9.9999999999999964 + z: -9.9999999999999964 error in location of the source: 3.55271368E-15 m @@ -279,48 +279,48 @@ in element 3156 in acoustic domain - using force point source: - xi coordinate of source in that element: -0.59999906249963331 - eta coordinate of source in that element: 1.0000000000000000 - gamma coordinate of source in that element: 0.92000000000000004 + using force point source: + xi coordinate of source in that element: -0.59999906249963331 + eta coordinate of source in that element: 1.0000000000000000 + gamma coordinate of source in that element: 0.92000000000000004 - component of direction vector in East direction: 0.00000000 - component of direction vector in North direction: 0.00000000 - component of direction vector in Vertical direction: -1.00000000 + component of direction vector in East direction: 0.00000000 + component of direction vector in North direction: 0.00000000 + component of direction vector in Vertical direction: -1.00000000 - nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 - nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 - nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 + nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 + nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 + nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 - at (x,y,z) coordinates = 4000.0000000000000 2000.0000000000000 -9.9999999999999964 + at (x,y,z) coordinates = 4000.0000000000000 2000.0000000000000 -9.9999999999999964 source time function: using Ricker source time function - using a source of dominant frequency 4.0000000000000000 - t0_ricker = 0.29999999999999999 + using a source of dominant frequency 4.0000000000000000 + t0_ricker = 0.29999999999999999 Ricker frequency: 4.0000000000000000 Hz time shift: 0.59999999999999998 seconds magnitude of the source: - factor = 10000000000.000000 + factor = 10000000000.000000 original (requested) position of the source: - latitude: 2000.0000000000000 - longitude: 4000.0000000000000 + latitude: 2000.0000000000000 + longitude: 4000.0000000000000 - x: 4000.0000000000000 - y: 2000.0000000000000 + x: 4000.0000000000000 + y: 2000.0000000000000 depth: 1.0000000000000000E-002 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 4000.0000000000000 - y: 2000.0000000000000 + x: 4000.0000000000000 + y: 2000.0000000000000 depth: 9.9999999999999967E-003 km - z: -9.9999999999999964 + z: -9.9999999999999964 error in location of the source: 3.55271368E-15 m @@ -331,48 +331,48 @@ in element 3158 in acoustic domain - using force point source: - xi coordinate of source in that element: 0.19999906250073243 - eta coordinate of source in that element: 1.0000000000000000 - gamma coordinate of source in that element: 0.92000000000000004 + using force point source: + xi coordinate of source in that element: 0.19999906250073243 + eta coordinate of source in that element: 1.0000000000000000 + gamma coordinate of source in that element: 0.92000000000000004 - component of direction vector in East direction: 0.00000000 - component of direction vector in North direction: 0.00000000 - component of direction vector in Vertical direction: -1.00000000 + component of direction vector in East direction: 0.00000000 + component of direction vector in North direction: 0.00000000 + component of direction vector in Vertical direction: -1.00000000 - nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 - nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 - nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 + nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 + nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 + nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 - at (x,y,z) coordinates = 4500.0000000000000 2000.0000000000000 -9.9999999999999964 + at (x,y,z) coordinates = 4500.0000000000000 2000.0000000000000 -9.9999999999999964 source time function: using Ricker source time function - using a source of dominant frequency 4.0000000000000000 - t0_ricker = 0.29999999999999999 + using a source of dominant frequency 4.0000000000000000 + t0_ricker = 0.29999999999999999 Ricker frequency: 4.0000000000000000 Hz time shift: 0.80000000000000004 seconds magnitude of the source: - factor = 10000000000.000000 + factor = 10000000000.000000 original (requested) position of the source: - latitude: 2000.0000000000000 - longitude: 4500.0000000000000 + latitude: 2000.0000000000000 + longitude: 4500.0000000000000 - x: 4500.0000000000000 - y: 2000.0000000000000 + x: 4500.0000000000000 + y: 2000.0000000000000 depth: 1.0000000000000000E-002 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 4500.0000000000000 - y: 2000.0000000000000 + x: 4500.0000000000000 + y: 2000.0000000000000 depth: 9.9999999999999967E-003 km - z: -9.9999999999999964 + z: -9.9999999999999964 error in location of the source: 3.55271368E-15 m @@ -383,48 +383,48 @@ in element 3160 in acoustic domain - using force point source: - xi coordinate of source in that element: 1.0000000000000000 - eta coordinate of source in that element: 1.0000000000000000 - gamma coordinate of source in that element: 0.92000000000000004 + using force point source: + xi coordinate of source in that element: 1.0000000000000000 + eta coordinate of source in that element: 1.0000000000000000 + gamma coordinate of source in that element: 0.92000000000000004 - component of direction vector in East direction: 0.00000000 - component of direction vector in North direction: 0.00000000 - component of direction vector in Vertical direction: -1.00000000 + component of direction vector in East direction: 0.00000000 + component of direction vector in North direction: 0.00000000 + component of direction vector in Vertical direction: -1.00000000 - nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 - nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 - nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 + nu1 = 1.0000000000000000 0.0000000000000000 0.0000000000000000 + nu2 = 0.0000000000000000 1.0000000000000000 0.0000000000000000 + nu3 = 0.0000000000000000 0.0000000000000000 1.0000000000000000 - at (x,y,z) coordinates = 5000.0000000000000 2000.0000000000000 -9.9999999999999947 + at (x,y,z) coordinates = 5000.0000000000000 2000.0000000000000 -9.9999999999999947 source time function: using Ricker source time function - using a source of dominant frequency 4.0000000000000000 - t0_ricker = 0.29999999999999999 + using a source of dominant frequency 4.0000000000000000 + t0_ricker = 0.29999999999999999 Ricker frequency: 4.0000000000000000 Hz time shift: 1.0000000000000000 seconds magnitude of the source: - factor = 10000000000.000000 + factor = 10000000000.000000 original (requested) position of the source: - latitude: 2000.0000000000000 - longitude: 5000.0000000000000 + latitude: 2000.0000000000000 + longitude: 5000.0000000000000 - x: 5000.0000000000000 - y: 2000.0000000000000 + x: 5000.0000000000000 + y: 2000.0000000000000 depth: 1.0000000000000000E-002 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 5000.0000000000000 - y: 2000.0000000000000 + x: 5000.0000000000000 + y: 2000.0000000000000 depth: 9.9999999999999950E-003 km - z: -9.9999999999999947 + z: -9.9999999999999947 error in location of the source: 5.32907052E-15 m @@ -461,254 +461,254 @@ station # 1 NETWORK STATION_3 - original latitude: 2000.00000 - original longitude: 1000.00000 - original x: 1000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 1000.00000 + original x: 1000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 1.50000000 - target x, y, z: 1000.00000 2000.00000 -500.000000 + horizontal distance: 1.50000000 + target x, y, z: 1000.00000 2000.00000 -500.000000 closest estimate found: 2.60489276E-13 m away receiver located in slice 0 in element 3029 in elastic domain - at coordinates: - xi = 0.60000035156253351 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 999.99999999999989 - y: 1999.9999999999998 + at coordinates: + xi = 0.60000035156253351 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 999.99999999999989 + y: 1999.9999999999998 depth: 499.99999999999994 m - z: -499.99999999999994 + z: -499.99999999999994 station # 2 NETWORK STATION_58 - original latitude: 2000.00000 - original longitude: 2000.00000 - original x: 2000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 2000.00000 + original x: 2000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 0.500000000 - target x, y, z: 2000.00000 2000.00000 -500.000000 + horizontal distance: 0.500000000 + target x, y, z: 2000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 0 in element 3034 in elastic domain - at coordinates: - xi = 0.20000046875018310 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 2000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = 0.20000046875018310 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 2000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 3 NETWORK STATION_114 - original latitude: 2000.00000 - original longitude: 3000.00000 - original x: 3000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 3000.00000 + original x: 3000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 0.500000000 - target x, y, z: 3000.00000 2000.00000 -500.000000 + horizontal distance: 0.500000000 + target x, y, z: 3000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 0 in element 3039 in elastic domain - at coordinates: - xi = -0.20000046875018310 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 3000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = -0.20000046875018310 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 3000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 4 NETWORK STATION_169 - original latitude: 2000.00000 - original longitude: 4000.00000 - original x: 4000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 4000.00000 + original x: 4000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 1.50000000 - target x, y, z: 4000.00000 2000.00000 -500.000000 + horizontal distance: 1.50000000 + target x, y, z: 4000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 1 in element 3028 in elastic domain - at coordinates: - xi = -0.59999906249963353 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 4000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = -0.59999906249963353 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 4000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 5 NETWORK STATION_225 - original latitude: 2000.00000 - original longitude: 5000.00000 - original x: 5000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 5000.00000 + original x: 5000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 2.50000000 - target x, y, z: 5000.00000 2000.00000 -500.000000 + horizontal distance: 2.50000000 + target x, y, z: 5000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 1 in element 3032 in elastic domain - at coordinates: - xi = 1.0000000000000000 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 5000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = 1.0000000000000000 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 5000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 6 NETWORK STATION_280 - original latitude: 2000.00000 - original longitude: 6000.00000 - original x: 6000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 6000.00000 + original x: 6000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 3.50000000 - target x, y, z: 6000.00000 2000.00000 -500.000000 + horizontal distance: 3.50000000 + target x, y, z: 6000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 1 in element 3053 in elastic domain - at coordinates: - xi = 0.60000093750145689 - eta = -1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 6000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = 0.60000093750145689 + eta = -1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 6000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 7 NETWORK STATION_336 - original latitude: 2000.00000 - original longitude: 7000.00000 - original x: 7000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 7000.00000 + original x: 7000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 4.50000000 - target x, y, z: 7000.00000 2000.00000 -500.000000 + horizontal distance: 4.50000000 + target x, y, z: 7000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 2 in element 3026 in elastic domain - at coordinates: - xi = 0.19999906250073243 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 7000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = 0.19999906250073243 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 7000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 8 NETWORK STATION_391 - original latitude: 2000.00000 - original longitude: 8000.00000 - original x: 8000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 8000.00000 + original x: 8000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 5.50000000 - target x, y, z: 8000.00000 2000.00000 -500.000000 + horizontal distance: 5.50000000 + target x, y, z: 8000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 2 in element 3031 in elastic domain - at coordinates: - xi = -0.19999906250073243 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 8000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = -0.19999906250073243 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 8000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 station # 9 NETWORK STATION_447 - original latitude: 2000.00000 - original longitude: 9000.00000 - original x: 9000.00000 - original y: 2000.00000 + original latitude: 2000.00000 + original longitude: 9000.00000 + original x: 9000.00000 + original y: 2000.00000 original depth: 500.000000 m - horizontal distance: 6.50000000 - target x, y, z: 9000.00000 2000.00000 -500.000000 + horizontal distance: 6.50000000 + target x, y, z: 9000.00000 2000.00000 -500.000000 closest estimate found: 0.00000000 m away receiver located in slice 2 in element 3036 in elastic domain - at coordinates: - xi = -0.59999812500585847 - eta = 1.0000000000000000 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 9000.0000000000000 - y: 2000.0000000000000 + at coordinates: + xi = -0.59999812500585847 + eta = 1.0000000000000000 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 9000.0000000000000 + y: 2000.0000000000000 depth: 500.00000000000000 m - z: -500.00000000000000 + z: -500.00000000000000 maximum error in location of all the receivers: 2.60489276E-13 m @@ -766,7 +766,7 @@ number of regular shaped elements : 0 number of irregular shaped elements: 3200 fused array done - bandwidth test (STREAM TRIAD): + bandwidth test (STREAM TRIAD): memory accesses = 7.62537003 MB timing min/max = 6.44000014E-04 s / 7.02999998E-04 s timing avg = 6.56299992E-04 s @@ -795,16 +795,16 @@ Time step # 5 Time: -0.292800009 seconds - Elapsed time in seconds = 0.10396699999999998 + Elapsed time in seconds = 0.10396699999999998 Elapsed time in hh:mm:ss = 0 h 00 m 00 s Mean elapsed time per time step in seconds = 2.07934007E-02 Max norm displacement vector U in all slices (m) = 7.94038840E-07 - Max norm pressure P in all slices (Pa) = 21.0943966 + Max norm pressure P in all slices (Pa) = 21.0943966 Time steps done = 5 out of 2000 Time steps remaining = 1995 - Estimated remaining time in seconds = 41.4828339 + Estimated remaining time in seconds = 41.4828339 Estimated remaining time in hh:mm:ss = 0 h 00 m 41 s - Estimated total run time in seconds = 41.5867996 + Estimated total run time in seconds = 41.5867996 Estimated total run time in hh:mm:ss = 0 h 00 m 41 s We have done 0.250000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 @@ -815,160 +815,160 @@ Time step # 200 Time: 5.82000017E-02 seconds - Elapsed time in seconds = 5.0051079999999999 + Elapsed time in seconds = 5.0051079999999999 Elapsed time in hh:mm:ss = 0 h 00 m 05 s Mean elapsed time per time step in seconds = 2.50255391E-02 - Max norm displacement vector U in all slices (m) = 0.655175209 - Max norm pressure P in all slices (Pa) = 237044.547 + Max norm displacement vector U in all slices (m) = 0.655175209 + Max norm pressure P in all slices (Pa) = 237044.547 Time steps done = 200 out of 2000 Time steps remaining = 1800 - Estimated remaining time in seconds = 45.0459709 + Estimated remaining time in seconds = 45.0459709 Estimated remaining time in hh:mm:ss = 0 h 00 m 45 s - Estimated total run time in seconds = 50.0510788 + Estimated total run time in seconds = 50.0510788 Estimated total run time in hh:mm:ss = 0 h 00 m 50 s We have done 10.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 400 Time: 0.418199986 seconds - Elapsed time in seconds = 10.318274000000001 + Elapsed time in seconds = 10.318274000000001 Elapsed time in hh:mm:ss = 0 h 00 m 10 s Mean elapsed time per time step in seconds = 2.57956851E-02 - Max norm displacement vector U in all slices (m) = 0.803838432 - Max norm pressure P in all slices (Pa) = 50770.8516 + Max norm displacement vector U in all slices (m) = 0.803838432 + Max norm pressure P in all slices (Pa) = 50770.8516 Time steps done = 400 out of 2000 Time steps remaining = 1600 - Estimated remaining time in seconds = 41.2730942 + Estimated remaining time in seconds = 41.2730942 Estimated remaining time in hh:mm:ss = 0 h 00 m 41 s - Estimated total run time in seconds = 51.5913696 + Estimated total run time in seconds = 51.5913696 Estimated total run time in hh:mm:ss = 0 h 00 m 51 s We have done 20.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 600 Time: 0.778199971 seconds - Elapsed time in seconds = 15.983793000000000 + Elapsed time in seconds = 15.983793000000000 Elapsed time in hh:mm:ss = 0 h 00 m 15 s Mean elapsed time per time step in seconds = 2.66396552E-02 - Max norm displacement vector U in all slices (m) = 0.671937823 - Max norm pressure P in all slices (Pa) = 55502.1211 + Max norm displacement vector U in all slices (m) = 0.671937823 + Max norm pressure P in all slices (Pa) = 55502.1211 Time steps done = 600 out of 2000 Time steps remaining = 1400 - Estimated remaining time in seconds = 37.2955170 + Estimated remaining time in seconds = 37.2955170 Estimated remaining time in hh:mm:ss = 0 h 00 m 37 s - Estimated total run time in seconds = 53.2793083 + Estimated total run time in seconds = 53.2793083 Estimated total run time in hh:mm:ss = 0 h 00 m 53 s We have done 30.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 800 Time: 1.13820004 seconds - Elapsed time in seconds = 21.545078999999998 + Elapsed time in seconds = 21.545078999999998 Elapsed time in hh:mm:ss = 0 h 00 m 21 s Mean elapsed time per time step in seconds = 2.69313492E-02 - Max norm displacement vector U in all slices (m) = 0.174651608 - Max norm pressure P in all slices (Pa) = 64505.8008 + Max norm displacement vector U in all slices (m) = 0.174651608 + Max norm pressure P in all slices (Pa) = 64505.8008 Time steps done = 800 out of 2000 Time steps remaining = 1200 - Estimated remaining time in seconds = 32.3176193 + Estimated remaining time in seconds = 32.3176193 Estimated remaining time in hh:mm:ss = 0 h 00 m 32 s - Estimated total run time in seconds = 53.8626976 + Estimated total run time in seconds = 53.8626976 Estimated total run time in hh:mm:ss = 0 h 00 m 53 s We have done 40.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 1000 Time: 1.49820006 seconds - Elapsed time in seconds = 27.196864999999999 + Elapsed time in seconds = 27.196864999999999 Elapsed time in hh:mm:ss = 0 h 00 m 27 s Mean elapsed time per time step in seconds = 2.71968655E-02 Max norm displacement vector U in all slices (m) = 2.76113432E-02 Max norm pressure P in all slices (Pa) = 1.15376422E-02 Time steps done = 1000 out of 2000 Time steps remaining = 1000 - Estimated remaining time in seconds = 27.1968651 + Estimated remaining time in seconds = 27.1968651 Estimated remaining time in hh:mm:ss = 0 h 00 m 27 s - Estimated total run time in seconds = 54.3937302 + Estimated total run time in seconds = 54.3937302 Estimated total run time in hh:mm:ss = 0 h 00 m 54 s We have done 50.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:07 Time step # 1200 Time: 1.85819995 seconds - Elapsed time in seconds = 32.842962999999997 + Elapsed time in seconds = 32.842962999999997 Elapsed time in hh:mm:ss = 0 h 00 m 32 s Mean elapsed time per time step in seconds = 2.73691360E-02 Max norm displacement vector U in all slices (m) = 1.57053489E-02 - Max norm pressure P in all slices (Pa) = 757.377686 + Max norm pressure P in all slices (Pa) = 757.377686 Time steps done = 1200 out of 2000 Time steps remaining = 800 - Estimated remaining time in seconds = 21.8953094 + Estimated remaining time in seconds = 21.8953094 Estimated remaining time in hh:mm:ss = 0 h 00 m 21 s - Estimated total run time in seconds = 54.7382698 + Estimated total run time in seconds = 54.7382698 Estimated total run time in hh:mm:ss = 0 h 00 m 54 s We have done 60.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:07 Time step # 1400 Time: 2.21819997 seconds - Elapsed time in seconds = 38.430399999999999 + Elapsed time in seconds = 38.430399999999999 Elapsed time in hh:mm:ss = 0 h 00 m 38 s Mean elapsed time per time step in seconds = 2.74502859E-02 Max norm displacement vector U in all slices (m) = 7.63012143E-03 - Max norm pressure P in all slices (Pa) = 1949.29749 + Max norm pressure P in all slices (Pa) = 1949.29749 Time steps done = 1400 out of 2000 Time steps remaining = 600 - Estimated remaining time in seconds = 16.4701710 + Estimated remaining time in seconds = 16.4701710 Estimated remaining time in hh:mm:ss = 0 h 00 m 16 s - Estimated total run time in seconds = 54.9005699 + Estimated total run time in seconds = 54.9005699 Estimated total run time in hh:mm:ss = 0 h 00 m 54 s We have done 70.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 1600 Time: 2.57820010 seconds - Elapsed time in seconds = 44.207546000000001 + Elapsed time in seconds = 44.207546000000001 Elapsed time in hh:mm:ss = 0 h 00 m 44 s Mean elapsed time per time step in seconds = 2.76297163E-02 Max norm displacement vector U in all slices (m) = 8.30366276E-03 - Max norm pressure P in all slices (Pa) = 2079.50928 + Max norm pressure P in all slices (Pa) = 2079.50928 Time steps done = 1600 out of 2000 Time steps remaining = 400 - Estimated remaining time in seconds = 11.0518866 + Estimated remaining time in seconds = 11.0518866 Estimated remaining time in hh:mm:ss = 0 h 00 m 11 s - Estimated total run time in seconds = 55.2594337 + Estimated total run time in seconds = 55.2594337 Estimated total run time in hh:mm:ss = 0 h 00 m 55 s We have done 80.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 1800 Time: 2.93820000 seconds - Elapsed time in seconds = 50.061988999999997 + Elapsed time in seconds = 50.061988999999997 Elapsed time in hh:mm:ss = 0 h 00 m 50 s Mean elapsed time per time step in seconds = 2.78122164E-02 Max norm displacement vector U in all slices (m) = 1.06844110E-02 - Max norm pressure P in all slices (Pa) = 1958.82959 + Max norm pressure P in all slices (Pa) = 1958.82959 Time steps done = 1800 out of 2000 Time steps remaining = 200 - Estimated remaining time in seconds = 5.56244326 + Estimated remaining time in seconds = 5.56244326 Estimated remaining time in hh:mm:ss = 0 h 00 m 05 s - Estimated total run time in seconds = 55.6244316 + Estimated total run time in seconds = 55.6244316 Estimated total run time in hh:mm:ss = 0 h 00 m 55 s We have done 90.0000000 % of that The run will finish approximately on (in local time): Sun Apr 25, 2021 16:08 Time step # 2000 Time: 3.29819989 seconds - Elapsed time in seconds = 56.035246000000001 + Elapsed time in seconds = 56.035246000000001 Elapsed time in hh:mm:ss = 0 h 00 m 56 s Mean elapsed time per time step in seconds = 2.80176234E-02 Max norm displacement vector U in all slices (m) = 7.10760430E-03 - Max norm pressure P in all slices (Pa) = 1644.60815 + Max norm pressure P in all slices (Pa) = 1644.60815 Time steps done = 2000 out of 2000 Time steps remaining = 0 - Estimated remaining time in seconds = 0.00000000 + Estimated remaining time in seconds = 0.00000000 Estimated remaining time in hh:mm:ss = 0 h 00 m 00 s - Estimated total run time in seconds = 56.0352478 + Estimated total run time in seconds = 56.0352478 Estimated total run time in hh:mm:ss = 0 h 00 m 56 s We have done 100.000000 % of that @@ -982,7 +982,7 @@ Writing the seismograms by main proc alone took 0.101990998 seconds Time loop finished. Timing info: - Total elapsed time in seconds = 56.166803999999999 + Total elapsed time in seconds = 56.166803999999999 Total elapsed time in hh:mm:ss = 0 h 00 m 56 s End of the simulation diff --git a/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_generate_databases.txt b/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_generate_databases.txt index 4a7eaec92..be9c93d5e 100644 --- a/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_generate_databases.txt +++ b/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_generate_databases.txt @@ -21,7 +21,7 @@ Surface shape functions defined by NGNOD2D = 4 control nodes Beware! Curvature (i.e. HEX27 elements) is not handled by our internal mesher - velocity model: default + velocity model: default suppressing UTM projection @@ -46,7 +46,7 @@ defined materials : 2 undefined materials : 0 total number of spectral elements: 23040 - absorbing boundaries: + absorbing boundaries: xmin,xmax : 480 480 ymin,ymax : 480 480 bottom,top: 2304 2304 @@ -60,14 +60,14 @@ create regions: - ...allocating arrays + ...allocating arrays NGLLX = 5 NGLLY = 5 NGLLZ = 5 NGNOD = 8 NGNOD2D = 4 - main process setup: + main process setup: nspec = 5760 separating regular/irregular element shapes @@ -84,16 +84,16 @@ File DATA/Par_file_faults not found: assuming that there are no faults - ...setting up jacobian + ...setting up jacobian ...indexing global points - creating ibool indexing : x min/max = 0.00000000 / 30000.0000 + creating ibool indexing : x min/max = 0.00000000 / 30000.0000 creating indirect addressing: nglob = 385769 creating unique point locations - ...preparing MPI interfaces + ...preparing MPI interfaces number of interfaces : 3 - creating MPI indexing : x min/max = 0.00000000 / 30000.0000 + creating MPI indexing : x min/max = 0.00000000 / 30000.0000 tolerance = 3.0000000000000001E-006 total MPI interface points: 31980 @@ -124,7 +124,7 @@ 90 % time remaining: 2.9842423744952404E-008 s 100 % time remaining: 0.0000000000000000 s - ...detecting acoustic-elastic-poroelastic surfaces + ...detecting acoustic-elastic-poroelastic surfaces total acoustic elements : 0 total elastic elements : 23040 total poroelastic elements: 0 @@ -134,18 +134,18 @@ elastic - poroelastic coupling : total number of faces = 0 - ...element inner/outer separation + ...element inner/outer separation for overlapping of communications with calculations: percentage of edge elements 8.15972137 % percentage of volume elements 91.8402786 % - ...element mesh coloring + ...element mesh coloring use coloring = F - ...external binary models - no external binary model used + ...external binary models + no external binary model used - ...creating mass matrix + ...creating mass matrix ...saving databases using binary file format @@ -178,29 +178,29 @@ ******** - Model: P velocity min,max = 4100.00000 5800.00000 - Model: S velocity min,max = 2400.00000 3350.00000 + Model: P velocity min,max = 4100.00000 5800.00000 + Model: S velocity min,max = 2400.00000 3350.00000 - Model: Poisson's ratio min,max = 0.239366516 0.249693319 + Model: Poisson's ratio min,max = 0.239366516 0.249693319 ******** ********************************************* *** Verification of simulation parameters *** ********************************************* - *** Xmin and Xmax of the model = 0.00000000 30000.0000 - *** Ymin and Ymax of the model = 0.00000000 25000.0000 - *** Zmin and Zmax of the model = -6000.00000 0.00000000 + *** Xmin and Xmax of the model = 0.00000000 30000.0000 + *** Ymin and Ymax of the model = 0.00000000 25000.0000 + *** Zmin and Zmax of the model = -6000.00000 0.00000000 - *** Max GLL point distance = 204.580078 - *** Min GLL point distance = 89.9335938 - *** Max/min ratio = 2.27479053 + *** Max GLL point distance = 204.580078 + *** Min GLL point distance = 89.9335938 + *** Max/min ratio = 2.27479053 - *** Max element size = 625.000000 - *** Min element size = 520.832031 - *** Max/min ratio = 1.20000303 + *** Max element size = 625.000000 + *** Min element size = 520.832031 + *** Max/min ratio = 1.20000303 - *** Minimum period resolved = 0.325520813 + *** Minimum period resolved = 0.325520813 *** Maximum suggested time step = 7.75289582E-03 Elapsed time for checking mesh resolution in seconds = 1.0793000000000497E-002 @@ -209,7 +209,7 @@ mesh regions done - min and max of elevation (i.e. height of the upper surface of the mesh) included in mesh in m is 0.0000000000000000 0.0000000000000000 + min and max of elevation (i.e. height of the upper surface of the mesh) included in mesh in m is 0.0000000000000000 0.0000000000000000 done mesh setup @@ -234,7 +234,7 @@ smallest and largest possible floating-point numbers are: 1.17549435E-38 3.40282347E+38 - Elapsed time for mesh generation and buffer creation in seconds = 6.6023550000000002 + Elapsed time for mesh generation and buffer creation in seconds = 6.6023550000000002 End of mesh generation done diff --git a/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_meshfem3D.txt b/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_meshfem3D.txt index 69542dd0d..f3fb8cf9a 100644 --- a/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_meshfem3D.txt +++ b/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_meshfem3D.txt @@ -41,8 +41,8 @@ Creating global slice addressing Spatial distribution of slice numbers: - 2 3 - 0 1 + 2 3 + 0 1 This is process 0 There are 4 MPI processes Processes are numbered from 0 to 3 @@ -63,19 +63,19 @@ region selected: - latitude min = 0.0000000000000000 - latitude max = 25000.000000000000 + latitude min = 0.0000000000000000 + latitude max = 25000.000000000000 - longitude min = 0.0000000000000000 - longitude max = 30000.000000000000 + longitude min = 0.0000000000000000 + longitude max = 30000.000000000000 this is given directly as UTM - UTM X min = 0.0000000000000000 - UTM X max = 30000.000000000000 + UTM X min = 0.0000000000000000 + UTM X max = 30000.000000000000 - UTM Y min = 0.0000000000000000 - UTM Y max = 25000.000000000000 + UTM Y min = 0.0000000000000000 + UTM Y max = 25000.000000000000 UTM size of model along X is 30.000000000000000 km UTM size of model along Y is 25.000000000000000 km @@ -95,27 +95,27 @@ number of interfaces: 1 mesh: - origin UTM minimum x/y (m) = 0.00000000 0.00000000 - origin UTM maximum x/y (m) = 30000.0000 25000.0000 + origin UTM minimum x/y (m) = 0.00000000 0.00000000 + origin UTM maximum x/y (m) = 30000.0000 25000.0000 reading interface 1 interface file : interface1.dat number of points x/y = 2 2 - origin x/y (m) = 0.00000000 0.00000000 - spacing x/y (m) = 20000.0000 20000.0000 + origin x/y (m) = 0.00000000 0.00000000 + spacing x/y (m) = 20000.0000 20000.0000 - dimension x-direction (m) = 0.00000000 / 20000.0000 - dimension y-direction (m) = 0.00000000 / 20000.0000 + dimension x-direction (m) = 0.00000000 / 20000.0000 + dimension y-direction (m) = 0.00000000 / 20000.0000 total number of file points = 4 should be 4 this point total is okay - original elevation min/max = 0.00000000 0.00000000 - interpolated mesh elevation min/max = 0.00000000 0.00000000 + original elevation min/max = 0.00000000 0.00000000 + interpolated mesh elevation min/max = 0.00000000 0.00000000 - interpolated mesh UTM minimum x/y (m) = 0.00000000 0.00000000 - interpolated mesh UTM maximum x/y (m) = 30000.0000 25000.0000 + interpolated mesh UTM minimum x/y (m) = 0.00000000 0.00000000 + interpolated mesh UTM maximum x/y (m) = 30000.0000 25000.0000 ************************** @@ -140,16 +140,16 @@ creating mesh region 1 (regular mesh) mesh dimensions: - Xmin and Xmax of the model = 0.00000000 30000.0000 - Ymin and Ymax of the model = 0.00000000 25000.0000 - Zmin and Zmax of the model = -6000.00000 0.00000000 + Xmin and Xmax of the model = 0.00000000 30000.0000 + Ymin and Ymax of the model = 0.00000000 25000.0000 + Zmin and Zmax of the model = -6000.00000 0.00000000 exact area = 750000000. (m^2) = 750.000000 (km^2) Max element size = 625.000000 (m) Min element size = 520.832031 (m) - Max/min ratio = 1.20000303 + Max/min ratio = 1.20000303 creating indirect addressing for unstructured mesh @@ -167,7 +167,7 @@ ************************** start computing the minimum and maximum edge size - done processing + done processing ------------ mesh quality parameter definitions: @@ -178,9 +178,9 @@ diagonal aspect ratio: 1. perfect, above 1. gives stretching factor ------------ - minimum length of an edge in the whole mesh (m) = 520.83333333332848 + minimum length of an edge in the whole mesh (m) = 520.83333333332848 - maximum length of an edge in the whole mesh (m) = 625.00000000000728 + maximum length of an edge in the whole mesh (m) = 625.00000000000728 *** *** max equiangle skewness = 2.2617277734851675E-015 in element 25 of slice 0 @@ -188,20 +188,20 @@ max deviation angle from a right angle (90 degrees) is therefore = 2.0355549961366507E-013 - worst angle in the mesh is therefore 89.999999999999801 + worst angle in the mesh is therefore 89.999999999999801 or 90.000000000000199 degrees - max edge aspect ratio = 1.2000000000000250 + max edge aspect ratio = 1.2000000000000250 - max diagonal aspect ratio = 1.0000000000000033 + max diagonal aspect ratio = 1.0000000000000033 *** *** Maximum suggested time step for simulation = 0.00729395 *** - *** Max CFL stability condition of the time scheme (must be below about 0.55 or so) = 0.47039999999999998 - *** computed using the maximum P wave velocity = 5800.0000000000000 + *** Max CFL stability condition of the time scheme (must be below about 0.55 or so) = 0.47039999999999998 + *** computed using the maximum P wave velocity = 5800.0000000000000 *** - that value is below the upper CFL limit of 0.55000000000000004 + that value is below the upper CFL limit of 0.55000000000000004 therefore the run should be stable creating histogram of mesh quality @@ -250,7 +250,7 @@ smallest and largest possible floating-point numbers are: 1.17549435E-38 3.40282347E+38 - Elapsed time for mesh generation and buffer creation in seconds = 0.18683500000000000 + Elapsed time for mesh generation and buffer creation in seconds = 0.18683500000000000 End of mesh generation done diff --git a/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_solver.txt b/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_solver.txt index 0d4232907..0cf54b155 100644 --- a/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_solver.txt +++ b/EXAMPLES/small_adjoint_multiple_sources/REF_SEIS/output_solver.txt @@ -25,7 +25,7 @@ smallest and largest possible floating-point numbers are: 1.17549435E-38 3.40282347E+38 - velocity model: default + velocity model: default Reading mesh databases... reads binary mesh files: proc***_external_mesh.bin @@ -60,33 +60,33 @@ ******** - Model: P velocity min,max = 4100.00000 5800.00000 - Model: S velocity min,max = 2400.00000 3350.00000 + Model: P velocity min,max = 4100.00000 5800.00000 + Model: S velocity min,max = 2400.00000 3350.00000 - Model: Poisson's ratio min,max = 0.239366516 0.249693319 + Model: Poisson's ratio min,max = 0.239366516 0.249693319 ******** ********************************************* *** Verification of simulation parameters *** ********************************************* - *** Xmin and Xmax of the model = 0.00000000 30000.0000 - *** Ymin and Ymax of the model = 0.00000000 25000.0000 - *** Zmin and Zmax of the model = -6000.00000 0.00000000 + *** Xmin and Xmax of the model = 0.00000000 30000.0000 + *** Ymin and Ymax of the model = 0.00000000 25000.0000 + *** Zmin and Zmax of the model = -6000.00000 0.00000000 - *** Max GLL point distance = 204.580078 - *** Min GLL point distance = 89.9335938 - *** Max/min ratio = 2.27479053 + *** Max GLL point distance = 204.580078 + *** Min GLL point distance = 89.9335938 + *** Max/min ratio = 2.27479053 - *** Max element size = 625.000000 - *** Min element size = 520.832031 - *** Max/min ratio = 1.20000303 + *** Max element size = 625.000000 + *** Min element size = 520.832031 + *** Max/min ratio = 1.20000303 - *** Minimum period resolved = 0.325520813 + *** Minimum period resolved = 0.325520813 *** Maximum suggested time step = 7.75289582E-03 *** for DT : 8.0000000000000002E-003 - *** Max stability for wave velocities = 0.515936255 + *** Max stability for wave velocities = 0.515936255 Elapsed time for checking mesh resolution in seconds = 7.9629999999999996E-003 saving VTK files for Courant number and minimum period @@ -123,10 +123,10 @@ in element 5582 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7599999999999973E-002 - eta coordinate of source in that element: 0.39263950920061519 - gamma coordinate of source in that element: 1.0000000000000000 + eta coordinate of source in that element: 0.39263950920061519 + gamma coordinate of source in that element: 1.0000000000000000 source time function: using Ricker source time function @@ -135,24 +135,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 0.0000000000000000 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 0.0000000000000000 km - z: 0.0000000000000000 + z: 0.0000000000000000 error in location of the source: 0.00000000 m @@ -163,10 +163,10 @@ in element 5006 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7600000000000022E-002 - eta coordinate of source in that element: 0.39263950920061547 - gamma coordinate of source in that element: -0.33333333333333320 + eta coordinate of source in that element: 0.39263950920061547 + gamma coordinate of source in that element: -0.33333333333333320 source time function: using Ricker source time function @@ -175,24 +175,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 1.0000000000000000 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 1.0000000000000000 km - z: -1000.0000000000000 + z: -1000.0000000000000 error in location of the source: 0.00000000 m @@ -203,10 +203,10 @@ in element 3854 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7599999999999973E-002 - eta coordinate of source in that element: 0.39263950920061513 - gamma coordinate of source in that element: 0.33333333333333298 + eta coordinate of source in that element: 0.39263950920061513 + gamma coordinate of source in that element: 0.33333333333333298 source time function: using Ricker source time function @@ -215,24 +215,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 2.0000000000000000 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 2.0000000000000000 km - z: -2000.0000000000000 + z: -2000.0000000000000 error in location of the source: 0.00000000 m @@ -243,10 +243,10 @@ in element 2702 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7599999999999973E-002 - eta coordinate of source in that element: 0.39263950920061519 - gamma coordinate of source in that element: 1.0000000000000000 + eta coordinate of source in that element: 0.39263950920061519 + gamma coordinate of source in that element: 1.0000000000000000 source time function: using Ricker source time function @@ -255,24 +255,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 3.0000000000000000 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 3.0000000000000000 km - z: -3000.0000000000000 + z: -3000.0000000000000 error in location of the source: 0.00000000 m @@ -283,10 +283,10 @@ in element 2126 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7600000000000022E-002 - eta coordinate of source in that element: 0.39263950920061541 - gamma coordinate of source in that element: -0.33333333333333459 + eta coordinate of source in that element: 0.39263950920061541 + gamma coordinate of source in that element: -0.33333333333333459 source time function: using Ricker source time function @@ -295,24 +295,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 4.0000000000000000 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 4.0000000000000009 km - z: -4000.0000000000009 + z: -4000.0000000000009 error in location of the source: 9.09494702E-13 m @@ -323,10 +323,10 @@ in element 974 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7599999999999973E-002 - eta coordinate of source in that element: 0.39263950920061513 - gamma coordinate of source in that element: 0.33333333333333476 + eta coordinate of source in that element: 0.39263950920061513 + gamma coordinate of source in that element: 0.33333333333333476 source time function: using Ricker source time function @@ -335,24 +335,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000018 + x: 8443.0000000000000 + y: 8696.0000000000018 depth: 5.0000000000000000 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000018 + x: 8443.0000000000000 + y: 8696.0000000000018 depth: 5.0000000000000000 km - z: -5000.0000000000000 + z: -5000.0000000000000 error in location of the source: 1.81898940E-12 m @@ -363,10 +363,10 @@ in element 398 in elastic domain - using moment tensor source: + using moment tensor source: xi coordinate of source in that element: 1.7600000000000022E-002 - eta coordinate of source in that element: 0.39263950920060847 - gamma coordinate of source in that element: -0.33333333333333165 + eta coordinate of source in that element: 0.39263950920060847 + gamma coordinate of source in that element: -0.33333333333333165 source time function: using Ricker source time function @@ -375,24 +375,24 @@ magnitude of the source: scalar moment M0 = 12247448.713915890 dyne-cm - moment magnitude Mw = -5.9746360562465766 + moment magnitude Mw = -5.9746360562465766 original (requested) position of the source: - latitude: 8696.0000000000000 - longitude: 8443.0000000000000 + latitude: 8696.0000000000000 + longitude: 8443.0000000000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 5.7999999999999998 km - topo elevation: 0.0000000000000000 + topo elevation: 0.0000000000000000 position of the source that will be used: - x: 8443.0000000000000 - y: 8696.0000000000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 5.7999999999999980 km - z: -5799.9999999999982 + z: -5799.9999999999982 error in location of the source: 1.81898940E-12 m @@ -429,30 +429,30 @@ station # 1 EX GRO1S - original latitude: 8696.00000 - original longitude: 8443.00000 - original x: 8443.00000 - original y: 8696.00000 + original latitude: 8696.00000 + original longitude: 8443.00000 + original x: 8443.00000 + original y: 8696.00000 original depth: 0.00000000 m - horizontal distance: 0.00000000 - target x, y, z: 8443.00000 8696.00000 0.00000000 + horizontal distance: 0.00000000 + target x, y, z: 8443.00000 8696.00000 0.00000000 closest estimate found: 0.00000000 m away receiver located in slice 0 in element 5582 in elastic domain - at coordinates: + at coordinates: xi = 1.7599999999999973E-002 - eta = 0.39263950920061519 - gamma = 1.0000000000000000 - rotation matrix: - nu1 = 1.00000000 0.00000000 0.00000000 - nu2 = 0.00000000 1.00000000 0.00000000 - nu3 = 0.00000000 0.00000000 1.00000000 - x: 8443.0000000000000 - y: 8696.0000000000000 + eta = 0.39263950920061519 + gamma = 1.0000000000000000 + rotation matrix: + nu1 = 1.00000000 0.00000000 0.00000000 + nu2 = 0.00000000 1.00000000 0.00000000 + nu3 = 0.00000000 0.00000000 1.00000000 + x: 8443.0000000000000 + y: 8696.0000000000000 depth: 0.0000000000000000 m - z: 0.0000000000000000 + z: 0.0000000000000000 maximum error in location of all the receivers: 0.00000000 m @@ -514,7 +514,7 @@ number of regular shaped elements : 0 number of irregular shaped elements: 5760 fused array done - bandwidth test (STREAM TRIAD): + bandwidth test (STREAM TRIAD): memory accesses = 13.2443275 MB timing min/max = 1.49599998E-03 s / 1.99500006E-03 s timing avg = 1.67779997E-03 s @@ -543,15 +543,15 @@ Time step # 5 Time: -0.208000004 seconds - Elapsed time in seconds = 0.21394099999999994 + Elapsed time in seconds = 0.21394099999999994 Elapsed time in hh:mm:ss = 0 h 00 m 00 s Mean elapsed time per time step in seconds = 4.27882001E-02 Max norm displacement vector U in all slices (m) = 1.94722578E-16 Time steps done = 5 out of 500 Time steps remaining = 495 - Estimated remaining time in seconds = 21.1801586 + Estimated remaining time in seconds = 21.1801586 Estimated remaining time in hh:mm:ss = 0 h 00 m 21 s - Estimated total run time in seconds = 21.3941002 + Estimated total run time in seconds = 21.3941002 Estimated total run time in hh:mm:ss = 0 h 00 m 21 s We have done 1.00000000 % of that The run will finish approximately on (in local time): Mon Apr 26, 2021 10:14 @@ -562,75 +562,75 @@ Time step # 100 Time: 0.551999986 seconds - Elapsed time in seconds = 4.8205130000000000 + Elapsed time in seconds = 4.8205130000000000 Elapsed time in hh:mm:ss = 0 h 00 m 04 s Mean elapsed time per time step in seconds = 4.82051298E-02 Max norm displacement vector U in all slices (m) = 3.71283425E-15 Time steps done = 100 out of 500 Time steps remaining = 400 - Estimated remaining time in seconds = 19.2820511 + Estimated remaining time in seconds = 19.2820511 Estimated remaining time in hh:mm:ss = 0 h 00 m 19 s - Estimated total run time in seconds = 24.1025658 + Estimated total run time in seconds = 24.1025658 Estimated total run time in hh:mm:ss = 0 h 00 m 24 s We have done 20.0000000 % of that The run will finish approximately on (in local time): Mon Apr 26, 2021 10:14 Time step # 200 Time: 1.35200000 seconds - Elapsed time in seconds = 10.022310999999998 + Elapsed time in seconds = 10.022310999999998 Elapsed time in hh:mm:ss = 0 h 00 m 10 s Mean elapsed time per time step in seconds = 5.01115546E-02 Max norm displacement vector U in all slices (m) = 2.08720205E-15 Time steps done = 200 out of 500 Time steps remaining = 300 - Estimated remaining time in seconds = 15.0334663 + Estimated remaining time in seconds = 15.0334663 Estimated remaining time in hh:mm:ss = 0 h 00 m 15 s - Estimated total run time in seconds = 25.0557766 + Estimated total run time in seconds = 25.0557766 Estimated total run time in hh:mm:ss = 0 h 00 m 25 s We have done 40.0000000 % of that The run will finish approximately on (in local time): Mon Apr 26, 2021 10:14 Time step # 300 Time: 2.15199995 seconds - Elapsed time in seconds = 15.335498999999999 + Elapsed time in seconds = 15.335498999999999 Elapsed time in hh:mm:ss = 0 h 00 m 15 s Mean elapsed time per time step in seconds = 5.11183292E-02 Max norm displacement vector U in all slices (m) = 1.67042986E-15 Time steps done = 300 out of 500 Time steps remaining = 200 - Estimated remaining time in seconds = 10.2236662 + Estimated remaining time in seconds = 10.2236662 Estimated remaining time in hh:mm:ss = 0 h 00 m 10 s - Estimated total run time in seconds = 25.5591640 + Estimated total run time in seconds = 25.5591640 Estimated total run time in hh:mm:ss = 0 h 00 m 25 s We have done 60.0000000 % of that The run will finish approximately on (in local time): Mon Apr 26, 2021 10:14 Time step # 400 Time: 2.95199990 seconds - Elapsed time in seconds = 20.803642000000000 + Elapsed time in seconds = 20.803642000000000 Elapsed time in hh:mm:ss = 0 h 00 m 20 s Mean elapsed time per time step in seconds = 5.20091057E-02 Max norm displacement vector U in all slices (m) = 1.37444267E-15 Time steps done = 400 out of 500 Time steps remaining = 100 - Estimated remaining time in seconds = 5.20091057 + Estimated remaining time in seconds = 5.20091057 Estimated remaining time in hh:mm:ss = 0 h 00 m 05 s - Estimated total run time in seconds = 26.0045528 + Estimated total run time in seconds = 26.0045528 Estimated total run time in hh:mm:ss = 0 h 00 m 26 s We have done 80.0000000 % of that The run will finish approximately on (in local time): Mon Apr 26, 2021 10:14 Time step # 500 Time: 3.75200009 seconds - Elapsed time in seconds = 26.348107000000002 + Elapsed time in seconds = 26.348107000000002 Elapsed time in hh:mm:ss = 0 h 00 m 26 s Mean elapsed time per time step in seconds = 5.26962131E-02 Max norm displacement vector U in all slices (m) = 1.19551480E-15 Time steps done = 500 out of 500 Time steps remaining = 0 - Estimated remaining time in seconds = 0.00000000 + Estimated remaining time in seconds = 0.00000000 Estimated remaining time in hh:mm:ss = 0 h 00 m 00 s - Estimated total run time in seconds = 26.3481064 + Estimated total run time in seconds = 26.3481064 Estimated total run time in hh:mm:ss = 0 h 00 m 26 s We have done 100.000000 % of that @@ -639,7 +639,7 @@ Writing the seismograms in parallel took 7.59299994E-02 seconds Time loop finished. Timing info: - Total elapsed time in seconds = 26.479069000000003 + Total elapsed time in seconds = 26.479069000000003 Total elapsed time in hh:mm:ss = 0 h 00 m 26 s End of the simulation diff --git a/src/shared/utm_geo.f90 b/src/shared/utm_geo.f90 index b59a7867a..979fd2b63 100644 --- a/src/shared/utm_geo.f90 +++ b/src/shared/utm_geo.f90 @@ -55,12 +55,12 @@ ! mistakes size in original routine: ! from UTM x/y to lon/lat and converted back again: ! 477415.5 5712313.5 -> 2.6741959317615298 51.561449486527074 (zone 31) -! 477415.50000327773 5712320.9415952200 <- 2.6741959317615298 51.561449486527074 +! 477415.50000327773 5712320.9415952200 <- 2.6741959317615298 51.561449486527074 ! has maximum UTM error ~ 7.44 m ! ! corrected version: ! 477415.5 5712313.5 -> 2.6741959317615298 51.561449479910003 (zone 31) -! 477415.49999999994 5712313.5002520066 <- 2.6741959317615298 51.561449479910003 +! 477415.49999999994 5712313.5002520066 <- 2.6741959317615298 51.561449479910003 ! has maximum UTM error ~ 0.00025 m ! ! error introduced in converting lon/lat to UTM: diff --git a/src/specfem3D/write_output_ASCII_or_binary.f90 b/src/specfem3D/write_output_ASCII_or_binary.f90 index 037f64b59..5744685fc 100644 --- a/src/specfem3D/write_output_ASCII_or_binary.f90 +++ b/src/specfem3D/write_output_ASCII_or_binary.f90 @@ -121,7 +121,7 @@ subroutine write_output_ASCII_or_binary(one_seismogram, & time_t_db = dble(NSTEP - it_current * subsamp_seismos) * DT - t0 endif - ! converts time to custom_real for output + ! converts time to CUSTOM_REAL for output time_t = real(time_t_db,kind=CUSTOM_REAL) if (USE_BINARY_FOR_SEISMOGRAMS) then diff --git a/src/specfem3D/write_seismograms.f90 b/src/specfem3D/write_seismograms.f90 index 77909fdb1..c2f0cfcc6 100644 --- a/src/specfem3D/write_seismograms.f90 +++ b/src/specfem3D/write_seismograms.f90 @@ -894,8 +894,8 @@ subroutine write_adj_seismograms2_to_file() ! SPECFEM reference frame convention: x = East, y = North, z = Up ! ! UTM convention: uses convention for N/E/Z - ! ‘longitude’ parameters simply refer to the x axis - ! ‘latitude’ parameters simply refer to the y axis + ! 'longitude' parameters simply refer to the x axis + ! 'latitude' parameters simply refer to the y axis ! ! SPECFEM3D Cartesian uses the traditional orientation code E/N/Z (East-West, North-South, Vertical) ! for three components when a UTM projection is used. @@ -1019,8 +1019,8 @@ subroutine write_channel_name(iorientation,channel) ! SPECFEM reference frame convention: x = East, y = North, z = Up ! ! UTM convention: uses convention for N/E/Z - ! ‘longitude’ parameters simply refer to the x axis - ! ‘latitude’ parameters simply refer to the y axis + ! 'longitude' parameters simply refer to the x axis + ! 'latitude' parameters simply refer to the y axis ! ! SPECFEM3D Cartesian uses the traditional orientation code E/N/Z (East-West, North-South, Vertical) ! for three components when a UTM projection is used. diff --git a/utils/convert_utm2geo.py b/utils/convert_utm2geo.py index 44d94c0ec..fbc28cf74 100755 --- a/utils/convert_utm2geo.py +++ b/utils/convert_utm2geo.py @@ -3,7 +3,7 @@ # converts between UTM and lon/lat coordinates # # usage example: -# convert from lon/lat to UTM: (zone 31) +# convert from lon/lat to UTM: (zone 31) # > ./convert_utm2geo.py 2.6 51.5 31 2 # convert from UTM x/y to lon/lat: # > ./convert_utm2geo.py 472234.9543355125 5705505.016786354 31 1