ceventLoop.f

Go to the documentation of this file.

!              eventLoop
      subroutine ceventloop
      implicit none
#include  "Zmanagerp.h"
#include  "Zmanager.h"
#include  "Ztrack.h"
#include  "Zprimary.h"
#include  "Zprimaryv.h"
#include  "Zcondc.h"
#include  "Ztrackp.h"


      type(track):: incident
      integer jold, twoir(2)
      integer fin, jeof


      do while (eventno .lt. abs( desteventno(1)) .and. 
     *      eventsintherun .lt. abs( desteventno(2)))
!             switch to the random number generator 1
         call rndsw(jold, 1)
         if(job .ne. 'newflesh') then
!              for newflesh, these are managed in chookBgEvent
!c            EventsInTheRun = EventsInTheRun + 1  ! moved to cmkIncident.f
!c            EventNo = EventNo + 1                ! from  v6.25
            if(refreshir) then
!               special case. init for ir of this event is
!               read from 14 which should have been opend by the
!               user in chookBrRun
               read(14, *, end=10) twoir
               call rnd1r(twoir)
            endif
            if(job .eq. 'flesh') then
!                get random seed at skelton making
               call creadseed(seedsave, eventno, jeof)
               if( jeof .ne. 0 ) goto 10
!                 reset the seed.
               call rnd1r(seedsave)
            endif
!                 save the initial seed of random number
            call rnd1s(seedsave)
            write(0,*) seedsave, eventno+1
            call cmkincident(incident, fin)
            if(fin .ne. 0) goto 10 ! all including cuoffed end
            if(job .eq. ' ' .and. seedfile .ne. ' ') then
               call cwriteseed
            endif

            call cinitracking( incident ) ! init for tracking each event
            call cinitstack
            call cpush(incident)
         endif
         if(eabsorb(1) .ne.  0 ) then
            call chookeabsorbi(0)
         endif
         if(job .eq. 'newflesh') then
!                  this may be used always; but next one is original
            call chookbgevent   ! user hook for each event beginning
            call csetmolunit    ! for newflesh, primary is fixed above
         else
            call csetmolunit
            call chookbgevent   ! user hook for each event beginning
         endif

         call ctrackingall
         if(job .ne. 'newflesh') then
!              next is managed by chookEnEven if new flesh
            call cfintracking   ! enf of tracking of each event
         endif
         call chookenevent    ! user hook for each event  end

      enddo
 10   continue
      end
      subroutine  csetmolunit
#include "Ztrack.h"
#include "Zobs.h"
#include "Zobsp.h"
#include "Zobsv.h"

!         set  Moliere unit at each observtion depth (2r.l above along primary)
      type(track)::inci
      type(coord)::angle
      integer i
      real*8 mu
      call cqincident(inci, angle)

      do  i = 1, noofsites
!         call cgetMoliereU(ObsSites(i).pos.depth, inci.vec.coszenith, ! bef v7.
          call cgetmoliereu(obssites(i)%pos%depth, -angle%r(3),
     *    obssites(i)%mu)
      enddo
      do i = 1, noofassites
!         call cgetMoliereU(ASObsSites(i).pos.depth, inci.vec.coszenith, 
          call cgetmoliereu(asobssites(i)%pos%depth, -angle%r(3),
     *       asobssites(i)%mu)
      enddo
      end