[cig-commits] r19200 - in seismo/2D/SPECFEM2D/trunk: DATA doc/USER_MANUAL src/meshfem2D src/specfem2D

dkomati1 at geodynamics.org dkomati1 at geodynamics.org
Mon Nov 14 16:06:37 PST 2011 

Author: dkomati1
Date: 2011-11-14 16:06:36 -0800 (Mon, 14 Nov 2011)
New Revision: 19200

Modified:
   seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in
   seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf
   seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
   seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
   seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
first step of making periodic conditions more flexible (not finished yet)


Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file.in	2011-11-15 00:06:36 UTC (rev 19200)
@@ -93,6 +93,15 @@
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 

Modified: seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf	2011-11-15 00:06:36 UTC (rev 19200)
@@ -22676,9 +22676,8 @@
 /Filter /FlateDecode
 >>
 stream
-xÚuMkƒ@†ïþŠ=*ÔÉÎìŽÇ$UH!iˆÚKÓƒIEÁ˜Bÿ}WW)
-ô4³óñ>ïÎ*w©VC d-ò« Ö Lªc	‘R"¿ˆw7Û'ë4ÙÒ³÷‘¿ŸbŒl¿¸UçKw[6÷òs3ÂZ b¦AØWbMÂ×ʳ]<£]Üwus®ú¶±Ï¢©=b÷ËCv«îV÷߶~”,‹li´W6—Ñ”O26¾PÌìë•ö:A‡Fy·…
-ëÝ!ûSØì›å/=íJct40pâàûæ!✿ Y—¨AG h:îë¹oOóý0x²	IN’;?rleÔ
+xÚuMkƒ@†ïþŠ=*ÔÉÎìNÔcš*¤	Q{ir0¡]Á$…þû®®RZèi¾ß÷Ù},½E¦•ÀY‹ò"ˆ5(›êDB¬”(ÏâÕ/vé:K·ôËg¡€B‚±›WצBbéoks¯ß‡5+¬Ú%f„C¥!Ñ$B­¬0»Ã1ºÃ]ßš·æÖWV¦
+ˆýÏ Ùoúk{ûrýƒdY+¦»ÚœG¨db¹PY›™ë•î2™ƒúë|_üjlòýfõãžõµß(‚xùÇö%@ÄÑp~Îh¢”@Ú†húÛ¼ âÇiþG䗐ÄÑÈKKgz
 endstream
 endobj
 2207 0 obj <<
@@ -32288,17 +32287,17 @@
 /ProcSet [ /PDF /Text ]
 >> endobj
 2684 0 obj <<
-/Length 924       
+/Length 698       
 /Filter /FlateDecode
 >>
 stream
-xÚuUKsÛ6¾ûWèVp&DRÅ[7ö8ãd22=4HBjPñˆ£Ÿ’-zÂ~À¾¹ßòÝæê÷Ûe¹(–´^×|±Ù.êb±*-ËÕbÓ/þ&×Ý£6gäi”ý0Ií]öÏæ˜-
-mVÅ*Z±E^qºfÉè¯l͈Ìò¢X¿úÅϯˆ#¢«œWZ!ÚZ3¡ô Zc…7ÊÎ*Ò'Oe'´ú/$(=žýoם	àêôb¶É@X'Õ8Ê77äÖ
-ÝÁÖXtÿ#f%ím:aå6Œc‚^ƇXr,’sÚT™ìœ2t«’¬(GaŽ§ßÅ„ª‚ìE÷(‰·B÷(ì丏¾gž£é„­ù±â'ºr2{ÎYŒø®¦½±^h°Í
-FÂàb¥Õ’‡¦OrLú"SÚ[Ӈûè	v(Õq¡w“ÛEF&évÒ¢¬4ž©xv,žÎ.ø’609g~6Q³,r'‚sùƒi…÷&±ÈAêÞÎïkÈ®µ&À÷
- Ö<‹«úÅ,òuAyQA¶³÷a¿ªµÂè¶a†·))+×G}aS6­p²Oq2N4Š{Û%“nœÙú§ÈqÌéc–óhÇâó¢a”³ú¼|lc0®Î‰nœôÞáÕ½†ñóÁ'y°àÜd5ð«Ûi3šá­*?5¤ð‡ŒùE–ž u€†ª›ð{=@¥ÅßXÅnÄÔZÕ29;Ë'Ý}ýr
-ªœ^š‡ù;|M´Ñ¹ó6t>ØØÂx7ŒðYÇô¦ö›×¯9{kâ üPýÑW{ÀóÅ>t—E—·Êu8‘ðk{N€7O¯¾ÂÚy¦‘SþÔĸ¦ž×EEþÜ™þœi×Ý»xúgÚ“$ñçµÍb+ŒÖ‰Ÿ„‡õ4f§ÍgOÙT¯:~‰>ŸG	“šå¼™È'¡F”&X”¶Rö-°0zd=ébX—pg¦y¯'8ohæÂ0$M¼ðÏÒ©QÞdeE`À÷c¯/p²,kÊšòȱ}´È³<ödƒç¢“Þh¬4/«hùŠ+²ÌkÚ”'ç±[Ø:“VRÌûl%ÝÜß%Òþì#|‘ÊMÈ«„†¦ý_~Y(ýO¡EQÓzuÊ¥SVçf¯oiúƒ“ê5v¸TKUцW/‹I
-áuA«¦XäËrÍÄwÞÄÇ«÷›«_“âê
+xÚuTMSÛ0½ó+|”gj×’#;¾(00ÐaŠ9µŽl+ŽŠ-¥’“äUBÒIOû¡Ý·»O+]ÔgŸ¯yDi¹,qT¯¢’DEž¥y^Duý@çí‹T1ÎÐÛÀ»~äÒšøW}çÒQ•V)|V%§Ë,$Õk'y¾@7l2&¹W
+³V%÷<Æõ\v:œ›©Ñj²Br!½(Ðñ§U^¹*3z?¿Ñh¦·éªÊN´S‘´ªÈ.žíª5Ìðn_g–¦-h*xŒZÙ7?ð>
+j‰]‹+­FЬӏR.\ƒ8J°k”R(ûÀŒaíz2ÜÚ9u‰n¥±ÂNv®Z¬ã2C¼]K5¨~}rE†¾rߟ'ý0a‰\0“¢eø¯dï(äZÈ~f4»dc£E×ó÷|ÏOç.§óÙqÿó%¼DRÉÄX=µvÒž?ïëw§C8ŸÆ&Ô&¸ü‡@ –c£•ß‚WÑí°š-ÈGïgS@½Ž—JòZ˜vÍ5øÙ.¡èM‚ëY
+o¾z ®°;	Y Ç‚²å³“¢ïkÕñSSßšÉÎGÑÌÓ4	õÀ`乎=JJŽoÌ
+%g^œuÏ¥™Uzß
+@~0ŽÝÚ¸Å=jàqànMãWn'’‘‰ô­š4h+놵/1#èïč/k‚ݪq~¦ÁtÓ@š™ú>D‚Ã*w\+.A¯ãœ"·à0=×'dž—iVå»öÇúŒ$N<'5ȍ[ˆ–[%aÒ$§$ÍòòxVÞM'^/.ÝcÚƒ{¶€:¥Aú¾g¹´_ÞÞŵ1Ù£KxâÂŒð®fûÁ
+«êoH˜p•ÿ”2-‹}/­èa:3£~é¹ê¶’¢5©Òý©Y(M+L‡	„à’¤Ô}RÉ¢<`Wþôìª>{Ý8…S
 endstream
 endobj
 2683 0 obj <<
@@ -32312,19 +32311,19 @@
 2680 0 obj <<
 /Type /Annot
 /Border[0 0 0]/H/I/C[0 1 1]
-/Rect [423.04 491.691 540.996 502.595]
+/Rect [423.04 527.557 540.996 538.461]
 /Subtype/Link/A<</Type/Action/S/URI/URI(jtromp-AT-princeton.edu)>>
 >> endobj
 2687 0 obj <<
 /Type /Annot
 /Border[0 0 0]/H/I/C[0 1 1]
-/Rect [71.004 479.736 90.929 490.64]
+/Rect [71.004 515.602 90.929 526.506]
 /Subtype/Link/A<</Type/Action/S/URI/URI(jtromp-AT-princeton.edu)>>
 >> endobj
 2681 0 obj <<
 /Type /Annot
 /Border[0 0 0]/H/I/C[0 1 1]
-/Rect [316.7 479.736 474.608 490.64]
+/Rect [316.7 515.602 474.608 526.506]
 /Subtype/Link/A<</Type/Action/S/URI/URI(cig-seismo at geodynamics.org)>>
 >> endobj
 2685 0 obj <<
@@ -34614,8 +34613,8 @@
 >> endobj
 2800 0 obj <<
 /Author()/Title()/Subject()/Creator(LaTeX with hyperref package)/Producer(pdfTeX-1.40.3)/Keywords()
-/CreationDate (D:20111016005609+02'00')
-/ModDate (D:20111016005609+02'00')
+/CreationDate (D:20111115003757+01'00')
+/ModDate (D:20111115003757+01'00')
 /Trapped /False
 /PTEX.Fullbanner (This is pdfTeX using libpoppler, Version 3.141592-1.40.3-2.2 (Web2C 7.5.6) kpathsea version 3.5.6)
 >> endobj
@@ -34628,67 +34627,67 @@
 0000000000 00000 f 
 0000000015 00000 n 
 0001207812 00000 n 
-0002281689 00000 n 
+0002281465 00000 n 
 0000000060 00000 n 
 0000000092 00000 n 
 0001228003 00000 n 
-0002281617 00000 n 
+0002281393 00000 n 
 0000000139 00000 n 
 0000000170 00000 n 
 0001228060 00000 n 
-0002281544 00000 n 
+0002281320 00000 n 
 0000000218 00000 n 
 0000000248 00000 n 
 0001231446 00000 n 
-0002281456 00000 n 
+0002281232 00000 n 
 0000000294 00000 n 
 0000000330 00000 n 
 0001594004 00000 n 
-0002281330 00000 n 
+0002281106 00000 n 
 0000000376 00000 n 
 0000000412 00000 n 
 0001594057 00000 n 
-0002281256 00000 n 
+0002281032 00000 n 
 0000000460 00000 n 
 0000000503 00000 n 
 0001663530 00000 n 
-0002281169 00000 n 
+0002280945 00000 n 
 0000000551 00000 n 
 0000000618 00000 n 
 0001734083 00000 n 
-0002281095 00000 n 
+0002280871 00000 n 
 0000000666 00000 n 
 0000000732 00000 n 
 0002088076 00000 n 
-0002280969 00000 n 
+0002280745 00000 n 
 0000000778 00000 n 
 0000000828 00000 n 
 0002103106 00000 n 
-0002280895 00000 n 
+0002280671 00000 n 
 0000000876 00000 n 
 0000000950 00000 n 
 0002106425 00000 n 
-0002280808 00000 n 
+0002280584 00000 n 
 0000000998 00000 n 
 0000001065 00000 n 
 0002106478 00000 n 
-0002280734 00000 n 
+0002280510 00000 n 
 0000001113 00000 n 
 0000001161 00000 n 
 0002115940 00000 n 
-0002280608 00000 n 
+0002280384 00000 n 
 0000001207 00000 n 
 0000001247 00000 n 
 0002115993 00000 n 
-0002280547 00000 n 
+0002280323 00000 n 
 0000001295 00000 n 
 0000001358 00000 n 
 0002118831 00000 n 
-0002280458 00000 n 
+0002280234 00000 n 
 0000001404 00000 n 
 0000001442 00000 n 
-0002145501 00000 n 
-0002280382 00000 n 
+0002145277 00000 n 
+0002280158 00000 n 
 0000001489 00000 n 
 0000001525 00000 n 
 0000001717 00000 n 
@@ -34697,7 +34696,7 @@
 0000001576 00000 n 
 0001179668 00000 n 
 0001179719 00000 n 
-0002279606 00000 n 
+0002279382 00000 n 
 0000005386 00000 n 
 0000005574 00000 n 
 0000005652 00000 n 
@@ -36830,8 +36829,8 @@
 0001180291 00000 n 
 0001180177 00000 n 
 0001179837 00000 n 
-0002278884 00000 n 
-0002279061 00000 n 
+0002278660 00000 n 
+0002278837 00000 n 
 0001181279 00000 n 
 0001181106 00000 n 
 0001180379 00000 n 
@@ -36979,62 +36978,62 @@
 0001188138 00000 n 
 0001185049 00000 n 
 0001207757 00000 n 
-0002278351 00000 n 
+0002278127 00000 n 
 0001189833 00000 n 
 0001190221 00000 n 
 0001193083 00000 n 
 0001199156 00000 n 
 0001200239 00000 n 
-0002277300 00000 n 
+0002277076 00000 n 
 0001202252 00000 n 
 0001205711 00000 n 
-0002125383 00000 n 
-0002125572 00000 n 
-0002135098 00000 n 
-0002125635 00000 n 
-0002135350 00000 n 
-0002135035 00000 n 
-0002143777 00000 n 
-0002124943 00000 n 
-0002140102 00000 n 
-0002125069 00000 n 
-0002130231 00000 n 
-0002130483 00000 n 
-0002129539 00000 n 
-0002130294 00000 n 
-0002140165 00000 n 
-0002139598 00000 n 
-0002130420 00000 n 
-0002125006 00000 n 
-0002125320 00000 n 
-0002130042 00000 n 
-0002129854 00000 n 
-0002129917 00000 n 
-0002129791 00000 n 
-0002130105 00000 n 
-0002139787 00000 n 
-0002139095 00000 n 
-0002143273 00000 n 
-0002129980 00000 n 
-0002129602 00000 n 
-0002129728 00000 n 
-0002143210 00000 n 
-0002130735 00000 n 
-0002134595 00000 n 
-0002143651 00000 n 
-0002134720 00000 n 
-0002125446 00000 n 
-0002134972 00000 n 
-0002134783 00000 n 
-0002134846 00000 n 
-0002130546 00000 n 
-0002143840 00000 n 
-0002143525 00000 n 
-0002134532 00000 n 
-0002139472 00000 n 
-0002125131 00000 n 
-0002125194 00000 n 
-0002143399 00000 n 
+0002125159 00000 n 
+0002125348 00000 n 
+0002134874 00000 n 
+0002125411 00000 n 
+0002135126 00000 n 
+0002134811 00000 n 
+0002143553 00000 n 
+0002124719 00000 n 
+0002139878 00000 n 
+0002124845 00000 n 
+0002130007 00000 n 
+0002130259 00000 n 
+0002129315 00000 n 
+0002130070 00000 n 
+0002139941 00000 n 
+0002139374 00000 n 
+0002130196 00000 n 
+0002124782 00000 n 
+0002125096 00000 n 
+0002129818 00000 n 
+0002129630 00000 n 
+0002129693 00000 n 
+0002129567 00000 n 
+0002129881 00000 n 
+0002139563 00000 n 
+0002138871 00000 n 
+0002143049 00000 n 
+0002129756 00000 n 
+0002129378 00000 n 
+0002129504 00000 n 
+0002142986 00000 n 
+0002130511 00000 n 
+0002134371 00000 n 
+0002143427 00000 n 
+0002134496 00000 n 
+0002125222 00000 n 
+0002134748 00000 n 
+0002134559 00000 n 
+0002134622 00000 n 
+0002130322 00000 n 
+0002143616 00000 n 
+0002143301 00000 n 
+0002134308 00000 n 
+0002139248 00000 n 
+0002124907 00000 n 
+0002124970 00000 n 
+0002143175 00000 n 
 0001213159 00000 n 
 0001213468 00000 n 
 0001213623 00000 n 
@@ -37130,41 +37129,41 @@
 0001210832 00000 n 
 0001207980 00000 n 
 0001227948 00000 n 
-0002279419 00000 n 
+0002279195 00000 n 
 0001213314 00000 n 
 0001214736 00000 n 
 0001226498 00000 n 
-0002134658 00000 n 
-0002135161 00000 n 
-0002135224 00000 n 
-0002135287 00000 n 
-0002134469 00000 n 
-0002134406 00000 n 
-0002125698 00000 n 
-0002143588 00000 n 
-0002143147 00000 n 
-0002139913 00000 n 
-0002129665 00000 n 
-0002130609 00000 n 
-0002130168 00000 n 
-0002125509 00000 n 
-0002140039 00000 n 
-0002143336 00000 n 
-0002130357 00000 n 
-0002139850 00000 n 
-0002139535 00000 n 
-0002130672 00000 n 
-0002139284 00000 n 
-0002139346 00000 n 
-0002139158 00000 n 
-0002139409 00000 n 
-0002143714 00000 n 
-0002134909 00000 n 
-0002139032 00000 n 
-0002139661 00000 n 
-0002135413 00000 n 
-0002139724 00000 n 
-0002139221 00000 n 
+0002134434 00000 n 
+0002134937 00000 n 
+0002135000 00000 n 
+0002135063 00000 n 
+0002134245 00000 n 
+0002134182 00000 n 
+0002125474 00000 n 
+0002143364 00000 n 
+0002142923 00000 n 
+0002139689 00000 n 
+0002129441 00000 n 
+0002130385 00000 n 
+0002129944 00000 n 
+0002125285 00000 n 
+0002139815 00000 n 
+0002143112 00000 n 
+0002130133 00000 n 
+0002139626 00000 n 
+0002139311 00000 n 
+0002130448 00000 n 
+0002139060 00000 n 
+0002139122 00000 n 
+0002138934 00000 n 
+0002139185 00000 n 
+0002143490 00000 n 
+0002134685 00000 n 
+0002138808 00000 n 
+0002139437 00000 n 
+0002135189 00000 n 
+0002139500 00000 n 
+0002138997 00000 n 
 0001230486 00000 n 
 0001230640 00000 n 
 0001230795 00000 n 
@@ -37174,8 +37173,8 @@
 0001228233 00000 n 
 0001231391 00000 n 
 0001231188 00000 n 
-0002279727 00000 n 
-0002139976 00000 n 
+0002279503 00000 n 
+0002139752 00000 n 
 0001232146 00000 n 
 0001231975 00000 n 
 0001231601 00000 n 
@@ -37187,7 +37186,7 @@
 0001232248 00000 n 
 0001593949 00000 n 
 0001594114 00000 n 
-0002278176 00000 n 
+0002277952 00000 n 
 0001597014 00000 n 
 0001663257 00000 n 
 0001663649 00000 n 
@@ -37195,8 +37194,8 @@
 0001594331 00000 n 
 0001663411 00000 n 
 0001663466 00000 n 
-0002278526 00000 n 
-0002278705 00000 n 
+0002278302 00000 n 
+0002278481 00000 n 
 0001663587 00000 n 
 0001661234 00000 n 
 0001661497 00000 n 
@@ -37232,7 +37231,7 @@
 0001735597 00000 n 
 0002088021 00000 n 
 0002088129 00000 n 
-0002279852 00000 n 
+0002279628 00000 n 
 0002091795 00000 n 
 0002091993 00000 n 
 0002092152 00000 n 
@@ -37247,8 +37246,8 @@
 0002091586 00000 n 
 0002088359 00000 n 
 0002093250 00000 n 
-0002277594 00000 n 
-0002125257 00000 n 
+0002277370 00000 n 
+0002125033 00000 n 
 0002102925 00000 n 
 0002102556 00000 n 
 0002102713 00000 n 
@@ -37256,7 +37255,7 @@
 0002095498 00000 n 
 0002093448 00000 n 
 0002102870 00000 n 
-0002277883 00000 n 
+0002277659 00000 n 
 0002102989 00000 n 
 0002103048 00000 n 
 0002099778 00000 n 
@@ -37265,7 +37264,7 @@
 0002100443 00000 n 
 0002100466 00000 n 
 0002100749 00000 n 
-0002143462 00000 n 
+0002143238 00000 n 
 0002106060 00000 n 
 0002106215 00000 n 
 0002108010 00000 n 
@@ -37292,141 +37291,141 @@
 0002116152 00000 n 
 0002117097 00000 n 
 0002117152 00000 n 
-0002279977 00000 n 
+0002279753 00000 n 
 0002118884 00000 n 
 0002118660 00000 n 
 0002117327 00000 n 
 0002118776 00000 n 
-0002277152 00000 n 
-0002277739 00000 n 
-0002277447 00000 n 
-0002278029 00000 n 
-0002120229 00000 n 
-0002120569 00000 n 
-0002120855 00000 n 
-0002120074 00000 n 
+0002276928 00000 n 
+0002277515 00000 n 
+0002277223 00000 n 
+0002277805 00000 n 
+0002120003 00000 n 
+0002120344 00000 n 
+0002120631 00000 n 
+0002119848 00000 n 
 0002119068 00000 n 
-0002120741 00000 n 
-0002120796 00000 n 
-0002120400 00000 n 
-0002121539 00000 n 
-0002121309 00000 n 
-0002120957 00000 n 
-0002121425 00000 n 
-0002121480 00000 n 
-0002125761 00000 n 
-0002124713 00000 n 
-0002121627 00000 n 
-0002124829 00000 n 
-0002124884 00000 n 
-0002279238 00000 n 
-0002130796 00000 n 
-0002129368 00000 n 
-0002125863 00000 n 
-0002129484 00000 n 
-0002135476 00000 n 
-0002134235 00000 n 
-0002130898 00000 n 
-0002134351 00000 n 
-0002280102 00000 n 
-0002140228 00000 n 
-0002138861 00000 n 
-0002135578 00000 n 
-0002138977 00000 n 
-0002143903 00000 n 
-0002142976 00000 n 
-0002140330 00000 n 
-0002143092 00000 n 
-0002145613 00000 n 
+0002120517 00000 n 
+0002120572 00000 n 
+0002120174 00000 n 
+0002121315 00000 n 
+0002121085 00000 n 
+0002120733 00000 n 
+0002121201 00000 n 
+0002121256 00000 n 
+0002125537 00000 n 
+0002124489 00000 n 
+0002121403 00000 n 
+0002124605 00000 n 
+0002124660 00000 n 
+0002279014 00000 n 
+0002130572 00000 n 
+0002129144 00000 n 
+0002125639 00000 n 
+0002129260 00000 n 
+0002135252 00000 n 
+0002134011 00000 n 
+0002130674 00000 n 
+0002134127 00000 n 
+0002279878 00000 n 
+0002140004 00000 n 
+0002138637 00000 n 
+0002135354 00000 n 
+0002138753 00000 n 
+0002143679 00000 n 
+0002142752 00000 n 
+0002140106 00000 n 
+0002142868 00000 n 
+0002145389 00000 n 
+0002145106 00000 n 
+0002143781 00000 n 
+0002145222 00000 n 
 0002145330 00000 n 
-0002144005 00000 n 
-0002145446 00000 n 
-0002145554 00000 n 
-0002276565 00000 n 
-0002145715 00000 n 
-0002146612 00000 n 
-0002146638 00000 n 
-0002146664 00000 n 
-0002146690 00000 n 
-0002146716 00000 n 
-0002147375 00000 n 
-0002147503 00000 n 
-0002147619 00000 n 
-0002148007 00000 n 
-0002148123 00000 n 
-0002148319 00000 n 
-0002148871 00000 n 
-0002149263 00000 n 
-0002150238 00000 n 
-0002150716 00000 n 
+0002276341 00000 n 
+0002145491 00000 n 
+0002146388 00000 n 
+0002146414 00000 n 
+0002146440 00000 n 
+0002146466 00000 n 
+0002146492 00000 n 
+0002147151 00000 n 
+0002147279 00000 n 
+0002147395 00000 n 
+0002147783 00000 n 
+0002147899 00000 n 
+0002148095 00000 n 
+0002148647 00000 n 
+0002149039 00000 n 
+0002150014 00000 n 
+0002150492 00000 n 
+0002153020 00000 n 
 0002153244 00000 n 
-0002153468 00000 n 
-0002156582 00000 n 
-0002156843 00000 n 
-0002158574 00000 n 
-0002158795 00000 n 
-0002161333 00000 n 
-0002161598 00000 n 
-0002163332 00000 n 
-0002163554 00000 n 
-0002165359 00000 n 
-0002165633 00000 n 
-0002167016 00000 n 
-0002167244 00000 n 
-0002175943 00000 n 
-0002176219 00000 n 
-0002193510 00000 n 
-0002193981 00000 n 
-0002205167 00000 n 
-0002205546 00000 n 
-0002210823 00000 n 
-0002211099 00000 n 
-0002227203 00000 n 
-0002227653 00000 n 
-0002248008 00000 n 
-0002248669 00000 n 
-0002259819 00000 n 
-0002260138 00000 n 
-0002276115 00000 n 
-0002280209 00000 n 
-0002280307 00000 n 
-0002281799 00000 n 
-0002281995 00000 n 
-0002282179 00000 n 
-0002282399 00000 n 
-0002282628 00000 n 
-0002282856 00000 n 
-0002283079 00000 n 
-0002283304 00000 n 
-0002283542 00000 n 
-0002283763 00000 n 
-0002283992 00000 n 
-0002284221 00000 n 
-0002284444 00000 n 
-0002284651 00000 n 
-0002284866 00000 n 
-0002285091 00000 n 
-0002285318 00000 n 
-0002285507 00000 n 
-0002285684 00000 n 
-0002285860 00000 n 
-0002286035 00000 n 
-0002286204 00000 n 
-0002286405 00000 n 
-0002286602 00000 n 
-0002286715 00000 n 
-0002286841 00000 n 
-0002286967 00000 n 
+0002156358 00000 n 
+0002156619 00000 n 
+0002158350 00000 n 
+0002158571 00000 n 
+0002161109 00000 n 
+0002161374 00000 n 
+0002163108 00000 n 
+0002163330 00000 n 
+0002165135 00000 n 
+0002165409 00000 n 
+0002166792 00000 n 
+0002167020 00000 n 
+0002175719 00000 n 
+0002175995 00000 n 
+0002193286 00000 n 
+0002193757 00000 n 
+0002204943 00000 n 
+0002205322 00000 n 
+0002210599 00000 n 
+0002210875 00000 n 
+0002226979 00000 n 
+0002227429 00000 n 
+0002247784 00000 n 
+0002248445 00000 n 
+0002259595 00000 n 
+0002259914 00000 n 
+0002275891 00000 n 
+0002279985 00000 n 
+0002280083 00000 n 
+0002281575 00000 n 
+0002281771 00000 n 
+0002281955 00000 n 
+0002282175 00000 n 
+0002282404 00000 n 
+0002282632 00000 n 
+0002282855 00000 n 
+0002283080 00000 n 
+0002283318 00000 n 
+0002283539 00000 n 
+0002283768 00000 n 
+0002283997 00000 n 
+0002284220 00000 n 
+0002284427 00000 n 
+0002284642 00000 n 
+0002284867 00000 n 
+0002285094 00000 n 
+0002285283 00000 n 
+0002285460 00000 n 
+0002285636 00000 n 
+0002285811 00000 n 
+0002285980 00000 n 
+0002286181 00000 n 
+0002286378 00000 n 
+0002286491 00000 n 
+0002286617 00000 n 
+0002286743 00000 n 
+0002286863 00000 n 
+0002286983 00000 n 
 0002287087 00000 n 
-0002287207 00000 n 
-0002287311 00000 n 
-0002287351 00000 n 
-0002287565 00000 n 
+0002287127 00000 n 
+0002287341 00000 n 
 trailer
 << /Size 2801
 /Root 2799 0 R
 /Info 2800 0 R
-/ID [<ED7E23133C711AF833A7FAE4AE0C8700> <ED7E23133C711AF833A7FAE4AE0C8700>] >>
+/ID [<1DE3040442FE9CC3EC2E009607173BC0> <1DE3040442FE9CC3EC2E009607173BC0>] >>
 startxref
-2287897
+2287673
 %%EOF

Modified: seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex	2011-11-15 00:06:36 UTC (rev 19200)
@@ -859,10 +859,6 @@
 
 %------------------------------------------------------------------------------------------------%
 
-We thank Paul Cristini from Laboratoire de Mecanique et d'Acoustique of Marseille, France, for very carefully testing version 6.1
-of the package and helping us locate and fix several important bugs, as well
-as introducing support for the Gmsh mesher in the package.
-
 The Gauss-Lobatto-Legendre subroutines in \texttt{gll\_library.f90}
 are based in part on software libraries from the Massachusetts Institute
 of Technology, Department of Mechanical Engineering (Cambridge, Massachusetts, USA).

Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90	2011-11-15 00:06:36 UTC (rev 19200)
@@ -142,6 +142,15 @@
 !  simulation will start at t = - t0)
   double precision :: USER_T0
 
+!! DK DK for horizontal periodic conditions: detect common points between left and right edges
+  logical :: ADD_PERIODIC_CONDITIONS
+
+!! DK DK horizontal periodicity distance for periodic conditions
+  double precision :: PERIODIC_horiz_dist
+
+!! DK DK grid point detection tolerance for periodic conditions
+  double precision :: PERIODIC_DETECT_TOL
+
 contains
 
   subroutine read_parameter_file()
@@ -400,9 +409,17 @@
 
   ! boolean defining whether to use any absorbing boundaries
   call read_value_logical_p(any_abs, 'solver.absorbing_conditions')
-  if(err_occurred() /= 0) stop 'error reading parameter 51 in Par_file'
+  if(err_occurred() /= 0) stop 'error reading parameter 51a in Par_file'
 
+  call read_value_logical_p(ADD_PERIODIC_CONDITIONS, 'solver.ADD_PERIODIC_CONDITIONS')
+  if(err_occurred() /= 0) stop 'error reading parameter 51b in Par_file'
 
+  call read_value_double_precision_p(PERIODIC_horiz_dist, 'solver.PERIODIC_horiz_dist')
+  if(err_occurred() /= 0) stop 'error reading parameter 51c in Par_file'
+
+  call read_value_double_precision_p(PERIODIC_DETECT_TOL, 'solver.PERIODIC_DETECT_TOL')
+  if(err_occurred() /= 0) stop 'error reading parameter 51d in Par_file'
+
   !-----------------
   ! external mesh parameters
 

Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/save_databases.f90	2011-11-15 00:06:36 UTC (rev 19200)
@@ -160,6 +160,15 @@
     write(15,*) 'USER_T0'
     write(15,*) USER_T0
 
+    write(15,*) 'ADD_PERIODIC_CONDITIONS'
+    write(15,*) ADD_PERIODIC_CONDITIONS
+
+    write(15,*) 'PERIODIC_horiz_dist'
+    write(15,*) PERIODIC_horiz_dist
+
+    write(15,*) 'PERIODIC_DETECT_TOL'
+    write(15,*) PERIODIC_DETECT_TOL
+
     write(15,*) 'nt deltat'
     write(15,*) nt,deltat
     write(15,*) 'NSOURCES'

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_databases.f90	2011-11-15 00:06:36 UTC (rev 19200)
@@ -216,6 +216,15 @@
   read(IIN,"(a80)") datlin
   read(IIN,*) USER_T0
 
+  read(IIN,"(a80)") datlin
+  read(IIN,*) ADD_PERIODIC_CONDITIONS
+
+  read(IIN,"(a80)") datlin
+  read(IIN,*) PERIODIC_horiz_dist
+
+  read(IIN,"(a80)") datlin
+  read(IIN,*) PERIODIC_DETECT_TOL
+
   !---- check parameters read
   if (myrank == 0 .and. ipass == 1) then
     write(IOUT,200) npgeo,NDIM

Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2011-11-14 23:34:42 UTC (rev 19199)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90	2011-11-15 00:06:36 UTC (rev 19200)
@@ -785,38 +785,30 @@
   double precision :: distmin, dist_current, angleforce_recv
   double precision, dimension(:), allocatable :: dist_tangential_detection_curve
   double precision :: x_final_receiver_dummy, z_final_receiver_dummy
-!!!!!!!!!!
+
   double precision, dimension(:,:,:),allocatable:: rho_local,vp_local,vs_local
 !!!! hessian
   real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: rhorho_el_hessian_final1, rhorho_el_hessian_final2
   real(kind=CUSTOM_REAL), dimension(:), allocatable :: rhorho_el_hessian_temp1, rhorho_el_hessian_temp2
   real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: rhorho_ac_hessian_final1, rhorho_ac_hessian_final2
-!$$  real(kind=CUSTOM_REAL), dimension(:), allocatable :: weight_line_x, weight_line_z, weight_surface!,weight_jacobian
-!$$  !integer, dimension(:), allocatable :: weight_gll
-!$$  real(kind=CUSTOM_REAL) :: zmin_yang, zmax_yang, xmin_yang, xmax_yang
 
 ! to help locate elements with a negative Jacobian using OpenDX
   logical :: found_a_negative_jacobian
 
-!! DK DK Feb 2010 for periodic conditions: detect common points between left and right edges
-  logical, parameter :: ADD_PERIODIC_CONDITIONS = .false.
+!! DK DK the horizontal periodicity distance for periodic conditions
+  double precision :: PERIODIC_horiz_dist
 
-!! DK DK the periodic conditions below are currently specific to a Gmsh model designed by Paul Cristini
+!! DK DK grid point detection tolerance for periodic conditions
+  double precision :: PERIODIC_DETECT_TOL
 
-!! DK DK the horizontal periodicity distance is:
-  double precision, parameter :: PERIODIC_horiz_dist =   0.3597d0
+  integer :: NSPEC_PERIO
 
-!! 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, dimension(:), allocatable :: numperio_left
+  integer, dimension(:), allocatable :: numperio_right
 
-  integer, parameter :: NSPEC_PERIO = 670 / 2  ! 414 / 2
+  logical, dimension(:,:), allocatable :: codeabs_perio_left
+  logical, dimension(:,:), allocatable :: codeabs_perio_right
 
-  integer, dimension(NSPEC_PERIO) :: numperio_left
-  integer, dimension(NSPEC_PERIO) :: numperio_right
-
-  logical, dimension(4,NSPEC_PERIO) :: codeabs_perio_left
-  logical, dimension(4,NSPEC_PERIO) :: codeabs_perio_right
-
   integer :: idummy1, idummy2, idummy3, idummy4, idummy5, idummy6, idummy7, idummy8
   integer :: ispecperio, ispecperio2, ispec2, i2, j2
   integer :: iglob_target_to_replace, ispec3, i3, j3
@@ -884,9 +876,8 @@
 
 !>NOISE_TOMOGRAPHY
 
+!! DK DK for periodic conditions: detect common points between left and right edges
 
-!! DK DK Feb 2010 for periodic conditions: detect common points between left and right edges
-
 !***********************************************************************
 !
 !             i n i t i a l i z a t i o n    p h a s e
@@ -1385,7 +1376,7 @@
     print *,'Zmin,Zmax of the local mesh for proc ',myrank,' = ',minval(coord(2,:)),maxval(coord(2,:))
     print *
 
-!! DK DK Feb 2010 for periodic conditions: detect common points between left and right edges
+!! DK DK for periodic conditions: detect common points between left and right edges
 
     if(ADD_PERIODIC_CONDITIONS) then
 
@@ -1401,6 +1392,10 @@
 
       print *
       open(unit=123,file='DATA/Database00000_left_edge_only',status='old')
+      read(123,*) NSPEC_PERIO
+      read(123,*) PERIODIC_horiz_dist
+      allocate(numperio_left(NSPEC_PERIO))
+      allocate(codeabs_perio_left(4,NSPEC_PERIO))
       do ispecperio = 1,NSPEC_PERIO
       read(123,*) numperio_left(ispecperio), &
          codeabs_perio_left(IBOTTOM,ispecperio), &
@@ -1411,8 +1406,13 @@
       enddo
       close(123)
       print *,'read ',NSPEC_PERIO,' elements for left periodic edge'
+      print *,'horizontal periodicity distance is ',PERIODIC_horiz_dist
 
       open(unit=123,file='DATA/Database00000_right_edge_only',status='old')
+      read(123,*) NSPEC_PERIO
+      read(123,*) PERIODIC_horiz_dist
+      allocate(numperio_right(NSPEC_PERIO))
+      allocate(codeabs_perio_right(4,NSPEC_PERIO))
       do ispecperio = 1,NSPEC_PERIO
       read(123,*) numperio_right(ispecperio), &
          codeabs_perio_right(IBOTTOM,ispecperio), &
@@ -1423,6 +1423,7 @@
       enddo
       close(123)
       print *,'read ',NSPEC_PERIO,' elements for right periodic edge'
+      print *,'horizontal periodicity distance is ',PERIODIC_horiz_dist
       print *
 
       print *,'because of periodic conditions, values computed by checkgrid() are not reliable'
@@ -1516,7 +1517,7 @@
 
     endif ! of if(ADD_PERIODIC_CONDITIONS)
 
-!! DK DK Feb 2010 end of periodic conditions: detect common points between left and right edges
+!! DK DK end of periodic conditions: detect common points between left and right edges
 
     ! reduces cache misses
     call get_global(nspec_outer,nspec,nglob,ibool)

More information about the CIG-COMMITS mailing list