source: JEM-EUSO/esaf_lal/tags/v1_r0/esaf/packages/simulation/externals/lowtran/src/esaflwtrans_v2.F @ 117

Last change on this file since 117 was 117, checked in by moretto, 11 years ago

ESAF version compilable on mac OS

File size: 1.7 KB
Line 
1      SUBROUTINE ESAFLWTRANS_V2(INIT,USR_H1,USR_H2,USR_ANGLE)
2c**********************************************************************
3c     Subroutine to calcul the radiative transfert throw the atmosphere
4c     Based on the LOWTRAN7 program
5c
6c     Input parameters :
7c     ------------------
8c     INIT :  flag for initialisation of lowtran
9C             0 initialisation of atmosphere's parameters and conditions
10c                 ( defined in the lowtran.input file )
11c             1 initialisation of geometrical conditions only when
12c               atmosphere's parameters have already been defined
13c
14c     USR_H1, USR_H2 : altitude of the 2 points in km
15c          USR_ANGLE : angle theta compare to Nadir (0=vertical)
16c                        origine on ground at Nadir (EUSO = 0.,0.,400.)
17c
18c     Output parameters :
19c     -------------------
20c     compare to LWTRANS, all the ouputs are available via the
21c     common block "tuple"
22c
23c     Author : Anne Stutz
24c     Last revised by Corinne Berat, April 2004
25c
26c***********************************************************************
27#include "lwtrn7.common"
28#include "userlun.common"
29#include "tuple.common"
30
31      INTEGER INIT
32      REAL USR_H1,USR_H2,USR_ANGLE
33      IFILE  = 8
34      ITAPE5 = 9
35      IPRNT  = 10
36      IFILE2 = 18
37     
38c     Atmosphere Model Initialisation
39      IF(INIT.EQ.0) THEN
40         NEVT=1
41         ILAMBDA=0
42         CALL INITCARD3_V2(0.,0.,0.)
43         CALL LWTRN7SUB
44      ENDIF   
45         
46c     Calcul the atmospheric transmission from xpos,ypos,zpos to EUSO
47      IF(INIT.EQ.1) THEN
48         NEVT = NEVT+1
49         ILAMBDA=0
50         CALL INITCARD3_V2(USR_H1,USR_H2,USR_ANGLE)
51         CALL LWTRN7SUB
52      ENDIF
53c
54      RETURN
55      END
56
57
58         
59
60
61
62
63
64
65
66
67
Note: See TracBrowser for help on using the repository browser.