<div dir="ltr"><br><div class="gmail_extra">HI, <br></div><div class="gmail_extra">Thanks to Paul and Nerses for the feedback.<br><br></div><div class="gmail_extra">I compare the same file in <br>
/group/clas/builds/centos62/trunk/reconstruction/ana/fill_seb_ntn.F<br><br></div><div class="gmail_extra">The way to extract bit 31 from iw(ind3+1) is exactly the same in both version.<br>
<br><br></div><div class="gmail_extra">ifarm1102> grep -i 'stbit' /group/clas/builds/centos62/trunk/reconstruction/ana/fill_seb_ntn.F ../ana/fill_seb_ntn.F<br>/group/clas/builds/centos62/trunk/reconstruction/ana/fill_seb_ntn.F: integer hbit,sbit,htbit,stbit,sbitold,hbitold,j,sect<br>
/group/clas/builds/centos62/trunk/reconstruction/ana/fill_seb_ntn.F: stbit=3<br>/group/clas/builds/centos62/trunk/reconstruction/ana/fill_seb_ntn.F: stbit=int(iw(ind3+1)/2**31)<br>/group/clas/builds/centos62/trunk/reconstruction/ana/fill_seb_ntn.F: htbit=int((iw(ind3+1)-stbit*2**31)/2**30)<br>
../ana/fill_seb_ntn.F: integer hbit,sbit,htbit,stbit,sbitold,hbitold,j,sect<br>../ana/fill_seb_ntn.F: stbit=3<br>../ana/fill_seb_ntn.F: stbit=int(iw(ind3+1)/2**31)<br>../ana/fill_seb_ntn.F: htbit=int((iw(ind3+1)-stbit*2**31)/2**30)<br>
<br><br>Nerse's suggestion of using preprocessor commands should solve the problem.<br>Here is my modification only for gfortran.<br> <br>#ifdef __GFORTRAN__<br> stbit=ibits(iw(ind3+1),31,1)<br> htbit=ibits(iw(ind3+1),30,1)<br>
#else<br> stbit=int(iw(ind3+1)/2**31)<br> htbit=int((iw(ind3+1)-stbit*2**31)/2**30)<br>#endif<br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Fri, Aug 30, 2013 at 2:11 PM, Nerses Gevorgyan <<a href="mailto:nerses@jlab.org">nerses@jlab.org</a>> wrote:<br><blockquote class="gmail_quote">
<div bgcolor="#FFFFFF" text="#000000">
<div>I agree with Paul it is very possible
somebody fixed that somewhere.<br>
<br>
Just in case if You'll not find it.<br>
-----------------------------------------<br>
stbit=int(iw(ind3+1)/2**31)<br>
htbit=int((iw(ind3+1)-stbit*2**31)/2**30)<br>
<br>
These type of commands are happening in more than one place.<br>
I guess it is just first one encountered.<br>
It is just a way of extracting some bits(#32 and #31) from 32 bit
from integer iw.<br>
<br>
I've used a replacement code using ibits intrinsic function
(available for FORTRAN95 and later versions).<br>
stbit=ibits(iw(ind3+1),31,1)<br>
htbit=ibits(iw(ind3+1),30,1)<br>
Then in the code using preprocessor commands use appropriate code
depending on g77 or gfortran compiler is used.<br>
<br>
HTH,<br>
Nerses<div><div class="h5"><br></div></div></div></div></blockquote></div>-- <br>With Best Regards,<br>Jixie Zhang<br>757-269-7735 <br>_______________________________________________
</div></div>