source: JEM-EUSO/esaf_lal/tags/v1_r0/esaf/packages/simulation/externals/slast/src/init/lowtran_init.F @ 117

      SUBROUTINE LOWTRAN_INIT
#include "detector.inc"
#include "paw.inc"
#include "userlun.common"
#include "lwtrn7.common"
#include "tuple.common"
*
      REAL H_step, theta_step, H1C,H2C
      REAL TRANS(100)
      INTEGER J,JINV
*
      IFILE  = 8
      ITAPE5 = 9
      IPRNT  = 10
      IFILE2 = 18
*     ... atmosphere model initialisation
      CALL LWTRANS(0,0.,0.,TRANS)
      CALL SYSTEM("ln -sf /dev/null TAPE6")
      CALL SYSTEM("ln -sf /dev/null TAPE7")
      CALL SYSTEM("ln -sf /dev/null TAPE8")

      IF(LWTR.EQ.0) GOTO 1000
      H_step = 1. !1.e-1           ! km
      theta_step = 10. !1.e0        ! deg.
      H1C = 0
      H2C = 100
      NH  = INT((H2C-H1C)/H_step)
      NA  = INT(90/theta_step)
      k = 0
*
      CALL BOOK_LOWTRAN_NTUPLE
      Ntot      = NA*NH
      i_portion = Ntot/10
      WRITE(6,*) ' '
      WRITE(6,98)
      WRITE(6,99) 'NTUPLE     '

      DO H = H1C, H2C, H_step 
        DO ANGL = 0., 90., theta_step
          k = k + 1      
          CALL LWTRANS(1,H,ANGL,TRANS)
          ZA   = H
          Teta = ANGL
          CALL HFNT(10)
          ntmp = k/i_portion
          ntmp = k - ntmp*i_portion
          IF(ntmp.EQ.0) WRITE(6,100)
        ENDDO
      ENDDO
      CALL HROUT(10,Icycle,' ')
      CALL HREND('LOWTRAN')
 98   FORMAT(1x,'COMPUTING THE LOWTRAN TABLE: ')
 99   FORMAT($,1x,'SAVE ',A,1x,'(breath deeply)')
 100  FORMAT($,1x,'.')
 101  FORMAT($,1x,'[done]')
 1000 RETURN
 1001 FORMAT(3x,f8.1,f12.2)
      END