Author: dkomati1
Date: 2011-02-14 16:12:10 -0800 (Mon, 14 Feb 2011)
New Revision: 17863
Modified:
seismo/2D/SPECFEM2D/trunk/Database00000_left_edge_only
seismo/2D/SPECFEM2D/trunk/Database00000_right_edge_only
seismo/2D/SPECFEM2D/trunk/include_for_periodic_conditions.f90
seismo/2D/SPECFEM2D/trunk/meshfem2D.F90
seismo/2D/SPECFEM2D/trunk/part_unstruct.F90
seismo/2D/SPECFEM2D/trunk/read_interfaces_file.f90
seismo/2D/SPECFEM2D/trunk/read_materials.f90
seismo/2D/SPECFEM2D/trunk/read_parameter_file.F90
seismo/2D/SPECFEM2D/trunk/read_regions.f90
seismo/2D/SPECFEM2D/trunk/read_source_file.f90
seismo/2D/SPECFEM2D/trunk/read_value_parameters.f90
seismo/2D/SPECFEM2D/trunk/save_databases.f90
seismo/2D/SPECFEM2D/trunk/save_gnuplot_file.f90
seismo/2D/SPECFEM2D/trunk/save_stations_file.f90
seismo/2D/SPECFEM2D/trunk/specfem2D.F90
Log:
- improved the routine I added last week to detect and implement periodic boundary conditions
- removed useless white spaces
Modified: seismo/2D/SPECFEM2D/trunk/Database00000_left_edge_only
===================================================================
--- seismo/2D/SPECFEM2D/trunk/Database00000_left_edge_only 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/Database00000_left_edge_only 2011-02-15 00:12:10 UTC (rev 17863)
@@ -1,670 +1,384 @@
- 19653 F F T F 1 3 1 3 1
+ 20260 F T F F 1 3 1 3 1
3 1 3
- 19654 F F T F 1 3 1 3 1
+ 20261 F T F F 1 3 1 3 1
3 1 3
- 20145 F F T F 1 3 1 3 1
+ 20733 F F T F 1 3 1 3 1
3 1 3
- 20146 F F T F 1 3 1 3 1
+ 20734 F F T F 1 3 1 3 1
3 1 3
- 21457 F F T F 1 3 1 3 1
+ 20829 F F T F 1 3 1 3 1
3 1 3
- 21458 F F T F 1 3 1 3 1
+ 20830 F F T F 1 3 1 3 1
3 1 3
- 22655 F F F T 1 3 1 3 1
+ 23351 F F F T 1 3 1 3 1
3 1 3
- 22658 F F F T 1 3 1 3 1
+ 23354 F F F T 1 3 1 3 1
3 1 3
- 24427 F F F T 1 3 1 3 1
+ 24169 F F T F 1 3 1 3 1
3 1 3
- 24430 F F F T 1 3 1 3 1
+ 24170 F F T F 1 3 1 3 1
3 1 3
- 26311 F F F T 1 3 1 3 1
+ 29305 F F T F 1 3 1 3 1
3 1 3
- 26314 F F F T 1 3 1 3 1
+ 29306 F F T F 1 3 1 3 1
3 1 3
- 26597 F F T F 1 3 1 3 1
+ 31255 F F F T 1 3 1 3 1
3 1 3
- 26598 F F T F 1 3 1 3 1
+ 31258 F F F T 1 3 1 3 1
3 1 3
- 26968 F T F F 1 3 1 3 1
+ 31365 F F T F 1 3 1 3 1
3 1 3
- 26969 F T F F 1 3 1 3 1
+ 31366 F F T F 1 3 1 3 1
3 1 3
- 28035 T F F F 1 3 1 3 1
+ 31727 F F F T 1 3 1 3 1
3 1 3
- 28036 T F F F 1 3 1 3 1
+ 31730 F F F T 1 3 1 3 1
3 1 3
- 29247 F F F T 1 3 1 3 1
+ 33899 F F F T 1 3 1 3 1
3 1 3
- 29250 F F F T 1 3 1 3 1
+ 33902 F F F T 1 3 1 3 1
3 1 3
- 31693 F F T F 1 3 1 3 1
+ 34565 F F T F 1 3 1 3 1
3 1 3
- 31694 F F T F 1 3 1 3 1
+ 34566 F F T F 1 3 1 3 1
3 1 3
- 32111 F F F T 1 3 1 3 1
+ 35379 F F F T 1 3 1 3 1
3 1 3
- 32114 F F F T 1 3 1 3 1
+ 35382 F F F T 1 3 1 3 1
3 1 3
- 33929 F F T F 1 3 1 3 1
+ 36379 F F F T 1 3 1 3 1
3 1 3
- 33930 F F T F 1 3 1 3 1
+ 36382 F F F T 1 3 1 3 1
3 1 3
- 34001 F F T F 1 3 1 3 1
+ 36475 F F F T 1 3 1 3 1
3 1 3
- 34002 F F T F 1 3 1 3 1
+ 36478 F F F T 1 3 1 3 1
3 1 3
- 34845 F F T F 1 3 1 3 1
+ 36589 F F T F 1 3 1 3 1
3 1 3
- 34846 F F T F 1 3 1 3 1
+ 36590 F F T F 1 3 1 3 1
3 1 3
- 35148 F T F F 1 3 1 3 1
+ 37207 F F F T 1 3 1 3 1
3 1 3
- 35149 F T T F 1 3 1 3 1
+ 37210 F F F T 1 3 1 3 1
3 1 3
- 35150 F F T F 1 3 1 3 1
+ 39156 F T F F 1 3 1 3 1
3 1 3
- 35577 F F T F 1 3 1 3 1
+ 39157 F T F F 1 3 1 3 1
3 1 3
- 35578 F F T F 1 3 1 3 1
+ 39159 F F F T 1 3 1 3 1
3 1 3
- 36337 F F T F 1 3 1 3 1
+ 39162 F F F T 1 3 1 3 1
3 1 3
- 36338 F F T F 1 3 1 3 1
+ 39421 F F T F 1 3 1 3 1
3 1 3
- 36715 F F F T 1 3 1 3 1
+ 39422 F F T F 1 3 1 3 1
3 1 3
- 36718 F F F T 1 3 1 3 1
+ 40883 F F F T 1 3 1 3 1
3 1 3
- 36925 F F T F 1 3 1 3 1
+ 40886 F F F T 1 3 1 3 1
3 1 3
- 36926 F F T F 1 3 1 3 1
+ 41944 F T F F 1 3 1 3 1
3 1 3
- 37455 F F F T 1 3 1 3 1
+ 41945 F T F F 1 3 1 3 1
3 1 3
- 37458 F F F T 1 3 1 3 1
+ 43259 F F F T 1 3 1 3 1
3 1 3
- 41015 T F F F 1 3 1 3 1
+ 43262 F F F T 1 3 1 3 1
3 1 3
- 41016 T F F F 1 3 1 3 1
+ 43709 F F T F 1 3 1 3 1
3 1 3
- 41729 F F T F 1 3 1 3 1
+ 43710 F F T F 1 3 1 3 1
3 1 3
- 41730 F F T F 1 3 1 3 1
+ 44409 F F T F 1 3 1 3 1
3 1 3
- 45457 F F T F 1 3 1 3 1
+ 44410 F F T F 1 3 1 3 1
3 1 3
- 45458 F F T F 1 3 1 3 1
+ 47347 T F F F 1 3 1 3 1
3 1 3
- 46101 F F T F 1 3 1 3 1
+ 47348 T F F F 1 3 1 3 1
3 1 3
- 46102 F F T F 1 3 1 3 1
+ 48103 F F F T 1 3 1 3 1
3 1 3
- 46779 F F F T 1 3 1 3 1
+ 48106 F F F T 1 3 1 3 1
3 1 3
- 46782 F F F T 1 3 1 3 1
+ 50005 F F T F 1 3 1 3 1
3 1 3
- 48083 F F F T 1 3 1 3 1
+ 50006 F F T F 1 3 1 3 1
3 1 3
- 48086 F F F T 1 3 1 3 1
+ 50009 F F T F 1 3 1 3 1
3 1 3
- 48265 F F T F 1 3 1 3 1
+ 50010 F F T F 1 3 1 3 1
3 1 3
- 48266 F F T F 1 3 1 3 1
+ 50047 F F F T 1 3 1 3 1
3 1 3
- 50063 T F F F 1 3 1 3 1
+ 50050 F F F T 1 3 1 3 1
3 1 3
- 50064 T F F F 1 3 1 3 1
+ 52369 F F T F 1 3 1 3 1
3 1 3
- 50069 F F T F 1 3 1 3 1
+ 52370 F F T F 1 3 1 3 1
3 1 3
- 50070 F F T F 1 3 1 3 1
+ 53259 F F F T 1 3 1 3 1
3 1 3
- 50813 F F T F 1 3 1 3 1
+ 53262 F F F T 1 3 1 3 1
3 1 3
- 50814 F F T F 1 3 1 3 1
+ 54285 F F T F 1 3 1 3 1
3 1 3
- 53659 F F F T 1 3 1 3 1
+ 54286 F F T F 1 3 1 3 1
3 1 3
- 53662 F F F T 1 3 1 3 1
+ 54819 T F F F 1 3 1 3 1
3 1 3
- 54501 F F T F 1 3 1 3 1
+ 54820 T F F F 1 3 1 3 1
3 1 3
- 54502 F F T F 1 3 1 3 1
+ 55361 F F T F 1 3 1 3 1
3 1 3
- 55036 F T F F 1 3 1 3 1
+ 55362 F F T F 1 3 1 3 1
3 1 3
- 55037 F T F F 1 3 1 3 1
+ 56353 F F T F 1 3 1 3 1
3 1 3
- 57207 F F F T 1 3 1 3 1
+ 56354 F F T F 1 3 1 3 1
3 1 3
- 57210 F F F T 1 3 1 3 1
+ 56641 F F T F 1 3 1 3 1
3 1 3
- 59505 F F T F 1 3 1 3 1
+ 56642 F F T F 1 3 1 3 1
3 1 3
- 59506 F F T F 1 3 1 3 1
+ 57181 F F T F 1 3 1 3 1
3 1 3
- 61392 F T F F 1 3 1 3 1
+ 57182 F F T F 1 3 1 3 1
3 1 3
- 61393 F T F F 1 3 1 3 1
+ 58703 F F F T 1 3 1 3 1
3 1 3
- 63137 F F T F 1 3 1 3 1
+ 58706 F F F T 1 3 1 3 1
3 1 3
- 63138 F F T F 1 3 1 3 1
+ 61807 F F F T 1 3 1 3 1
3 1 3
- 64447 F F F T 1 3 1 3 1
+ 61810 F F F T 1 3 1 3 1
3 1 3
- 64450 F F F T 1 3 1 3 1
+ 62261 F F T F 1 3 1 3 1
3 1 3
- 64559 T F F F 1 3 1 3 1
+ 62262 F F T F 1 3 1 3 1
3 1 3
- 64560 T F F F 1 3 1 3 1
+ 62265 F F T F 1 3 1 3 1
3 1 3
- 65097 F F T F 1 3 1 3 1
+ 62266 F F T F 1 3 1 3 1
3 1 3
- 65098 F F T F 1 3 1 3 1
+ 62701 F F T F 1 3 1 3 1
3 1 3
- 65363 F F F T 1 3 1 3 1
+ 62702 F F T F 1 3 1 3 1
3 1 3
- 65366 F F F T 1 3 1 3 1
+ 63907 F F F T 1 3 1 3 1
3 1 3
- 66587 F F F T 1 3 1 3 1
+ 63910 F F F T 1 3 1 3 1
3 1 3
- 66590 F F F T 1 3 1 3 1
+ 64217 F F T F 1 3 1 3 1
3 1 3
- 67189 F F T F 1 3 1 3 1
+ 64218 F F T F 1 3 1 3 1
3 1 3
- 67190 F F T F 1 3 1 3 1
+ 65631 F F F T 1 3 1 3 1
3 1 3
- 69643 F F F T 1 3 1 3 1
+ 65634 F F F T 1 3 1 3 1
3 1 3
- 69646 F F F T 1 3 1 3 1
+ 67413 F F T F 1 3 1 3 1
3 1 3
- 71489 F F T F 1 3 1 3 1
+ 67414 F F T F 1 3 1 3 1
3 1 3
- 71490 F F T F 1 3 1 3 1
+ 67693 F F T F 1 3 1 3 1
3 1 3
- 71493 F F T F 1 3 1 3 1
+ 67694 F F T F 1 3 1 3 1
3 1 3
- 71494 F F T F 1 3 1 3 1
+ 68825 F F T F 1 3 1 3 1
3 1 3
- 71959 F F F T 1 3 1 3 1
+ 68826 F F T F 1 3 1 3 1
3 1 3
- 71962 F F F T 1 3 1 3 1
+ 68829 F F T F 1 3 1 3 1
3 1 3
- 71963 F F F T 1 3 1 3 1
+ 68830 F F T F 1 3 1 3 1
3 1 3
- 71966 F F F T 1 3 1 3 1
+ 71244 F T F F 1 3 1 3 1
3 1 3
- 72693 F F T F 1 3 1 3 1
+ 71245 F T F F 1 3 1 3 1
3 1 3
- 72694 F F T F 1 3 1 3 1
+ 71581 F F T F 1 3 1 3 1
3 1 3
- 74163 F F F T 1 3 1 3 1
+ 71582 F F T F 1 3 1 3 1
3 1 3
- 74166 F F F T 1 3 1 3 1
+ 73317 F F T F 1 3 1 3 1
3 1 3
- 76021 F F T F 1 3 1 3 1
+ 73318 F F T F 1 3 1 3 1
3 1 3
- 76022 F F T F 1 3 1 3 1
+ 73729 F F T F 1 3 1 3 1
3 1 3
- 77155 T F F F 1 3 1 3 1
+ 73730 F F T F 1 3 1 3 1
3 1 3
- 77156 T F F F 1 3 1 3 1
+ 76349 F F T F 1 3 1 3 1
3 1 3
- 77709 F F T F 1 3 1 3 1
+ 76350 F F T F 1 3 1 3 1
3 1 3
- 77710 F F T F 1 3 1 3 1
+ 77419 F F F T 1 3 1 3 1
3 1 3
- 77879 F F F T 1 3 1 3 1
+ 77422 F F F T 1 3 1 3 1
3 1 3
- 77882 F F F T 1 3 1 3 1
+ 77739 F F F T 1 3 1 3 1
3 1 3
- 77884 F T F F 1 3 1 3 1
+ 77742 F F F T 1 3 1 3 1
3 1 3
- 77885 F T F F 1 3 1 3 1
+ 77860 F T F F 1 3 1 3 1
3 1 3
- 79405 F F T F 1 3 1 3 1
+ 77861 F T F F 1 3 1 3 1
3 1 3
- 79406 F F T F 1 3 1 3 1
+ 79429 F F T F 1 3 1 3 1
3 1 3
- 80793 F F T F 1 3 1 3 1
+ 79430 F F T F 1 3 1 3 1
3 1 3
- 80794 F F T F 1 3 1 3 1
+ 79855 F F F T 1 3 1 3 1
3 1 3
- 81689 F F T F 1 3 1 3 1
+ 79858 F F F T 1 3 1 3 1
3 1 3
- 81690 F F T F 1 3 1 3 1
+ 79891 F F F T 1 3 1 3 1
3 1 3
- 82947 T F F F 1 3 1 3 1
+ 79894 F F F T 1 3 1 3 1
3 1 3
- 82948 T F F F 1 3 1 3 1
+ 80663 F F F T 1 3 1 3 1
3 1 3
- 84955 F F F T 1 3 1 3 1
+ 80666 F F F T 1 3 1 3 1
3 1 3
- 84958 F F F T 1 3 1 3 1
+ 81227 F F F T 1 3 1 3 1
3 1 3
- 88368 F T F F 1 3 1 3 1
+ 81230 F F F T 1 3 1 3 1
3 1 3
- 88369 F T F F 1 3 1 3 1
+ 81761 F F T F 1 3 1 3 1
3 1 3
- 88857 F F T F 1 3 1 3 1
+ 81762 F F T F 1 3 1 3 1
3 1 3
- 88858 F F T F 1 3 1 3 1
+ 82675 T F F F 1 3 1 3 1
3 1 3
- 88861 F F T F 1 3 1 3 1
+ 82676 T F F F 1 3 1 3 1
3 1 3
- 88862 F F T F 1 3 1 3 1
+ 83315 F F F T 1 3 1 3 1
3 1 3
- 88941 F F T F 1 3 1 3 1
+ 83318 F F F T 1 3 1 3 1
3 1 3
- 88942 F F T F 1 3 1 3 1
+ 84127 F F F T 1 3 1 3 1
3 1 3
- 88947 F F F T 1 3 1 3 1
+ 84130 F F F T 1 3 1 3 1
3 1 3
- 88950 F F F T 1 3 1 3 1
+ 84988 F T F F 1 3 1 3 1
3 1 3
- 89521 F F T F 1 3 1 3 1
+ 84989 F T F F 1 3 1 3 1
3 1 3
- 89522 F F T F 1 3 1 3 1
+ 85069 F F T F 1 3 1 3 1
3 1 3
- 90421 F F T F 1 3 1 3 1
+ 85070 F F T F 1 3 1 3 1
3 1 3
- 90422 F F T F 1 3 1 3 1
+ 85512 T F F F 1 3 1 3 1
3 1 3
- 91627 F F F T 1 3 1 3 1
+ 85513 T F F F 1 3 1 3 1
3 1 3
- 91630 F F F T 1 3 1 3 1
+ 85537 F T F F 1 3 1 3 1
3 1 3
- 92235 F F F T 1 3 1 3 1
+ 85538 F F T F 1 3 1 3 1
3 1 3
- 92238 F F F T 1 3 1 3 1
+ 85669 F T F F 1 3 1 3 1
3 1 3
- 93705 F F T F 1 3 1 3 1
+ 85670 F F T F 1 3 1 3 1
3 1 3
- 93706 F F T F 1 3 1 3 1
+ 85836 T F F F 1 3 1 3 1
3 1 3
- 96521 F F T F 1 3 1 3 1
+ 85837 T F F F 1 3 1 3 1
3 1 3
- 96522 F F T F 1 3 1 3 1
+ 86358 T F F F 1 3 1 3 1
3 1 3
- 96925 F F T F 1 3 1 3 1
+ 86359 T F F F 1 3 1 3 1
3 1 3
- 96926 F F T F 1 3 1 3 1
+ 86536 F T F F 1 3 1 3 1
3 1 3
- 97607 F F F T 1 3 1 3 1
+ 86537 F F T F 1 3 1 3 1
3 1 3
- 97610 F F F T 1 3 1 3 1
+ 86782 F T F F 1 3 1 3 1
3 1 3
- 98645 F F T F 1 3 1 3 1
+ 86783 F F T F 1 3 1 3 1
3 1 3
- 98646 F F T F 1 3 1 3 1
+ 87535 F T F F 1 3 1 3 1
3 1 3
- 100437 F F T F 1 3 1 3 1
+ 87536 F F T F 1 3 1 3 1
3 1 3
- 100438 F F T F 1 3 1 3 1
+ 87613 F T F F 1 3 1 3 1
3 1 3
- 100603 F F F T 1 3 1 3 1
+ 87614 F F T F 1 3 1 3 1
3 1 3
- 100606 F F F T 1 3 1 3 1
+ 87762 T F F F 1 3 1 3 1
3 1 3
- 100775 F F F T 1 3 1 3 1
+ 87763 T F F F 1 3 1 3 1
3 1 3
- 100778 F F F T 1 3 1 3 1
+ 88975 F T F F 1 3 1 3 1
3 1 3
- 101465 F F T F 1 3 1 3 1
+ 88976 F F T F 1 3 1 3 1
3 1 3
- 101466 F F T F 1 3 1 3 1
+ 89838 T F F F 1 3 1 3 1
3 1 3
- 101469 F F T F 1 3 1 3 1
+ 89839 T F F F 1 3 1 3 1
3 1 3
- 101470 F F T F 1 3 1 3 1
+ 90433 F T F F 1 3 1 3 1
3 1 3
- 101845 F F T F 1 3 1 3 1
+ 90434 F F T F 1 3 1 3 1
3 1 3
- 101846 F F T F 1 3 1 3 1
+ 91852 F T F F 1 3 1 3 1
3 1 3
- 102097 F F T F 1 3 1 3 1
+ 91853 F F T F 1 3 1 3 1
3 1 3
- 102098 F F T F 1 3 1 3 1
+ 92244 T F F F 1 3 1 3 1
3 1 3
- 102811 F F F T 1 3 1 3 1
+ 92245 T F F F 1 3 1 3 1
3 1 3
- 102814 F F F T 1 3 1 3 1
+ 92676 T F F F 1 3 1 3 1
3 1 3
- 105885 F F T F 1 3 1 3 1
+ 92677 T F F F 1 3 1 3 1
3 1 3
- 105886 F F T F 1 3 1 3 1
+ 95265 T F F F 1 3 1 3 1
3 1 3
- 106449 F F T F 1 3 1 3 1
+ 95266 T F F F 1 3 1 3 1
3 1 3
- 106450 F F T F 1 3 1 3 1
+ 97152 T F F F 1 3 1 3 1
3 1 3
- 106835 F F F T 1 3 1 3 1
+ 97153 T F F F 1 3 1 3 1
3 1 3
- 106838 F F F T 1 3 1 3 1
+ 99090 T F F F 1 3 1 3 1
3 1 3
- 108441 F F T F 1 3 1 3 1
+ 99091 T F F F 1 3 1 3 1
3 1 3
- 108442 F F T F 1 3 1 3 1
+ 99222 T F F F 1 3 1 3 1
3 1 3
- 108633 F F T F 1 3 1 3 1
+ 99223 T F F F 1 3 1 3 1
3 1 3
- 108634 F F T F 1 3 1 3 1
+ 100542 F F F T 1 3 1 3 1
3 1 3
- 109133 F F T F 1 3 1 3 1
+ 100544 F F F T 1 3 1 3 1
3 1 3
- 109134 F F T F 1 3 1 3 1
+ 101289 T F F F 1 3 1 3 1
3 1 3
- 109379 F F F T 1 3 1 3 1
+ 101290 T F F F 1 3 1 3 1
3 1 3
- 109382 F F F T 1 3 1 3 1
+ 103386 T F F F 1 3 1 3 1
3 1 3
- 110491 T F F F 1 3 1 3 1
+ 103387 T F F F 1 3 1 3 1
3 1 3
- 110492 T F F F 1 3 1 3 1
+ 103398 T F F F 1 3 1 3 1
3 1 3
- 110563 F F F T 1 3 1 3 1
+ 103399 T F F F 1 3 1 3 1
3 1 3
- 110566 F F F T 1 3 1 3 1
+ 103425 T F F F 1 3 1 3 1
3 1 3
- 111156 F T F F 1 3 1 3 1
+ 103426 T F F F 1 3 1 3 1
3 1 3
- 111157 F T F F 1 3 1 3 1
+ 103536 T F F F 1 3 1 3 1
3 1 3
- 111159 F F F T 1 3 1 3 1
+ 103537 T F F F 1 3 1 3 1
3 1 3
- 111162 F F F T 1 3 1 3 1
+ 103588 F T F F 1 3 1 3 1
3 1 3
- 111197 F F T F 1 3 1 3 1
+ 103589 F F T F 1 3 1 3 1
3 1 3
- 111198 F F T F 1 3 1 3 1
+ 103671 F F F T 1 3 1 3 1
3 1 3
- 111440 F T F F 1 3 1 3 1
+ 103673 F F F T 1 3 1 3 1
3 1 3
- 111441 F T F F 1 3 1 3 1
- 3 1 3
- 111665 F F T F 1 3 1 3 1
- 3 1 3
- 111666 F F T F 1 3 1 3 1
- 3 1 3
- 113201 F F T F 1 3 1 3 1
- 3 1 3
- 113202 F F T F 1 3 1 3 1
- 3 1 3
- 113623 F F F T 1 3 1 3 1
- 3 1 3
- 113626 F F F T 1 3 1 3 1
- 3 1 3
- 113743 F F F T 1 3 1 3 1
- 3 1 3
- 113746 F F F T 1 3 1 3 1
- 3 1 3
- 114600 F T F F 1 3 1 3 1
- 3 1 3
- 114601 F T F F 1 3 1 3 1
- 3 1 3
- 115371 F F F T 1 3 1 3 1
- 3 1 3
- 115374 F F F T 1 3 1 3 1
- 3 1 3
- 115641 F F T F 1 3 1 3 1
- 3 1 3
- 115642 F F T F 1 3 1 3 1
- 3 1 3
- 115655 F F F T 1 3 1 3 1
- 3 1 3
- 115658 F F F T 1 3 1 3 1
- 3 1 3
- 117055 T F F F 1 3 1 3 1
- 3 1 3
- 117056 T F F F 1 3 1 3 1
- 3 1 3
- 118248 F T F F 1 3 1 3 1
- 3 1 3
- 118249 F T F F 1 3 1 3 1
- 3 1 3
- 118744 F T F F 1 3 1 3 1
- 3 1 3
- 118745 F T F F 1 3 1 3 1
- 3 1 3
- 119851 F F F T 1 3 1 3 1
- 3 1 3
- 119854 F F F T 1 3 1 3 1
- 3 1 3
- 120107 F F F T 1 3 1 3 1
- 3 1 3
- 120110 F F F T 1 3 1 3 1
- 3 1 3
- 120567 F F F T 1 3 1 3 1
- 3 1 3
- 120570 F F F T 1 3 1 3 1
- 3 1 3
- 120909 F F T F 1 3 1 3 1
- 3 1 3
- 120910 F F T F 1 3 1 3 1
- 3 1 3
- 120968 F T F F 1 3 1 3 1
- 3 1 3
- 120969 F T F F 1 3 1 3 1
- 3 1 3
- 121456 F T F F 1 3 1 3 1
- 3 1 3
- 121457 F T F F 1 3 1 3 1
- 3 1 3
- 121572 T F F F 1 3 1 3 1
- 3 1 3
- 121573 T F F F 1 3 1 3 1
- 3 1 3
- 121600 F T F F 1 3 1 3 1
- 3 1 3
- 121601 F F T F 1 3 1 3 1
- 3 1 3
- 121893 T F F F 1 3 1 3 1
- 3 1 3
- 121894 T F F F 1 3 1 3 1
- 3 1 3
- 122055 T F F F 1 3 1 3 1
- 3 1 3
- 122056 T F F F 1 3 1 3 1
- 3 1 3
- 122154 T F F F 1 3 1 3 1
- 3 1 3
- 122155 T F F F 1 3 1 3 1
- 3 1 3
- 122614 F T F F 1 3 1 3 1
- 3 1 3
- 122615 F F T F 1 3 1 3 1
- 3 1 3
- 123685 F T F F 1 3 1 3 1
- 3 1 3
- 123686 F F T F 1 3 1 3 1
- 3 1 3
- 123709 F T F F 1 3 1 3 1
- 3 1 3
- 123710 F F T F 1 3 1 3 1
- 3 1 3
- 123849 T F F F 1 3 1 3 1
- 3 1 3
- 123850 T F F F 1 3 1 3 1
- 3 1 3
- 125632 F T F F 1 3 1 3 1
- 3 1 3
- 125633 F F T F 1 3 1 3 1
- 3 1 3
- 125748 T F F F 1 3 1 3 1
- 3 1 3
- 125749 T F F F 1 3 1 3 1
- 3 1 3
- 126009 T F F F 1 3 1 3 1
- 3 1 3
- 126010 T F F F 1 3 1 3 1
- 3 1 3
- 126349 F T F F 1 3 1 3 1
- 3 1 3
- 126350 F F T F 1 3 1 3 1
- 3 1 3
- 127488 T F F F 1 3 1 3 1
- 3 1 3
- 127489 T F F F 1 3 1 3 1
- 3 1 3
- 127494 T F F F 1 3 1 3 1
- 3 1 3
- 127495 T F F F 1 3 1 3 1
- 3 1 3
- 127648 F T F F 1 3 1 3 1
- 3 1 3
- 127649 F F T F 1 3 1 3 1
- 3 1 3
- 128163 T F F F 1 3 1 3 1
- 3 1 3
- 128164 T F F F 1 3 1 3 1
- 3 1 3
- 129844 F T F F 1 3 1 3 1
- 3 1 3
- 129845 F F T F 1 3 1 3 1
- 3 1 3
- 129871 F T F F 1 3 1 3 1
- 3 1 3
- 129872 F F T F 1 3 1 3 1
- 3 1 3
- 129942 T F F F 1 3 1 3 1
- 3 1 3
- 129943 T F F F 1 3 1 3 1
- 3 1 3
- 130075 F T F F 1 3 1 3 1
- 3 1 3
- 130076 F F T F 1 3 1 3 1
- 3 1 3
- 130081 F T F F 1 3 1 3 1
- 3 1 3
- 130082 F F T F 1 3 1 3 1
- 3 1 3
- 130225 F T F F 1 3 1 3 1
- 3 1 3
- 130226 F F T F 1 3 1 3 1
- 3 1 3
- 130290 T F F F 1 3 1 3 1
- 3 1 3
- 130291 T F F F 1 3 1 3 1
- 3 1 3
- 130774 F T F F 1 3 1 3 1
- 3 1 3
- 130775 F F T F 1 3 1 3 1
- 3 1 3
- 130890 T F F F 1 3 1 3 1
- 3 1 3
- 130891 T F F F 1 3 1 3 1
- 3 1 3
- 130972 F T F F 1 3 1 3 1
- 3 1 3
- 130973 F F T F 1 3 1 3 1
- 3 1 3
- 131017 F T F F 1 3 1 3 1
- 3 1 3
- 131018 F F T F 1 3 1 3 1
- 3 1 3
- 131038 F T F F 1 3 1 3 1
- 3 1 3
- 131039 F F T F 1 3 1 3 1
- 3 1 3
- 131097 T F F F 1 3 1 3 1
- 3 1 3
- 131098 T F F F 1 3 1 3 1
- 3 1 3
- 131380 F T F F 1 3 1 3 1
- 3 1 3
- 131381 F F T F 1 3 1 3 1
- 3 1 3
- 133200 T F F F 1 3 1 3 1
- 3 1 3
- 133201 T F F F 1 3 1 3 1
- 3 1 3
- 133216 F T F F 1 3 1 3 1
- 3 1 3
- 133217 F F T F 1 3 1 3 1
- 3 1 3
- 133317 T F F F 1 3 1 3 1
- 3 1 3
- 133318 T F F F 1 3 1 3 1
- 3 1 3
- 134562 T F F F 1 3 1 3 1
- 3 1 3
- 134563 T F F F 1 3 1 3 1
- 3 1 3
- 137743 F T F F 1 3 1 3 1
- 3 1 3
- 137744 F F T F 1 3 1 3 1
- 3 1 3
- 138208 F T F F 1 3 1 3 1
- 3 1 3
- 138209 F F T F 1 3 1 3 1
- 3 1 3
- 138214 F T F F 1 3 1 3 1
- 3 1 3
- 138215 F F T F 1 3 1 3 1
- 3 1 3
- 140356 F T F F 1 3 1 3 1
- 3 1 3
- 140357 F F T F 1 3 1 3 1
- 3 1 3
- 141912 T F F F 1 3 1 3 1
- 3 1 3
- 141913 T F F F 1 3 1 3 1
- 3 1 3
- 142369 F T F F 1 3 1 3 1
- 3 1 3
- 142370 F F T F 1 3 1 3 1
- 3 1 3
- 142488 T F F F 1 3 1 3 1
- 3 1 3
- 142489 T F F F 1 3 1 3 1
- 3 1 3
- 143719 F T F F 1 3 1 3 1
- 3 1 3
- 143720 F F T F 1 3 1 3 1
- 3 1 3
- 145321 F T F F 1 3 1 3 1
- 3 1 3
- 145322 F F T F 1 3 1 3 1
- 3 1 3
- 145329 T F F F 1 3 1 3 1
- 3 1 3
- 145330 T F F F 1 3 1 3 1
- 3 1 3
- 145336 F T F F 1 3 1 3 1
- 3 1 3
- 145337 F F T F 1 3 1 3 1
- 3 1 3
- 145668 T F F F 1 3 1 3 1
- 3 1 3
- 145669 T F F F 1 3 1 3 1
- 3 1 3
- 146469 T F F F 1 3 1 3 1
- 3 1 3
- 146470 T F F F 1 3 1 3 1
- 3 1 3
- 146884 F T F F 1 3 1 3 1
- 3 1 3
- 146885 F F T F 1 3 1 3 1
- 3 1 3
- 146961 T F F F 1 3 1 3 1
- 3 1 3
- 146962 T F F F 1 3 1 3 1
- 3 1 3
- 147219 F F F T 1 3 1 3 1
- 3 1 3
- 147221 F F F T 1 3 1 3 1
- 3 1 3
- 147334 F T F F 1 3 1 3 1
- 3 1 3
- 147335 F F T F 1 3 1 3 1
- 3 1 3
- 147433 F T F F 1 3 1 3 1
- 3 1 3
- 147434 F F T F 1 3 1 3 1
- 3 1 3
- 147450 F F F T 1 3 1 3 1
- 3 1 3
- 147452 F F F T 1 3 1 3 1
- 3 1 3
Modified: seismo/2D/SPECFEM2D/trunk/Database00000_right_edge_only
===================================================================
--- seismo/2D/SPECFEM2D/trunk/Database00000_right_edge_only 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/Database00000_right_edge_only 2011-02-15 00:12:10 UTC (rev 17863)
@@ -1,670 +1,384 @@
- 21695 F F F T 1 3 1 3 1
+ 19451 F F F T 1 3 1 3 1
3 1 3
- 21698 F F F T 1 3 1 3 1
+ 19451 F F F T 1 3 1 3 1
3 1 3
- 22883 F F F T 1 3 1 3 1
+ 19454 F F F T 1 3 1 3 1
3 1 3
- 22885 F F T F 1 3 1 3 1
+ 20963 F F F T 1 3 1 3 1
3 1 3
- 22886 F F T T 1 3 1 3 1
+ 20966 F F F T 1 3 1 3 1
3 1 3
- 23145 F F T F 1 3 1 3 1
+ 21121 F F T F 1 3 1 3 1
3 1 3
- 23146 F F T F 1 3 1 3 1
+ 21122 F F T F 1 3 1 3 1
3 1 3
- 23563 F F F T 1 3 1 3 1
+ 21887 F F F T 1 3 1 3 1
3 1 3
- 23566 F F F T 1 3 1 3 1
+ 21889 F F T F 1 3 1 3 1
3 1 3
- 23799 F F F T 1 3 1 3 1
+ 21890 F F T T 1 3 1 3 1
3 1 3
- 23802 F F F T 1 3 1 3 1
+ 22527 F F F T 1 3 1 3 1
3 1 3
- 23931 F F F T 1 3 1 3 1
+ 22530 F F F T 1 3 1 3 1
3 1 3
- 23934 F F F T 1 3 1 3 1
+ 23347 F F F T 1 3 1 3 1
3 1 3
- 25091 F F F T 1 3 1 3 1
+ 23350 F F F T 1 3 1 3 1
3 1 3
- 25094 F F F T 1 3 1 3 1
+ 26429 F F T F 1 3 1 3 1
3 1 3
- 25519 F F F T 1 3 1 3 1
+ 26430 F F T F 1 3 1 3 1
3 1 3
- 25522 F F F T 1 3 1 3 1
+ 26995 F F F T 1 3 1 3 1
3 1 3
- 26729 F F T F 1 3 1 3 1
+ 26998 F F F T 1 3 1 3 1
3 1 3
- 26730 F F T F 1 3 1 3 1
+ 27085 F F T F 1 3 1 3 1
3 1 3
- 29601 F F T F 1 3 1 3 1
+ 27086 F F T F 1 3 1 3 1
3 1 3
- 29602 F F T F 1 3 1 3 1
+ 28529 F F T F 1 3 1 3 1
3 1 3
- 30729 F F T F 1 3 1 3 1
+ 28530 F F T F 1 3 1 3 1
3 1 3
- 30730 F F T F 1 3 1 3 1
+ 29785 F F T F 1 3 1 3 1
3 1 3
- 30781 F F T F 1 3 1 3 1
+ 29786 F F T F 1 3 1 3 1
3 1 3
- 30782 F F T F 1 3 1 3 1
+ 29995 F F F T 1 3 1 3 1
3 1 3
- 32705 F F T F 1 3 1 3 1
+ 29998 F F F T 1 3 1 3 1
3 1 3
- 32706 F F T F 1 3 1 3 1
+ 30265 F F T F 1 3 1 3 1
3 1 3
- 33353 F F T F 1 3 1 3 1
+ 30266 F F T F 1 3 1 3 1
3 1 3
- 33354 F F T F 1 3 1 3 1
+ 30811 F F F T 1 3 1 3 1
3 1 3
- 33515 F F F T 1 3 1 3 1
+ 30814 F F F T 1 3 1 3 1
3 1 3
- 33518 F F F T 1 3 1 3 1
+ 30959 F F F T 1 3 1 3 1
3 1 3
- 34613 F F T F 1 3 1 3 1
+ 30962 F F F T 1 3 1 3 1
3 1 3
- 34614 F F T F 1 3 1 3 1
+ 32157 F F T F 1 3 1 3 1
3 1 3
- 34867 F F F T 1 3 1 3 1
+ 32158 F F T F 1 3 1 3 1
3 1 3
- 34870 F F F T 1 3 1 3 1
+ 33103 F F F T 1 3 1 3 1
3 1 3
- 36299 F F F T 1 3 1 3 1
+ 33106 F F F T 1 3 1 3 1
3 1 3
- 36302 F F F T 1 3 1 3 1
+ 33619 F F F T 1 3 1 3 1
3 1 3
- 36947 F F F T 1 3 1 3 1
+ 33622 F F F T 1 3 1 3 1
3 1 3
- 36950 F F F T 1 3 1 3 1
+ 33963 F F F T 1 3 1 3 1
3 1 3
- 37977 F F T F 1 3 1 3 1
+ 33966 F F F T 1 3 1 3 1
3 1 3
- 37978 F F T F 1 3 1 3 1
+ 36401 F F T F 1 3 1 3 1
3 1 3
- 38027 F F F T 1 3 1 3 1
+ 36402 F F T F 1 3 1 3 1
3 1 3
- 38030 F F F T 1 3 1 3 1
+ 38157 F F T F 1 3 1 3 1
3 1 3
- 38907 F F F T 1 3 1 3 1
+ 38158 F F T F 1 3 1 3 1
3 1 3
- 38910 F F F T 1 3 1 3 1
+ 39117 F F T F 1 3 1 3 1
3 1 3
- 39539 T F F F 1 3 1 3 1
+ 39118 F F T F 1 3 1 3 1
3 1 3
- 39540 T F F F 1 3 1 3 1
+ 39256 F T F F 1 3 1 3 1
3 1 3
- 39663 F F F T 1 3 1 3 1
+ 39257 F T F F 1 3 1 3 1
3 1 3
- 39666 F F F T 1 3 1 3 1
+ 40909 F F T F 1 3 1 3 1
3 1 3
- 40359 F F F T 1 3 1 3 1
+ 40910 F F T F 1 3 1 3 1
3 1 3
- 40362 F F F T 1 3 1 3 1
+ 41079 F F F T 1 3 1 3 1
3 1 3
- 41695 F F F T 1 3 1 3 1
+ 41082 F F F T 1 3 1 3 1
3 1 3
- 41698 F F F T 1 3 1 3 1
+ 41543 F F F T 1 3 1 3 1
3 1 3
- 43887 F F F T 1 3 1 3 1
+ 41546 F F F T 1 3 1 3 1
3 1 3
- 43890 F F F T 1 3 1 3 1
+ 45471 F F F T 1 3 1 3 1
3 1 3
- 43989 F F T F 1 3 1 3 1
+ 45474 F F F T 1 3 1 3 1
3 1 3
- 43990 F F T F 1 3 1 3 1
+ 46299 F F F T 1 3 1 3 1
3 1 3
- 44348 F T F F 1 3 1 3 1
+ 46302 F F F T 1 3 1 3 1
3 1 3
- 44349 F T F F 1 3 1 3 1
+ 47677 F F T F 1 3 1 3 1
3 1 3
- 45169 F F T F 1 3 1 3 1
+ 47678 F F T F 1 3 1 3 1
3 1 3
- 45170 F F T F 1 3 1 3 1
+ 48473 F F T F 1 3 1 3 1
3 1 3
- 45947 F F F T 1 3 1 3 1
+ 48474 F F T F 1 3 1 3 1
3 1 3
- 45950 F F F T 1 3 1 3 1
+ 48477 F F T F 1 3 1 3 1
3 1 3
- 46555 F F F T 1 3 1 3 1
+ 48478 F F T F 1 3 1 3 1
3 1 3
- 46558 F F F T 1 3 1 3 1
+ 50485 F F T F 1 3 1 3 1
3 1 3
- 47423 F F F T 1 3 1 3 1
+ 50486 F F T F 1 3 1 3 1
3 1 3
- 47426 F F F T 1 3 1 3 1
+ 52596 F T F F 1 3 1 3 1
3 1 3
- 48449 F F T F 1 3 1 3 1
+ 52597 F T F F 1 3 1 3 1
3 1 3
- 48450 F F T F 1 3 1 3 1
+ 53865 F F T F 1 3 1 3 1
3 1 3
- 48457 F F T F 1 3 1 3 1
+ 53866 F F T F 1 3 1 3 1
3 1 3
- 48458 F F T F 1 3 1 3 1
+ 53869 F F T F 1 3 1 3 1
3 1 3
- 50317 F F T F 1 3 1 3 1
+ 53870 F F T F 1 3 1 3 1
3 1 3
- 50318 F F T F 1 3 1 3 1
+ 54047 F F F T 1 3 1 3 1
3 1 3
- 53221 F F T F 1 3 1 3 1
+ 54050 F F F T 1 3 1 3 1
3 1 3
- 53222 F F T F 1 3 1 3 1
+ 54339 F F F T 1 3 1 3 1
3 1 3
- 53499 F F F T 1 3 1 3 1
+ 54342 F F F T 1 3 1 3 1
3 1 3
- 53502 F F F T 1 3 1 3 1
+ 55541 F F T F 1 3 1 3 1
3 1 3
- 54391 F F F T 1 3 1 3 1
+ 55542 F F T F 1 3 1 3 1
3 1 3
- 54394 F F F T 1 3 1 3 1
+ 55675 F F F T 1 3 1 3 1
3 1 3
- 56153 F F T F 1 3 1 3 1
+ 55678 F F F T 1 3 1 3 1
3 1 3
- 56154 F F T F 1 3 1 3 1
+ 56493 F F T F 1 3 1 3 1
3 1 3
- 56789 F F T F 1 3 1 3 1
+ 56494 F F T F 1 3 1 3 1
3 1 3
- 56790 F F T F 1 3 1 3 1
+ 56845 F F T F 1 3 1 3 1
3 1 3
- 56964 F T F F 1 3 1 3 1
+ 56846 F F T F 1 3 1 3 1
3 1 3
- 56965 F T F F 1 3 1 3 1
+ 57883 F F F T 1 3 1 3 1
3 1 3
- 57847 F F F T 1 3 1 3 1
+ 57886 F F F T 1 3 1 3 1
3 1 3
- 57850 F F F T 1 3 1 3 1
+ 60429 F F T F 1 3 1 3 1
3 1 3
- 61423 F F F T 1 3 1 3 1
+ 60430 F F T F 1 3 1 3 1
3 1 3
- 61426 F F F T 1 3 1 3 1
+ 62187 F F F T 1 3 1 3 1
3 1 3
- 61861 F F T F 1 3 1 3 1
+ 62190 F F F T 1 3 1 3 1
3 1 3
- 61862 F F T F 1 3 1 3 1
+ 63323 F F F T 1 3 1 3 1
3 1 3
- 62819 F F F T 1 3 1 3 1
+ 63326 F F F T 1 3 1 3 1
3 1 3
- 62822 F F F T 1 3 1 3 1
+ 64221 F F T F 1 3 1 3 1
3 1 3
- 63693 F F T F 1 3 1 3 1
+ 64222 F F T F 1 3 1 3 1
3 1 3
- 63694 F F T F 1 3 1 3 1
+ 65079 F F F T 1 3 1 3 1
3 1 3
- 64589 F F T F 1 3 1 3 1
+ 65082 F F F T 1 3 1 3 1
3 1 3
- 64590 F F T F 1 3 1 3 1
+ 66679 F F F T 1 3 1 3 1
3 1 3
- 64735 F F F T 1 3 1 3 1
+ 66682 F F F T 1 3 1 3 1
3 1 3
- 64738 F F F T 1 3 1 3 1
+ 66715 F F F T 1 3 1 3 1
3 1 3
- 65081 F F T F 1 3 1 3 1
+ 66718 F F F T 1 3 1 3 1
3 1 3
- 65082 F F T F 1 3 1 3 1
+ 67787 F F F T 1 3 1 3 1
3 1 3
- 66993 F F T F 1 3 1 3 1
+ 67790 F F F T 1 3 1 3 1
3 1 3
- 66994 F F T F 1 3 1 3 1
+ 67981 F F T F 1 3 1 3 1
3 1 3
- 67203 F F F T 1 3 1 3 1
+ 67982 F F T F 1 3 1 3 1
3 1 3
- 67206 F F F T 1 3 1 3 1
+ 70021 F F T F 1 3 1 3 1
3 1 3
- 67805 F F T F 1 3 1 3 1
+ 70022 F F T F 1 3 1 3 1
3 1 3
- 67806 F F T F 1 3 1 3 1
+ 70239 T F F F 1 3 1 3 1
3 1 3
- 68207 T F F F 1 3 1 3 1
+ 70240 T F F F 1 3 1 3 1
3 1 3
- 68208 T F F F 1 3 1 3 1
+ 70245 F F T F 1 3 1 3 1
3 1 3
- 68861 F F T F 1 3 1 3 1
+ 70246 F F T F 1 3 1 3 1
3 1 3
- 68862 F F T F 1 3 1 3 1
+ 70437 F F T F 1 3 1 3 1
3 1 3
- 69541 F F T F 1 3 1 3 1
+ 70438 F F T F 1 3 1 3 1
3 1 3
- 69542 F F T F 1 3 1 3 1
+ 70525 F F T F 1 3 1 3 1
3 1 3
- 71071 F F F T 1 3 1 3 1
+ 70526 F F T F 1 3 1 3 1
3 1 3
- 71074 F F F T 1 3 1 3 1
+ 71649 F F T F 1 3 1 3 1
3 1 3
- 71603 T F F F 1 3 1 3 1
+ 71650 F F T F 1 3 1 3 1
3 1 3
- 71604 T F F F 1 3 1 3 1
+ 72359 F F F T 1 3 1 3 1
3 1 3
- 72027 F F F T 1 3 1 3 1
+ 72362 F F F T 1 3 1 3 1
3 1 3
- 72030 F F F T 1 3 1 3 1
+ 72655 F F F T 1 3 1 3 1
3 1 3
- 73012 F T F F 1 3 1 3 1
+ 72658 F F F T 1 3 1 3 1
3 1 3
- 73013 F T F F 1 3 1 3 1
+ 73511 F F F T 1 3 1 3 1
3 1 3
- 73545 F F T F 1 3 1 3 1
+ 73514 F F F T 1 3 1 3 1
3 1 3
- 73546 F F T F 1 3 1 3 1
+ 74111 F F F T 1 3 1 3 1
3 1 3
- 74511 T F F F 1 3 1 3 1
+ 74114 F F F T 1 3 1 3 1
3 1 3
- 74512 T F F F 1 3 1 3 1
+ 74635 F F F T 1 3 1 3 1
3 1 3
- 74691 F F F T 1 3 1 3 1
+ 74638 F F F T 1 3 1 3 1
3 1 3
- 74694 F F F T 1 3 1 3 1
+ 76563 F F F T 1 3 1 3 1
3 1 3
- 74983 F F F T 1 3 1 3 1
+ 76566 F F F T 1 3 1 3 1
3 1 3
- 74986 F F F T 1 3 1 3 1
+ 77651 F F F T 1 3 1 3 1
3 1 3
- 75201 F F T F 1 3 1 3 1
+ 77654 F F F T 1 3 1 3 1
3 1 3
- 75202 F F T F 1 3 1 3 1
+ 77803 F F F T 1 3 1 3 1
3 1 3
- 75205 F F T F 1 3 1 3 1
+ 77806 F F F T 1 3 1 3 1
3 1 3
- 75206 F F T F 1 3 1 3 1
+ 80359 T F F F 1 3 1 3 1
3 1 3
- 75663 F F F T 1 3 1 3 1
+ 80360 T F F F 1 3 1 3 1
3 1 3
- 75666 F F F T 1 3 1 3 1
+ 80527 F F F T 1 3 1 3 1
3 1 3
- 76103 F F F T 1 3 1 3 1
+ 80530 F F F T 1 3 1 3 1
3 1 3
- 76106 F F F T 1 3 1 3 1
+ 80741 F F T F 1 3 1 3 1
3 1 3
- 76327 F F F T 1 3 1 3 1
+ 80742 F F T F 1 3 1 3 1
3 1 3
- 76330 F F F T 1 3 1 3 1
+ 80813 F F T F 1 3 1 3 1
3 1 3
- 77967 F F F T 1 3 1 3 1
+ 80814 F F T F 1 3 1 3 1
3 1 3
- 77970 F F F T 1 3 1 3 1
+ 81115 F F F T 1 3 1 3 1
3 1 3
- 77997 F F T F 1 3 1 3 1
+ 81118 F F F T 1 3 1 3 1
3 1 3
- 77998 F F T F 1 3 1 3 1
+ 82251 T F F F 1 3 1 3 1
3 1 3
- 78069 F F T F 1 3 1 3 1
+ 82252 T F F F 1 3 1 3 1
3 1 3
- 78070 F F T F 1 3 1 3 1
+ 83501 F F T F 1 3 1 3 1
3 1 3
- 79713 F F T F 1 3 1 3 1
+ 83502 F F T F 1 3 1 3 1
3 1 3
- 79714 F F T F 1 3 1 3 1
+ 84237 F F T F 1 3 1 3 1
3 1 3
- 80129 F F T F 1 3 1 3 1
+ 84238 F F T F 1 3 1 3 1
3 1 3
- 80130 F F T F 1 3 1 3 1
+ 85223 T F F F 1 3 1 3 1
3 1 3
- 81763 F F F T 1 3 1 3 1
+ 85224 T F F F 1 3 1 3 1
3 1 3
- 81766 F F F T 1 3 1 3 1
+ 85437 T F F F 1 3 1 3 1
3 1 3
- 81771 F F F T 1 3 1 3 1
+ 85438 T F F F 1 3 1 3 1
3 1 3
- 81774 F F F T 1 3 1 3 1
+ 85443 T F F F 1 3 1 3 1
3 1 3
- 81876 F T F F 1 3 1 3 1
+ 85444 T F F F 1 3 1 3 1
3 1 3
- 81877 F T F F 1 3 1 3 1
+ 85464 T F F F 1 3 1 3 1
3 1 3
- 83631 T F F F 1 3 1 3 1
+ 85465 T F F F 1 3 1 3 1
3 1 3
- 83632 T F F F 1 3 1 3 1
+ 85555 F T F F 1 3 1 3 1
3 1 3
- 84057 F F T F 1 3 1 3 1
+ 85556 F F T F 1 3 1 3 1
3 1 3
- 84058 F F T F 1 3 1 3 1
+ 85563 T F F F 1 3 1 3 1
3 1 3
- 84497 F F T F 1 3 1 3 1
+ 85564 T F F F 1 3 1 3 1
3 1 3
- 84498 F F T F 1 3 1 3 1
+ 85681 F T F F 1 3 1 3 1
3 1 3
- 84501 F F T F 1 3 1 3 1
+ 85682 F F T F 1 3 1 3 1
3 1 3
- 84502 F F T F 1 3 1 3 1
+ 85717 F T F F 1 3 1 3 1
3 1 3
- 85403 F F F T 1 3 1 3 1
+ 85718 F F T F 1 3 1 3 1
3 1 3
- 85406 F F F T 1 3 1 3 1
+ 86245 F T F F 1 3 1 3 1
3 1 3
- 85789 F F T F 1 3 1 3 1
+ 86246 F F T F 1 3 1 3 1
3 1 3
- 85790 F F T F 1 3 1 3 1
+ 86773 F T F F 1 3 1 3 1
3 1 3
- 90001 F F T F 1 3 1 3 1
+ 86774 F F T F 1 3 1 3 1
3 1 3
- 90002 F F T F 1 3 1 3 1
+ 87628 F T F F 1 3 1 3 1
3 1 3
- 91001 F F T F 1 3 1 3 1
+ 87629 F F T F 1 3 1 3 1
3 1 3
- 91002 F F T F 1 3 1 3 1
+ 88435 F T F F 1 3 1 3 1
3 1 3
- 92115 F F F T 1 3 1 3 1
+ 88436 F F T F 1 3 1 3 1
3 1 3
- 92118 F F F T 1 3 1 3 1
+ 89532 T F F F 1 3 1 3 1
3 1 3
- 93151 F F F T 1 3 1 3 1
+ 89533 T F F F 1 3 1 3 1
3 1 3
- 93154 F F F T 1 3 1 3 1
+ 89608 F T F F 1 3 1 3 1
3 1 3
- 95037 F F T F 1 3 1 3 1
+ 89609 F F T F 1 3 1 3 1
3 1 3
- 95038 F F T F 1 3 1 3 1
+ 91861 F T F F 1 3 1 3 1
3 1 3
- 97001 F F T F 1 3 1 3 1
+ 91862 F F T F 1 3 1 3 1
3 1 3
- 97002 F F T F 1 3 1 3 1
+ 91938 T F F F 1 3 1 3 1
3 1 3
- 97361 F F T F 1 3 1 3 1
+ 91939 T F F F 1 3 1 3 1
3 1 3
- 97362 F F T F 1 3 1 3 1
+ 92050 F T F F 1 3 1 3 1
3 1 3
- 100609 F F T F 1 3 1 3 1
+ 92051 F F T F 1 3 1 3 1
3 1 3
- 100610 F F T F 1 3 1 3 1
+ 92146 F T F F 1 3 1 3 1
3 1 3
- 101277 F F T F 1 3 1 3 1
+ 92147 F F T F 1 3 1 3 1
3 1 3
- 101278 F F T F 1 3 1 3 1
+ 92181 F F F T 1 3 1 3 1
3 1 3
- 102543 F F F T 1 3 1 3 1
+ 92183 F F F T 1 3 1 3 1
3 1 3
- 102546 F F F T 1 3 1 3 1
+ 94545 T F F F 1 3 1 3 1
3 1 3
- 102897 F F T F 1 3 1 3 1
+ 94546 T F F F 1 3 1 3 1
3 1 3
- 102898 F F T F 1 3 1 3 1
+ 97086 T F F F 1 3 1 3 1
3 1 3
- 103671 F F F T 1 3 1 3 1
+ 97087 T F F F 1 3 1 3 1
3 1 3
- 103674 F F F T 1 3 1 3 1
+ 100210 F T F F 1 3 1 3 1
3 1 3
- 103989 F F T F 1 3 1 3 1
+ 100211 F F T F 1 3 1 3 1
3 1 3
- 103990 F F T F 1 3 1 3 1
- 3 1 3
- 105461 F F T F 1 3 1 3 1
- 3 1 3
- 105462 F F T F 1 3 1 3 1
- 3 1 3
- 105903 F F F T 1 3 1 3 1
- 3 1 3
- 105906 F F F T 1 3 1 3 1
- 3 1 3
- 107629 F F T F 1 3 1 3 1
- 3 1 3
- 107630 F F T F 1 3 1 3 1
- 3 1 3
- 107755 F F F T 1 3 1 3 1
- 3 1 3
- 107758 F F F T 1 3 1 3 1
- 3 1 3
- 107997 F F T F 1 3 1 3 1
- 3 1 3
- 107998 F F T F 1 3 1 3 1
- 3 1 3
- 108629 F F T F 1 3 1 3 1
- 3 1 3
- 108630 F F T F 1 3 1 3 1
- 3 1 3
- 109521 F F T F 1 3 1 3 1
- 3 1 3
- 109522 F F T F 1 3 1 3 1
- 3 1 3
- 109979 F F F T 1 3 1 3 1
- 3 1 3
- 109982 F F F T 1 3 1 3 1
- 3 1 3
- 111359 F F F T 1 3 1 3 1
- 3 1 3
- 111362 F F F T 1 3 1 3 1
- 3 1 3
- 111569 F F T F 1 3 1 3 1
- 3 1 3
- 111570 F F T F 1 3 1 3 1
- 3 1 3
- 112599 F F F T 1 3 1 3 1
- 3 1 3
- 112602 F F F T 1 3 1 3 1
- 3 1 3
- 113001 F F T F 1 3 1 3 1
- 3 1 3
- 113002 F F T F 1 3 1 3 1
- 3 1 3
- 113013 F F T F 1 3 1 3 1
- 3 1 3
- 113014 F F T F 1 3 1 3 1
- 3 1 3
- 114455 F F F T 1 3 1 3 1
- 3 1 3
- 114458 F F F T 1 3 1 3 1
- 3 1 3
- 114535 T F F F 1 3 1 3 1
- 3 1 3
- 114536 T F F F 1 3 1 3 1
- 3 1 3
- 115089 F F T F 1 3 1 3 1
- 3 1 3
- 115090 F F T F 1 3 1 3 1
- 3 1 3
- 115109 F F T F 1 3 1 3 1
- 3 1 3
- 115110 F F T F 1 3 1 3 1
- 3 1 3
- 115117 F F T F 1 3 1 3 1
- 3 1 3
- 115118 F F T F 1 3 1 3 1
- 3 1 3
- 115165 F F T F 1 3 1 3 1
- 3 1 3
- 115166 F F T F 1 3 1 3 1
- 3 1 3
- 115625 F F T F 1 3 1 3 1
- 3 1 3
- 115626 F F T F 1 3 1 3 1
- 3 1 3
- 115967 F F F T 1 3 1 3 1
- 3 1 3
- 115970 F F F T 1 3 1 3 1
- 3 1 3
- 117047 T F F F 1 3 1 3 1
- 3 1 3
- 117048 T F F F 1 3 1 3 1
- 3 1 3
- 117115 T F F F 1 3 1 3 1
- 3 1 3
- 117116 T F F F 1 3 1 3 1
- 3 1 3
- 117819 F F F T 1 3 1 3 1
- 3 1 3
- 117822 F F F T 1 3 1 3 1
- 3 1 3
- 118359 T F F F 1 3 1 3 1
- 3 1 3
- 118360 T F F F 1 3 1 3 1
- 3 1 3
- 119624 F T F F 1 3 1 3 1
- 3 1 3
- 119625 F T F F 1 3 1 3 1
- 3 1 3
- 120007 F F F T 1 3 1 3 1
- 3 1 3
- 120010 F F F T 1 3 1 3 1
- 3 1 3
- 120025 F F T F 1 3 1 3 1
- 3 1 3
- 120026 F F T F 1 3 1 3 1
- 3 1 3
- 121235 T F F F 1 3 1 3 1
- 3 1 3
- 121236 T F F F 1 3 1 3 1
- 3 1 3
- 121263 F F F T 1 3 1 3 1
- 3 1 3
- 121266 F F F T 1 3 1 3 1
- 3 1 3
- 121497 T F F F 1 3 1 3 1
- 3 1 3
- 121498 T F F F 1 3 1 3 1
- 3 1 3
- 121503 T F F F 1 3 1 3 1
- 3 1 3
- 121504 T F F F 1 3 1 3 1
- 3 1 3
- 121521 T F F F 1 3 1 3 1
- 3 1 3
- 121522 T F F F 1 3 1 3 1
- 3 1 3
- 121621 F T F F 1 3 1 3 1
- 3 1 3
- 121622 F F T F 1 3 1 3 1
- 3 1 3
- 121626 T F F F 1 3 1 3 1
- 3 1 3
- 121627 T F F F 1 3 1 3 1
- 3 1 3
- 121735 F T F F 1 3 1 3 1
- 3 1 3
- 121736 F F T F 1 3 1 3 1
- 3 1 3
- 122338 F T F F 1 3 1 3 1
- 3 1 3
- 122339 F F T F 1 3 1 3 1
- 3 1 3
- 122857 F T F F 1 3 1 3 1
- 3 1 3
- 122858 F F T F 1 3 1 3 1
- 3 1 3
- 126580 F T F F 1 3 1 3 1
- 3 1 3
- 126581 F F T F 1 3 1 3 1
- 3 1 3
- 127486 F T F F 1 3 1 3 1
- 3 1 3
- 127487 F F T F 1 3 1 3 1
- 3 1 3
- 127680 T F F F 1 3 1 3 1
- 3 1 3
- 127681 T F F F 1 3 1 3 1
- 3 1 3
- 127902 T F F F 1 3 1 3 1
- 3 1 3
- 127903 T F F F 1 3 1 3 1
- 3 1 3
- 128077 F T F F 1 3 1 3 1
- 3 1 3
- 128078 F F T F 1 3 1 3 1
- 3 1 3
- 129691 F T F F 1 3 1 3 1
- 3 1 3
- 129692 F F T F 1 3 1 3 1
- 3 1 3
- 130441 F T F F 1 3 1 3 1
- 3 1 3
- 130442 F F T F 1 3 1 3 1
- 3 1 3
- 130794 T F F F 1 3 1 3 1
- 3 1 3
- 130795 T F F F 1 3 1 3 1
- 3 1 3
- 130965 T F F F 1 3 1 3 1
- 3 1 3
- 130966 T F F F 1 3 1 3 1
- 3 1 3
- 131031 T F F F 1 3 1 3 1
- 3 1 3
- 131032 T F F F 1 3 1 3 1
- 3 1 3
- 131095 F T F F 1 3 1 3 1
- 3 1 3
- 131096 F F T F 1 3 1 3 1
- 3 1 3
- 131335 F T F F 1 3 1 3 1
- 3 1 3
- 131336 F F T F 1 3 1 3 1
- 3 1 3
- 133168 F T F F 1 3 1 3 1
- 3 1 3
- 133169 F F T F 1 3 1 3 1
- 3 1 3
- 133198 F T F F 1 3 1 3 1
- 3 1 3
- 133199 F F T F 1 3 1 3 1
- 3 1 3
- 133203 T F F F 1 3 1 3 1
- 3 1 3
- 133204 T F F F 1 3 1 3 1
- 3 1 3
- 133263 T F F F 1 3 1 3 1
- 3 1 3
- 133264 T F F F 1 3 1 3 1
- 3 1 3
- 133269 T F F F 1 3 1 3 1
- 3 1 3
- 133270 T F F F 1 3 1 3 1
- 3 1 3
- 134574 T F F F 1 3 1 3 1
- 3 1 3
- 134575 T F F F 1 3 1 3 1
- 3 1 3
- 136371 T F F F 1 3 1 3 1
- 3 1 3
- 136372 T F F F 1 3 1 3 1
- 3 1 3
- 137533 F T F F 1 3 1 3 1
- 3 1 3
- 137534 F F T F 1 3 1 3 1
- 3 1 3
- 138216 T F F F 1 3 1 3 1
- 3 1 3
- 138217 T F F F 1 3 1 3 1
- 3 1 3
- 138313 F T F F 1 3 1 3 1
- 3 1 3
- 138314 F F T F 1 3 1 3 1
- 3 1 3
- 140307 T F F F 1 3 1 3 1
- 3 1 3
- 140308 T F F F 1 3 1 3 1
- 3 1 3
- 140383 F T F F 1 3 1 3 1
- 3 1 3
- 140384 F F T F 1 3 1 3 1
- 3 1 3
- 140476 F T F F 1 3 1 3 1
- 3 1 3
- 140477 F F T F 1 3 1 3 1
- 3 1 3
- 140833 F T F F 1 3 1 3 1
- 3 1 3
- 140834 F F T F 1 3 1 3 1
- 3 1 3
- 145278 T F F F 1 3 1 3 1
- 3 1 3
- 145279 T F F F 1 3 1 3 1
- 3 1 3
- 145297 F T F F 1 3 1 3 1
- 3 1 3
- 145298 F F T F 1 3 1 3 1
- 3 1 3
- 145410 F F F T 1 3 1 3 1
- 3 1 3
- 145412 F F F T 1 3 1 3 1
- 3 1 3
- 146145 T F F F 1 3 1 3 1
- 3 1 3
- 146146 T F F F 1 3 1 3 1
- 3 1 3
- 146238 T F F F 1 3 1 3 1
- 3 1 3
- 146239 T F F F 1 3 1 3 1
- 3 1 3
- 147148 F T F F 1 3 1 3 1
- 3 1 3
- 147149 F F T F 1 3 1 3 1
- 3 1 3
- 147175 F T F F 1 3 1 3 1
- 3 1 3
- 147176 F F T F 1 3 1 3 1
- 3 1 3
Modified: seismo/2D/SPECFEM2D/trunk/include_for_periodic_conditions.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/include_for_periodic_conditions.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/include_for_periodic_conditions.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -13,7 +13,17 @@
print *,iglob,' and ',iglob2,' are the same periodic point, merging them'
! print *,'horiz dist is = ',abs(coord(1,iglob) - coord(1,iglob2))
! print *,ispec,i,j,ispec2,i2,j2
- ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+! ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+!--------------------------------------------------------------------------------
+ iglob_target_to_replace = ibool(i2,j2,ispec2)
+ do ispec3 = 1,nspec
+ do j3 = 1,NGLLZ
+ do i3 = 1,NGLLX
+ if(ibool(i3,j3,ispec3) == iglob_target_to_replace) ibool(i3,j3,ispec3) = ibool(i,j,ispec)
+ enddo
+ enddo
+ enddo
+!--------------------------------------------------------------------------------
endif
enddo
endif
@@ -27,7 +37,17 @@
print *,iglob,' and ',iglob2,' are the same periodic point, merging them'
! print *,'horiz dist is = ',abs(coord(1,iglob) - coord(1,iglob2))
! print *,ispec,i,j,ispec2,i2,j2
- ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+! ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+!--------------------------------------------------------------------------------
+ iglob_target_to_replace = ibool(i2,j2,ispec2)
+ do ispec3 = 1,nspec
+ do j3 = 1,NGLLZ
+ do i3 = 1,NGLLX
+ if(ibool(i3,j3,ispec3) == iglob_target_to_replace) ibool(i3,j3,ispec3) = ibool(i,j,ispec)
+ enddo
+ enddo
+ enddo
+!--------------------------------------------------------------------------------
endif
enddo
endif
@@ -41,7 +61,17 @@
print *,iglob,' and ',iglob2,' are the same periodic point, merging them'
! print *,'horiz dist is = ',abs(coord(1,iglob) - coord(1,iglob2))
! print *,ispec,i,j,ispec2,i2,j2
- ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+! ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+!--------------------------------------------------------------------------------
+ iglob_target_to_replace = ibool(i2,j2,ispec2)
+ do ispec3 = 1,nspec
+ do j3 = 1,NGLLZ
+ do i3 = 1,NGLLX
+ if(ibool(i3,j3,ispec3) == iglob_target_to_replace) ibool(i3,j3,ispec3) = ibool(i,j,ispec)
+ enddo
+ enddo
+ enddo
+!--------------------------------------------------------------------------------
endif
enddo
endif
@@ -55,7 +85,17 @@
print *,iglob,' and ',iglob2,' are the same periodic point, merging them'
! print *,'horiz dist is = ',abs(coord(1,iglob) - coord(1,iglob2))
! print *,ispec,i,j,ispec2,i2,j2
- ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+! ibool(i2,j2,ispec2) = ibool(i,j,ispec)
+!--------------------------------------------------------------------------------
+ iglob_target_to_replace = ibool(i2,j2,ispec2)
+ do ispec3 = 1,nspec
+ do j3 = 1,NGLLZ
+ do i3 = 1,NGLLX
+ if(ibool(i3,j3,ispec3) == iglob_target_to_replace) ibool(i3,j3,ispec3) = ibool(i,j,ispec)
+ enddo
+ enddo
+ enddo
+!--------------------------------------------------------------------------------
endif
enddo
endif
Modified: seismo/2D/SPECFEM2D/trunk/meshfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/meshfem2D.F90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/meshfem2D.F90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -341,16 +341,16 @@
npoints_interface_top
integer :: number_of_layers
integer :: nz,nxread,nzread
-
- integer :: ilayer,ipoint_current
+
+ integer :: ilayer,ipoint_current
integer, dimension(:), pointer :: nz_layer
double precision, dimension(:), allocatable :: &
xinterface_bottom,zinterface_bottom,coefs_interface_bottom, &
xinterface_top,zinterface_top,coefs_interface_top
- integer :: nspec
- integer :: nbregion
-
+ integer :: nspec
+ integer :: nbregion
+
! external functions
integer, external :: num_4, num_9
double precision, external :: value_spline
@@ -387,29 +387,29 @@
if( ios /= 0 ) stop 'error opening DATA/Par_file file'
! reads in parameters in DATA/Par_file
- call read_parameter_file()
+ call read_parameter_file()
! reads in mesh elements
if ( read_external_mesh ) then
call read_external_mesh_file(mesh_file, num_start, ngnod)
- else
+ else
call read_interfaces_file(interfacesfile,max_npoints_interface, &
number_of_interfaces,npoints_interface_bottom, &
number_of_layers,nz_layer,nx,nz,nxread,nzread,ngnod, &
- nelmnts,elmnts)
+ nelmnts,elmnts)
endif
allocate(num_material(nelmnts))
num_material(:) = 0
-
+
! assigns materials to mesh elements
if ( read_external_mesh ) then
call read_mat(materials_file, num_material)
- else
+ else
call read_regions(nbregion,nb_materials,icodemat,cp,cs, &
rho_s,Qp,Qs,aniso3,aniso4,aniso5,aniso6,aniso7,aniso8, &
- nelmnts,num_material,nxread,nzread)
+ nelmnts,num_material,nxread,nzread)
endif
close(IIN)
@@ -908,7 +908,7 @@
! *** generate the databases for the solver
call save_databases(nspec,num_material, &
my_interfaces,my_nb_interfaces, &
- nnodes_tangential_curve,nodes_tangential_curve)
+ nnodes_tangential_curve,nodes_tangential_curve)
! print position of the source
do i_source=1,NSOURCE
@@ -922,7 +922,7 @@
if (generate_STATIONS) then
call save_stations_file(nreceiverlines,nrec,xdeb,zdeb,xfin,zfin,enreg_surf, &
xinterface_top,zinterface_top,coefs_interface_top, &
- npoints_interface_top,max_npoints_interface)
+ npoints_interface_top,max_npoints_interface)
endif
print *
Modified: seismo/2D/SPECFEM2D/trunk/part_unstruct.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/part_unstruct.F90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/part_unstruct.F90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -63,7 +63,7 @@
integer, dimension(:), pointer :: adjwgt
integer, dimension(:), pointer :: adjncy
-
+
integer :: nnodes
double precision, dimension(:,:), pointer :: nodes_coords
integer, dimension(:), pointer :: nnodes_elmnts
@@ -75,29 +75,29 @@
! interface data
integer :: ninterfaces
integer, dimension(:), pointer :: tab_size_interfaces, tab_interfaces
-
+
integer :: nelem_acoustic_surface
integer, dimension(:,:), pointer :: acoustic_surface
integer :: nelem_acoustic_surface_loc
-
+
integer :: nelemabs
integer, dimension(:,:), pointer :: abs_surface
logical, dimension(:,:), pointer :: abs_surface_char
integer, dimension(:), pointer :: abs_surface_merge
integer :: nelemabs_loc
-
+
integer :: nelemabs_merge
integer, dimension(:), pointer :: ibegin_bottom,iend_bottom,ibegin_top,iend_top, &
jbegin_left,jend_left,jbegin_right,jend_right
-
+
! for acoustic/elastic coupled elements
integer :: nedges_coupled
integer, dimension(:,:), pointer :: edges_coupled
-
+
! for acoustic/poroelastic coupled elements
integer :: nedges_acporo_coupled
integer, dimension(:,:), pointer :: edges_acporo_coupled
-
+
! for poroelastic/elastic coupled elements
integer :: nedges_elporo_coupled
integer, dimension(:,:), pointer :: edges_elporo_coupled
@@ -893,7 +893,7 @@
write(IIN_database,*) local_elmnt, surface(2,i), local_nodes(1), -1
endif
-
+
if ( surface(2,i) == 2 ) then
do l = glob2loc_nodes_nparts(surface(3,i)), &
glob2loc_nodes_nparts(surface(3,i)+1)-1
@@ -1579,5 +1579,5 @@
endif
end subroutine write_fluidsolid_edges_database
-
+
end module part_unstruct
Modified: seismo/2D/SPECFEM2D/trunk/read_interfaces_file.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/read_interfaces_file.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/read_interfaces_file.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -47,10 +47,10 @@
! note: we use this module definition only to be able to allocate
! arrays for receiverlines and materials in this subroutine rather than in the main
! routine in meshfem2D.F90
-
- ! note 2: the filename ending is .F90 to have pre-compilation with pragmas
+
+ ! note 2: the filename ending is .F90 to have pre-compilation with pragmas
! (like #ifndef USE_MPI) working properly
-
+
implicit none
contains
@@ -61,7 +61,7 @@
nelmnts,elmnts)
implicit none
include "constants.h"
-
+
character(len=100) :: interfacesfile
integer :: max_npoints_interface,number_of_interfaces,npoints_interface_bottom
@@ -142,7 +142,7 @@
nelmnts = nxread * nzread
allocate(elmnts(0:ngnod*nelmnts-1))
-
+
if ( ngnod == 4 ) then
num_elmnt = 0
do j = 1, nzread
@@ -174,6 +174,6 @@
endif
- end subroutine read_interfaces_file
-
+ end subroutine read_interfaces_file
+
end module interfaces_file
Modified: seismo/2D/SPECFEM2D/trunk/read_materials.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/read_materials.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/read_materials.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -52,9 +52,9 @@
implicit none
include "constants.h"
-
+
integer :: nb_materials
-
+
integer, dimension(nb_materials) :: icodemat
double precision, dimension(nb_materials) :: rho_s,cp,cs, &
@@ -66,7 +66,7 @@
integer :: imaterial,i,icodematread
double precision :: val0read,val1read,val2read,val3read,val4read, &
val5read,val6read,val7read,val8read,val9read,val10read,val11read,val12read
-
+
! initializes material properties
icodemat(:) = 0
cp(:) = 0.d0
@@ -99,16 +99,16 @@
val4read,val5read,val6read,val7read, &
val8read,val9read,val10read,val11read,val12read)
- ! checks material id
+ ! checks material id
if(i < 1 .or. i > nb_materials) stop 'Wrong material number!'
icodemat(i) = icodematread
! sets material properties
if(icodemat(i) == ISOTROPIC_MATERIAL) then
-
+
! isotropic materials
-
+
rho_s(i) = val0read
cp(i) = val1read
cs(i) = val2read
@@ -126,9 +126,9 @@
phi(i) = 1.d0 ! acoustic
endif
elseif (icodemat(i) == ANISOTROPIC_MATERIAL) then
-
+
! anisotropic materials
-
+
rho_s(i) = val0read
cp(i) = val1read
cs(i) = val2read
@@ -140,10 +140,10 @@
aniso8(i) = val8read
Qp(i) = val9read
Qs(i) = val10read
- else
-
+ else
+
! poroelastic materials
-
+
rho_s(i) = val0read
rho_f(i) = val1read
phi(i) = val2read
@@ -195,5 +195,5 @@
endif
print *
enddo
-
- end subroutine read_materials
\ No newline at end of file
+
+ end subroutine read_materials
Modified: seismo/2D/SPECFEM2D/trunk/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/read_parameter_file.F90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/read_parameter_file.F90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -47,13 +47,13 @@
! note: we use this module definition only to be able to allocate
! arrays for receiverlines and materials in this subroutine rather than in the main
! routine in meshfem2D.F90
-
- ! note 2: the filename ending is .F90 to have pre-compilation with pragmas
+
+ ! note 2: the filename ending is .F90 to have pre-compilation with pragmas
! (like #ifndef USE_MPI) working properly
-
+
implicit none
character(len=100) :: interfacesfile,title
-
+
integer :: SIMULATION_TYPE
logical :: SAVE_FORWARD,read_external_mesh
@@ -74,7 +74,7 @@
logical :: p_sv
logical :: any_abs,absbottom,absright,abstop,absleft
-
+
integer :: nt
double precision :: deltat
@@ -87,15 +87,15 @@
integer :: seismotype
logical :: generate_STATIONS
-
+
integer :: nreceiverlines
- double precision :: anglerec
+ double precision :: anglerec
logical :: rec_normal_to_surface
integer, dimension(:), pointer :: nrec
double precision, dimension(:), pointer :: xdeb,zdeb,xfin,zfin
logical, dimension(:), pointer :: enreg_surf
-
+
integer :: NTSTEP_BETWEEN_OUTPUT_INFO
logical :: output_postscript_snapshot,output_color_image
integer :: imagetype
@@ -117,12 +117,12 @@
contains
subroutine read_parameter_file()
-
+
! reads in DATA/Par_file
-
+
implicit none
include "constants.h"
-
+
! local parameters
integer :: ios,ireceiverlines
@@ -178,7 +178,7 @@
! read time step parameters
call read_value_integer(IIN,IGNORE_JUNK,nt)
call read_value_double_precision(IIN,IGNORE_JUNK,deltat)
-
+
! read source infos
call read_value_integer(IIN,IGNORE_JUNK,NSOURCE)
call read_value_logical(IIN,IGNORE_JUNK,force_normal_to_surface)
@@ -269,21 +269,21 @@
Qp,Qs,rho_s,rho_f,phi,tortuosity, &
permxx,permxz,permzz,kappa_s,kappa_f,kappa_fr, &
eta_f,mu_fr)
-
+
! checks input parameters
call check_parameters()
-
+
end subroutine read_parameter_file
-
+
!
!-------------------------------------------------------------------------------------------------
!
-
+
subroutine check_parameters()
-
+
implicit none
-
+
! checks partitioning
if ( nproc <= 0 ) then
print *, 'Number of processes (nproc) must be greater than or equal to one.'
@@ -320,8 +320,8 @@
else
plot_lowerleft_corner_only = .false.
endif
-
+
end subroutine check_parameters
-
+
end module parameter_file
-
+
Modified: seismo/2D/SPECFEM2D/trunk/read_regions.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/read_regions.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/read_regions.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -55,7 +55,7 @@
integer, dimension(nb_materials) :: icodemat
double precision, dimension(nb_materials) :: rho_s,cp,cs, &
aniso3,aniso4,aniso5,aniso6,aniso7,aniso8,Qp,Qs
-
+
integer :: nelmnts
integer,dimension(nelmnts) :: num_material
integer :: nxread,nzread
@@ -64,7 +64,7 @@
integer :: iregion,ixdebregion,ixfinregion,izdebregion,izfinregion,imaterial_number
integer :: i,j
double precision :: vpregion,vsregion,poisson_ratio
-
+
! read the material numbers for each region
call read_value_integer(IIN,IGNORE_JUNK,nbregion)
@@ -90,8 +90,8 @@
print *,'IZ from ',izdebregion,' to ',izfinregion
if(icodemat(imaterial_number) /= ANISOTROPIC_MATERIAL .and. icodemat(imaterial_number) /= POROELASTIC_MATERIAL) then
-
- ! isotropic material
+
+ ! isotropic material
vpregion = cp(imaterial_number)
vsregion = cs(imaterial_number)
print *,'Material # ',imaterial_number,' isotropic'
@@ -110,7 +110,7 @@
print *,'Qs = ',Qs(imaterial_number)
elseif(icodemat(imaterial_number) == POROELASTIC_MATERIAL) then
- ! poroelastic material
+ ! poroelastic material
print *,'Material # ',imaterial_number,' isotropic'
print *,'Material is poroelastic'
else
@@ -140,6 +140,6 @@
enddo
- if(minval(num_material) <= 0) stop 'Velocity model not entirely set...'
-
- end subroutine read_regions
\ No newline at end of file
+ if(minval(num_material) <= 0) stop 'Velocity model not entirely set...'
+
+ end subroutine read_regions
Modified: seismo/2D/SPECFEM2D/trunk/read_source_file.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/read_source_file.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/read_source_file.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -43,15 +43,15 @@
!========================================================================
module source_file
-
+
implicit none
! source parameters
- integer, dimension(:),pointer :: source_type,time_function_type
+ integer, dimension(:),pointer :: source_type,time_function_type
double precision, dimension(:),pointer :: xs,zs,f0,t0,angleforce, &
- Mxx,Mzz,Mxz,factor
+ Mxx,Mzz,Mxz,factor
logical, dimension(:),pointer :: source_surf
-
+
contains
subroutine read_source_file(NSOURCE,deltat,f0_attenuation)
@@ -60,7 +60,7 @@
implicit none
include "constants.h"
-
+
integer :: NSOURCE
double precision :: deltat,f0_attenuation
@@ -83,22 +83,22 @@
allocate(Mzz(NSOURCE))
allocate(factor(NSOURCE))
- ! counts lines
+ ! counts lines
open(unit=IIN_SOURCE,file='DATA/SOURCE',iostat=ios,status='old',action='read')
if(ios /= 0) stop 'error opening DATA/SOURCE file'
-
+
icounter = 0
do while(ios == 0)
read(IIN_SOURCE,"(a)",iostat=ios) dummystring
if(ios == 0) icounter = icounter + 1
enddo
close(IIN_SOURCE)
-
+
if(mod(icounter,NLINES_PER_SOURCE) /= 0) &
stop 'total number of lines in SOURCE file should be a multiple of NLINES_PER_SOURCE'
nsources = icounter / NLINES_PER_SOURCE
-
+
if(nsources < 1) stop 'need at least one source in SOURCE file'
if(nsources /= NSOURCE) &
stop 'total number of sources read is different than declared in Par_file'
@@ -122,13 +122,13 @@
! note: this is slightly different than in specfem2D.f90,
! t0 will be set outside of this next if statement, i.e. it will be set for all sources
! regardless of their type (it just makes a distinction between type 5 sources and the rest)
-
+
! if Dirac source time function, use a very thin Gaussian instead
! if Heaviside source time function, use a very thin error function instead
if(time_function_type(i_source) == 4 .or. time_function_type(i_source) == 5) then
f0(i_source) = 1.d0 / (10.d0 * deltat)
endif
-
+
! time delay of the source in seconds, use a 20 % security margin (use 2 / f0 if error function)
if(time_function_type(i_source)== 5) then
t0(i_source) = 2.0d0 / f0(i_source) + t0(i_source)
@@ -160,4 +160,4 @@
end subroutine read_source_file
end module source_file
-
\ No newline at end of file
+
Modified: seismo/2D/SPECFEM2D/trunk/read_value_parameters.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/read_value_parameters.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/read_value_parameters.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -100,7 +100,7 @@
logical ignore_junk
character(len=*) value_to_read
character(len=100) string_read
-
+
call read_next_line(iin,ignore_junk,string_read)
value_to_read = string_read
Modified: seismo/2D/SPECFEM2D/trunk/save_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/save_databases.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/save_databases.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -46,8 +46,8 @@
subroutine save_databases(nspec,num_material, &
my_interfaces,my_nb_interfaces, &
nnodes_tangential_curve,nodes_tangential_curve )
-
+
! generates the databases for the solver
use part_unstruct
@@ -55,7 +55,7 @@
use source_file
implicit none
include "constants.h"
-
+
integer :: nspec
integer, dimension(nelmnts) :: num_material
@@ -64,7 +64,7 @@
integer :: nnodes_tangential_curve
double precision, dimension(2,nnodes_tangential_curve) :: nodes_tangential_curve
-
+
! local parameters
integer :: iproc,i_source,i,ios
integer :: npgeo
@@ -72,17 +72,17 @@
integer :: nedges_coupled_loc
integer :: nedges_acporo_coupled_loc
integer :: nedges_elporo_coupled_loc
-
+
character(len=256) :: prname
-
-
+
+
do iproc = 0, nproc-1
! opens Database file
write(prname, "('./OUTPUT_FILES/Database',i5.5)") iproc
open(unit=15,file=trim(prname),status='unknown',iostat=ios)
if( ios /= 0 ) stop 'error saving databases'
-
+
write(15,*) '#'
write(15,*) '# Database for SPECFEM2D'
write(15,*) '# Dimitri Komatitsch, (c) University of Pau, France'
@@ -246,18 +246,18 @@
write(15,*) 'List of tangential detection curve nodes:'
!write(15,*) nnodes_tangential_curve
- write(15,*) force_normal_to_surface,rec_normal_to_surface
+ write(15,*) force_normal_to_surface,rec_normal_to_surface
if (force_normal_to_surface .or. rec_normal_to_surface) then
do i = 1, nnodes_tangential_curve
write(15,*) nodes_tangential_curve(1,i),nodes_tangential_curve(2,i)
enddo
endif
-
+
! closes Database file
close(15)
-
+
enddo
end subroutine save_databases
-
\ No newline at end of file
+
Modified: seismo/2D/SPECFEM2D/trunk/save_gnuplot_file.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/save_gnuplot_file.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/save_gnuplot_file.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -48,10 +48,10 @@
! creates a Gnuplot file that displays the grid
implicit none
-
+
integer :: ngnod,nx,nz
double precision, dimension(0:nx,0:nz) :: x,z
-
+
! local parameters
integer :: ios,istepx,istepz,ili,icol
@@ -60,7 +60,7 @@
open(unit=20,file='OUTPUT_FILES/gridfile.gnu',status='unknown',iostat=ios)
if( ios /= 0 ) stop 'error saving gnuplot file'
-
+
! draw horizontal lines of the grid
print *,'drawing horizontal lines of the grid'
istepx = 1
@@ -100,7 +100,7 @@
! create a Gnuplot script to display the grid
open(unit=20,file='OUTPUT_FILES/plotgnu',status='unknown',iostat=ios)
if( ios /= 0 ) stop 'error saving plotgnu file'
-
+
write(20,*) '#set term X11'
write(20,*) 'set term postscript landscape monochrome solid "Helvetica" 22'
write(20,*) 'set output "grid.ps"'
@@ -114,5 +114,5 @@
print *,'Grid saved in Gnuplot format...'
print *
-
- end subroutine save_gnuplot_file
\ No newline at end of file
+
+ end subroutine save_gnuplot_file
Modified: seismo/2D/SPECFEM2D/trunk/save_stations_file.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/save_stations_file.f90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/save_stations_file.f90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -45,7 +45,7 @@
subroutine save_stations_file(nreceiverlines,nrec,xdeb,zdeb,xfin,zfin,enreg_surf, &
xinterface_top,zinterface_top,coefs_interface_top, &
npoints_interface_top,max_npoints_interface)
-
+
implicit none
integer :: nreceiverlines
@@ -57,13 +57,13 @@
double precision, dimension(max_npoints_interface) :: xinterface_top, &
zinterface_top,coefs_interface_top
integer :: npoints_interface_top
-
+
!local parameters
integer :: ireceiverlines,irec,irec_global_number,ios
integer :: nrec_total
double precision :: xrec,zrec
double precision, external :: value_spline
-
+
print *
print *,'writing the DATA/STATIONS_target file'
print *
@@ -80,7 +80,7 @@
open(unit=15,file='DATA/STATIONS_target',status='unknown',iostat=ios)
if( ios /= 0 ) stop 'error saving STATIONS file'
-
+
irec_global_number = 0
! loop on all the receiver lines
@@ -119,4 +119,4 @@
close(15)
end subroutine save_stations_file
-
\ No newline at end of file
+
Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.F90 2011-02-14 23:01:54 UTC (rev 17862)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.F90 2011-02-15 00:12:10 UTC (rev 17863)
@@ -334,7 +334,7 @@
Mxx,Mzz,Mxz,f0,t0,factor,angleforce,hdur,hdur_gauss
real(kind=CUSTOM_REAL), dimension(:,:,:,:),allocatable :: sourcearray
double precision :: t0_start
-
+
double precision, dimension(:,:), allocatable :: coorg
double precision, dimension(:), allocatable :: coorgread
@@ -446,7 +446,7 @@
double precision :: vpImin,vpImax,vpIImin,vpIImax
integer :: colors,numbers,subsamp,imagetype,NTSTEP_BETWEEN_OUTPUT_INFO,NTSTEP_BETWEEN_OUTPUT_SEISMO,seismotype
- integer :: numat,ngnod,nspec,pointsdisp,nelemabs,nelem_acoustic_surface,ispecabs,UPPER_LIMIT_DISPLAY
+ integer :: numat,numat_local,ngnod,nspec,pointsdisp,nelemabs,nelem_acoustic_surface,ispecabs,UPPER_LIMIT_DISPLAY
logical interpol,meshvect,modelvect,boundvect,assign_external_model,initialfield, &
outputgrid,gnuplot,TURN_ATTENUATION_ON,output_postscript_snapshot,output_color_image, &
@@ -786,7 +786,7 @@
!! DK DK the length of an edge is about 1d-003, thus use e.g. 1/300 of that
double precision, parameter :: PERIODIC_DETECT_TOL = 1d-003 / 300.d0
- integer, parameter :: NSPEC_PERIO = 670 / 2
+ integer, parameter :: NSPEC_PERIO = 384 / 2 ! 670 / 2
integer, dimension(NSPEC_PERIO) :: numperio_left
integer, dimension(NSPEC_PERIO) :: numperio_right
@@ -796,6 +796,7 @@
integer :: idummy1, idummy2, idummy3, idummy4, idummy5, idummy6, idummy7, idummy8
integer :: ispecperio, ispecperio2, ispec2, i2, j2
+ integer :: iglob_target_to_replace, ispec3, i3, j3
!! DK DK Feb 2010 for periodic conditions: detect common points between left and right edges
@@ -1025,7 +1026,7 @@
! note: this is slightly different than in meshfem2D.f90,
! t0 will only be set within this if statement, i.e. only for type 4 or 5 sources
! (since we set f0 to a new values for these two types of sources)
-
+
! if Dirac source time function, use a very thin Gaussian instead
! if Heaviside source time function, use a very thin error function instead
if(time_function_type(i_source) == 4 .or. time_function_type(i_source) == 5) then
@@ -1038,7 +1039,7 @@
t0(i_source) = 1.20d0 / f0(i_source)+t0(i_source)
endif
endif
-
+
! for the source time function
aval(i_source) = pi*pi*f0(i_source)*f0(i_source)
@@ -1064,7 +1065,7 @@
endif
! checks if automatically set t0 is too small
- ! note: times in seismograms are shifted by t0(1)
+ ! note: times in seismograms are shifted by t0(1)
if( t0_start <= USER_T0 ) then
! sets new simulation start time such that
! simulation starts at t = - t0 = - USER_T0
@@ -1075,25 +1076,25 @@
write(IOUT,*) ' fix new simulation start time. . . . . = ', - t0_start
write(IOUT,*)
endif
-
+
! loops over all sources
do i_source=1,NSOURCE
! gets the given, initial time shifts
if( time_function_type(i_source) == 5 ) then
- t0(i_source) = t0(i_source) - 2.0d0 / f0(i_source)
+ t0(i_source) = t0(i_source) - 2.0d0 / f0(i_source)
else
- t0(i_source) = t0(i_source) - 1.20d0 / f0(i_source)
+ t0(i_source) = t0(i_source) - 1.20d0 / f0(i_source)
endif
-
+
! sets new t0 according to simulation start time,
! using absolute time shifts for each source such that
! a zero time shift would have the maximum gaussian at time t = (it-1)*DT - t0_start = 0
t0(i_source) = USER_T0 + t0(i_source)
-
+
if( myrank == 0 .and. ipass == 1) then
write(IOUT,*) ' source ',i_source,'uses t0. . . . . . = ',t0(i_source)
endif
-
+
enddo
else
@@ -1766,7 +1767,11 @@
enddo
close(123)
print *,'read ',NSPEC_PERIO,' elements for right periodic edge'
+print *
+print *,'because of periodic conditions, values computed by checkgrid() are not reliable'
+print *
+
!---------------------------------------------------------------------------
do ispecperio = 1,NSPEC_PERIO
@@ -3944,7 +3949,7 @@
write(IOUT,*)
write(IOUT,*) '*** calculation of the initial plane wave ***'
write(IOUT,*)
- write(IOUT,*) 'To change the initial plane wave, change source_type in DATA/Par_file'
+ write(IOUT,*) 'To change the initial plane wave, change source_type in DATA/SOURCE'
write(IOUT,*) 'and use 1 for a plane P wave, 2 for a plane SV wave, 3 for a Rayleigh wave'
write(IOUT,*)
@@ -3966,11 +3971,15 @@
endif
endif
! only implemented for homogeneous media therefore only 1 material supported
- if (numat==1) then
+ numat_local = numat
+ if (numat /= 1) then
+ if (myrank == 0) write(IOUT,*) 'not possible to have several materials with a plane wave, using the first material'
+ numat_local = 1
+ endif
- mu = poroelastcoef(2,1,numat)
- lambdaplus2mu = poroelastcoef(3,1,numat)
- denst = density(1,numat)
+ mu = poroelastcoef(2,1,numat_local)
+ lambdaplus2mu = poroelastcoef(3,1,numat_local)
+ denst = density(1,numat_local)
cploc = sqrt(lambdaplus2mu/denst)
csloc = sqrt(mu/denst)
@@ -4049,9 +4058,6 @@
B_plane(1)=0.d0; B_plane(2)=0.d0
C_plane(1)=0.d0; C_plane(2)=0.d0
endif
- else
- call exit_MPI('not possible to have several materials with a plane wave')
- endif
! get minimum and maximum values of mesh coordinates
xmin = minval(coord(1,:))
@@ -4259,7 +4265,7 @@
! compute current time
time = (it-1)*deltat
- stf_used = 0.0
+ stf_used = 0._CUSTOM_REAL
! loop on all the sources
do i_source=1,NSOURCE