<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Jordi,<div><br></div><div>To print the station code of the intensities as they are flagged, edit the program 'grind' and go to around line 2906 (this may vary somewhat depending on what revision of ShakeMap you are running -- if you have trouble locating the relevant lines, search for 'sub flag_stations_by_ipe' and scroll down a few lines). You should see something that looks like this:</div><div><br></div><div><div> for($i = 0; $i < $n; $i++) {</div><div> if (defined $obs->[$i] and $obs->[$i] !~ /NaN/i) {</div><div> #---------------------------------------------------------------------</div><div> # Beyond the GMPE's fitted range, we want to toss large data points</div><div> # but not necessarily small ones</div><div> #---------------------------------------------------------------------</div><div> $d = $obs->[$i] - $ests->[$i];</div><div> $mi_sd = $misda->[$i] * $outlier_level;</div><div> if ($d > $mi_sd or ($d < -$mi_sd and $dists->[$i] < $bias_max_range)) {</div><div> $sta = $stahp->{$cc->[$i]};</div><div> $flagged{$sta->code}{intensity} = 1;</div><div> $obs->[$i] = 'NaN';</div><div> $nflagged++;</div><div> }</div><div> }</div><div> }</div><div><br></div><div>Add a line something like this:</div><div><br></div><div><div> for($i = 0; $i < $n; $i++) {</div><div> if (defined $obs->[$i] and $obs->[$i] !~ /NaN/i) {</div><div> #---------------------------------------------------------------------</div><div> # Beyond the GMPE's fitted range, we want to toss large data points</div><div> # but not necessarily small ones</div><div> #---------------------------------------------------------------------</div><div> $d = $obs->[$i] - $ests->[$i];</div><div> $mi_sd = $misda->[$i] * $outlier_level;</div><div> if ($d > $mi_sd or ($d < -$mi_sd and $dists->[$i] < $bias_max_range)) {</div><div> $sta = $stahp->{$cc->[$i]};</div><div> $flagged{$sta->code}{intensity} = 1;</div><div>printf "Flagging station '%s'\n", $sta->code;</div><div> $obs->[$i] = 'NaN';</div><div> $nflagged++;</div><div> }</div><div> }</div></div><div><br></div><div>That will print a line like:</div><div><br></div><div>Flagging station 'JZM'</div><div><br></div><div>Each time a station is flagged. Modify the printf line as you see fit.</div><div><br></div><div>Bruce</div><div><br></div><div><div>On Oct 18, 2012, at 5:56 AM, Roviró Pou, Jordi wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div lang="CA" link="blue" vlink="purple"><div class="WordSection1" style="page: WordSection1; "><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB">Hello,<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB">I need some help to know how can I do to show the flagged stations in every iteration, only for IPE, in the grind module.<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB">I try to use the Print option into the grind.pm code, but it’s so difficult for me (I've never worked with perl language.)<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB">My idea is to know which stations are flagged in each iteration to understand how it works the bias methodology.<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; "><o:p> </o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">The –verbose option explains to me:<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; "><o:p> </o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><i><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">Using IPE::Isard2008 for IPE<o:p></o:p></span></i></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><i><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">Intensity Bias=-0.08<o:p></o:p></span></i></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><i><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">24 stations flagged out this iteration<o:p></o:p></span></i></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><i><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">Intensity Bias=-0.05<o:p></o:p></span></i></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><i><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">0 stations flagged out this iteration<o:p></o:p></span></i></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><i><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">Uncertainty computation skipped (no Intensity > 6.0)<o:p></o:p></span></i></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; "><o:p> </o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">I want to know which stations are.<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; ">Thank you very much!<o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB" style="font-family: 'Times New Roman', serif; "><o:p> </o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span style="font-size: 10pt; font-family: Verdana, sans-serif; color: black; "><o:p> </o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span style="font-size: 10pt; font-family: Verdana, sans-serif; color: black; "><br></span><span style="font-family: Verdana, sans-serif; color: black; "><span><image002.jpg></span><o:p></o:p></span></div><div class="MsoNormal" align="center" style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; text-align: center; "><span style="font-family: Verdana, sans-serif; color: black; "><hr size="1" width="100%" noshade="" align="center" style="color: rgb(192, 0, 0); "></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 37.5pt; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><b><span style="font-size: 8pt; font-family: Verdana, sans-serif; color: black; ">Jordi Roviró i Pou</span></b><span style="font-size: 8pt; font-family: Verdana, sans-serif; color: black; "><br>Unitat de Risc Sísmic<br>Institut Geològic de Catalunya<br><br><a href="http://mercuri.icc.cat/website/mob_nf/mob1/mob2/inici2.htm?&zoom=3&XY=429136:4583504&ident=%0d%0aInstitut%20Geologic%20de%20Catalunya%3cbr%3e%20c/Balmes,%20209-211%3cbr%3e%20E-08006%20Barcelona%3cbr%3e%20Telefon%2093%20553%2084%2030&titol=IGC&eina=2&base=m&adreca=Balmes,%20209" target="_blank" style="color: blue; text-decoration: underline; "><span style="color: black; text-decoration: none; ">Balmes, 209-211 | E-08006 Barcelona</span></a><span class="Apple-converted-space"> </span>| Tel. 34 93 553 84 30 (2418)<br></span><span style="font-size: 8pt; font-family: Verdana, sans-serif; color: navy; "><a href="mailto:mail@igc.cat" target="_blank" style="color: blue; text-decoration: underline; "><span style="color: navy; text-decoration: none; ">jroviro@igc.cat</span></a></span><span style="font-size: 8pt; font-family: Verdana, sans-serif; color: black; "><span class="Apple-converted-space"> </span>|<span class="Apple-converted-space"> </span></span><span style="font-size: 8pt; font-family: Verdana, sans-serif; color: navy; "><a href="http://www.igc.cat/" target="_blank" style="color: blue; text-decoration: underline; "><span style="color: navy; text-decoration: none; ">www.igc.cat</span></a></span><span style="font-size: 8pt; font-family: Verdana, sans-serif; color: black; "><o:p></o:p></span></div><div style="margin-top: 0cm; margin-right: 0cm; margin-left: 0cm; margin-bottom: 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; "><span lang="EN-GB"><o:p> </o:p></span></div></div>_______________________________________________<br>Shake-dev mailing list<br><a href="mailto:Shake-dev@geohazards.usgs.gov" style="color: blue; text-decoration: underline; ">Shake-dev@geohazards.usgs.gov</a><br><a href="https://geohazards.usgs.gov/mailman/listinfo/shake-dev" style="color: blue; text-decoration: underline; ">https://geohazards.usgs.gov/mailman/listinfo/shake-dev</a><br></div></blockquote></div><br></div></body></html>