showSpec.f

Go to the documentation of this file.

#include  "BlockData/cblkGene.h"
      implicit none
#include "Zglobalc.h"
#include "Zmanagerp.h"
#include "ZrigCut.h"
#include "Zptcl.h"
#include "Zobs.h"
#include "Zobsp.h"
#include "Zprimary.h"
#include "Zprimaryc.h"
#include "Zprimaryv.h"
#include "Zincidentp.h"
#ifdef NEXT486
#define IMAG_P dimag
#else
#define IMAG_P imag
#endif


      integer  i, j

      real*8 cosmin, cosmax, rigc
      real*8 flux, E, azmmin, azmmax, dE, zen1, zen2


      call creadparam(5)
      call cbeginrun
      call cprintprim(errorout)
      cosmax = imag_p(coszenith)
      cosmin = real(coszenith)
      azmmax = imag_p(azimuth) + xaxisfromsouth
      azmmin = real(Azimuth) + XaxisFromSouth
      de=10.**0.05
      if(cutofffile .eq. ' ') then
         rigc = 0.
         call csetcosdeg(.false., .false.)
      else
         rigc = 100.
         if(zenvalue .eq. 'deg') then
            call csetcosdeg(.false., .true.)
         else
            call csetcosdeg(.false., .false.)
         endif
      endif
      zen1 = cosmin
      zen2 = cosmax
      do i = 1, prim%no_of_comps
         e = prim%each(i)%energy(1)
         do j = 1, 10000
            call cprimflux(prim%each(i), e, rigc, zen1, zen2,
     *      azmmin, azmmax,  flux)
!            if(flux .gt. 0.) then
               write(*,*)
     *         prim%each(i)%label, sngl(e), sngl(flux)
!            endif
            e = e*de
            if(e .gt. prim%each(i)%energy(prim%each(i)%no_of_seg+1))
     *          goto 100            
         enddo
 100     continue
!         write(*,*)
      enddo
      end
      subroutine chooktrace
      end
      subroutine chookceren
      end
      subroutine chookcerens
      end
      subroutine chookcerene
      end
      subroutine chookbgrun
      end