cMCSconstForMedia.f

Go to the documentation of this file.

      subroutine cmcsconstformedia(md, pm, aMCS)
!      for this media, sampling table is read and 
!      stored in aMCS.
!      It is assumed that pm = -1 comes first for a given md
!      and next call is for pm=1. So the file opened for
!      pm=-1 is kept open for pm=1. After it is read for
!     pm=1, the file will be closed.
      use modmcscontrol
      use modxsecmedia
      use modmcs
      implicit none
#include "Zmanager.h"
#include "Zmanagerp.h"
      !  CosOrEpi is used from above
      type(xsmedia),intent(in):: md  ! media 
      integer,intent(in):: pm   ! >0 for e+, < 0 for e-
      type(mcsconst),intent(out):: aMCS

      integer::icon, ie

      if( cosorepi == "epics") then
!     filename is declared in modTPXS
         if(mcsdir == " " ) then
            mcsdir ="$EPICSTOP/Data/MCS/"
         endif
!     filename="$EPICSTOP/Data/MCS/"//trim(md%name)
      else
         if(mcsdir == " " ) then
            mcsdir ="$COSMOSTOP/Data/MCS/"
         endif
!        filename="$COSMOSTOP/Data/MCS/"//trim(md%name)
      endif
      filename=trim(mcsdir)//trim(mcsparam)//"/"//trim(md%name)

      if( pm == -1 ) then
         call copenf(tempdev, filename, icon)
         if( icon /= 0 ) then
            write(0,*) trim(filename), ' cannot be opened for MCS'
            if( cosorepi == "epics") then
               write(0,*) 'If the path: '
               write(0,*)  trim(mcsdir)//trim(mcsparam)
               write(0,*)
     *         ' is correct, ',trim(md%name),' may be missing'  
               write(0,*) ' If so, you have to do next:'
               write(0,*) 
     *         '1: goto $EPICSTOP/Util/Elemag/MixedMCS'
               write(0,*) 
     *              '2: Edit paramdata file there: fix the 1st line'
               write(0,*)
     *         '   Normally only 1st term may be given; (0.05~1.9)'
               write(0,*) '3: then do'
               write(0,'(a,a,a,a)') '   ./ForManyMedia.sh ',
     *             ' $EPICSTOP/Data/MCS/',trim(mcsparam),
     *             " "//trim(md%name)
            endif
            stop
         endif
         read(tempdev, *)  ! skip 1 line  which is  C1forHardScat
         call creadmcstab(tempdev, amcs)
         read(tempdev, *)  
      elseif(pm == 1 ) then
         call creadmcstab(tempdev, amcs)
         close(tempdev)
      else
         write(0,*)' pm =',pm, ' invalid for  cMCSconstForMedia'
         stop
      endif

      amcs%loglambdah(:) =log( amcs%lambdah(:) )

      ie = amcs%minNon0mucEindex  ! for ie, 0 is stored
      amcs%loglambdas1(ie:) =log( amcs%lambdas1(ie:) )
      amcs%loglambdas2(ie:) =log( amcs%lambdas2(ie:) )
      amcs%logmuc(ie:) =log( amcs%muc(ie:) )

      end  subroutine cmcsconstformedia