Hi everyone!
I’m a new user of SW4 and I’m trying to test a grid I created in UTM coordinates.
I paste here the first lines of my ppfile to have an easier look of the issue :
EjidoLamb
250
9 4068522.12 4070522.12
9 515985.17 517985.17
10 0 40000
-99 -99 -99 -99
.FALSE.
4068522.12 515985.17 10
1 0 1060 620 1975
2 5.4 1060 620 1975
3 149.4 1480 795 2200
4 317 2080 1060 2300
5 367 2770 1500 2400
6 1500 4450 2500 2700
7 3300 4750 2800 2750
8 11000 6100 3500 2750
9 20000 6500 3750 2850
10 40000 7900 4450 3400
4068772.12 515985.17 10
1 0 1060 620 1975
2 8 1060 620 1975
3 150.4 1480 795 2200
4 336.3 2080 1060 2300
5 407 2770 1500 2400
6 1500 4450 2500 2700
7 3300 4750 2800 2750
8 11000 6100 3500 2750
9 20000 6500 3750 2850
10 40000 7900 4450 3400
The input file is the next one :
fileio path=Dalias_lamb1.out
grid x=2000 y=2000 z=40000 nz=3200
time t=15.0
pfile filename=Ejido1bien_2000x2000.ppmod vsmin=620 vpmin=1060 smoothingsize=2 style=cartesian
source type=C6SmoothBump x=1000 y=1000 depth=10000 fz=1e13 freq=1 t0=0
rec x=1500 y=1800 z=0 file=Ejido2000lamb
Running them on the sw4 I got the next error :
Input file: Dalias_Lamb2000.in
- Processing the grid command…
- Setting h to 1.25039075e+01 from z/(nz-1) (z=4.00000000e+04, nz=3200)
- Setting nx to 161 to be consistent with h=1.25039075e+01
- Setting ny to 161 to be consistent with h=1.25039075e+01
- Changing x from 2.00000000e+03 to 2.00062520e+03 to be consistent with h=1.25039075e+01
- Changing y from 2.00000000e+03 to 2.00062520e+03 to be consistent with h=1.25039075e+01
*** No topography command found in input file. Using z=0 as free surface boundary ***
Global grid sizes (without ghost points)
Grid h Nx Ny Nz Points
0 12.5039 161 161 3200 82947200
Total number of grid points (without ghost points): 8.29472e+07
*** Reading data from Pfile Ejido1bien_2000x2000.ppmod in directory ./
Pfile model name (string): ‘EjidoLamb’
Step size in x and y: 2.50000000e+02
Number of x-direction points: 9
Min x: 4.06852212e+06 Max x: 4.07052212e+06
Number of y-direction points: 9
Min y: 5.15985170e+05 Max y: 5.17985170e+05
Number of depth points: 10
Min depth: 0.00000000e+00 Max depth: 4.00000000e+04
Optional indices: Sediment: -99 MoHo: -99 410: -99 660: -99
Attenuation Q-factors available: no
ppmod: NOT allocating arrays for Qp and Qs
******* Done reading Pfile **********
Execution time, reading input file 8.47671604e+00 seconds
Making Output Directory: Dalias_lamb1.out/
… Done!
Geographic and Cartesian coordinates of the corners of the computational grid:
0: Lon= 1.180000e+02, Lat=3.700000e+01, x=0.000000e+00, y=0.000000e+00
1: Lon= 1.180000e+02, Lat=3.701797e+01, x=2.000625e+03, y=0.000000e+00
2: Lon= 1.180225e+02, Lat=3.701797e+01, x=2.000625e+03, y=2.000625e+03
3: Lon= 1.180225e+02, Lat=3.700000e+01, x=0.000000e+00, y=2.000625e+03
Considering all material blocks
Assigning material properties from pfile data…
outside = 82999042, material = 0
----------- Material properties ranges ---------------
-1.00000000e+00 kg/m^3 <= Density <= -1.00000000e+00 kg/m^3
1.00000000e+00 m/s <= Vp <= 1.00000000e+00 m/s
1.00000000e+00 m/s <= Vs <= 1.00000000e+00 m/s
1.00000000e+00 <= Vp/Vs <= 1.00000000e+00
-1.00000000e+00 Pa <= mu <= -1.00000000e+00 Pa
1.00000000e+00 Pa <= lambda <= 1.00000000e+00 Pa
------------------------------------------------------
Fatal input error: Density= -1.00000000e+00 in grid g= 0 at point (-1,-1,-1)
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 1.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
To sum up my thoughts : It seems that it cannot read the materials on the ppfile and possibly there is an error on the ppfile linked to the coordinate system used and how the software works with it. I’m struggling to understand what I’m doing wrong and came here to ask for some help. Someone here to shed me some light?
Thanks in advance,
Helena.