From e789a2b7af803f01e6c314febb39cf7f41795443 Mon Sep 17 00:00:00 2001 From: Piotr Date: Thu, 28 Mar 2019 14:32:31 +0100 Subject: [PATCH] Fixed bug in ptc_twiss that did not add deltap to the user orbit, but just overwrote it with deltap --- src/madx_ptc_twiss.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/madx_ptc_twiss.f90 b/src/madx_ptc_twiss.f90 index a0aae6401..107614296 100644 --- a/src/madx_ptc_twiss.f90 +++ b/src/madx_ptc_twiss.f90 @@ -716,14 +716,14 @@ subroutine ptc_twiss(tab_name,summary_tab_name) orbit(3)=get_value('ptc_twiss ','y ') orbit(4)=get_value('ptc_twiss ','py ') orbit(6)=-get_value('ptc_twiss ','t ') ! swap of t sign - orbit(5)=orbit(5)+get_value('ptc_twiss ','pt ') + orbit(5)=get_value('ptc_twiss ','pt ') if(mytime) then call Convert_dp_to_dt (deltap, dt) else dt=deltap endif - if(icase.eq.5 .or. icase.eq.56) orbit(5)=dt + if(icase.eq.5 .or. icase.eq.56) orbit(5) = dt + orbit(5) closed_orbit = get_value('ptc_twiss ','closed_orbit ') .ne. 0