cg Version
cg Version
cg Version
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
! ----------------------------------------------------------------! Print out stack information (for limited number cases and lines)! ----------------------------------------------------------------if (ncount .le. nwrite .and. ilines .le. nlines) thenilines = ilines + 1write(6,100) e(np),x(np),y(np),z(np),u(np),v(np),w(np),* iql,irl,iarg100 FORMAT(7G15.7,3I5)end if! ------------------------------------! Output particle information for plot! ------------------------------------if (ncount.le.maxpict) thencall plotxyz(iarg,np,iq(np),x(np),y(np),z(np),e(np),ir(np),* wt(np),time(np))end ifreturnend!--------------------------last line of ausgab.f------------------------!-------------------------------howfar.f--------------------------------! <strong>Version</strong>: 070627-1600! Reference: T. Torii and T. Sugita, "Development of PRESTA-CG! Incorporating Combinatorial Geometry in EGS4/PRESTA", JNC TN1410 2002-201,! Japan Nuclear Cycle Development Institute (2002).! Improved version is provided by T. Sugita. 7/27/2004!-----------------------------------------------------------------------!23456789|123456789|123456789|123456789|123456789|123456789|123456789|12! ----------------------------------------------------------------------! Required (geometry) subroutine for use with the EGS5 Code System! ----------------------------------------------------------------------! This is a CG-HOWFAR.! ----------------------------------------------------------------------ccccccccccsubroutine howfarimplicit noneinclude 'include/egs5_h.f' ! Main EGS "header" fileinclude 'include/egs5_epcont.f' ! COMMONs required by EGS5 codeinclude 'include/egs5_stack.f'include 'auxcommons/geom_common.f' ! geom-common fileinteger i,j,jjj,ir_np,nozone,jty,knointeger irnear,irnext,irlold,irlfg,itvlfg,ihit<strong>cg</strong>double precision xidd,yidd,zidd,x_np,y_np,z_np,u_np,v_np,w_npdouble precision tval,tval0,tval00,tval10,tvalmn,delhowdouble precision atvaltmpinteger iq_npir_np = ir(np)iq_np = iq(np) + 2if(ir_np.le.0) thenwrite(6,*) 'Stopped in howfar with ir(np) izonin'stopend ifif(ir_np.EQ.izonin) thenidisc=1returnend iftval=1.d+30itvalm=0body checku_np=u(np)v_np=v(np)w_np=w(np)60
Change language