<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">HI,<div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 25 Oct 2015, at 10:44, Rick Walsh <<a href="mailto:rickmwalsh@gmail.com" class="">rickmwalsh@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">I like this concept. I'm not sure how to turn it into something meaningful (other than the single statistic of "was the free gas volume more or less than permitted"), but it's worth discussing.</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""></div></blockquote><div><br class=""></div></div>as a first step in this direction, I added a number of printf’s to the VPM-B part. In particular, this displays the allowable (bottom) gradients for all compartments and for each iteration of the CVA how the volume of free gas compares to the critical value.</div><div class=""><br class=""></div><div class="">This is of course not production code, it should not go into master. But for all those who want to understand VPM-B a but better, I pushed it to github. You can pull it form </div><div class=""><br class=""></div><div class=""><a href="https://github.com/atdotde/subsurface.git" class="">https://github.com/atdotde/subsurface.git</a> analyze_vpm</div><div class=""><br class=""></div><div class="">Have fun.</div><div class=""><br class=""></div><div class="">What I learned so far: There are typically three rounds of the CVA: one with the vpmb_start_gradient, one with an updated gradient and then finally one with almost the same gradient again that results in the same run time. So it converges very fast in essentially one step of updating.</div><div class=""><br class=""></div><div class="">Best</div><div class="">Robert</div></body></html>