[mesa-users] pgstar SIGSEGV using 7184 and latest SDK

[Bill Paxton]

On Nov 28, 2014, at 12:12 PM, Marcin Mackiewicz wrote:

> Now we have another SIGSEGV very similar to the previous one. Terminal output is attached also with the modified pgstar_abundance.f file.

damn!  that's not supposed to happen!    I didn't let me test run long enough for that to happen.   Shows I need to slow down.

let's try it again.  same routine.

         integer function abundance_line(cnt, j, str)
            integer, intent(in) :: cnt, j
            character (len=*), intent(in) :: str
            integer :: i, ii
            real :: x, frac, y, dx
            include 'formats'
            call set_line_style(cnt)
            call pgslw(lw)
            call pgline(npts, xvec(grid_min:grid_max), yvec(grid_min:grid_max))
            call pgslw(lw_sav)
            call pgsch(txt_scale*s% Abundance_line_txt_scale_factor)
            abundance_line = cnt + 1
            do i=1,num_labels
               ii = iloc_abundance_label(i)
               if (ii > grid_max .or. ii < grid_min) cycle
               x = xloc_abundance_label(i)
               if (ii > grid_min) then
                  dx = xvec(ii)-xvec(ii-1)
                  if (abs(dx) > 1e-20) then
                     frac = (x-xvec(ii-1))/dx
                  else
                     frac = 1.0
                  end if
                  y = frac*yvec(ii) + (1.0-frac)*yvec(ii-1) + dylbl
               else
                  y = yvec(ii) + dylbl
               end if
               if (y < ymin .or. y > ymax) cycle
               call pgptxt(x, y, 0.0, 0.5, str)
            end do
         end function abundance_line


-Bill



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20141128/e6198638/attachment.html>