[cig-commits] [commit] devel, master: added utils/Roland_Sylvain_gravity/output_file_of_the_code.txt and utils/Roland_Sylvain_gravity/saved_observation_grid_real_x_y_z_used_by_the_code.txt.bz2 and updated utils/Roland_Sylvain_gravity/constants.h.in (a63cd7b)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:19:11 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit a63cd7b969fe0cc9f4445f8251234d95569f8cda
Author: Dimitri Komatitsch <see_my_GitHub_web_page>
Date: Tue May 27 03:11:00 2014 +0200
added utils/Roland_Sylvain_gravity/output_file_of_the_code.txt and utils/Roland_Sylvain_gravity/saved_observation_grid_real_x_y_z_used_by_the_code.txt.bz2 and updated utils/Roland_Sylvain_gravity/constants.h.in
>---------------------------------------------------------------
a63cd7b969fe0cc9f4445f8251234d95569f8cda
utils/Roland_Sylvain_gravity/constants.h.in | 10 +
.../output_file_of_the_code.txt | 976 +++++++++++++++++++++
...vation_grid_real_x_y_z_used_by_the_code.txt.bz2 | Bin 24399082 -> 25158686 bytes
3 files changed, 986 insertions(+)
diff --git a/utils/Roland_Sylvain_gravity/constants.h.in b/utils/Roland_Sylvain_gravity/constants.h.in
index edfd73e..199c7e0 100644
--- a/utils/Roland_Sylvain_gravity/constants.h.in
+++ b/utils/Roland_Sylvain_gravity/constants.h.in
@@ -416,6 +416,16 @@
! and results for a 3D Earth with topography)
logical, parameter :: REUSE_EXISTING_OBSERVATION_SURF = .true.
+! we may want to shift the reference frame to a pre-computed center of mass
+ logical, parameter :: SHIFT_TO_THIS_CENTER_OF_MASS = .true.
+
+! position of the center of mass of the Earth for this density model and mesh,
+! converted to meters and then made non-dimensional
+ double precision, parameter :: x_shift = 0.606220633681674d0 * 1000.d0 / R_EARTH ! km
+ double precision, parameter :: y_shift = 0.433103991863316d0 * 1000.d0 / R_EARTH ! km
+ double precision, parameter :: z_shift = 0.520078637496872d0 * 1000.d0 / R_EARTH ! km
+! distance to center = 0.908605697566306 km
+
! altitude of the observation points in km
double precision, parameter :: altitude_of_observation_points = 255.d0
diff --git a/utils/Roland_Sylvain_gravity/output_file_of_the_code.txt b/utils/Roland_Sylvain_gravity/output_file_of_the_code.txt
new file mode 100644
index 0000000..eb2bf18
--- /dev/null
+++ b/utils/Roland_Sylvain_gravity/output_file_of_the_code.txt
@@ -0,0 +1,976 @@
+
+ ****************************
+ *** Specfem3D MPI Mesher ***
+ ****************************
+
+
+ There are 600 MPI processes
+ Processes are numbered from 0 to 599
+
+ There are 160 elements along xi in each chunk
+ There are 160 elements along eta in each chunk
+
+ There are 10 slices along xi in each chunk
+ There are 10 slices along eta in each chunk
+ There is a total of 100 slices in each chunk
+ There are 6 chunks in the global mesh
+ There is a total of 600 slices in the global mesh
+
+ NGLLX = 5
+ NGLLY = 5
+ NGLLZ = 5
+
+ Shape functions defined by NGNOD = 27 control nodes
+ Surface shape functions defined by NGNOD2D = 9 control nodes
+
+ model: s40rts
+ no oceans
+ incorporating ellipticity
+ incorporating surface topography
+ no self-gravitation
+ no rotation
+ no attenuation
+
+ incorporating 3-D lateral variations
+ no heterogeneities in the mantle
+ incorporating crustal variations
+ using one layer only in PREM crust
+ incorporating anisotropy
+ no inner-core anisotropy
+ no general mantle anisotropy
+
+ Reference radius of the Earth used is 6371.00000000000 km
+
+ Central cube is at a radius of 950.000000000000 km
+ creating global slice addressing
+
+ Spatial distribution of the slices
+ 99 89 79 69 59 49 39 29 19 9
+ 98 88 78 68 58 48 38 28 18 8
+ 97 87 77 67 57 47 37 27 17 7
+ 96 86 76 66 56 46 36 26 16 6
+ 95 85 75 65 55 45 35 25 15 5
+ 94 84 74 64 54 44 34 24 14 4
+ 93 83 73 63 53 43 33 23 13 3
+ 92 82 72 62 52 42 32 22 12 2
+ 91 81 71 61 51 41 31 21 11 1
+ 90 80 70 60 50 40 30 20 10 0
+
+ 299 289 279 269 259 249 239 229 219 209 199 189 179 169 159 149 139 129 119 109 499 489 479 469 459 449 439 429 419 409
+ 298 288 278 268 258 248 238 228 218 208 198 188 178 168 158 148 138 128 118 108 498 488 478 468 458 448 438 428 418 408
+ 297 287 277 267 257 247 237 227 217 207 197 187 177 167 157 147 137 127 117 107 497 487 477 467 457 447 437 427 417 407
+ 296 286 276 266 256 246 236 226 216 206 196 186 176 166 156 146 136 126 116 106 496 486 476 466 456 446 436 426 416 406
+ 295 285 275 265 255 245 235 225 215 205 195 185 175 165 155 145 135 125 115 105 495 485 475 465 455 445 435 425 415 405
+ 294 284 274 264 254 244 234 224 214 204 194 184 174 164 154 144 134 124 114 104 494 484 474 464 454 444 434 424 414 404
+ 293 283 273 263 253 243 233 223 213 203 193 183 173 163 153 143 133 123 113 103 493 483 473 463 453 443 433 423 413 403
+ 292 282 272 262 252 242 232 222 212 202 192 182 172 162 152 142 132 122 112 102 492 482 472 462 452 442 432 422 412 402
+ 291 281 271 261 251 241 231 221 211 201 191 181 171 161 151 141 131 121 111 101 491 481 471 461 451 441 431 421 411 401
+ 290 280 270 260 250 240 230 220 210 200 190 180 170 160 150 140 130 120 110 100 490 480 470 460 450 440 430 420 410 400
+
+ 599 589 579 569 559 549 539 529 519 509
+ 598 588 578 568 558 548 538 528 518 508
+ 597 587 577 567 557 547 537 527 517 507
+ 596 586 576 566 556 546 536 526 516 506
+ 595 585 575 565 555 545 535 525 515 505
+ 594 584 574 564 554 544 534 524 514 504
+ 593 583 573 563 553 543 533 523 513 503
+ 592 582 572 562 552 542 532 522 512 502
+ 591 581 571 561 551 541 531 521 511 501
+ 590 580 570 560 550 540 530 520 510 500
+
+ 399 389 379 369 359 349 339 329 319 309
+ 398 388 378 368 358 348 338 328 318 308
+ 397 387 377 367 357 347 337 327 317 307
+ 396 386 376 366 356 346 336 326 316 306
+ 395 385 375 365 355 345 335 325 315 305
+ 394 384 374 364 354 344 334 324 314 304
+ 393 383 373 363 353 343 333 323 313 303
+ 392 382 372 362 352 342 332 322 312 302
+ 391 381 371 361 351 341 331 321 311 301
+ 390 380 370 360 350 340 330 320 310 300
+
+
+ incorporating topography
+ topography/bathymetry: min/max = -10803 8333
+
+ incorporating crustal model: CRUST1.0
+
+
+
+
+ *******************************************
+ creating mesh in region 1
+ this region is the crust and mantle
+ *******************************************
+
+
+ first pass
+
+ ...allocating arrays
+
+ ...setting up layers
+
+ ...creating mesh elements
+ creating layer 1 out of 10
+ 10.0% current clock (NOT elapsed) time is: 02h 24min 30sec
+ creating layer 2 out of 10
+ 20.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 3 out of 10
+ 30.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 4 out of 10
+ 40.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 5 out of 10
+ 50.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 6 out of 10
+ 60.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 7 out of 10
+ 70.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 8 out of 10
+ 80.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 9 out of 10
+ 90.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+ creating layer 10 out of 10
+ 100.0% current clock (NOT elapsed) time is: 02h 24min 31sec
+
+
+ ...creating global addressing
+
+ ...creating MPI buffers
+
+ second pass
+
+ ...allocating arrays
+
+ ...setting up layers
+
+ ...creating mesh elements
+ creating layer 1 out of 10
+ 10.0% current clock (NOT elapsed) time is: 02h 24min 39sec
+ creating layer 2 out of 10
+ 20.0% current clock (NOT elapsed) time is: 02h 24min 42sec
+ creating layer 3 out of 10
+ 30.0% current clock (NOT elapsed) time is: 02h 24min 42sec
+ creating layer 4 out of 10
+ 40.0% current clock (NOT elapsed) time is: 02h 24min 46sec
+ creating layer 5 out of 10
+ 50.0% current clock (NOT elapsed) time is: 02h 24min 54sec
+ creating layer 6 out of 10
+ 60.0% current clock (NOT elapsed) time is: 02h 24min 55sec
+ creating layer 7 out of 10
+ 70.0% current clock (NOT elapsed) time is: 02h 24min 56sec
+ creating layer 8 out of 10
+ 80.0% current clock (NOT elapsed) time is: 02h 24min 59sec
+ creating layer 9 out of 10
+ 90.0% current clock (NOT elapsed) time is: 02h 25min 01sec
+ creating layer 10 out of 10
+ 100.0% current clock (NOT elapsed) time is: 02h 25min 11sec
+
+
+ ...precomputing Jacobian
+
+ ...creating chunk buffers
+
+ ----- creating chunk buffers -----
+
+ There are 10 slices along xi in each chunk
+ There are 10 slices along eta in each chunk
+ There is a total of 100 slices in each chunk
+ There are 6 chunks
+ There is a total of 600 slices in all the chunks
+
+ There is a total of 120 messages to assemble faces between chunks
+
+ Generating message 1 for faces out of 120
+ Generating message 2 for faces out of 120
+ Generating message 3 for faces out of 120
+ Generating message 4 for faces out of 120
+ Generating message 5 for faces out of 120
+ Generating message 6 for faces out of 120
+ Generating message 7 for faces out of 120
+ Generating message 8 for faces out of 120
+ Generating message 9 for faces out of 120
+ Generating message 10 for faces out of 120
+ Generating message 11 for faces out of 120
+ Generating message 12 for faces out of 120
+ Generating message 13 for faces out of 120
+ Generating message 14 for faces out of 120
+ Generating message 15 for faces out of 120
+ Generating message 16 for faces out of 120
+ Generating message 17 for faces out of 120
+ Generating message 18 for faces out of 120
+ Generating message 19 for faces out of 120
+ Generating message 20 for faces out of 120
+ Generating message 21 for faces out of 120
+ Generating message 22 for faces out of 120
+ Generating message 23 for faces out of 120
+ Generating message 24 for faces out of 120
+ Generating message 25 for faces out of 120
+ Generating message 26 for faces out of 120
+ Generating message 27 for faces out of 120
+ Generating message 28 for faces out of 120
+ Generating message 29 for faces out of 120
+ Generating message 30 for faces out of 120
+ Generating message 31 for faces out of 120
+ Generating message 32 for faces out of 120
+ Generating message 33 for faces out of 120
+ Generating message 34 for faces out of 120
+ Generating message 35 for faces out of 120
+ Generating message 36 for faces out of 120
+ Generating message 37 for faces out of 120
+ Generating message 38 for faces out of 120
+ Generating message 39 for faces out of 120
+ Generating message 40 for faces out of 120
+ Generating message 41 for faces out of 120
+ Generating message 42 for faces out of 120
+ Generating message 43 for faces out of 120
+ Generating message 44 for faces out of 120
+ Generating message 45 for faces out of 120
+ Generating message 46 for faces out of 120
+ Generating message 47 for faces out of 120
+ Generating message 48 for faces out of 120
+ Generating message 49 for faces out of 120
+ Generating message 50 for faces out of 120
+ Generating message 51 for faces out of 120
+ Generating message 52 for faces out of 120
+ Generating message 53 for faces out of 120
+ Generating message 54 for faces out of 120
+ Generating message 55 for faces out of 120
+ Generating message 56 for faces out of 120
+ Generating message 57 for faces out of 120
+ Generating message 58 for faces out of 120
+ Generating message 59 for faces out of 120
+ Generating message 60 for faces out of 120
+ Generating message 61 for faces out of 120
+ Generating message 62 for faces out of 120
+ Generating message 63 for faces out of 120
+ Generating message 64 for faces out of 120
+ Generating message 65 for faces out of 120
+ Generating message 66 for faces out of 120
+ Generating message 67 for faces out of 120
+ Generating message 68 for faces out of 120
+ Generating message 69 for faces out of 120
+ Generating message 70 for faces out of 120
+ Generating message 71 for faces out of 120
+ Generating message 72 for faces out of 120
+ Generating message 73 for faces out of 120
+ Generating message 74 for faces out of 120
+ Generating message 75 for faces out of 120
+ Generating message 76 for faces out of 120
+ Generating message 77 for faces out of 120
+ Generating message 78 for faces out of 120
+ Generating message 79 for faces out of 120
+ Generating message 80 for faces out of 120
+ Generating message 81 for faces out of 120
+ Generating message 82 for faces out of 120
+ Generating message 83 for faces out of 120
+ Generating message 84 for faces out of 120
+ Generating message 85 for faces out of 120
+ Generating message 86 for faces out of 120
+ Generating message 87 for faces out of 120
+ Generating message 88 for faces out of 120
+ Generating message 89 for faces out of 120
+ Generating message 90 for faces out of 120
+ Generating message 91 for faces out of 120
+ Generating message 92 for faces out of 120
+ Generating message 93 for faces out of 120
+ Generating message 94 for faces out of 120
+ Generating message 95 for faces out of 120
+ Generating message 96 for faces out of 120
+ Generating message 97 for faces out of 120
+ Generating message 98 for faces out of 120
+ Generating message 99 for faces out of 120
+ Generating message 100 for faces out of 120
+ Generating message 101 for faces out of 120
+ Generating message 102 for faces out of 120
+ Generating message 103 for faces out of 120
+ Generating message 104 for faces out of 120
+ Generating message 105 for faces out of 120
+ Generating message 106 for faces out of 120
+ Generating message 107 for faces out of 120
+ Generating message 108 for faces out of 120
+ Generating message 109 for faces out of 120
+ Generating message 110 for faces out of 120
+ Generating message 111 for faces out of 120
+ Generating message 112 for faces out of 120
+ Generating message 113 for faces out of 120
+ Generating message 114 for faces out of 120
+ Generating message 115 for faces out of 120
+ Generating message 116 for faces out of 120
+ Generating message 117 for faces out of 120
+ Generating message 118 for faces out of 120
+ Generating message 119 for faces out of 120
+ Generating message 120 for faces out of 120
+
+ all the messages for chunk faces have the right size
+
+ Generating message 1 for corners out of 8
+ Generating message 2 for corners out of 8
+ Generating message 3 for corners out of 8
+ Generating message 4 for corners out of 8
+ Generating message 5 for corners out of 8
+ Generating message 6 for corners out of 8
+ Generating message 7 for corners out of 8
+ Generating message 8 for corners out of 8
+
+ ...preparing MPI interfaces
+
+ crust/mantle region:
+ #max of points in MPI buffers along xi npoin2D_xi = 3669
+ #max of array elements transferred npoin2D_xi*NDIM = 11007
+
+ #max of points in MPI buffers along eta npoin2D_eta = 3669
+ #max of array elements transferred npoin2D_eta*NDIM = 11007
+
+ crust mantle MPI:
+ maximum interfaces: 8
+ MPI addressing maximum interfaces: 8
+ MPI addressing : all interfaces okay
+
+ total MPI interface points : 9064584
+ unique MPI interface points: 8544000
+ maximum valence : 3
+ total unique MPI interface points: 8544000
+
+
+ ...element inner/outer separation
+
+ for overlapping of communications with calculations:
+
+ percentage of edge elements in crust/mantle 38.33333 %
+ percentage of volume elements in crust/mantle 61.66667 %
+
+
+ ...element mesh coloring
+ mesh coloring: F
+
+ ...creating mass matrix
+
+ calculated top area: 12.5600450357667
+ calculated bottom area: 3.74933399372606
+
+ ...computing Roland_Sylvain integrals
+ for Roland_Sylvain integrals 1 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 100 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 200 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 300 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 400 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 500 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 600 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 700 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 800 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 900 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1000 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1100 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1200 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1300 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1400 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1500 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1600 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1700 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1800 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 1900 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 2000 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 2100 elements computed out of
+ 2160
+ for Roland_Sylvain integrals 2160 elements computed out of
+ 2160
+
+ *******************************************
+ creating mesh in region 2
+ this region is the outer core
+ *******************************************
+
+
+ first pass
+
+ ...allocating arrays
+
+ ...setting up layers
+
+ ...creating mesh elements
+ creating layer 1 out of 2
+ 50.0% current clock (NOT elapsed) time is: 06h 39min 17sec
+ creating layer 2 out of 2
+ 100.0% current clock (NOT elapsed) time is: 06h 39min 17sec
+
+
+ ...creating global addressing
+
+ ...creating MPI buffers
+
+ second pass
+
+ ...allocating arrays
+
+ ...setting up layers
+
+ ...creating mesh elements
+ creating layer 1 out of 2
+ 50.0% current clock (NOT elapsed) time is: 06h 39min 17sec
+ creating layer 2 out of 2
+ 100.0% current clock (NOT elapsed) time is: 06h 39min 17sec
+
+
+ ...precomputing Jacobian
+
+ ...creating chunk buffers
+
+ ----- creating chunk buffers -----
+
+ There are 10 slices along xi in each chunk
+ There are 10 slices along eta in each chunk
+ There is a total of 100 slices in each chunk
+ There are 6 chunks
+ There is a total of 600 slices in all the chunks
+
+ There is a total of 120 messages to assemble faces between chunks
+
+ Generating message 1 for faces out of 120
+ Generating message 2 for faces out of 120
+ Generating message 3 for faces out of 120
+ Generating message 4 for faces out of 120
+ Generating message 5 for faces out of 120
+ Generating message 6 for faces out of 120
+ Generating message 7 for faces out of 120
+ Generating message 8 for faces out of 120
+ Generating message 9 for faces out of 120
+ Generating message 10 for faces out of 120
+ Generating message 11 for faces out of 120
+ Generating message 12 for faces out of 120
+ Generating message 13 for faces out of 120
+ Generating message 14 for faces out of 120
+ Generating message 15 for faces out of 120
+ Generating message 16 for faces out of 120
+ Generating message 17 for faces out of 120
+ Generating message 18 for faces out of 120
+ Generating message 19 for faces out of 120
+ Generating message 20 for faces out of 120
+ Generating message 21 for faces out of 120
+ Generating message 22 for faces out of 120
+ Generating message 23 for faces out of 120
+ Generating message 24 for faces out of 120
+ Generating message 25 for faces out of 120
+ Generating message 26 for faces out of 120
+ Generating message 27 for faces out of 120
+ Generating message 28 for faces out of 120
+ Generating message 29 for faces out of 120
+ Generating message 30 for faces out of 120
+ Generating message 31 for faces out of 120
+ Generating message 32 for faces out of 120
+ Generating message 33 for faces out of 120
+ Generating message 34 for faces out of 120
+ Generating message 35 for faces out of 120
+ Generating message 36 for faces out of 120
+ Generating message 37 for faces out of 120
+ Generating message 38 for faces out of 120
+ Generating message 39 for faces out of 120
+ Generating message 40 for faces out of 120
+ Generating message 41 for faces out of 120
+ Generating message 42 for faces out of 120
+ Generating message 43 for faces out of 120
+ Generating message 44 for faces out of 120
+ Generating message 45 for faces out of 120
+ Generating message 46 for faces out of 120
+ Generating message 47 for faces out of 120
+ Generating message 48 for faces out of 120
+ Generating message 49 for faces out of 120
+ Generating message 50 for faces out of 120
+ Generating message 51 for faces out of 120
+ Generating message 52 for faces out of 120
+ Generating message 53 for faces out of 120
+ Generating message 54 for faces out of 120
+ Generating message 55 for faces out of 120
+ Generating message 56 for faces out of 120
+ Generating message 57 for faces out of 120
+ Generating message 58 for faces out of 120
+ Generating message 59 for faces out of 120
+ Generating message 60 for faces out of 120
+ Generating message 61 for faces out of 120
+ Generating message 62 for faces out of 120
+ Generating message 63 for faces out of 120
+ Generating message 64 for faces out of 120
+ Generating message 65 for faces out of 120
+ Generating message 66 for faces out of 120
+ Generating message 67 for faces out of 120
+ Generating message 68 for faces out of 120
+ Generating message 69 for faces out of 120
+ Generating message 70 for faces out of 120
+ Generating message 71 for faces out of 120
+ Generating message 72 for faces out of 120
+ Generating message 73 for faces out of 120
+ Generating message 74 for faces out of 120
+ Generating message 75 for faces out of 120
+ Generating message 76 for faces out of 120
+ Generating message 77 for faces out of 120
+ Generating message 78 for faces out of 120
+ Generating message 79 for faces out of 120
+ Generating message 80 for faces out of 120
+ Generating message 81 for faces out of 120
+ Generating message 82 for faces out of 120
+ Generating message 83 for faces out of 120
+ Generating message 84 for faces out of 120
+ Generating message 85 for faces out of 120
+ Generating message 86 for faces out of 120
+ Generating message 87 for faces out of 120
+ Generating message 88 for faces out of 120
+ Generating message 89 for faces out of 120
+ Generating message 90 for faces out of 120
+ Generating message 91 for faces out of 120
+ Generating message 92 for faces out of 120
+ Generating message 93 for faces out of 120
+ Generating message 94 for faces out of 120
+ Generating message 95 for faces out of 120
+ Generating message 96 for faces out of 120
+ Generating message 97 for faces out of 120
+ Generating message 98 for faces out of 120
+ Generating message 99 for faces out of 120
+ Generating message 100 for faces out of 120
+ Generating message 101 for faces out of 120
+ Generating message 102 for faces out of 120
+ Generating message 103 for faces out of 120
+ Generating message 104 for faces out of 120
+ Generating message 105 for faces out of 120
+ Generating message 106 for faces out of 120
+ Generating message 107 for faces out of 120
+ Generating message 108 for faces out of 120
+ Generating message 109 for faces out of 120
+ Generating message 110 for faces out of 120
+ Generating message 111 for faces out of 120
+ Generating message 112 for faces out of 120
+ Generating message 113 for faces out of 120
+ Generating message 114 for faces out of 120
+ Generating message 115 for faces out of 120
+ Generating message 116 for faces out of 120
+ Generating message 117 for faces out of 120
+ Generating message 118 for faces out of 120
+ Generating message 119 for faces out of 120
+ Generating message 120 for faces out of 120
+
+ all the messages for chunk faces have the right size
+
+ Generating message 1 for corners out of 8
+ Generating message 2 for corners out of 8
+ Generating message 3 for corners out of 8
+ Generating message 4 for corners out of 8
+ Generating message 5 for corners out of 8
+ Generating message 6 for corners out of 8
+ Generating message 7 for corners out of 8
+ Generating message 8 for corners out of 8
+
+ ...preparing MPI interfaces
+
+ outer core region:
+ #max of points in MPI buffers along xi npoin2D_xi = 877
+ #max of array elements transferred npoin2D_xi*NDIM = 2631
+
+ #max of points in MPI buffers along eta npoin2D_eta = 877
+ #max of array elements transferred npoin2D_eta*NDIM = 2631
+
+ outer core MPI:
+ maximum interfaces: 8
+ MPI addressing maximum interfaces: 8
+ MPI addressing : all interfaces okay
+
+ total MPI interface points : 2259240
+ unique MPI interface points: 1948800
+ maximum valence : 3
+ total assembled MPI interface points: 1948800
+
+
+ ...element inner/outer separation
+ percentage of edge elements in outer core 81.39535 %
+ percentage of volume elements in outer core 18.60465 %
+
+
+ ...element mesh coloring
+ mesh coloring: F
+
+ ...creating mass matrix
+
+ calculated top area: 3.74933399372606
+ calculated bottom area: 0.461558583668424
+
+ ...computing Roland_Sylvain integrals
+ for Roland_Sylvain integrals 1 elements computed out of
+ 172
+ for Roland_Sylvain integrals 100 elements computed out of
+ 172
+ for Roland_Sylvain integrals 172 elements computed out of
+ 172
+
+ *******************************************
+ creating mesh in region 3
+ this region is the inner core
+ *******************************************
+
+
+ first pass
+
+ ...allocating arrays
+
+ ...setting up layers
+
+ ...creating mesh elements
+ creating layer 1 out of 1
+ 100.0% current clock (NOT elapsed) time is: 06h 59min 31sec
+
+ creating central cube
+
+ ...creating global addressing
+
+ ...creating MPI buffers
+
+ second pass
+
+ ...allocating arrays
+
+ ...setting up layers
+
+ ...creating mesh elements
+ creating layer 1 out of 1
+ 100.0% current clock (NOT elapsed) time is: 06h 59min 31sec
+
+ creating central cube
+
+ ...precomputing Jacobian
+
+ ...creating chunk buffers
+
+ ----- creating chunk buffers -----
+
+ There are 10 slices along xi in each chunk
+ There are 10 slices along eta in each chunk
+ There is a total of 100 slices in each chunk
+ There are 6 chunks
+ There is a total of 600 slices in all the chunks
+
+ There is a total of 120 messages to assemble faces between chunks
+
+ Generating message 1 for faces out of 120
+ Generating message 2 for faces out of 120
+ Generating message 3 for faces out of 120
+ Generating message 4 for faces out of 120
+ Generating message 5 for faces out of 120
+ Generating message 6 for faces out of 120
+ Generating message 7 for faces out of 120
+ Generating message 8 for faces out of 120
+ Generating message 9 for faces out of 120
+ Generating message 10 for faces out of 120
+ Generating message 11 for faces out of 120
+ Generating message 12 for faces out of 120
+ Generating message 13 for faces out of 120
+ Generating message 14 for faces out of 120
+ Generating message 15 for faces out of 120
+ Generating message 16 for faces out of 120
+ Generating message 17 for faces out of 120
+ Generating message 18 for faces out of 120
+ Generating message 19 for faces out of 120
+ Generating message 20 for faces out of 120
+ Generating message 21 for faces out of 120
+ Generating message 22 for faces out of 120
+ Generating message 23 for faces out of 120
+ Generating message 24 for faces out of 120
+ Generating message 25 for faces out of 120
+ Generating message 26 for faces out of 120
+ Generating message 27 for faces out of 120
+ Generating message 28 for faces out of 120
+ Generating message 29 for faces out of 120
+ Generating message 30 for faces out of 120
+ Generating message 31 for faces out of 120
+ Generating message 32 for faces out of 120
+ Generating message 33 for faces out of 120
+ Generating message 34 for faces out of 120
+ Generating message 35 for faces out of 120
+ Generating message 36 for faces out of 120
+ Generating message 37 for faces out of 120
+ Generating message 38 for faces out of 120
+ Generating message 39 for faces out of 120
+ Generating message 40 for faces out of 120
+ Generating message 41 for faces out of 120
+ Generating message 42 for faces out of 120
+ Generating message 43 for faces out of 120
+ Generating message 44 for faces out of 120
+ Generating message 45 for faces out of 120
+ Generating message 46 for faces out of 120
+ Generating message 47 for faces out of 120
+ Generating message 48 for faces out of 120
+ Generating message 49 for faces out of 120
+ Generating message 50 for faces out of 120
+ Generating message 51 for faces out of 120
+ Generating message 52 for faces out of 120
+ Generating message 53 for faces out of 120
+ Generating message 54 for faces out of 120
+ Generating message 55 for faces out of 120
+ Generating message 56 for faces out of 120
+ Generating message 57 for faces out of 120
+ Generating message 58 for faces out of 120
+ Generating message 59 for faces out of 120
+ Generating message 60 for faces out of 120
+ Generating message 61 for faces out of 120
+ Generating message 62 for faces out of 120
+ Generating message 63 for faces out of 120
+ Generating message 64 for faces out of 120
+ Generating message 65 for faces out of 120
+ Generating message 66 for faces out of 120
+ Generating message 67 for faces out of 120
+ Generating message 68 for faces out of 120
+ Generating message 69 for faces out of 120
+ Generating message 70 for faces out of 120
+ Generating message 71 for faces out of 120
+ Generating message 72 for faces out of 120
+ Generating message 73 for faces out of 120
+ Generating message 74 for faces out of 120
+ Generating message 75 for faces out of 120
+ Generating message 76 for faces out of 120
+ Generating message 77 for faces out of 120
+ Generating message 78 for faces out of 120
+ Generating message 79 for faces out of 120
+ Generating message 80 for faces out of 120
+ Generating message 81 for faces out of 120
+ Generating message 82 for faces out of 120
+ Generating message 83 for faces out of 120
+ Generating message 84 for faces out of 120
+ Generating message 85 for faces out of 120
+ Generating message 86 for faces out of 120
+ Generating message 87 for faces out of 120
+ Generating message 88 for faces out of 120
+ Generating message 89 for faces out of 120
+ Generating message 90 for faces out of 120
+ Generating message 91 for faces out of 120
+ Generating message 92 for faces out of 120
+ Generating message 93 for faces out of 120
+ Generating message 94 for faces out of 120
+ Generating message 95 for faces out of 120
+ Generating message 96 for faces out of 120
+ Generating message 97 for faces out of 120
+ Generating message 98 for faces out of 120
+ Generating message 99 for faces out of 120
+ Generating message 100 for faces out of 120
+ Generating message 101 for faces out of 120
+ Generating message 102 for faces out of 120
+ Generating message 103 for faces out of 120
+ Generating message 104 for faces out of 120
+ Generating message 105 for faces out of 120
+ Generating message 106 for faces out of 120
+ Generating message 107 for faces out of 120
+ Generating message 108 for faces out of 120
+ Generating message 109 for faces out of 120
+ Generating message 110 for faces out of 120
+ Generating message 111 for faces out of 120
+ Generating message 112 for faces out of 120
+ Generating message 113 for faces out of 120
+ Generating message 114 for faces out of 120
+ Generating message 115 for faces out of 120
+ Generating message 116 for faces out of 120
+ Generating message 117 for faces out of 120
+ Generating message 118 for faces out of 120
+ Generating message 119 for faces out of 120
+ Generating message 120 for faces out of 120
+
+ all the messages for chunk faces have the right size
+
+ Generating message 1 for corners out of 8
+ Generating message 2 for corners out of 8
+ Generating message 3 for corners out of 8
+ Generating message 4 for corners out of 8
+ Generating message 5 for corners out of 8
+ Generating message 6 for corners out of 8
+ Generating message 7 for corners out of 8
+ Generating message 8 for corners out of 8
+
+ ...preparing MPI interfaces
+
+ inner core region:
+ #max of points in MPI buffers along xi npoin2D_xi = 801
+ #max of array elements transferred npoin2D_xi*NDIM = 2403
+
+ #max of points in MPI buffers along eta npoin2D_eta = 801
+ #max of array elements transferred npoin2D_eta*NDIM = 2403
+
+
+ including central cube
+
+ inner core MPI:
+ inner core with central cube MPI:
+ maximum interfaces: 31
+ MPI addressing maximum interfaces: 31
+ MPI addressing : all interfaces okay
+
+ total MPI interface points : 601992
+ unique MPI interface points: 458200
+ maximum valence : 7
+ total assembled MPI interface points: 458200
+
+
+ ...element inner/outer separation
+ percentage of edge elements in inner core 54.54546 %
+ percentage of volume elements in inner core 45.45454 %
+
+
+ ...element mesh coloring
+ mesh coloring: F
+
+ ...creating mass matrix
+
+ calculated top area: 0.461558583668424
+ calculated bottom area: 0.211696672374717
+
+ ...computing Roland_Sylvain integrals
+ for Roland_Sylvain integrals 1 elements computed out of
+ 88
+ for Roland_Sylvain integrals 88 elements computed out of
+ 88
+
+ calculated volume: 4.18410575573324
+
+ computed total Earth mass for this density model and mesh:
+ 5.973928520702748E+024 kg
+ (should be not too far from 5.974E+24 kg)
+
+ average density for this density model and mesh: 5521.21365725308 kg/m
+ 3
+ (should be not too far from 5514 kg/m3)
+
+ position of the center of mass of the Earth for this density model and mesh:
+ x = -4.252974378716729E-013 km
+ y = 1.437719760222687E-012 km
+ z = 2.087146141064491E-013 km
+ distance to center = 1.513762798372517E-012 km
+
+
+ minval of norm of g vector on whole observation surface =
+ 9.07016969348020 m.s-2
+ maxval of norm of g vector on whole observation surface =
+ 9.09813353875059 m.s-2
+
+ minval of G_xx on whole observation surface = -1374.31671074432 Eotvos
+ maxval of G_xx on whole observation surface = 2738.87372493946 Eotvos
+
+ minval of G_yy on whole observation surface = -1373.67994555197 Eotvos
+ maxval of G_yy on whole observation surface = 2742.06663652309 Eotvos
+
+ minval of G_zz on whole observation surface = -1374.48230626873 Eotvos
+ maxval of G_zz on whole observation surface = 2746.20413470664 Eotvos
+
+ minval of G_xy on whole observation surface = -2056.96533920110 Eotvos
+ maxval of G_xy on whole observation surface = 2053.69297111018 Eotvos
+
+ minval of G_xz on whole observation surface = -2056.13544781585 Eotvos
+ maxval of G_xz on whole observation surface = 2057.99063009606 Eotvos
+
+ minval of G_yz on whole observation surface = -2058.48841667374 Eotvos
+ maxval of G_yz on whole observation surface = 2056.65339660025 Eotvos
+
+ Minval and maxval of trace of G, which in principle should be zero:
+
+ minval of abs(G_xx + G_yy + G_zz) on whole observation surface =
+ 0.000000000000000E+000 Eotvos
+ maxval of abs(G_xx + G_yy + G_zz) on whole observation surface =
+ 7.725020623183809E-011 Eotvos
+
+ -----------------------------
+
+ displaying the fields computed at:
+ ix_observation = 166 out of 500
+ iy_observation = 125 out of 500
+ of mesh chunk 3
+
+ computed g_x = 8.11917857575563 m.s-2
+ computed g_y = 3.38656178001491 m.s-2
+ computed g_z = 2.21918985060483 m.s-2
+
+ computed norm of g vector = 9.07274297155466 m.s-2
+ (should be not too far from 9.049929 m.s-2)
+
+ computed G_xx = 1919.19497952206 Eotvos
+ computed G_yy = -794.837027501745 Eotvos
+ computed G_zz = -1124.35795202034 Eotvos
+
+ G tensor should be traceless, G_xx + G_yy + G_zz = 0.
+ Actual sum obtained = -2.410160959698260E-011
+ i.e., 1.2558187E-12 % of max(abs(G_xx),abs(G_yy),abs(G_zz))
+
+ computed G_xy = 1370.96465913758 Eotvos
+ computed G_xz = 900.313902022877 Eotvos
+ computed G_yz = 376.174174998562 Eotvos
+
+ Repartition of elements in regions:
+ ----------------------------------
+
+ total number of elements in each slice: 2420
+
+ - crust and mantle: 89.25620 %
+ - outer core: 7.107438 %
+ - inner core: 3.636364 %
+
+ for some mesh statistics, see comments in file OUTPUT_FILES/values_from_mesher.
+ h
+
+ Load balancing = 100 % by definition
+
+
+ total number of time steps in the solver will be: 1100
+
+
+ time-stepping of the solver will be: 0.142500000000000
+
+ using double precision for the calculations
+
+ smallest and largest possible floating-point numbers are:
+ 2.225073858507201E-308 1.797693134862316E+308
+
+
+ Elapsed time for mesh generation and buffer creation in seconds =
+ 16986.0214519501
+ Elapsed time for mesh generation and buffer creation in hh:mm:ss = 4 h 43 m 06 s
+
+ End of mesh generation
+
+########## ########## ########## ########## ########## ########## ########## ##########
+Execution Sum Up
+########## ########## ########## ########## ########## ########## ########## ##########
+Jobid : 1723478
+Jobname : SPECFEM128
+User : komatits
+Account : gen6351 at standard
+Limits : time = 1-00:00:00 , memory/task = 4000 Mo
+Date : submit = 25/05/2014 02:22:48 , start = 25/05/2014 02:22:48
+Execution : partition = standard , QoS = normal , Comment = (null)
+Resources : ntasks = 608 , cpus/task = 1 , ncpus = 608 , nodes = 38
+ Nodes=curie[3060-3064,3102-3115,3117-3119,5864,5866,5983,5986-5990,5992-5999] CPU_IDs=0-15 Mem=64000
+
+Memory / step
+--------------
+ Resident Size (Mo) Virtual Size (Go)
+JobID Max (Node:Task) AveTask Max (Node:Task) AveTask
+----------- ------------------------ ------- -------------------------- -------
+
+Accounting / step
+------------------
+
+JobID JobName Ntasks Ncpus Nnodes Layout Elapsed Ratio CPusage Eff State
+------------ ------------ ------ ----- ------ ------- ------- ----- ------- --- -----
+1723478 SPECFEM128 - 608 38 - 04:43:19 100 - - -
+1723478.0 xmeshfem3D 600 600 38 BBlock 04:43:11 99.9 04:40:04 98.8 COMPLETED
+########## ########## ########## ########## ########## ########## ########## ##########
diff --git a/doc/USER_MANUAL/manual_SPECFEM3D_GLOBE.pdf b/utils/Roland_Sylvain_gravity/saved_observation_grid_real_x_y_z_used_by_the_code.txt.bz2
similarity index 69%
copy from doc/USER_MANUAL/manual_SPECFEM3D_GLOBE.pdf
copy to utils/Roland_Sylvain_gravity/saved_observation_grid_real_x_y_z_used_by_the_code.txt.bz2
index 249ff7c..b21b621 100644
Binary files a/doc/USER_MANUAL/manual_SPECFEM3D_GLOBE.pdf and b/utils/Roland_Sylvain_gravity/saved_observation_grid_real_x_y_z_used_by_the_code.txt.bz2 differ
More information about the CIG-COMMITS
mailing list