[cig-commits] r21173 - in seismo/2D/SPECFEM2D/trunk: doc/USER_MANUAL src/specfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sat Dec 15 16:53:58 PST 2012
Author: dkomati1
Date: 2012-12-15 16:53:57 -0800 (Sat, 15 Dec 2012)
New Revision: 21173
Modified:
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/specfem2D/define_external_model.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90
seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
Log:
better implementation of define_external_model(): now call the routine only once for the whole mesh instead of once for each mesh point;
this is much faster for large meshes and avoids fragmenting memory by calling the routine millions of times when there are a few million grid points
Modified: seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf 2012-12-15 16:09:50 UTC (rev 21172)
+++ seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.pdf 2012-12-16 00:53:57 UTC (rev 21173)
@@ -5069,12 +5069,13 @@
/ProcSet [ /PDF ]
>> endobj
221 0 obj <<
-/Length 270
+/Length 271
/Filter /FlateDecode
>>
stream
-xÚuQßOƒ0~篸GšH×»¶B ç3“‘e€/ÎtÌ($l3ñ¿·¥,fš=}×»û~´ð qïÊ`–) sB ÜiÅ¥-•<‘Ê <‡Å:]déŠîÙKù’k!7˜øyuh‘áªîNõ‡[³Â
-Ð.iMN8’ŠE)i…µ'nI£'®‡¶{kŽ}çU×2ÒáC 6á=~ûþVhQs ¯îvc¨ˆˆ cs¡´6S®E¾)üÚÌÃ2ß,çkFý~Êåõé×=jt à|ã˜'·lŸ"ŽjçëÄ\L)'e!á4½mÞ;ïÿùz~G47¾ ä¨ÿþ.-ƒ ”£pN
+xÚu‘Mo‚@†ïüŠ9BRÖÙ Bcè¥z@ņ¤ÅµIÿ}ÓØÆÓÌÎÇû>»+á $<{òN|ʽQª`$YC~ b-”Mµ‘"V
+ò=¼ùÙ2™¤Éœ¦Á&%Ø@H( Æ®_œª6‰¥?/›KùÑYa
+h‡˜©•F„ZYav‹kbt‹Ë¶nvÕùظcÑÔ±ÿ ûU{ªÏß®¾–,‹llðW6û *$ÒX.TÖfàš,V™¹0[¬fã›Ê=£ãaàêåå×=mK Út¾Q$âÇ?¶¯"öj×ëDB:J4"b ÞvZíªÏíõý \B©[ù÷gIîý !–nÃ
endstream
endobj
220 0 obj <<
@@ -10904,31 +10905,31 @@
/ProcSet [ /PDF /Text ]
>> endobj
662 0 obj <<
-/Length 2187
+/Length 2284
/Filter /FlateDecode
>>
stream
-xÚ¥Xmoã¸þž_!Ü'»X3")QÒý°½½)îÚ´ëE?ÜbÓ6±z;½l6ýõá²%ëri‹ M‘ÃáðáóÌ( ŽAüù&tÏ?monï"ˆ)%â`{(3ÉÓ`»~YýpÊ›^·ëˆÃU´þ×ö/4#bIšpœ›(c)hÂ?†ª2Õ‘&ô'MuñuÍÕ·ôkôî KñÁ™Œ‚ŒeJ(g1æ,
-c²¸]§áª¦yíPÍ,wÖrŒ–£põν}iô÷hyܪÞÞ©äb)ØG&Á†§,ŽÝ–Ùí“©n_õp3™—м |”¡óÖGl”¹q]ÞÚߟÃ8 :MKL<ã B 8m—A»×®È”¥qì‡Õήþu0k¯¾â¿¼ÐUOýæ@Ï—z Ʈ˅ÞÓ¯gÓŸ¨Õäm^ºpnMS·=xËÙz©pµ=™ÎÏ)ܨz蛡Ÿí½Ó¦Û ߬8d Ǩq–ÅîtËúØæ%ØJ¬òjO®Ê›îT÷®¿>ГN Ϧ|-Â…“zm]ù)yOÏ=lbpÐ
-‚µbJ7¡ëuãÆã ÑøVïúº}Y8 Ž‡ÒÇúoŸ¶Ÿ¶w÷?ýøñvád’„%aêGCÜd¤Î`†¥¾šnÈ óo=n®|G‘’iÊ™N#… †‘‚¯¾[ðMHAlÝjGØ•‚èO|ü
-³ö9Ï ¶Ì³$r÷ööM̲$
-.ºÄçzk<a™ ±ù y…øOܹI§=u×w»Ö4=½ø †¢Ð ½Ä+`ƒ/žOfw¢n7™/ µ
-]ìiˆq¦ZÝ´ºƒsÖîáÙ:Uœ •bi’Mƒ›·m8yîB¢ƒÙßÂ5h
-ήÌû݉¸
-~X¦=Xçvv/ÐçÖƒVN 6ô…š©\«0•v«àqüSW·O£q½?¾Åæ±Õº"Khï§êí.A%Š"f#TLÙP`\?ÅE-€Ü”ùQ‚„ÅYB’s%ÁE7;Â
-¼†à5ãÅA±Gò ˆ¦
-¤6ûÒQƒûu¢c£È!ÙÑ£1ß,ZðzíÉyoßZ>Ág5ëp ÂÊvµ_½l`(r !H&,J£)‚¬‰KeZ n| ܘ°_;IDèvC«„íu£«ý¨¤µÓ½ç“%¶ØQ86NWktº€«®Ý"¿Âè™'jl›j‰èdÆD4ªÊ‡÷Û÷·yûx ©ø-¦E*q)¾?noð¨Ã€Hl1‘)&cìÊ›_oB&âò qÙ¶/ce»ƒ(;oïK|¨oþÞ.‹2?é°˜,[äÿ‡Étnx .†T?œ÷OIØ$œ"IÄ"éugŽ€‹õ&Q|1þý (åêŸ(%g?" Ç]hâ»-°±×D]CÑÓ6¾åeS¸G]é6·(ŸO‹‚(9S#L^M™TV~èó y íb¢QjL1¯Íod–±0T³‹fÑûh)áq¤¤¹BÁŒXŒc“•,‚›ÑS¨fa; rÔPÖ·ƒfK)W ™™ºö$£°£ñ]«!pnE ›îeÒÞ÷nxêl˜÷D×8Èî¹H©ZÉ+×8åÕÑ™{5\<˜3^:ä˜Î<ŽËºÙôàôàKŒ3ÆEäwi“ÛD ½âãdŽ'»)h» sÐôÓ¦ºð$‘_ÊlC M:
- ÉÐÜîêªës`XvZr&„’DŒ!·ê†k@ªFù+ý¤¤Ï§š–Úsɧª¶«÷Ž
-‡)vA#aÁd¤´2ÿâ"沈™QÂÒ3Ø0Õã>EÆEìIbãi\U †b$à);ïòÎó´“W2¼Ð¯"o×<] Ý öRtSÂ
-6™dßÈ
-áöüPèÑ0FÛìõYü‹C]—‚ŠºTC¦Ó½óò£«±[¨óÎõÍ´ü
-8©Mh/
-¦Ï ÓQ¿lbW;,&¸@{o•¤8úœáP)„~y¥¤jö$ýí—î¬Ø•—åQuÇ3 |ìL¶]{y…Q°ìËxTíòÞxAï]eÛë®_òhz¬<” ¸ (Ðü÷ã†2[ƒù:ë:° Þéo5^Ü—k…dR+tº´olnǯ§ÌÏ ]È[ï°! Âöû?ÜßSó\7Àw„P‹žs5èYʪ‘”"õ¥ ôÓ5ËÎAAc-ÅÜÍrÉ]³x¬RÕãØóÇL?x¨\]Þ©Ú¿{[
-é—>:ƒáŠÆrp?ÇèD}]3,æ‹ý¸òn–ðË©;g©gò¢¨†¦¨û%†ã,ä9_ŸÂ2‚×Jð·Á¿Að,#÷y–’(cÖ\ù+(nÏZ”w9VúÈ`•–—TC@‘&vm9c‰È~_ñ!k=ß‘Å”ybè¿K˜ñØB‡LŒ1UQïLÿB±ðÃW{
-Äí{Ç}ƒÊ¯F\Ì>PFT Ãs΄ÐEC¤C0th uhÔ‡¥mf¶î{ë>!*2»LÂÐpò‹Ã1ÂípëåãÂ蟓‘çÖô=
-—˜ƒzÎW›´ ¨g»‰Y ;ºåàš”jþ-
-êîêUµ×¤cw uÖá•*ŠËÈ–‰nûrC>ƒ-;dáä’Ì¿4òH€ƒžñÅëWËɳJÜw=h•Žd™¬î«Þ¢GðAÆD*K,¹ \q¿‘
-\À¡t ÑÓ lGï–î_3Mj|°_7NûÀ+«}è¬Ë—±ýdªÜf:Ð~øëÏ4ÐçÇع7]Sä/.mÃO±Š).çÐS>e1¸B*{³´zVÏdžc_mó§ëD¾€užÒ Ëß¡[p—ƒrÀI7˜×>KŠX±$Ë°OÛûŸ?HÚº%>g©³Ý‡,įb"fe\Ϋøÿ ”K/
+xÚ¥XÛŽÛÈ}Ÿ¯ öI
+¬ ö…7y˜¬íÅkgËÈ"ë@ Ä–Ô0Eryñxòõ©êjR¤DÏÎ& 6ûR]]}ªê}ïàùÞO7¾{þe}sûNIOø, Eà÷^$¼PLòØ[gÞ¯‹ iÕêz¹¿P˯ÿJ+‹âˆã
+ß[©„Å\Ñ‚tEaŠ-hš Ëüë’û‹^Ò·¦Ò»½>‰7N¤ò–„"tÎb? ‰ëeì/JZWwÅ…äÆJP²ò¯ÜèS¥_£äá¨ýóö]¶‚s$~äxÄ"?¡ýØíÖ·Ïj¸š¬s¦jAGé; ±auÄÆ)5®k¤}ÿì~×hÚb¢ç`‚P€ÒVø©2(÷Z³8úi¥“«ëÌR‹¯ø—æºh©ßìéùTvÔØ• 9×½=šöH*Ó<×¹S»«ª²nA[Ζ+ú‹õÑ4ýšÜÍ*»¶êÚ‹³7Ú4+ÔŒø,âh5Î’ÀÝî©<Ôé d‰P,Ò"£FS¤Us,[×_îéI·GSºþ‚ÌI½ûº,ú%iKÏ Ž 6ØëÚÁJ1'· iuåæã
+ÑüZïÚ²~š¹ Šû²·õß> >7ïî~ûñvæf"†›
+v“*<ƒ¶újš.ÍÍôp¸Ó+²”ŒcÉxj)5Ì|ñÃŒnB2¶u» àT!X¢ãW8˜•Ïy‚@€x)ç·²#K¢Øu3°ÏõÑ ö‰ °ùi…øoÜ©I·=eÓ6»ÚT-
+|ö}‘ë†Ñ¬ñaàñhvGêv‹ùÜUÛU¾ÐyFSŒUëªÖ
+ܳv„g«Ô àì`†,ŽgÜÎÖu‰8yl0„¨7ÉnÁc
+Jƒ»;¥íîH±
+Þ lSïr;{èsûA+¥G úBÍÊ®•›B»]ð8Û”õv®³ÃKd j ’„BÑ?ÝÔ^î T”R ‚]sªÈ0®º €Ì)=hB ± ‰AòA‘7êfðÀk^>APàô
+{¦r mvÐ…÷v¤k#Ëa°£Ge¾Y´ {e$Œ´·£6žà³¸èpÂÎve¿û©‚©[A2b*VSKYãÌ4Ü` Ü€°_º”ˆÐmºZ»–éJÙIK—÷ 6°.„cãxµG£spuí6Ù÷; šõÛ¦˜ t2aB
+YåÍÝúîö!7{Hß‹t (KŒ“ïÛõ
+^µïqXL$1SQàíN7¿ÝøLðœ1nÛÁ ´ÝžJ†ÎÛûSè½)oþ¿^.W!SI¿èÍKäÿ‡ÈøR$Äup 8ŸŸHØÄœÂQ:"ƒ¾3ÀÅr…|¡½ J¹ø'¦’ÑÝÈÂy£¹ŠÝØÈ4…®.oi†FßÒS•»]è:µ(Â×ílB”œ…Lž¥LaÄ$[7õñˆq å"Ñ8i¤˜×âW2I˜ï‡ŽfÑ»±!a3„¤Ë …©: ²¯%+*Ï@~¢ 53Ç 1 9”µu§Ù å‚ia8Jì¡ ‰ …ïj
+†s;bâ±t^Ðß›nÛX3g®q=…S‘¨JI êÙ Óâ ©3+»mîDCÄè™U,€ó©v¦16l$âY+‡ SÁ`©}Šìg3hé,lÙÄŸé–åàñ`Ëiã1(p4‡£%ÿÐv¤´Ó4Ã8 ‰Ì(˜„,,i—’‘£«¸›µ*n·Õ.ëqÅdtÁJ“ï;¤µ˜7
+z¥…Æ.m\ Ý
+ŸŽÓ?Ñ[žÖK /앉n¦ó d Š¸¥EwÚZXŒÄ I«™ö;»»k¥SdLØ"æj¯ÜnÛ5-
+lÝ’3ÌÔF~¬¶Æ')J0
+¦ýðËæÿ6÷ïï~z»y÷ËÌU$‚ù|Hò™P…MK’ dðýk€¦#ÆÇ·ëO·»²hÚ1;2SÌú=Q[·
+“©Léa ˆëw´Ð”a ;᯴†¬çRÒ… ÇΊsh_NÜw hÄ#‡ K ±ÙÖƒ3Hðm|íA2N›x5 Ûç$û!׃P4„Éô9³÷û2ÏÇ”™E \µyÕ] µôL¥~®P§´Ç! Ú’dT*v4õë*pÕßl‰‰ë¥~ C89sT*fQ¯ž- «¾Ž¨WYiÎœ«è‰ÕÀ›(@£ëmgÚ³ì²'H0 ¶}®É ´]JŽf÷wß&ZÝ´sMì³â¾d½1fCÈŽ~ߎ i«è¾R¾6,„<ÎÅwëÒžŸÍ“õxBÖ}²#DϯV\Þ\ÑhJZ÷êâؾûøãý=5GuŸ j39æÚ
+™~w*š×¤Ã¾T.KàäÞ$(¬¶‡ZÕ'CÅÔ9§zo3m×Å\}NÒàYÙ«¡·Ö@ŸõyÒ
+|¶nx(ç³K„N Ô–e~ÄÜ|± Ç^]vcƒÝÎÉŸ+åd™ ÕrEWåe;sSŠ³XÉs½5¥‚áPð—¿!ñ$!õyâHö`Í ¡v†bH¦„êwO›Æ ŠþY¦Hó}s±\@‘éG#F1Ã
+X$’ßglPuÄÃégKž‰ ?Vðàµù aíAq*/w¦}"ÃXøáP¦!l÷½Ã¹1—@=ÀÅÅ ªõ|E 8ày¡‹¦ô~ ÚB å~ÍY/='XE94ÊmýÞp‹à n;r1ФDõ\y¬MÛöÓ·Oƒ §ª9Ô35ÝH¬LG>šIyõ)íé\õªLîéÇÐ]C™¼¦æRÙ* /{Sè Ô²}âO|äòC1Δ/ ÷.7‡‘û,KÙ9ŒûP÷EkÏMq¢EìEÂ¥0èœ < †¯ ¶ÿëNŽÒoŸà8z7ç~IÀB5ùDòËÊ%>ÐÊ&>TÖ•;ØÞš"µÚ ÞÓÄ=Agfš*OŸðhö&”Ù!——ÈÑ7”‰ #6  &zÊÄå凈ÿž-d¾
endstream
endobj
661 0 obj <<
@@ -14375,7 +14376,7 @@
/D [661 0 R /XYZ 72 720 null]
>> endobj
664 0 obj <<
-/D [661 0 R /XYZ 134.097 286.39 null]
+/D [661 0 R /XYZ 134.097 288.512 null]
>> endobj
660 0 obj <<
/Font << /F43 222 0 R /F44 223 0 R /F67 376 0 R >>
@@ -14383,31 +14384,33 @@
/ProcSet [ /PDF /Text /ImageC ]
>> endobj
677 0 obj <<
-/Length 3250
+/Length 3303
/Filter /FlateDecode
>>
stream
-xÚZY“Û6~÷¯˜òËR©‡ ooíƒcsìÚqÙònªœ”CAÊ©d4ò¯ß Š#O’}ˆ³Ñç×
-WÛ«àê»'ß®žÜ¼J²«ÜÏš\6W)½JÃìjµ¾úè½øþùÛÕí»Å’Æù‹eœÞ»Ex Þ¼ùáÍw8°úþïúׂ÷o»äç÷oo_¼º}M_.~]ýxó*ŠFEqæa
-dè³H¤æ< aö÷jEÄóüjISXâì_(ùW —1 -6‹edÞ©îT#õ~ Z±ónA¯ÛJ;uƒ#¸,µË2oÇÔÜß$ö8õ~Ï«Vâ¸](»í–ËV,ˆWW¶$Ø(9“ÜÌáH@ê±
-{øž‰»Ú»
-±]J.ä¾Æþç Bˆ·ÂÁ-¯×§Ší—êÊW™ %ÀBü<Ž‘ ¢>D`*rž©#U£Ýql¬ù‚ÆæFe}àÄþ£(K3µ9™F¿Šw÷ý6ûÉÎx-ÕÒ\Ã1üíkÖr3ŠG6RÔ•+±¥½Hî§áWäÚ ¦@ æŽ`þ¼R'ÕMÛØ9E½62èÙ\ï]k;w 籊•§V¬4«Ë®Uôê¯M/l±¹â’=_ f·a-¶Žf°³¢gÒ!°áÞðª0ÓDõ<ï`ÒnϚϰžäš©÷p3¸ã ¿êùíl°.Qr³HQ¥~•¢£ý%éÈþH f™µ¿]Û žÝÜ G_ DÁ—`í×´(Öö¦w
-kN7²Þ´GÖð›ÛŸéË—?¼›±kš‡~dwUáë[‚ Ç~ç9]iÞ)5AY©–Vceш;«G¢zl#÷- k:7ëYécã#ÉóèW<“F~ ¤ZÅâˆØ#ù 7À¹dÔ’ýþ§ï^ÜÎ\.Ìü0N‡Û´4T—u1輨ë;?+%~˜»ÉËç«ç7sGQ?Ls;k-^´µ5Q¹«»Ò vg( œ%@76uYÖÊ Ee8wgí±Ó3ÇÞ€ˆpDÄ’$O‚Üe¡T"(ø'Ù5›éã[Ä‘ŸqÏW%0R6£¤ FÚA‡lkzën%¡Ÿ¤½BùmÓq.œ?K{¦j{´íjkƒƒÐÚ£LSßwšywf¬—¢d=±f¶¾¬bZa¦+ töÞƒ—×ÐN©WCos’»h¥$aìSà²c–(×%[Ñ J8²ÖTzñ*wÔí/Jrd ½»½——dGB? {1¦P³I;HõuïêÜž·*º\Ô*ÜÙ¡åË×h¡ åËŸ!…N(©w{:ð‹ê û)!ŽzÓ$6êM“¢ Ø +;ŽŸü·N=gNÙsâ¤Wv2£ç™ŸÐÐNСIï[2Ù~k½7$à ¢¨;€*fÆ¡áô™_«Ï ¢ùé³’<Ð@ Ia –¡à“Ù™WÒRŽ²”@»‰aù$I]ó0¡]I5TÑTý€QT ›èÊÆMõŠª(»µñw!D˜JAWƒGŒ« - ¡•//kG¦»J³ Å_c|¡·i¸i „âk yo/øü!×øûþs©Ù÷ ;•\lwšgˬó
-HNÕ$½ðï°@K-ÕÏ°ßpSõ̯?Š‚N‹ 3¼T›TvS‡—s{J\ÀÌ~š—ª¼T /ƒÈ>W1>'Ãt18êÅÝæ ¯«ò´
-4»¼uX–<̲ Y™^I•|Ô珱¦Aט¸``p0Ï ßÌxÛ52À~wô Xô[›1Nš ø`mOT¢¬ÜPáÜNS¬4&Ϲ#?‰"¨¼eͧ(ùÜ~©ŸÙò{ÐóFIêiâúÞVìù§MWʧÕýf©é¡=òÅ®®1Új#b¼ËRíŠývkÀ’hŒD at Bz@BÇ [}ÎzA:'½à§sˆ2¨É à4†ßïDñY©Œézü¤~e—¼:g€=;CµB¦4w[óF`ò§nÖŠ!œ;2U@¯‡¸á׎dÝ2ž:íïX'¥`ÕÜ)*©zÕŠfÃH”dÃ):ž¿
-+°YÏ@˜™ö¤ÅÜõ{Ž¶(`åèÆÏFo.«4¡±Ÿ$8±
-l §‰UÇ4AÝIc#ZèX×{Q±ÊLù[Ù¢I U—¥ãxpLI–‘T jQµË «¹bß•P™ó;©¨ÚaÏ„]à È.¥ hƒå Mßã`ÆÖõØÓø©ò
-F«„ÅôBþ£K0küËm‘2øÀ«’`ÀöÐ9e
-t]rž!ØZ%=p f´‡Äã’•sš€Ïgø¹³‘yf ’Ùì÷.æƒÅq'Š 43X f‡ºZKìÐÅøµ‚Ñh&ÓNw+ÔÇæÃÑÚlðvæ @MÓ,GT¶lX;YaÇô ZB:ÇB–Ó²9ìá؉ ŒÁƒEO½w°@SÈÿ»‘ýŒCR:¸_¥ðº0¶]·–PßµLT6…‚X:CFùq6Ô#of5€†~> ~ ¥þÃf’Ž
-ÅŒ!—ÓÎD¦õŽ_q8¦Ap¯¹ Z ²!µüiø«‘XÈ0bd?â“ØV¬´d¶;wpÇìV y¢d[4Óþ+÷µÊsmµ± mKþ7{’a?W6xh)~MÔ矋¼Ó äÞÊb7©8®Ô`Sï¦ÐQ˜þ:Jî°ÊQ¹‰à+±µUlj|‹#°‡Ùoe'{+-] ÅÚrÞÌy”<*aZ¯Ñ%ènçÚT=@'úöŸEx¶4 ‹Š®pš®:»_i°úHÞ36 ¹W …Á¥ÇYZ¾ð¦žÈ¨O ŠXº&ô‚ìJíìPýªÎ›”LNÌÇ„V×£Uõ| ^º²?Ô²] š²Lٿئí¹Ö ÚU׶îR°Þ^MF<TÌÌöëb8““úùÔâÑí7‚·§±ÖLmàA]”FüÖHÖ*òílûâkÏ‚ª–è-¯xà ®(HaìT§—éF,7ëZ~8wÔYÞÞ²Û–Â$Ûðò4<P”p'K£NO{颴ëó™´ÓºƒÝÍ)d7œ>¸E[½›äk.Eücîrº©[ºÒ®jy`RºžYSÓ0lݶŽ3çf<‘\ïw°3þ ü¤hëãô&Ÿ:yè§iïΟ“TiþÐÎœ '>IzwpñèBîfŽ^Ò>¢‰±oêfj‡®Óæ×}>`Þ€•“åXuéí#ß™OШ4Ö¹²â3ÛZß/»ÞÚLÏûçJo^Ì)þ¯.¡ ’ c¨Ç–²¬J[ÿ~#!³™Éµ ¡ˆ¨›TkÅX©0M‡ Ã9t ©y y²“–ìÙgŽÎý’f×$“$Gû¼ Ø„~I ,ÇÞµéBM±
-Méq4ýL'ä§CÙb‘i¦%€å⋵
-øMþ‚¢*Éc'É×· @j¨ÜdÆÊ'€ã74˜°ÖÿÃEüöñù娌ª¬@Aµ±ë’Æ–Ž}îS.0Ê>®
-1~wæíäNlZk<*UÐÉPZjg}PäGÃ[9Çøªgÿ¶ À/;Ùl@}¤ÆYÆ³ó ¨Ÿ&éd‰ÁSä¿z$ø ßíêý¾¬3ÕÈcÓ Í¿"¯VÂ3(Ž•%všÊ´Ì¹j¾vGö; ’ÄoåœX¹æÚ¤Ð3kfwÖµ…ñ¬¡Ð¤&j-²¹;I¡®¿³0nâ!
-üC|Ô§a„A :Öµ®ìC† ¡-¶Cc‚æ” ‚ä £‚šˆ‚{²õ¢‘ æΫêji€*ô Ñ5A¿E/ÐÜ+<ªŸÿõÞh/ÜXuz4<0}åm‰$£ŒÆšá»?M¿ÌFÀ#ÅÂ
-Ë ñ
-ô™WP½¶*OJ@ê ¼éLû¦’©•oÈRa.–¬P Xï¶GiŠ•|ò|äx`^ÿ¤ãSDàÓ8™Ãò,ö»–e»»
-büúKåõýýõë/_®_ß¹$‘$?Ü¿¶×ÈÓáE+×ÿ€‰ÿRßJ"¦â#ÖØ¥Ÿ$tkcWõLÉg·Á3몗‰N˜çÃðSúÔ€‚(&Þ›º5ýƸòtxwŽˆ e“ij«D+í?š³*Æè"q0³ôÇÁ@ñ÷Hj-6Šúÿ/›Q s3™Éxøòâ?W¢/ršùA:~qœ>²?pá)ŸúKÛ
-קŠS\~ž7¤¨zYeÎ%è_¸ý¿]âÍ7—/S ÝÒ7SÿY¹D|è 0P¨2~Ó7¯€ÓìïíêÉÿ <‡¡d
+xÚZY“Û6~÷¯Pùe¥ÔC€··öÁñ‘x7v\ y7UIÊ‹‘ eŠT2ù×o
+€ÇpäI²Oq6úüº¡`¶›³ïž|»zrý:Éf9É–ÌVÛYÊfi˜ÍV›ÙÏóß?¿zõa±dq0Èb'ÁüÃ"
+æ ß½{óî; X}ÿ
+7?þ° A0ÿ·[òÓÍûW/^¿zË^.~]ýóúuõŠâŒa
+d˜³h¤ç< ,aîw¶Œ"Jâ<Ÿ-Y
+kCœý cñ`þìçe G7{Q <œûþlÅ‚Åó~´JlÛÛj]Ëc£ìärƒ× Ï«º©y‰ uÕ6²vž´½Y‹uSÕg¼a’önÈℤyìnøqõ懛ë N„1É? ïµd9#I2[Òœ$4¹|}¹],à ›Ÿ«V7Òù/AÀôtçí‚ÁÅwÊÄAUã.Kݲl¾çzîÖÕá ÊFá¸[¨ÚÝN¨F.è¼*lI±ÿX®„#€t®9©{ÄË»š
+»Ör·TBªC…ýÏ”Òù
+w¢ÚœK~Xê+Ï2%
+”’<Ž‘ rDaêHŠ VNFúx£¢:ŠÚÊñ$ «
+ÈþjÞÝùm£ñZºe¸†cV¹ŽÞ;ŠGÖJVåH°î"Z®ée¹zÁ¬‘‡ù@° «ç¬«•gsu8¶ëÜ œÇK^œ¹æ…•XU€ªWVZ[/ly ¹â’ƒØHî¶á
+¶Nv°u¢çj@`-¶`åÚN“åCò¼…Iû¯?Ãzš£±éßÜ îxƯjkÛ¼†KÂ.ÒTé_èÆIc80ËœÕí›æøìúút:u”k±g6ß” X»ëBÞÖ¼>_«jÛœÀ±\¿ú‰½|ùæÄ1³<$‘ÛU[„±gíÆbÂyƒ‹® ï´š ¬ŒgQNcµcºuz$ËgÀ6Ì¿}Mv=/6~¦y ýŠg²ˆäAjT øïŽSÜ ßš1GöÍ ?¼x5å©2Æiw»€–ê¢Zw:ïœc3yVJI˜S·ÉËç«ç“N‘‘0ÍݬÎÍ¢ñí«¶°‡ÝZ
+@#' @ç_Eeü¿,-çN«æüül`o at DØ#bI“€Ð ·,´Vª´ÆÖâ“jëíxƒþ-âˆ$wí7ZZa¤mFK9ŒŒƒùÎöVÜJB’¤^¡HS·‚LESJ²Ô3ÕØ{dl×X „Ö e†zßiçÝÚ1/E=È=±v¶¹¬fÚÚN×8Øû Z\\A;eóJÇâ“T¢ÛÅ(%…pÇ€ËCpD
+]²í±€#½h•t¦âÅ«ÝQ{¸(Éž-xw{§.ÉŽ†$ ½ƒèShØd ¤þºêÜA4:º\Ô*ÜyÒ¿|öR¾üRX6äJO½›óQ\T”Ô›%±Uo–$ýÀæxÑ
+ü¿µ:Ș9SÊžSP¯ìtBÏ3’°ÐM0¡Éì[pÕ ‘kü6zoIÀ9@źjªØÇZ(Ðgq¥?sˆæ÷HŸ$æ1øO {$…1X††OvgQ*G9ÊRí6N„aDh’ÍÆv-Õñl¨¢ØDWÖ‡lºW–ë¢ÝXB„)5r·xĺÊÐÑYû"ð²nd¼«² üµÆηµ°žƒÐ|£ùûKþèþÞ<b.³û~àçBÈÝÞð,2ëþ4gzš‡óðï°@K#ÝϱßrS÷L¯?É‚Nƒ ¼Ô›”nÓ/§öT¸€Ûý /uy©[#^Ñ\¡c|N»é bp†eƒ»M^•ÅyhwyÿX–È2HüÌJ¦å£?/Î
+…ºúÄ ƒƒifØ„M7š¶.‘q î,ô»µ ïÈknÜÚT¾•Q|p¶'KÙH^l¥ ¨pßNS¬,ê&O¹#’DÑ ¨¼çõ§,ÄÔ~)Iƒì ò{ÐóFIJ‚4úÞF ħm[®µOÿªûÍRÙ±;òžª0Ú#j½ËRïŠýnkÀ’hŒB at B= a}?'½ “’_ðÓ —²!a6/€Ó8~ëÏZ@˜`LWSà'%Y”]ò¶úœöì-Õ™ÖÊÞm#j‰ÉŸ¾Y#»DpêÈ4'qä!nøµ#¹E÷€ŒÇNû;Þ*%y9uŠ®i^µ¢É0%YwŠ‰ç/eÍ×ج¦ ÌÌ<iñ#wý^ -ÊX
+ôÆàg«7—Uš²˜$ÉØlƒGÄÇiâÔ1MPwÒØŠ:6ÕA–¼´Á_þÖB¶hÓBÝe at i^ S’¥]$Õ#ÇJ–ÍÒb=W Ú *{~«L Õr챫;Ù¥Õ´ aN ,]ú ¶nÆžÆOµ÷Ð0Z',¶ò S‚Ùà'Xnƒ”Á^• ¶‡Î1k ë’ó ÁÖÂ(ñÀ%˜Ð ÷+vÌ*<ŸãçÞEæe˜Jf#°ï]Ì‹Ó^®÷ÐÌt`˜ «r£°Ãcà× Æ ™Ì88ÓQŸ8ÚÖf·³çè j›v9¢²eÍ›Ñ
+{8¦çÐ’jp,d9
+Ÿ
+¤;1ƒ€³äñ©÷ hj ù[#²Ÿ°bHJ;÷«¾Wçìg®–PÝ6\–.…º*VÒ "qÖ•c¯'5€…$ï ?A©¿ÙŽÒ±®˜ÑårÆù÷ÈtÞñ+Ç2XîuT%zȆôò§áS¬Fb!Àˆ•}OrWò‘Ù쇃{î¶Í“¯‡µA;í¿êPé<×Uk¸Ð®s'Yö mƒw€€:Øâ×H}þµÁ; @îZïGÇ• ¬«ÃÑz
+³_ÇUÉ=V9Êa"øZî\ŸYß2ØÃìw²Sƒ{k-]žäÆqÚNy”<é*aF¯Ñ%˜îÁµ/¨z€ôí?‹,˜»Ò€)*…S·å½û«÷ä=a{åQølQ¹zœ£å‹¨«‘Œ|ò QÄÑ5¢|`[g‡êŸ0}Þ¨dr.a>&´¦ «ç;pð£Šc¥šå±® ËT¾øÅ·çš?´-¯\Ýeͽ½ÚŒ¸«˜Ù훡‹\êçc‹G·_KÑœûZãm éW7&tQYñÿY#ÙèÈwèL°ñÅWÏ‚²’Žè(EÍ ®(vH¡ïTÇÝöXn×5âxßQ?dyÇnW
+S|+Šs÷@QÀ &=õÒEiÕÎçs妵G·;ŸRȶ;½s‹®z7ºÞF(éÌ}@N×UíJWÆU-\©á…€gVÖ ÕØ4 ;·mâÌ}3 ×û ìL|B?iÚ|œÞæ“@'IšzGpÿ9 ß 'ÎŽBï. ½Vû‰£—,‚hdìÛª [ã±Á´Å•Ïì; ð¯ d9N]¼=bä»çÌã*s®|ý™ïœïW·6Ûsó\ëÍ‹)ÍxàÕ%$!õC?¶EY]õûµ‚Ìf"ׂ„"bäÚ(ÆJÇ€qzÜåƒC—š'2íÿ,й_Ò¤€àšt”䟗à›4o"‰…åØ»±]¨I8Vâ -=ö¦_‚é4ƒü´+[ì!2MÁ´°\|±¶A±ä/(ª¶‘< $ùæ ÈAM •Û̘ÁAyxÁo0á",ùÃEüæñùe¯Œª@Cµ¾ëRÖ–N>
+û” ŒÒǵ.Æïïy;µ—ÛÆÏ…JtÚ•–šI ‘¨›â*ç_Íìþ¿6 ¸@âå&Ûí ¨÷€T?ËxvFÒ$ m0Db°Áù¯ É>˜w;Ät hôû²ÉT£9?Úfhÿ5Í+-<‹âxQ`§lAËž«Wàkwä ¸£*)üÖΉ¯!×ÜØzbÍäΦ¶ÐŸÕšôD£E.w§)tÀõ÷Æ<¤… aˆú,Œ0ÈCǦ2•}è1°#tÅvhŒÐœ–Cƒ LTÐQ¡'Û,ê bê¼²*—¨Bï ]ô;ô̓ƣæùßìöòÇÕ¤GÝÓWÞ–hÒ Áh¬¾û³Ôùˬ<R,¬`°Ìß@Ÿ}5kË⬤_òh>žéÞT²ûEjí²ÔÅ_˜‹%çtÖ»ÝQ†b-Ÿ<ï9 ˜çŸ‚ôb|
+Bƒ‹“©wQu/ö-Ëu·%ÄøÍ Š^oïî®Þ~ùrõöîË%‰$Iÿáþ»Fžv/Z¹ù§ Lü¯þÖ± ¹Á.ó$aZ[·Ê3%O»užÙT½ltÂ< †Ÿ²§D1¿«Ûo+O»w爒(‹†¹F[J÷ǯ ÂòUŒÞEâ`\gñÇAGñ{$µ‘[Íóÿ—f˜Qtÿ‘f2£ñðåÅ*C‰>¼h ÐŒiÿ dxàø‘ý ùä/íÿ#w.Å7ëËÏó–]/KB:¸û —`ÿ·K¼ûæ0‰5_žü}!Â
endstream
endobj
676 0 obj <<
@@ -14421,62 +14424,62 @@
665 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 1]
-/Rect [247.899 663.241 543.291 674.145]
+/Rect [247.899 643.741 543.291 654.645]
/Subtype/Link/A<</Type/Action/S/URI/URI(http://www.spice-rtn.org/library/software/EX2DDIR)>>
>> endobj
666 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [201.374 651.286 244.412 662.19]
+/Rect [201.374 631.786 244.412 642.69]
/Subtype /Link
/A << /S /GoTo /D (cite.BeIfNiSk94) >>
>> endobj
667 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [248.227 651.286 270.145 662.19]
+/Rect [248.227 631.786 270.145 642.69]
/Subtype /Link
/A << /S /GoTo /D (cite.BeIfNiSk94) >>
>> endobj
668 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [95.91 378.143 189.628 389.047]
+/Rect [95.91 359.751 189.628 370.655]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
669 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [193.444 378.143 215.361 389.047]
+/Rect [193.444 359.751 215.361 370.655]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
670 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[1 0 0]
-/Rect [290.33 378.143 304.776 389.047]
+/Rect [290.33 359.751 304.776 370.655]
/Subtype /Link
/A << /S /GoTo /D (figure.4.2) >>
>> endobj
671 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [188.094 342.278 281.811 353.182]
+/Rect [188.094 323.885 281.811 334.789]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
672 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [285.627 342.278 307.545 353.182]
+/Rect [285.627 323.885 307.545 334.789]
/Subtype /Link
/A << /S /GoTo /D (cite.KoTr02a) >>
>> endobj
673 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 1]
-/Rect [99.228 282.502 196.862 293.406]
+/Rect [99.228 264.11 196.862 275.014]
/Subtype/Link/A<</Type/Action/S/URI/URI(www.llnl.gov/sac)>>
>> endobj
678 0 obj <<
@@ -14487,24 +14490,19 @@
/ProcSet [ /PDF /Text ]
>> endobj
692 0 obj <<
-/Length 1908
+/Length 1940
/Filter /FlateDecode
>>
stream
-xÚ…XKsã6¾ûWè¶Tà{R9Ì:öx¶2^ïHI¥*›rÁd±B´ÇùõÛ@)q2¥ñht7º¿~@tñ¼ ‹ÿÚ\\Ý&Ù"'y“Åf·Hù"
-³Åf»ø=¸¾{ÿ°¹ù¼\ñ˜Y®â„Ÿ—
-~¹¿ÿxÿ76w78Xÿçç%£4øÕ ùmýps}{ó‰ÿ´ücóï«Û( ŠâŒÐ05¬,š êóߛ̓]°Ës£xóœ$,ZõÅ— J¢,äÜ’ŒÇvÓŸt W ë|ñ“ºø/ü<O²òÜW#öÖ:¥I¥‹˜¦$ÉTý¶|î[éÄß-W!£Áµª¢-;ÕàŽÚáWï i·9ÝŽ¤-Eó\¹MÑlON~}× Îv}SèÒ‹…îEU½á¬ïä–Ì.Ú«¬BN(\Ü^å^iy ú‡!J Yü%[…#]ÖnÍ(lW<QWÖ}%P
-3/TÛÊî šmç(ÕɉB6Z¶cvéqÓ[àÊ_Õ*ÅÕÎpiU¹uçg„ëÕWŒ‘<v÷”¢Õû/½ø H at v”$ÁÆ2¥<è4l›!ó aÑr„¯ ¨ÆW†yÙ!Ê&g$Œ<È5RLâ díþŽÃrŒDéÀá”Ò9&i”
-‘4£FDòlòîœ'i2„"ŸÑ"%qš”‘MirO³,t°í[´Ò™ÔŒ’ÆþÄ?gX2 ˆ’Á|t&°È ªÅtðÐÎx¢æŒðÁUèK…^® œcM4Ь@6‰h8ÅŽ Ã,¨U
-è3ã<€+ºÕq8æV’Y}Y²8í®•ê ¹QÁ,a”™•×½lpdqX6ÏH1È=ïXHÂàNŠ%§ÁKÙ•Û-:<Û{qçc˘eH`9viø±o
-º”çÉÑ¥vfï_õ¤EÙH“ò8xzÃUvµ£&À¢,ø莿© ‡Vî @A1ñMÓãN! sÀÿŠC½÷’\€ÚE cÔÎh¹KV”hE-MÆ9GYDfá ³_Ö7Ÿ 7sX‹"’æÌÓ•ÍQÆ9S¸qLO
-](À%Z4º#ûþIF’˜ÂŒÛ"3„ |ÿr TWêrÉc¬KÄq£š•Ð¾˜}Qõ6¿AòÜì ³qæ2óÒ}õ^h Â$½¤¨¾Úº£Þò Š ©dv–0öxÍÙàHÈ´ÕpS¹{Oжú7¤$ベˆ£xœ£õ‘fÅ33ßψþ˜ˆ&
-ÒAG·–®Ð:œkljf>ä:¶ê±*;ýx¨$Xîq«ýØÊZ½H¢¿ê¹¬•‡$=¦õ âÐfÓÖàÏFw¾º2ÿd¾î.qzhÕS%k?õÜêŠ'æX<Ø„0o.n›4<±—êáB®ûX é¾k@¬¤±i¶nÔõ‡ƒju‡3Q¨¾Óeq%+a¾¿Yx\R5Þâ›í„Ø ™^ÎÐ:%ì VsÉíh³ÃCqn7n昃Æ:ƒu÷Þ
-?hL&zÞ–I \Žýæ|Û(7Z˦Ç.nÀÃÐóMD×Pm CªF¾ŒL¡äÎaH ¡þaË´å4¸SF›WÇXá·íȇÕú×)zÖw>
-bZËú©4UW_…K?¶nZ|!œžø¦Ùä-/iör ÌX`°Ç*gÞ`ÎwÉÛmiDÛDsDû»IókM5-̸SÉZ çÁÙ™ñ5àu³!ò7Ë DXRœÔæБò$…¨Šc Æ(>ˆç%$m=àªlN€ðuõ—£®À ÎFPóf²O8É£¡7:<v/Ø’ýˆ ¢[S Î}Yš§“:çT˜ý jçƒh w¥Ï™ßËÚßrŒy;ÿ¸]ßý½g‰³sÏ@†@ϘT løÎ Vç¼ËSoe¼QFD. ]éd¡Ã`€ƒu˜“QêýD¥·UïEÕ`C 6ï÷VF¢q3 sÞÍÄULÙÄͶ}Àz,;i_½ïwa³Î_Y¦SÛo°÷‚÷ÿN
-Ï]˜à" Tc_œf¥†’aúb[7ŸÊ’q'£ ÇóøÆM ›¡³)µã‰ÍyÛÈÊé€Ý! ž$ž÷ÎñBg;é”´PÀ·¶QCã×dO{öwNiúǤ&$'Uy+¡³\ØúhõÈñù¿uMµüFЄa1O0ùÄ-ïM—˜F}×W8~xÓ{|\¦AW´åAã:\ÖTå™×§–¡ÂúæÓc'^ )[rx›‹óœ$ÉЄ˜z
-꾶¥)G(Ë4ÕFþƒ¹ƒðš]·ÐŸ”M‰³]«j ý,ž¼ TÙJ\°Õ Ž’à²òKï–ÆØuÄÉ @·ÿx… È
-îÃ0¸¾ÿ¼6CžŸDÛÉ `·¢à\Px: G ‚Djô}çЖ‘aðQãØ)Gª¬a` ¢!èµt–…û–ºiáA.F¢OËF¤_8¾ˆh´²ë+ÝÍõ 'û<úæ´¤’oý=ö‹'¸
+xÚ…ÛrÛ¶òÝ_¡·Cͱ`¼§Ó‡Ôµãt×ÔNgÚŽ¦ S’PHÐŽûõ]`Š”™dôÀÅb±»Ø;D ºxwöÃæìâ:Í)R–.6»EÆY”/6ÛÅŸÁåÍÛ»ÍÕÇåŠ%4ˆÉr•¤4ø¸ŒiðÛííûÛw¸±¹¹B`ýËÏËÒàwäõÝÕåõÕöãòïÍO×q<'9¡QjXYabhΨS̯6g! t.¢ ,N +HÆ‹²>ûtFIœGŒY’1l7ýI‡¸x_‹ ÕÙ¯ðóü=ÉÊs_Ø[ëL”YBò8[$4#i‘¢ê×ò±o…·{³\E!
+.U}àìTƒ;j‡_½w¤Ýž Ät;’Vòæ±r›¼Ùžœ|Çû®\íú¦ÔÒ‹á¥îyU½àªïÄ–ÌÖá‚Ñè‚«$²\EŒP0‡½à˜ªvÂæÄ.Ø çK^ *ù¥¬äGÅvÁÀà0<)˜t«´8óE^6
+Óà_Ñ*„´¬ ÎØËb<Q'ë¾rÂ̺Tm+ºƒj¶£T''JÑhÑŽÙeÇMï€ oi« T;Ã¥UrëÎÏ*hnKáž!)wOÁ[½ÿÔóŒM$Vœ¦ÁÆ2¥,è4l0ô i9ÂWTã+ÃZv˜da8Ž×"$Qìý¢‘b’î3Æîï8L">$q6p §Ò9&Yœ‰<£FLŠ|òæ5F²t¨lF‹Œ$Y>R"B84¥)<ÍžW`°ªÑ¶oÑJ¯¤æ”D4ñ'þ?Ã2ŠxkMgÊXð¨ZB"Þ®CÂW¡/z¹và kŠ ¦Hœ’˜F.v\µqU jUCô¸àŠ;®…•d°OË0 Dû‚زR#7*f™Á<ïEƒCÙ<"Å ÷xß"‚Á—ŒO²“Û-:<Û{q ×cË4Ô°\xnø…_Jt)+Ò£KíÊÞ¾êAsÙS,‹$xxAlx±£&Áâ<x¨ C+vP Ÿ™ü¦Ùq§äŽ¹©bŸÔ{/É%¨EÚ ì„vF+\¡°¢xËka*Îë(‹)„Y4„Ùoë«÷›¹X‹c’¡§“ÍQÆk¦pã„„,›Ú°—hÞèŽìgø§9I“p”fÌö "Ä0¯ñ_Aƒƒê¤–K–¸À:Ç8nT³r´OfŸW½oP<7{Çl\¹ÌZº¯ÞsPÉmÀ€¤g)ª¯¶î¨·<¨¢G*™%À>^‹pðL¨´ÕpS¹{O,µúŠ2’³ÁLÄQD <ÎÐúH³b9IB?NñþXˆ& M¤ƒŽn-\Ÿw x±iÎÅPë´ÚªûJvúþP °ÜýV5ú¾µzDÖsU«ˆHv,ë;LÄa 1̦“É?2î|vÍüÁ$|ÿØãòЪ‡JÔ~é¹Õš'Ö@ lA˜73sUØKõp!7ü¬‡ é¾i@줉¶êúÃAµºÃ/UßiY^ˆŠ›¯‰ß<:¢T·XàV¯v" ÏÆLÏghö«¹âv´Ùa<·WsÌŒAJcÁºÇûáÈs' &“=/Ë4 ÎÇ~s¾m”¸Ö¢édzš÷à0rND×Ðm`BªF¾ŒM£dÎQDR£þ1 —+xÐàFmž c…߶w"ïVëߧѳ¾ñYÐZÔ-o„éˆ}æ®üؾiMðaŒ0H#œÙMÝò’f/çƒ{ {ìr¶à
+Öéü¾ÝJ#Úšc´¿9 š£“
+¿Ö Þ|mІÇI™¿Yæ êÄ’üÔ ¶†Ž”wARòª<ÖaÌâ\BÑÖC\Éæ$>¯þuÔøÃÙzÞL•b)#E<ÌF‡ûî G²ïñCtkšÃk_Á`–Ù¤Ï9f„ÀüE‡ÞyÇÛûô5ó[Uû‹ Ó;‹<™úg}óuÏ„$É_{*zÆ”
+üØÀ†ïL`vÎ[€žz+á2"r… 褓… €u˜“!õ~¢Ò˪„çªjp „ûð•‘hÜ ‘9ïæò*¡áÄÍv|Áza M ÷ý2lÖù+ëÑÌÙ>õÝÎ^)”&Á5¾¶aÈ0P}ð <=…™‹mßH})KÇIžŽ’ Ïã;5Ü “ÔŽ'çm#*§N‡ <<ïã…ÏvÂ)iCŸúF
+_S=íÙ?¥Ùßs=a¨Â[“måÒÖg« _ÿ[7T‹/$½‰0l0æÉ&ŸÈzk&¸Ô 껾BøîEïñq™]ÙʃF<\Ötå™×£–¡ „õÕ‡ûŽ?Á8"[rx™Ëó‚¤é0„˜~
+ê>·Ò´#”e†j#ÿÎÜ{Í.[˜Od#qµkUÐÏüAÁË@ÉV Âv8þA€Ëä§Þ¡1ŒqêHÒ“ Ýþï-6r 7q EÁåíǵY|àm'$ÀnÅÁ5¸ ôt6
+©Ñ÷ žC[bDFÁ{[`; ¤ ©²À6DCÐkÑ)K÷' ºiãA.F¢/ËF¤_8¾Ñ ´¢ë+ÝÍÅ×É¿N£¾9m©äKÿÎýø–Fù
endstream
endobj
691 0 obj <<
@@ -14636,14 +14634,14 @@
682 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [468.671 154.839 511.843 165.743]
+/Rect [468.671 134.914 511.843 145.818]
/Subtype /Link
/A << /S /GoTo /D (cite.TaLiTr07) >>
>> endobj
683 0 obj <<
/Type /Annot
/Border[0 0 0]/H/I/C[0 1 0]
-/Rect [515.761 154.839 537.679 165.743]
+/Rect [515.761 134.914 537.679 145.818]
/Subtype /Link
/A << /S /GoTo /D (cite.TaLiTr07) >>
>> endobj
@@ -14654,16 +14652,16 @@
/D [691 0 R /XYZ 172.661 519.651 null]
>> endobj
695 0 obj <<
-/D [691 0 R /XYZ 72 380.489 null]
+/D [691 0 R /XYZ 72 360.564 null]
>> endobj
696 0 obj <<
-/D [691 0 R /XYZ 72 331.008 null]
+/D [691 0 R /XYZ 72 311.082 null]
>> endobj
50 0 obj <<
-/D [691 0 R /XYZ 72 267.766 null]
+/D [691 0 R /XYZ 72 247.841 null]
>> endobj
690 0 obj <<
-/Font << /F68 554 0 R /F44 223 0 R /F11 382 0 R /F8 679 0 R /F14 694 0 R /F67 376 0 R /F43 222 0 R >>
+/Font << /F68 554 0 R /F44 223 0 R /F43 222 0 R /F11 382 0 R /F8 679 0 R /F14 694 0 R /F67 376 0 R >>
/XObject << /Im9 674 0 R >>
/ProcSet [ /PDF /Text ]
>> endobj
@@ -18884,8 +18882,8 @@
>> endobj
900 0 obj <<
/Author()/Title()/Subject()/Creator(LaTeX with hyperref package)/Producer(pdfTeX-1.40.3)/Keywords()
-/CreationDate (D:20121119005933+01'00')
-/ModDate (D:20121119005933+01'00')
+/CreationDate (D:20121215211022+01'00')
+/ModDate (D:20121215211022+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
@@ -18897,84 +18895,84 @@
0000000004 00000 f
0000000000 00000 f
0000000015 00000 n
-0000462923 00000 n
-0001761299 00000 n
+0000462924 00000 n
+0001761482 00000 n
0000000060 00000 n
0000000092 00000 n
-0000485648 00000 n
-0001761227 00000 n
+0000485649 00000 n
+0001761410 00000 n
0000000139 00000 n
0000000170 00000 n
-0000485704 00000 n
-0001761154 00000 n
+0000485705 00000 n
+0001761337 00000 n
0000000218 00000 n
0000000248 00000 n
-0000489168 00000 n
-0001761030 00000 n
+0000489169 00000 n
+0001761213 00000 n
0000000294 00000 n
0000000330 00000 n
-0000491749 00000 n
-0001760969 00000 n
+0000491750 00000 n
+0001761152 00000 n
0000000378 00000 n
0000000459 00000 n
-0000853687 00000 n
-0001760844 00000 n
+0000853688 00000 n
+0001761027 00000 n
0000000505 00000 n
0000000541 00000 n
-0000853739 00000 n
-0001760770 00000 n
+0000853740 00000 n
+0001760953 00000 n
0000000589 00000 n
0000000632 00000 n
-0000926501 00000 n
-0001760683 00000 n
+0000926502 00000 n
+0001760866 00000 n
0000000680 00000 n
0000000747 00000 n
-0000996868 00000 n
-0001760596 00000 n
+0000996869 00000 n
+0001760779 00000 n
0000000795 00000 n
0000000861 00000 n
-0000998504 00000 n
-0001760522 00000 n
+0000998505 00000 n
+0001760705 00000 n
0000000909 00000 n
0000000979 00000 n
-0001351240 00000 n
-0001760397 00000 n
+0001351338 00000 n
+0001760580 00000 n
0000001025 00000 n
0000001075 00000 n
-0001366242 00000 n
-0001760323 00000 n
+0001366425 00000 n
+0001760506 00000 n
0000001123 00000 n
0000001197 00000 n
-0001568905 00000 n
-0001760236 00000 n
+0001569088 00000 n
+0001760419 00000 n
0000001245 00000 n
0000001289 00000 n
-0001568957 00000 n
-0001760149 00000 n
+0001569140 00000 n
+0001760332 00000 n
0000001337 00000 n
0000001385 00000 n
-0001576885 00000 n
-0001760062 00000 n
+0001577068 00000 n
+0001760245 00000 n
0000001433 00000 n
0000001500 00000 n
-0001580751 00000 n
-0001759988 00000 n
+0001580934 00000 n
+0001760171 00000 n
0000001548 00000 n
0000001598 00000 n
-0001584070 00000 n
-0001759863 00000 n
+0001584253 00000 n
+0001760046 00000 n
0000001644 00000 n
0000001684 00000 n
-0001584122 00000 n
-0001759802 00000 n
+0001584305 00000 n
+0001759985 00000 n
0000001732 00000 n
0000001795 00000 n
-0001588142 00000 n
-0001759714 00000 n
+0001588325 00000 n
+0001759897 00000 n
0000001841 00000 n
0000001894 00000 n
-0001616929 00000 n
-0001759639 00000 n
+0001617112 00000 n
+0001759822 00000 n
0000001941 00000 n
0000001977 00000 n
0000002182 00000 n
@@ -18983,7 +18981,7 @@
0000002028 00000 n
0000432641 00000 n
0000432692 00000 n
-0001758788 00000 n
+0001758971 00000 n
0000005854 00000 n
0000006045 00000 n
0000006093 00000 n
@@ -19110,693 +19108,693 @@
0000431642 00000 n
0000432110 00000 n
0000432193 00000 n
-0000433272 00000 n
-0000433161 00000 n
+0000433273 00000 n
+0000433162 00000 n
0000432810 00000 n
-0001758082 00000 n
-0001758255 00000 n
-0000434250 00000 n
-0000434082 00000 n
-0000433357 00000 n
-0000434193 00000 n
-0000435552 00000 n
-0000435703 00000 n
-0000435856 00000 n
-0000436009 00000 n
-0000436158 00000 n
-0000436311 00000 n
-0000436462 00000 n
-0000436615 00000 n
-0000436768 00000 n
-0000436920 00000 n
-0000437073 00000 n
-0000437223 00000 n
-0000437376 00000 n
-0000437529 00000 n
-0000437681 00000 n
-0000437834 00000 n
-0000437986 00000 n
-0000438137 00000 n
-0000438288 00000 n
-0000438439 00000 n
-0000438701 00000 n
-0000435269 00000 n
-0000434335 00000 n
-0000438591 00000 n
-0000438644 00000 n
-0000443042 00000 n
-0000443382 00000 n
-0000443549 00000 n
-0000443724 00000 n
-0000443878 00000 n
-0000444032 00000 n
-0000444342 00000 n
-0000444497 00000 n
-0000444653 00000 n
-0000444807 00000 n
-0000444960 00000 n
-0000445118 00000 n
-0000445276 00000 n
-0000445432 00000 n
-0000445588 00000 n
-0000445938 00000 n
-0000446113 00000 n
-0000446538 00000 n
-0000446694 00000 n
-0000446850 00000 n
-0000447025 00000 n
-0000447201 00000 n
-0000447357 00000 n
-0000447513 00000 n
-0000447665 00000 n
-0000447817 00000 n
-0000447973 00000 n
-0000448129 00000 n
-0000448285 00000 n
-0000448441 00000 n
-0000448597 00000 n
-0000448752 00000 n
-0000448905 00000 n
-0000449056 00000 n
-0000449212 00000 n
-0000449370 00000 n
-0000449524 00000 n
-0000449678 00000 n
-0000449832 00000 n
-0000449985 00000 n
-0000450138 00000 n
-0000450291 00000 n
-0000450445 00000 n
-0000450599 00000 n
-0000450755 00000 n
-0000450911 00000 n
-0000451067 00000 n
-0000451223 00000 n
-0000451377 00000 n
-0000451533 00000 n
-0000451687 00000 n
-0000451840 00000 n
-0000451994 00000 n
-0000452148 00000 n
-0000452304 00000 n
-0000452460 00000 n
-0000452616 00000 n
-0000452772 00000 n
-0000452927 00000 n
-0000453083 00000 n
-0000453237 00000 n
-0000453390 00000 n
-0000453546 00000 n
-0000453702 00000 n
-0000453858 00000 n
-0000454013 00000 n
-0000454171 00000 n
-0000454329 00000 n
-0000454482 00000 n
-0000454636 00000 n
-0000454794 00000 n
-0000454952 00000 n
-0000455108 00000 n
-0000455265 00000 n
-0000455418 00000 n
-0000455571 00000 n
-0000455725 00000 n
-0000455879 00000 n
-0000456033 00000 n
-0000456187 00000 n
-0000456341 00000 n
-0000456495 00000 n
-0000456804 00000 n
-0000456958 00000 n
-0000457106 00000 n
-0000457258 00000 n
-0000457411 00000 n
-0000457563 00000 n
-0000457716 00000 n
-0000457869 00000 n
-0000458027 00000 n
-0000458185 00000 n
-0000458343 00000 n
-0000458501 00000 n
-0000458658 00000 n
-0000458815 00000 n
-0000458979 00000 n
-0000459144 00000 n
-0000459299 00000 n
-0000459454 00000 n
-0000459763 00000 n
-0000459917 00000 n
-0000460074 00000 n
-0000460232 00000 n
-0000460385 00000 n
-0000460538 00000 n
-0000460693 00000 n
-0000460848 00000 n
-0000461005 00000 n
-0000461161 00000 n
-0000461318 00000 n
-0000461475 00000 n
-0000461629 00000 n
-0000461783 00000 n
-0000461943 00000 n
-0000462103 00000 n
-0000462256 00000 n
-0000462409 00000 n
-0000462717 00000 n
-0000462974 00000 n
-0000441911 00000 n
-0000438786 00000 n
-0000462870 00000 n
-0001757561 00000 n
-0000443212 00000 n
-0000444188 00000 n
-0000445764 00000 n
-0000446326 00000 n
-0000456651 00000 n
-0001756533 00000 n
-0000459609 00000 n
-0000462565 00000 n
-0001603903 00000 n
-0001608872 00000 n
+0001758265 00000 n
+0001758438 00000 n
+0000434251 00000 n
+0000434083 00000 n
+0000433358 00000 n
+0000434194 00000 n
+0000435553 00000 n
+0000435704 00000 n
+0000435857 00000 n
+0000436010 00000 n
+0000436159 00000 n
+0000436312 00000 n
+0000436463 00000 n
+0000436616 00000 n
+0000436769 00000 n
+0000436921 00000 n
+0000437074 00000 n
+0000437224 00000 n
+0000437377 00000 n
+0000437530 00000 n
+0000437682 00000 n
+0000437835 00000 n
+0000437987 00000 n
+0000438138 00000 n
+0000438289 00000 n
+0000438440 00000 n
+0000438702 00000 n
+0000435270 00000 n
+0000434336 00000 n
+0000438592 00000 n
+0000438645 00000 n
+0000443043 00000 n
+0000443383 00000 n
+0000443550 00000 n
+0000443725 00000 n
+0000443879 00000 n
+0000444033 00000 n
+0000444343 00000 n
+0000444498 00000 n
+0000444654 00000 n
+0000444808 00000 n
+0000444961 00000 n
+0000445119 00000 n
+0000445277 00000 n
+0000445433 00000 n
+0000445589 00000 n
+0000445939 00000 n
+0000446114 00000 n
+0000446539 00000 n
+0000446695 00000 n
+0000446851 00000 n
+0000447026 00000 n
+0000447202 00000 n
+0000447358 00000 n
+0000447514 00000 n
+0000447666 00000 n
+0000447818 00000 n
+0000447974 00000 n
+0000448130 00000 n
+0000448286 00000 n
+0000448442 00000 n
+0000448598 00000 n
+0000448753 00000 n
+0000448906 00000 n
+0000449057 00000 n
+0000449213 00000 n
+0000449371 00000 n
+0000449525 00000 n
+0000449679 00000 n
+0000449833 00000 n
+0000449986 00000 n
+0000450139 00000 n
+0000450292 00000 n
+0000450446 00000 n
+0000450600 00000 n
+0000450756 00000 n
+0000450912 00000 n
+0000451068 00000 n
+0000451224 00000 n
+0000451378 00000 n
+0000451534 00000 n
+0000451688 00000 n
+0000451841 00000 n
+0000451995 00000 n
+0000452149 00000 n
+0000452305 00000 n
+0000452461 00000 n
+0000452617 00000 n
+0000452773 00000 n
+0000452928 00000 n
+0000453084 00000 n
+0000453238 00000 n
+0000453391 00000 n
+0000453547 00000 n
+0000453703 00000 n
+0000453859 00000 n
+0000454014 00000 n
+0000454172 00000 n
+0000454330 00000 n
+0000454483 00000 n
+0000454637 00000 n
+0000454795 00000 n
+0000454953 00000 n
+0000455109 00000 n
+0000455266 00000 n
+0000455419 00000 n
+0000455572 00000 n
+0000455726 00000 n
+0000455880 00000 n
+0000456034 00000 n
+0000456188 00000 n
+0000456342 00000 n
+0000456496 00000 n
+0000456805 00000 n
+0000456959 00000 n
+0000457107 00000 n
+0000457259 00000 n
+0000457412 00000 n
+0000457564 00000 n
+0000457717 00000 n
+0000457870 00000 n
+0000458028 00000 n
+0000458186 00000 n
+0000458344 00000 n
+0000458502 00000 n
+0000458659 00000 n
+0000458816 00000 n
+0000458980 00000 n
+0000459145 00000 n
+0000459300 00000 n
+0000459455 00000 n
+0000459764 00000 n
+0000459918 00000 n
+0000460075 00000 n
+0000460233 00000 n
+0000460386 00000 n
+0000460539 00000 n
+0000460694 00000 n
+0000460849 00000 n
+0000461006 00000 n
+0000461162 00000 n
+0000461319 00000 n
+0000461476 00000 n
+0000461630 00000 n
+0000461784 00000 n
+0000461944 00000 n
+0000462104 00000 n
+0000462257 00000 n
+0000462410 00000 n
+0000462718 00000 n
+0000462975 00000 n
+0000441912 00000 n
+0000438787 00000 n
+0000462871 00000 n
+0001757744 00000 n
+0000443213 00000 n
+0000444189 00000 n
+0000445765 00000 n
+0000446327 00000 n
+0000456652 00000 n
+0001756716 00000 n
+0000459610 00000 n
+0000462566 00000 n
+0001604086 00000 n
+0001609055 00000 n
+0001609116 00000 n
0001608933 00000 n
-0001608750 00000 n
-0001608994 00000 n
-0001614167 00000 n
-0001613496 00000 n
-0001594404 00000 n
-0001598640 00000 n
-0001604694 00000 n
-0001598701 00000 n
-0001608508 00000 n
-0001604633 00000 n
-0001614411 00000 n
-0001593917 00000 n
-0001613618 00000 n
-0001594039 00000 n
-0001599493 00000 n
-0001603720 00000 n
-0001598823 00000 n
-0001599554 00000 n
-0001613679 00000 n
0001609177 00000 n
-0001599735 00000 n
-0001593978 00000 n
-0001594343 00000 n
+0001614350 00000 n
+0001613679 00000 n
+0001594587 00000 n
+0001598823 00000 n
+0001604877 00000 n
+0001598884 00000 n
+0001608691 00000 n
+0001604816 00000 n
+0001614594 00000 n
+0001594100 00000 n
+0001613801 00000 n
+0001594222 00000 n
+0001599676 00000 n
+0001603903 00000 n
+0001599006 00000 n
+0001599737 00000 n
+0001613862 00000 n
+0001609360 00000 n
+0001599918 00000 n
+0001594161 00000 n
+0001594526 00000 n
+0001599493 00000 n
0001599310 00000 n
-0001599127 00000 n
-0001599188 00000 n
-0001599066 00000 n
0001599371 00000 n
-0001609357 00000 n
-0001608689 00000 n
-0001613862 00000 n
0001599249 00000 n
-0001598884 00000 n
-0001599005 00000 n
-0001613801 00000 n
-0001604208 00000 n
-0001614228 00000 n
-0001604330 00000 n
-0001594587 00000 n
-0001604574 00000 n
+0001599554 00000 n
+0001609540 00000 n
+0001608872 00000 n
+0001614045 00000 n
+0001599432 00000 n
+0001599067 00000 n
+0001599188 00000 n
+0001613984 00000 n
0001604391 00000 n
+0001614411 00000 n
0001604513 00000 n
-0001604452 00000 n
-0001603781 00000 n
-0001615264 00000 n
-0001614106 00000 n
-0001604147 00000 n
-0001609055 00000 n
-0000467097 00000 n
-0000467257 00000 n
-0000467417 00000 n
-0000467575 00000 n
-0000467731 00000 n
-0000467886 00000 n
-0000468048 00000 n
-0000468209 00000 n
-0000468364 00000 n
-0000468519 00000 n
-0000468681 00000 n
-0000468842 00000 n
-0000469004 00000 n
-0000469166 00000 n
-0000469326 00000 n
-0000469484 00000 n
-0000469660 00000 n
-0000469834 00000 n
-0000469990 00000 n
-0000470146 00000 n
-0000470300 00000 n
-0000470454 00000 n
-0000470608 00000 n
-0000470762 00000 n
-0000470924 00000 n
-0000471086 00000 n
-0000471244 00000 n
-0000471402 00000 n
-0000471562 00000 n
-0000471724 00000 n
-0000471890 00000 n
-0000472056 00000 n
-0000472213 00000 n
-0000472371 00000 n
-0000472533 00000 n
-0000472695 00000 n
-0000472848 00000 n
-0000473001 00000 n
-0000473308 00000 n
-0000473462 00000 n
-0000473616 00000 n
-0000473770 00000 n
-0000473924 00000 n
-0000474078 00000 n
-0000474231 00000 n
-0000474384 00000 n
-0000474538 00000 n
-0000474691 00000 n
-0000474847 00000 n
-0000475001 00000 n
-0000475153 00000 n
-0000475305 00000 n
-0000475463 00000 n
-0000475621 00000 n
-0000475775 00000 n
-0000475928 00000 n
-0000476082 00000 n
-0000476236 00000 n
-0000476388 00000 n
-0000476540 00000 n
-0000476694 00000 n
-0000476848 00000 n
-0000477004 00000 n
-0000477159 00000 n
-0000477315 00000 n
-0000477470 00000 n
-0000477625 00000 n
-0000477781 00000 n
-0000477934 00000 n
-0000478087 00000 n
-0000478241 00000 n
-0000478396 00000 n
-0000478552 00000 n
-0000478709 00000 n
-0000478863 00000 n
-0000479016 00000 n
-0000479169 00000 n
-0000479322 00000 n
-0000479477 00000 n
-0000479633 00000 n
-0000479809 00000 n
-0000479985 00000 n
-0000480144 00000 n
-0000480304 00000 n
-0000480464 00000 n
-0000480624 00000 n
-0000480785 00000 n
-0000480943 00000 n
-0000481100 00000 n
-0000481258 00000 n
-0000481415 00000 n
-0000481569 00000 n
-0000481727 00000 n
-0000481884 00000 n
-0000482046 00000 n
-0000482208 00000 n
-0000482362 00000 n
-0000482516 00000 n
-0000482670 00000 n
-0000482824 00000 n
-0000483147 00000 n
-0000483309 00000 n
-0000483463 00000 n
-0000483617 00000 n
-0000483998 00000 n
-0000484153 00000 n
-0000484309 00000 n
-0000484485 00000 n
-0000484660 00000 n
-0000484814 00000 n
-0000484968 00000 n
-0000485124 00000 n
-0000485280 00000 n
-0000485437 00000 n
-0000485760 00000 n
-0000466038 00000 n
-0000463085 00000 n
-0000485595 00000 n
-0001758605 00000 n
-0000473155 00000 n
-0000482986 00000 n
-0000483808 00000 n
-0001594099 00000 n
-0001594160 00000 n
-0001613984 00000 n
+0001594770 00000 n
+0001604757 00000 n
+0001604574 00000 n
+0001604696 00000 n
+0001604635 00000 n
0001603964 00000 n
+0001615447 00000 n
+0001614289 00000 n
+0001604330 00000 n
+0001609238 00000 n
+0000467098 00000 n
+0000467258 00000 n
+0000467418 00000 n
+0000467576 00000 n
+0000467732 00000 n
+0000467887 00000 n
+0000468049 00000 n
+0000468210 00000 n
+0000468365 00000 n
+0000468520 00000 n
+0000468682 00000 n
+0000468843 00000 n
+0000469005 00000 n
+0000469167 00000 n
+0000469327 00000 n
+0000469485 00000 n
+0000469661 00000 n
+0000469835 00000 n
+0000469991 00000 n
+0000470147 00000 n
+0000470301 00000 n
+0000470455 00000 n
+0000470609 00000 n
+0000470763 00000 n
+0000470925 00000 n
+0000471087 00000 n
+0000471245 00000 n
+0000471403 00000 n
+0000471563 00000 n
+0000471725 00000 n
+0000471891 00000 n
+0000472057 00000 n
+0000472214 00000 n
+0000472372 00000 n
+0000472534 00000 n
+0000472696 00000 n
+0000472849 00000 n
+0000473002 00000 n
+0000473309 00000 n
+0000473463 00000 n
+0000473617 00000 n
+0000473771 00000 n
+0000473925 00000 n
+0000474079 00000 n
+0000474232 00000 n
+0000474385 00000 n
+0000474539 00000 n
+0000474692 00000 n
+0000474848 00000 n
+0000475002 00000 n
+0000475154 00000 n
+0000475306 00000 n
+0000475464 00000 n
+0000475622 00000 n
+0000475776 00000 n
+0000475929 00000 n
+0000476083 00000 n
+0000476237 00000 n
+0000476389 00000 n
+0000476541 00000 n
+0000476695 00000 n
+0000476849 00000 n
+0000477005 00000 n
+0000477160 00000 n
+0000477316 00000 n
+0000477471 00000 n
+0000477626 00000 n
+0000477782 00000 n
+0000477935 00000 n
+0000478088 00000 n
+0000478242 00000 n
+0000478397 00000 n
+0000478553 00000 n
+0000478710 00000 n
+0000478864 00000 n
+0000479017 00000 n
+0000479170 00000 n
+0000479323 00000 n
+0000479478 00000 n
+0000479634 00000 n
+0000479810 00000 n
+0000479986 00000 n
+0000480145 00000 n
+0000480305 00000 n
+0000480465 00000 n
+0000480625 00000 n
+0000480786 00000 n
+0000480944 00000 n
+0000481101 00000 n
+0000481259 00000 n
+0000481416 00000 n
+0000481570 00000 n
+0000481728 00000 n
+0000481885 00000 n
+0000482047 00000 n
+0000482209 00000 n
+0000482363 00000 n
+0000482517 00000 n
+0000482671 00000 n
+0000482825 00000 n
+0000483148 00000 n
+0000483310 00000 n
+0000483464 00000 n
+0000483618 00000 n
+0000483999 00000 n
+0000484154 00000 n
+0000484310 00000 n
+0000484486 00000 n
+0000484661 00000 n
+0000484815 00000 n
+0000484969 00000 n
+0000485125 00000 n
+0000485281 00000 n
+0000485438 00000 n
+0000485761 00000 n
+0000466039 00000 n
+0000463086 00000 n
+0000485596 00000 n
+0001758788 00000 n
+0000473156 00000 n
+0000482987 00000 n
+0000483809 00000 n
+0001594282 00000 n
+0001594343 00000 n
+0001614167 00000 n
+0001604147 00000 n
+0001604452 00000 n
+0001608508 00000 n
+0001608569 00000 n
+0001608630 00000 n
0001604269 00000 n
-0001608325 00000 n
-0001608386 00000 n
-0001608447 00000 n
-0001604086 00000 n
+0001604208 00000 n
+0001598945 00000 n
+0001599797 00000 n
+0001614533 00000 n
+0001613923 00000 n
+0001613557 00000 n
+0001599127 00000 n
0001604025 00000 n
-0001598762 00000 n
-0001599614 00000 n
-0001614350 00000 n
+0001599615 00000 n
+0001594831 00000 n
0001613740 00000 n
-0001613374 00000 n
-0001598944 00000 n
-0001603842 00000 n
-0001599432 00000 n
-0001594648 00000 n
-0001613557 00000 n
-0001613923 00000 n
-0001599674 00000 n
-0001613313 00000 n
-0001609116 00000 n
-0001614289 00000 n
-0001608629 00000 n
-0001609237 00000 n
-0001608569 00000 n
-0001609296 00000 n
-0001608811 00000 n
-0000486396 00000 n
-0000486231 00000 n
-0000485871 00000 n
-0000486343 00000 n
-0001758903 00000 n
-0000488620 00000 n
-0000488773 00000 n
-0000488927 00000 n
-0000491289 00000 n
-0000489220 00000 n
-0000488472 00000 n
-0000486481 00000 n
-0000489115 00000 n
-0001757390 00000 n
-0001613435 00000 n
-0000491805 00000 n
-0000491149 00000 n
-0000489331 00000 n
-0000491696 00000 n
-0000491493 00000 n
-0000493717 00000 n
-0000853445 00000 n
-0000853857 00000 n
-0000493585 00000 n
-0000491916 00000 n
-0000853634 00000 n
-0000853795 00000 n
-0000856791 00000 n
-0000923138 00000 n
-0000856679 00000 n
-0000853991 00000 n
-0000923023 00000 n
-0000923076 00000 n
-0000921007 00000 n
-0000921269 00000 n
-0000922447 00000 n
-0000922470 00000 n
-0000922658 00000 n
-0000926295 00000 n
-0000928564 00000 n
-0000937075 00000 n
-0000926617 00000 n
-0000926163 00000 n
-0000923290 00000 n
-0000926448 00000 n
-0001757732 00000 n
-0001757907 00000 n
-0000926557 00000 n
-0000996924 00000 n
-0000928452 00000 n
-0000926780 00000 n
-0000996691 00000 n
-0000996744 00000 n
-0000996806 00000 n
-0001759020 00000 n
-0000928871 00000 n
-0000929051 00000 n
-0000937005 00000 n
-0000937025 00000 n
-0000937053 00000 n
-0000937382 00000 n
-0000937562 00000 n
-0000996620 00000 n
-0000996640 00000 n
-0000996668 00000 n
-0000998560 00000 n
-0000998339 00000 n
-0000997089 00000 n
-0000998451 00000 n
-0001001051 00000 n
-0001079810 00000 n
-0001192833 00000 n
-0001351353 00000 n
-0001000939 00000 n
-0000998671 00000 n
-0001351187 00000 n
-0001351292 00000 n
-0001355040 00000 n
-0001355237 00000 n
-0001355394 00000 n
-0001355551 00000 n
-0001355704 00000 n
-0001355859 00000 n
+0001614106 00000 n
+0001599857 00000 n
+0001613496 00000 n
+0001609299 00000 n
+0001614472 00000 n
+0001608812 00000 n
+0001609420 00000 n
+0001608752 00000 n
+0001609479 00000 n
+0001608994 00000 n
+0000486397 00000 n
+0000486232 00000 n
+0000485872 00000 n
+0000486344 00000 n
+0001759086 00000 n
+0000488621 00000 n
+0000488774 00000 n
+0000488928 00000 n
+0000491290 00000 n
+0000489221 00000 n
+0000488473 00000 n
+0000486482 00000 n
+0000489116 00000 n
+0001757573 00000 n
+0001613618 00000 n
+0000491806 00000 n
+0000491150 00000 n
+0000489332 00000 n
+0000491697 00000 n
+0000491494 00000 n
+0000493718 00000 n
+0000853446 00000 n
+0000853858 00000 n
+0000493586 00000 n
+0000491917 00000 n
+0000853635 00000 n
+0000853796 00000 n
+0000856792 00000 n
+0000923139 00000 n
+0000856680 00000 n
+0000853992 00000 n
+0000923024 00000 n
+0000923077 00000 n
+0000921008 00000 n
+0000921270 00000 n
+0000922448 00000 n
+0000922471 00000 n
+0000922659 00000 n
+0000926296 00000 n
+0000928565 00000 n
+0000937076 00000 n
+0000926618 00000 n
+0000926164 00000 n
+0000923291 00000 n
+0000926449 00000 n
+0001757915 00000 n
+0001758090 00000 n
+0000926558 00000 n
+0000996925 00000 n
+0000928453 00000 n
+0000926781 00000 n
+0000996692 00000 n
+0000996745 00000 n
+0000996807 00000 n
+0001759203 00000 n
+0000928872 00000 n
+0000929052 00000 n
+0000937006 00000 n
+0000937026 00000 n
+0000937054 00000 n
+0000937383 00000 n
+0000937563 00000 n
+0000996621 00000 n
+0000996641 00000 n
+0000996669 00000 n
+0000998561 00000 n
+0000998340 00000 n
+0000997090 00000 n
+0000998452 00000 n
+0001001149 00000 n
+0001079908 00000 n
+0001192931 00000 n
+0001351452 00000 n
+0001001037 00000 n
+0000998672 00000 n
+0001351285 00000 n
+0001351390 00000 n
+0001355192 00000 n
+0001355389 00000 n
+0001355546 00000 n
+0001355703 00000 n
+0001355856 00000 n
0001356011 00000 n
-0001356166 00000 n
-0001356321 00000 n
-0001358802 00000 n
-0001356537 00000 n
-0001354844 00000 n
-0001351513 00000 n
-0001356484 00000 n
-0001756821 00000 n
-0001594221 00000 n
-0001366066 00000 n
-0001365701 00000 n
-0001365857 00000 n
+0001356163 00000 n
+0001356318 00000 n
+0001356473 00000 n
+0001358985 00000 n
+0001356688 00000 n
+0001354996 00000 n
+0001351612 00000 n
+0001356635 00000 n
+0001757004 00000 n
+0001594404 00000 n
+0001366249 00000 n
+0001365884 00000 n
+0001366040 00000 n
+0001567807 00000 n
+0001567963 00000 n
+0001568118 00000 n
+0001568270 00000 n
+0001568422 00000 n
+0001568575 00000 n
+0001366481 00000 n
+0001358845 00000 n
+0001356824 00000 n
+0001366196 00000 n
+0001757286 00000 n
+0001366311 00000 n
+0001366368 00000 n
+0001363114 00000 n
+0001363325 00000 n
+0001363381 00000 n
+0001363775 00000 n
+0001363797 00000 n
+0001364078 00000 n
+0001614228 00000 n
+0001369107 00000 n
+0001568727 00000 n
+0001568881 00000 n
+0001571255 00000 n
+0001569196 00000 n
+0001368919 00000 n
+0001366658 00000 n
+0001569035 00000 n
+0001369732 00000 n
+0001369965 00000 n
+0001370021 00000 n
+0001370070 00000 n
+0001370135 00000 n
+0001370207 00000 n
+0001567284 00000 n
+0001567305 00000 n
+0001567464 00000 n
0001567624 00000 n
-0001567780 00000 n
-0001567935 00000 n
-0001568087 00000 n
-0001568239 00000 n
-0001568392 00000 n
-0001366298 00000 n
-0001358662 00000 n
-0001356673 00000 n
-0001366013 00000 n
-0001757103 00000 n
-0001366128 00000 n
-0001366185 00000 n
-0001362931 00000 n
-0001363142 00000 n
-0001363198 00000 n
-0001363592 00000 n
-0001363614 00000 n
-0001363895 00000 n
-0001614045 00000 n
-0001368924 00000 n
-0001568544 00000 n
-0001568698 00000 n
-0001571072 00000 n
-0001569013 00000 n
-0001368736 00000 n
-0001366475 00000 n
-0001568852 00000 n
-0001369549 00000 n
-0001369782 00000 n
-0001369838 00000 n
-0001369887 00000 n
-0001369952 00000 n
-0001370024 00000 n
-0001567101 00000 n
-0001567122 00000 n
-0001567281 00000 n
-0001567441 00000 n
-0001594526 00000 n
-0001603659 00000 n
-0001594465 00000 n
-0001576463 00000 n
-0001576617 00000 n
-0001576941 00000 n
-0001570932 00000 n
-0001569153 00000 n
-0001576770 00000 n
-0001576823 00000 n
-0001759137 00000 n
-0001576194 00000 n
-0001576394 00000 n
-0001576441 00000 n
-0001579773 00000 n
-0001579928 00000 n
-0001580082 00000 n
-0001580236 00000 n
-0001580388 00000 n
-0001580543 00000 n
-0001580869 00000 n
-0001579601 00000 n
-0001577081 00000 n
-0001580698 00000 n
-0001756677 00000 n
-0001580807 00000 n
-0001594282 00000 n
-0001581938 00000 n
-0001581711 00000 n
-0001581031 00000 n
-0001581823 00000 n
-0001581876 00000 n
-0001584178 00000 n
-0001583905 00000 n
-0001582074 00000 n
-0001584017 00000 n
-0001586448 00000 n
-0001586169 00000 n
-0001584276 00000 n
-0001586281 00000 n
-0001586334 00000 n
-0001756388 00000 n
-0001757246 00000 n
-0001586391 00000 n
-0001588194 00000 n
-0001587977 00000 n
-0001586623 00000 n
-0001588089 00000 n
-0001756962 00000 n
-0001589306 00000 n
-0001589645 00000 n
-0001589927 00000 n
-0001589158 00000 n
-0001588368 00000 n
-0001589817 00000 n
-0001589870 00000 n
-0001589476 00000 n
-0001759254 00000 n
-0001590619 00000 n
-0001590397 00000 n
-0001590025 00000 n
-0001590509 00000 n
-0001590562 00000 n
0001594709 00000 n
-0001593695 00000 n
-0001590704 00000 n
-0001593807 00000 n
-0001593860 00000 n
-0001758428 00000 n
-0001599796 00000 n
-0001598475 00000 n
-0001594807 00000 n
-0001598587 00000 n
-0001604754 00000 n
-0001603494 00000 n
-0001599894 00000 n
-0001603606 00000 n
-0001609417 00000 n
-0001608160 00000 n
-0001604852 00000 n
-0001608272 00000 n
-0001614472 00000 n
-0001613148 00000 n
-0001609515 00000 n
-0001613260 00000 n
-0001759371 00000 n
-0001615325 00000 n
-0001615099 00000 n
-0001614570 00000 n
-0001615211 00000 n
-0001617038 00000 n
-0001616764 00000 n
-0001615423 00000 n
-0001616876 00000 n
-0001616981 00000 n
-0001755767 00000 n
-0001617136 00000 n
-0001618032 00000 n
-0001618057 00000 n
-0001618082 00000 n
-0001618107 00000 n
-0001618132 00000 n
-0001618790 00000 n
-0001619130 00000 n
-0001619245 00000 n
-0001619632 00000 n
-0001619863 00000 n
-0001620058 00000 n
-0001620663 00000 n
-0001621128 00000 n
-0001622102 00000 n
-0001622579 00000 n
-0001625106 00000 n
-0001625328 00000 n
-0001630894 00000 n
-0001631180 00000 n
-0001632910 00000 n
-0001633129 00000 n
-0001637257 00000 n
-0001637555 00000 n
-0001639288 00000 n
-0001639508 00000 n
-0001641466 00000 n
-0001641746 00000 n
-0001643128 00000 n
-0001643354 00000 n
-0001653299 00000 n
-0001653585 00000 n
-0001671582 00000 n
-0001672095 00000 n
-0001683280 00000 n
-0001683657 00000 n
-0001688933 00000 n
-0001689207 00000 n
-0001705444 00000 n
-0001705897 00000 n
-0001726495 00000 n
-0001727170 00000 n
-0001738319 00000 n
-0001738636 00000 n
-0001755307 00000 n
-0001759464 00000 n
-0001759565 00000 n
-0001761408 00000 n
-0001761599 00000 n
-0001761781 00000 n
-0001761995 00000 n
-0001762235 00000 n
-0001762449 00000 n
-0001762653 00000 n
-0001762885 00000 n
-0001763106 00000 n
-0001763333 00000 n
-0001763544 00000 n
-0001763772 00000 n
-0001763990 00000 n
-0001764196 00000 n
-0001764422 00000 n
-0001764628 00000 n
-0001764854 00000 n
-0001765078 00000 n
-0001765263 00000 n
-0001765433 00000 n
-0001765602 00000 n
-0001765772 00000 n
-0001765941 00000 n
-0001766103 00000 n
-0001766298 00000 n
-0001766494 00000 n
-0001766690 00000 n
-0001766839 00000 n
-0001766954 00000 n
-0001767079 00000 n
-0001767191 00000 n
-0001767303 00000 n
-0001767394 00000 n
-0001767499 00000 n
-0001767537 00000 n
-0001767747 00000 n
+0001603842 00000 n
+0001594648 00000 n
+0001576646 00000 n
+0001576800 00000 n
+0001577124 00000 n
+0001571115 00000 n
+0001569336 00000 n
+0001576953 00000 n
+0001577006 00000 n
+0001759320 00000 n
+0001576377 00000 n
+0001576577 00000 n
+0001576624 00000 n
+0001579956 00000 n
+0001580111 00000 n
+0001580265 00000 n
+0001580419 00000 n
+0001580571 00000 n
+0001580726 00000 n
+0001581052 00000 n
+0001579784 00000 n
+0001577264 00000 n
+0001580881 00000 n
+0001756860 00000 n
+0001580990 00000 n
+0001594465 00000 n
+0001582121 00000 n
+0001581894 00000 n
+0001581214 00000 n
+0001582006 00000 n
+0001582059 00000 n
+0001584361 00000 n
+0001584088 00000 n
+0001582257 00000 n
+0001584200 00000 n
+0001586631 00000 n
+0001586352 00000 n
+0001584459 00000 n
+0001586464 00000 n
+0001586517 00000 n
+0001756571 00000 n
+0001757429 00000 n
+0001586574 00000 n
+0001588377 00000 n
+0001588160 00000 n
+0001586806 00000 n
+0001588272 00000 n
+0001757145 00000 n
+0001589489 00000 n
+0001589828 00000 n
+0001590110 00000 n
+0001589341 00000 n
+0001588551 00000 n
+0001590000 00000 n
+0001590053 00000 n
+0001589659 00000 n
+0001759437 00000 n
+0001590802 00000 n
+0001590580 00000 n
+0001590208 00000 n
+0001590692 00000 n
+0001590745 00000 n
+0001594892 00000 n
+0001593878 00000 n
+0001590887 00000 n
+0001593990 00000 n
+0001594043 00000 n
+0001758611 00000 n
+0001599979 00000 n
+0001598658 00000 n
+0001594990 00000 n
+0001598770 00000 n
+0001604937 00000 n
+0001603677 00000 n
+0001600077 00000 n
+0001603789 00000 n
+0001609600 00000 n
+0001608343 00000 n
+0001605035 00000 n
+0001608455 00000 n
+0001614655 00000 n
+0001613331 00000 n
+0001609698 00000 n
+0001613443 00000 n
+0001759554 00000 n
+0001615508 00000 n
+0001615282 00000 n
+0001614753 00000 n
+0001615394 00000 n
+0001617221 00000 n
+0001616947 00000 n
+0001615606 00000 n
+0001617059 00000 n
+0001617164 00000 n
+0001755950 00000 n
+0001617319 00000 n
+0001618215 00000 n
+0001618240 00000 n
+0001618265 00000 n
+0001618290 00000 n
+0001618315 00000 n
+0001618973 00000 n
+0001619313 00000 n
+0001619428 00000 n
+0001619815 00000 n
+0001620046 00000 n
+0001620241 00000 n
+0001620846 00000 n
+0001621311 00000 n
+0001622285 00000 n
+0001622762 00000 n
+0001625289 00000 n
+0001625511 00000 n
+0001631077 00000 n
+0001631363 00000 n
+0001633093 00000 n
+0001633312 00000 n
+0001637440 00000 n
+0001637738 00000 n
+0001639471 00000 n
+0001639691 00000 n
+0001641649 00000 n
+0001641929 00000 n
+0001643311 00000 n
+0001643537 00000 n
+0001653482 00000 n
+0001653768 00000 n
+0001671765 00000 n
+0001672278 00000 n
+0001683463 00000 n
+0001683840 00000 n
+0001689116 00000 n
+0001689390 00000 n
+0001705627 00000 n
+0001706080 00000 n
+0001726678 00000 n
+0001727353 00000 n
+0001738502 00000 n
+0001738819 00000 n
+0001755490 00000 n
+0001759647 00000 n
+0001759748 00000 n
+0001761591 00000 n
+0001761782 00000 n
+0001761964 00000 n
+0001762178 00000 n
+0001762418 00000 n
+0001762632 00000 n
+0001762836 00000 n
+0001763068 00000 n
+0001763289 00000 n
+0001763516 00000 n
+0001763727 00000 n
+0001763955 00000 n
+0001764173 00000 n
+0001764379 00000 n
+0001764605 00000 n
+0001764811 00000 n
+0001765037 00000 n
+0001765261 00000 n
+0001765446 00000 n
+0001765616 00000 n
+0001765785 00000 n
+0001765955 00000 n
+0001766124 00000 n
+0001766286 00000 n
+0001766481 00000 n
+0001766677 00000 n
+0001766873 00000 n
+0001767022 00000 n
+0001767137 00000 n
+0001767262 00000 n
+0001767374 00000 n
+0001767486 00000 n
+0001767577 00000 n
+0001767682 00000 n
+0001767720 00000 n
+0001767930 00000 n
trailer
<< /Size 901
/Root 899 0 R
/Info 900 0 R
-/ID [<C0CED3F414CF57D43D6D066A52144D54> <C0CED3F414CF57D43D6D066A52144D54>] >>
+/ID [<6EEB80CBFA3043184897E8BB0BCCB99A> <6EEB80CBFA3043184897E8BB0BCCB99A>] >>
startxref
-1768078
+1768261
%%EOF
Modified: seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex
===================================================================
--- seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex 2012-12-15 16:09:50 UTC (rev 21172)
+++ seismo/2D/SPECFEM2D/trunk/doc/USER_MANUAL/manual_SPECFEM2D.tex 2012-12-16 00:53:57 UTC (rev 21173)
@@ -631,7 +631,7 @@
\includegraphics[width=2in]{figures/image0000400.jpg}
\includegraphics[width=2in]{figures/image0000500.jpg}
\par\end{centering}
-\caption{Wavefield snapshots of the default example generated by \texttt{xspecfem2D} when parameter \texttt{output\_color\_image} is set to \texttt{.true.}. To create smaller (subsampled) images you can change parameter \texttt{factor\_subsample\_image = 1} to a higher value in file \texttt{setup/constants.h} and recompile the whole code using \texttt{make all}. This can be useful in the case of very large models.}
+\caption{Wavefield snapshots of the default example generated by \texttt{xspecfem2D} when parameter \texttt{output\_color\_image} is set to \texttt{.true.}. To create smaller (subsampled) images you can change double precision parameter \texttt{factor\_subsample\_image = 1.0} to a higher value in file \texttt{DATA/Par\_file}. This can be useful in the case of very large models. The number of pixels of the image in each direction must be smaller than parameter \texttt{NX\_NZ\_IMAGE\_MAX} defined in file \texttt{SETUP/constants.h.in}, again to avoid creating huge images in the case of very large models.}
\label{fig:example.solver}
\end{figure}
%%
@@ -641,7 +641,7 @@
\begin{itemize}
\item the \texttt{DATA/Par\_file} given with the code works fine, you can use it without any modification to test the code
-\item the seismograms \texttt{OUTPUT\_FILES/*.sem*} are simple ASCII files with two columns: time in the first colum and amplitude in the second, therefore they can be visualized with any tool you like, for instance "\texttt{gnuplot}"
+\item the seismograms \texttt{OUTPUT\_FILES/*.sem*} are simple ASCII files with two columns: time in the first column and amplitude in the second, therefore they can be visualized with any tool you like, for instance "\texttt{gnuplot}"
\item if you set flag \texttt{assign\_external\_model} to \texttt{.true.} in \texttt{DATA/Par\_file}, the velocity and density model that is given at the end of \texttt{DATA/Par\_file} is then ignored and overwritten by the external velocity and density model that you define yourself in \texttt{define\_external\_model.f90}
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90 2012-12-15 16:09:50 UTC (rev 21172)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/define_external_model.f90 2012-12-16 00:53:57 UTC (rev 21173)
@@ -42,58 +42,78 @@
!
!========================================================================
+ subroutine define_external_model(coord,material_element,ibool,rho,vp,vs,QKappa_attenuation,Qmu_attenuation, &
+ c11,c13,c15,c33,c35,c55,nspec,nglob)
- subroutine define_external_model(x,y,iflag_element,rho,vp,vs,QKappa_attenuation, &
- Qmu_attenuation,c11,c13,c15,c33,c35,c55)
-
implicit none
include "constants.h"
! user can modify this routine to assign any different external Earth model (rho, vp, vs)
-! based on the x and y coordinates of that grid point and the flag of the region it belongs to
+! based on the x and y coordinates of that grid point and the material number of the region it belongs to
- integer, intent(in) :: iflag_element
+ integer, intent(in) :: nspec,nglob
- double precision, intent(in) :: x,y
+ double precision, dimension(NDIM,nglob), intent(in) :: coord
- double precision, intent(out) :: rho,vp,vs
- double precision, intent(out) :: QKappa_attenuation,Qmu_attenuation
- double precision, intent(out) :: c11,c15,c13,c33,c35,c55
+ integer, dimension(nspec), intent(in) :: material_element
+ integer, dimension(NGLLX,NGLLZ,nspec), intent(in) :: ibool
+
+ double precision, dimension(NGLLX,NGLLZ,nspec), intent(out) :: rho,vp,vs,QKappa_attenuation,Qmu_attenuation, &
+ c11,c15,c13,c33,c35,c55
+
+ integer :: i,j,ispec,iglob
+
+ double precision :: x,z
+
! completely dummy routine here, just to demonstrate how the model can be assigned
! and how such a routine can be written
- if(iflag_element == 1 .or. x < 1700.d0 .or. y >= 2300.d0) then
- rho = 2000.d0
- vp = 3000.d0
- vs = vp / sqrt(3.d0)
- QKappa_attenuation = 9999. ! this means no attenuation
- Qmu_attenuation = 9999. ! this means no attenuation
- c11 = 169.d9
- c13 = 122.d9
- c15 = 0.d0
- c33 = c11
- c35 = 0.d0
- c55 = 75.3d9
+! loop on all the elements of the mesh, and inside each element loop on all the GLL points
+ do ispec = 1,nspec
+ do j = 1,NGLLZ
+ do i = 1,NGLLX
- else if(iflag_element == 2) then
- rho = 2500.d0
- vp = 3600.d0
- vs = vp / 2.d0
- QKappa_attenuation = 120.
- Qmu_attenuation = 120.
- c11 = 0.d0 ! this means no anisotropy
- c13 = 0.d0
- c15 = 0.d0
- c33 = 0.d0
- c35 = 0.d0
- c55 = 0.d0
+ iglob = ibool(i,j,ispec)
+ x = coord(1,iglob)
+ z = coord(2,iglob)
+
+ if(material_element(ispec) == 1 .or. x < 1700.d0 .or. z >= 2300.d0) then
+ rho(i,j,ispec) = 2000.d0
+ vp(i,j,ispec) = 3000.d0
+ vs(i,j,ispec) = vp(i,j,ispec) / sqrt(3.d0)
+ QKappa_attenuation(i,j,ispec) = 9999. ! this means no attenuation
+ Qmu_attenuation(i,j,ispec) = 9999. ! this means no attenuation
+ c11(i,j,ispec) = 169.d9
+ c13(i,j,ispec) = 122.d9
+ c15(i,j,ispec) = 0.d0
+ c33(i,j,ispec) = c11(i,j,ispec)
+ c35(i,j,ispec) = 0.d0
+ c55(i,j,ispec) = 75.3d9
+
+ else if(material_element(ispec) == 2) then
+ rho(i,j,ispec) = 2500.d0
+ vp(i,j,ispec) = 3600.d0
+ vs(i,j,ispec) = vp(i,j,ispec) / 2.d0
+ QKappa_attenuation(i,j,ispec) = 120.
+ Qmu_attenuation(i,j,ispec) = 120.
+ c11(i,j,ispec) = 0.d0 ! this means no anisotropy
+ c13(i,j,ispec) = 0.d0
+ c15(i,j,ispec) = 0.d0
+ c33(i,j,ispec) = 0.d0
+ c35(i,j,ispec) = 0.d0
+ c55(i,j,ispec) = 0.d0
+
else
- stop 'wrong flag number in external model'
+ stop 'wrong material number in external model'
endif
+ enddo
+ enddo
+ enddo
+
end subroutine define_external_model
@@ -103,31 +123,13 @@
!
!========================================================================
- subroutine define_external_model_ak135f(x,y,iflag_element,rho,vp,vs,QKappa_attenuation, &
- Qmu_attenuation,c11,c13,c15,c33,c35,c55)
+ subroutine define_external_model_ak135f(coord,material_element,ibool,rho,vp,vs,QKappa_attenuation,Qmu_attenuation, &
+ c11,c13,c15,c33,c35,c55,nspec,nglob)
implicit none
include "constants.h"
-! user can modify this routine to assign any different external Earth model (rho, vp, vs)
-! based on the x and y coordinates of that grid point and the flag of the region it belongs to
-
- integer, intent(in) :: iflag_element
-
- double precision, intent(in) :: x,y
-
- double precision, intent(out) :: rho,vp,vs
- double precision, intent(out) :: QKappa_attenuation,Qmu_attenuation
- double precision, intent(out) :: c11,c15,c13,c33,c35,c55
-
-! routine to assign the AK135F_NO_MUD Earth model
-
- integer, parameter :: IREGION_MANTLE_CRUST_ABOVE_d670 = 1
- integer, parameter :: IREGION_MANTLE_BELOW_d670 = 2
- integer, parameter :: IREGION_OUTER_CORE = 3
- integer, parameter :: IREGION_INNER_CORE = 4
-
!--------------------------------------------------------------------------------------------------
!
! taken from S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
@@ -156,25 +158,39 @@
!--------------------------------------------------------------------------------------------------
-! number of layers in ak135-F
+ integer, intent(in) :: nspec,nglob
+
+ double precision, dimension(NDIM,nglob), intent(in) :: coord
+
+ integer, dimension(nspec), intent(in) :: material_element
+
+ integer, dimension(NGLLX,NGLLZ,nspec), intent(in) :: ibool
+
+ double precision, dimension(NGLLX,NGLLZ,nspec), intent(out) :: rho,vp,vs,QKappa_attenuation,Qmu_attenuation, &
+ c11,c15,c13,c33,c35,c55
+
+! number of layers in ak135-f
integer, parameter :: NR_AK135F_NO_MUD = 136
- double precision, dimension(NR_AK135F_NO_MUD), save :: radius_ak135
- double precision, dimension(NR_AK135F_NO_MUD), save :: density_ak135
- double precision, dimension(NR_AK135F_NO_MUD), save :: vp_ak135
- double precision, dimension(NR_AK135F_NO_MUD), save :: vs_ak135
- double precision, dimension(NR_AK135F_NO_MUD), save :: Qkappa_ak135
- double precision, dimension(NR_AK135F_NO_MUD), save :: Qmu_ak135
+ double precision, dimension(NR_AK135F_NO_MUD) :: radius_ak135
+ double precision, dimension(NR_AK135F_NO_MUD) :: density_ak135
+ double precision, dimension(NR_AK135F_NO_MUD) :: vp_ak135
+ double precision, dimension(NR_AK135F_NO_MUD) :: vs_ak135
+ double precision, dimension(NR_AK135F_NO_MUD) :: Qkappa_ak135
+ double precision, dimension(NR_AK135F_NO_MUD) :: Qmu_ak135
- logical, save :: first_call_to_this_routine = .true.
+! region flags to assign the AK135F_NO_MUD Earth model
+ integer, parameter :: IREGION_MANTLE_CRUST_ABOVE_d670 = 1
+ integer, parameter :: IREGION_MANTLE_BELOW_d670 = 2
+ integer, parameter :: IREGION_OUTER_CORE = 3
+ integer, parameter :: IREGION_INNER_CORE = 4
- integer :: i
- double precision :: r,frac
+ integer :: i,j,ispec,iglob,ii
+ double precision :: x,z,r,frac
+
! define all the values in the model once and for all
- if(first_call_to_this_routine) then
-
radius_ak135( 1) = 0.000000000000000E+000
radius_ak135( 2) = 50710.0000000000
radius_ak135( 3) = 101430.000000000
@@ -1006,72 +1022,84 @@
! Qmu_ak135(133:136) = Qmu_ak135(132)
! endif
- first_call_to_this_routine = .false.
+! loop on all the elements of the mesh, and inside each element loop on all the GLL points
+ do ispec = 1,nspec
- endif ! of definition of the model the first time this routine is called
+ if(material_element(ispec) /= IREGION_MANTLE_CRUST_ABOVE_d670 .and. &
+ material_element(ispec) /= IREGION_MANTLE_BELOW_d670 .and. &
+ material_element(ispec) /= IREGION_OUTER_CORE .and. &
+ material_element(ispec) /= IREGION_INNER_CORE) stop 'wrong flag number in external model'
- if(iflag_element /= IREGION_MANTLE_CRUST_ABOVE_d670 .and. &
- iflag_element /= IREGION_MANTLE_BELOW_d670 .and. &
- iflag_element /= IREGION_OUTER_CORE .and. &
- iflag_element /= IREGION_INNER_CORE) stop 'wrong flag number in external model'
+ do j = 1,NGLLZ
+ do i = 1,NGLLX
+ iglob = ibool(i,j,ispec)
+
+ x = coord(1,iglob)
+ z = coord(2,iglob)
+
! compute the radius
- r = sqrt(x**2 + y**2)
+ r = sqrt(x**2 + z**2)
- i = 1
- do while(r >= radius_ak135(i) .and. i /= NR_AK135F_NO_MUD)
- i = i + 1
+ ii = 1
+ do while(r >= radius_ak135(ii) .and. ii /= NR_AK135F_NO_MUD)
+ ii = ii + 1
enddo
! make sure we stay in the right region and never take a point above
! and a point below the ICB or the CMB and interpolate between them,
! which would lead to a wrong value (keeping in mind that we interpolate
! between points i-1 and i below)
- if(iflag_element == IREGION_INNER_CORE .and. i > 24) i = 24
+ if(material_element(ispec) == IREGION_INNER_CORE .and. ii > 24) ii = 24
- if(iflag_element == IREGION_OUTER_CORE .and. i < 26) i = 26
- if(iflag_element == IREGION_OUTER_CORE .and. i > 69) i = 69
+ if(material_element(ispec) == IREGION_OUTER_CORE .and. ii < 26) ii = 26
+ if(material_element(ispec) == IREGION_OUTER_CORE .and. ii > 69) ii = 69
- if((iflag_element == IREGION_MANTLE_CRUST_ABOVE_d670 .or. iflag_element == IREGION_MANTLE_BELOW_d670) .and. i < 71) i = 71
+ if((material_element(ispec) == IREGION_MANTLE_CRUST_ABOVE_d670 .or. &
+ material_element(ispec) == IREGION_MANTLE_BELOW_d670) .and. ii < 71) ii = 71
- if(i == 1) then
- rho = density_ak135(i)
- vp = vp_ak135(i)
- vs = vs_ak135(i)
- Qmu_attenuation = Qmu_ak135(i)
- Qkappa_attenuation = Qkappa_ak135(i)
+ if(ii == 1) then
+ rho(i,j,ispec) = density_ak135(1)
+ vp(i,j,ispec) = vp_ak135(1)
+ vs(i,j,ispec) = vs_ak135(1)
+ Qmu_attenuation(i,j,ispec) = Qmu_ak135(1)
+ Qkappa_attenuation(i,j,ispec) = Qkappa_ak135(1)
else
-! interpolate from radius_ak135(i-1) to r using the values at i-1 and i
- frac = (r-radius_ak135(i-1))/(radius_ak135(i)-radius_ak135(i-1))
+! interpolate from radius_ak135(ii-1) to r using the values at ii-1 and ii
+ frac = (r-radius_ak135(ii-1))/(radius_ak135(ii)-radius_ak135(ii-1))
- rho = density_ak135(i-1) + frac * (density_ak135(i)-density_ak135(i-1))
- vp = vp_ak135(i-1) + frac * (vp_ak135(i)-vp_ak135(i-1))
- vs = vs_ak135(i-1) + frac * (vs_ak135(i)-vs_ak135(i-1))
- Qmu_attenuation = Qmu_ak135(i-1) + frac * (Qmu_ak135(i)-Qmu_ak135(i-1))
- Qkappa_attenuation = Qkappa_ak135(i-1) + frac * (Qkappa_ak135(i)-Qkappa_ak135(i-1))
+ rho(i,j,ispec) = density_ak135(ii-1) + frac * (density_ak135(ii)-density_ak135(ii-1))
+ vp(i,j,ispec) = vp_ak135(ii-1) + frac * (vp_ak135(ii)-vp_ak135(ii-1))
+ vs(i,j,ispec) = vs_ak135(ii-1) + frac * (vs_ak135(ii)-vs_ak135(ii-1))
+ Qmu_attenuation(i,j,ispec) = Qmu_ak135(ii-1) + frac * (Qmu_ak135(ii)-Qmu_ak135(ii-1))
+ Qkappa_attenuation(i,j,ispec) = Qkappa_ak135(ii-1) + frac * (Qkappa_ak135(ii)-Qkappa_ak135(ii-1))
endif
! make sure Vs is zero in the outer core even if roundoff errors on depth
! also set fictitious attenuation to a very high value (attenuation is not used in the fluid)
- if(iflag_element == IREGION_OUTER_CORE) then
- vs = 0.d0
- Qkappa_attenuation = 9999.d0
- Qmu_attenuation = 9999.d0
+ if(material_element(ispec) == IREGION_OUTER_CORE) then
+ vs(i,j,ispec) = 0.d0
+ Qkappa_attenuation(i,j,ispec) = 9999.d0
+ Qmu_attenuation(i,j,ispec) = 9999.d0
endif
+ enddo
+ enddo
+ enddo
+
! convert to m/s
- vp=vp*1000.0d0
- vs=vs*1000.0d0
+ vp(:,:,:)=vp(:,:,:)*1000.0d0
+ vs(:,:,:)=vs(:,:,:)*1000.0d0
! no anisotropy
- c11 = 0.d0
- c13 = 0.d0
- c15 = 0.d0
- c33 = 0.d0
- c35 = 0.d0
- c55 = 0.d0
+ c11(:,:,:) = 0.d0
+ c13(:,:,:) = 0.d0
+ c15(:,:,:) = 0.d0
+ c33(:,:,:) = 0.d0
+ c35(:,:,:) = 0.d0
+ c55(:,:,:) = 0.d0
end subroutine define_external_model_ak135f
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90 2012-12-15 16:09:50 UTC (rev 21172)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/read_external_model.f90 2012-12-16 00:53:57 UTC (rev 21173)
@@ -106,49 +106,44 @@
! need to be modified to provide QKappa & Qmu values
QKappa_attenuationext(i,j,ispec) = 10.d0
Qmu_attenuationext(i,j,ispec) = 10.d0
- end do
- end do
- end do
+ enddo
+ enddo
+ enddo
close(1001)
else
+
+ call define_external_model(coord,kmato,ibool,rhoext,vpext,vsext,QKappa_attenuationext,Qmu_attenuationext, &
+ c11ext,c13ext,c15ext,c33ext,c35ext,c55ext,nspec,nglob)
+
+! check that the external model that has just been defined makes sense
do ispec = 1,nspec
do j = 1,NGLLZ
do i = 1,NGLLX
- iglob = ibool(i,j,ispec)
- call define_external_model(coord(1,iglob),coord(2,iglob),kmato(ispec),&
- rhoext(i,j,ispec),vpext(i,j,ispec),vsext(i,j,ispec), &
- QKappa_attenuationext(i,j,ispec),Qmu_attenuationext(i,j,ispec),&
- c11ext(i,j,ispec),c13ext(i,j,ispec),c15ext(i,j,ispec), &
- c33ext(i,j,ispec),c35ext(i,j,ispec),c55ext(i,j,ispec))
-
if(c11ext(i,j,ispec) > TINYVAL .or. c13ext(i,j,ispec) > TINYVAL .or. c15ext(i,j,ispec) > TINYVAL .or. &
c33ext(i,j,ispec) > TINYVAL .or. c35ext(i,j,ispec) > TINYVAL .or. c55ext(i,j,ispec) > TINYVAL) then
! vp, vs : assign dummy values, trick to avoid floating point errors in the case of an anisotropic medium
vpext(i,j,ispec) = 20.d0
vsext(i,j,ispec) = 10.d0
- end if
+ endif
! check that the element type is not redefined compared to what is defined initially in DATA/Par_file
if((c11ext(i,j,ispec) > TINYVAL .or. c13ext(i,j,ispec) > TINYVAL .or. c15ext(i,j,ispec) > TINYVAL .or. &
c33ext(i,j,ispec) > TINYVAL .or. c35ext(i,j,ispec) > TINYVAL .or. c55ext(i,j,ispec) > TINYVAL) &
- .and. .not. anisotropic(ispec)) then
+ .and. .not. anisotropic(ispec)) &
stop 'error: non anisotropic material in DATA/Par_file or external mesh redefined as anisotropic in define_external_model()'
- endif
- if(vsext(i,j,ispec) < TINYVAL .and. (elastic(ispec) .or. anisotropic(ispec))) then
+ if(vsext(i,j,ispec) < TINYVAL .and. (elastic(ispec) .or. anisotropic(ispec))) &
stop 'error: non acoustic material in DATA/Par_file or external mesh redefined as acoustic in define_external_model()'
- endif
- if(vsext(i,j,ispec) > TINYVAL .and. .not. elastic(ispec)) then
+ if(vsext(i,j,ispec) > TINYVAL .and. .not. elastic(ispec)) &
stop 'error: acoustic material in DATA/Par_file or external mesh redefined as non acoustic in define_external_model()'
- endif
- end do
- end do
- end do
- end if
+ enddo
+ enddo
+ enddo
+ endif
! initializes
any_acoustic = .false.
Modified: seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-12-15 16:09:50 UTC (rev 21172)
+++ seismo/2D/SPECFEM2D/trunk/src/specfem2D/specfem2D.F90 2012-12-16 00:53:57 UTC (rev 21173)
@@ -2169,10 +2169,10 @@
call plotgll(knods,ibool,coorg,coord,nglob,npgeo,ngnod,nspec)
if(myrank == 0 .and. ipass == 1) &
- write(IOUT,*) 'assign_external_model = ', assign_external_model
+ write(IOUT,*) 'Assigning an external velocity and density model...'
-!if ( assign_external_model .and. ipass == 1) then
- if ( assign_external_model) then
+! if (assign_external_model .and. ipass == 1) then
+ if (assign_external_model) then
call read_external_model(any_acoustic,any_elastic,any_poroelastic, &
elastic,poroelastic,anisotropic,nspec,nglob,N_SLS,ibool, &
f0_attenuation,inv_tau_sigma_nu1_sent,phi_nu1_sent, &
More information about the CIG-COMMITS
mailing list