[mesa-users] pgstar SIGSEGV using 7184 and latest SDK
Marcin Mackiewicz
astrophysics at o2.pl
Sat Nov 29 03:22:50 EST 2014
Bill, it worked!
Thanks for the patch.
M.
Dnia 28 listopada 2014 22:44 Bill Paxton <paxton at kitp.ucsb.edu> napisał(a):
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/20141129/822cf4f8/attachment.html>
More information about the Mesa-users
mailing list