csmpColInA2.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	csmpcolina2 (pj, ia, nc)
subroutine	csampcolln (pj, prob, np, nx, nc)

Function/Subroutine Documentation

csampcolln()

subroutine csampcolln	(	type(ptcl)	pj,
		real*4, dimension(np, nx)	prob,
		integer	np,
		integer	nx,
		integer	nc
	)

Definition at line 118 of file csmpColInA2.f.

References cinelx(), d0, parameter(), and rndc().

Referenced by csmpcolina2().

      implicit none

#include  "Zptcl.h"

        type(ptcl):: pj  !input projectile particle
        integer np       ! input table size for prob.
        integer nx       ! input. table size for xsection.
        integer nc       ! output number of collisions.
        real*4  prob(np, nx)  ! table for n-coll. prob.

        real*8  xs
!
        integer   idxxs, i

        real*8 u
        real*8  xs1, xs2, dxs
        integer nn
        parameter(nn = 46)  ! this should be nx
!         parameter ( xs1 =log10(15.), xs2 =log10(320.) )
        parameter( xs1 = 1.176091259d0, xs2 =  2.505149978d0,
     *  dxs =( xs2 - xs1)/(nn-1))
!
        if(nx .ne. nn) then
           write(*,*) ' error input to csampCollN'
           stop 9999
        endif
!             get cross-section for proton target.        
!        call cxpXsec(pj, xs)
        call cinelx(pj, 1.d0, 1.d0,  xs)
        xs = log10(xs)
!
        idxxs = (xs - xs1)/dxs + 1
        if( (xs - idxxs * dxs - xs1) .gt.
     *        (idxxs * dxs + dxs + xs1 - xs) ) then
            idxxs = idxxs +1
         endif
         idxxs =max(1, min(idxxs, nn))

         call rndc(u)
         do   i=1,  np
            if(u .le. prob(i,  idxxs) ) then
               nc=i
               goto 100
            endif
         enddo
 100     continue

Here is the call graph for this function:

Here is the caller graph for this function:

csmpcolina2()

subroutine csmpcolina2	(	type(ptcl)	pj,
		integer	ia,
		integer	nc
	)

Definition at line 3 of file csmpColInA2.f.

References cerrormsg(), cinelx(), csampcolln(), d0, parameter(), and rndc().

Referenced by chacoladhoc().

        implicit none
#include  "Zptcl.h"
        type(ptcl):: pj  !input projectile particle
        integer ia       ! target mass no.
        integer nc       ! output number of collisions.
        call cerrormsg('*** IBM AIX must use SucInt=0 ***', 0)

Here is the call graph for this function:

Here is the caller graph for this function: