source: trunk/MagicSoft/Simulation/Detector/Reflector/attenu.f@ 428

Last change on this file since 428 was 428, checked in by harald, 24 years ago
During the Monte Carlo Software meeting, a new version of attenu.f was promised by Jose Carlos, Frank Schroeder and Aitor Ibarra. I got this version for the repository. A small change in the Makefile was also neccessary. Thanks to Jose Carlos for his support to run it on alpha!
File size: 22.1 KB
Line 
1*********************************************************************
2* *
3* Atmospheric Absorption for Rayleigh, Mie and Ozone. *
4* *
5* Created: May-98 *
6* Author: Aitor Ibarra Ibaibarriaga *
7* Jose Carlos Gonzalez *
8* *
9*********************************************************************
10
11c @T \newpage
12
13
14c @section Source code of {\tt attenu.f}
15
16* @text
17* Copyright $\copyright$ 1998, Aitor Ibarra Ibaibarriaga
18* @endtext
19c @code
20
21 SUBROUTINE attenu(wavelength, height, theta, tr_atmos)
22
23C----------------------------------------------------------------------C
24C RHO (DENSITY) F(UNCTION) C
25C C
26C CALCULATES DENSITY (G/CM**3) OF ATMOSPHERE DEPENDING ON HEIGHT (CM) C
27C (US STANDARD ATMOSPHERE) C
28C THIS FUNCTION IS CALLED FROM ININKG, UPDATE, CERENE, CERENH C
29C ARGUMENT: C
30C ARG = HEIGHT IN CM C
31C----------------------------------------------------------------------C
32
33 COMMON /ATMOS/ AATM,BATM,CATM,LAHG,RHOF,LONG,OZ_ABI,AE_ABI
34 DOUBLE PRECISION AATM(5),BATM(5),CATM(5),LAHG(5),RHOF(5),Lm(12)
35 DOUBLE PRECISION H,RT, m,OZ_ABI(48,12),AE_ABI(28,12)
36 DOUBLE PRECISION XR, TOT_OZ, TOT_AE, T_Ray,T_Mie, T_Oz,
37 + RHO_TOT, RHO_FI, RHOFP
38
39 DOUBLE PRECISION Rcos2, Rsin2
40 DOUBLE PRECISION x1, x2, x3, x4
41 DOUBLE PRECISION e1, e2, e3, e4
42
43 REAL wavelength, height, theta, tr_atmos
44 real trr, trm, tro
45 REAL LONG(12)
46c fs: define obervation level
47 double precision obslev
48 INTEGER I,CON_OZ,CON_MI J, ROW
49
50 DATA OZ_ABI /
51 + 0.2880000D+00,0.5405000D+00,0.7775000D+00,0.1009000D+01,
52 + 0.1241500D+01,0.1480500D+01,0.1750500D+01,0.2085000D+01,
53 + 0.2514500D+01,0.3087500D+01,0.3864500D+01,0.4817500D+01,
54 + 0.5847500D+01,0.6917500D+01,0.8052500D+01,0.9287499D+01,
55 + 0.1068750D+02,0.1231250D+02,0.1415750D+02,0.1617750D+02,
56 + 0.1827250D+02,0.2034750D+02,0.2232750D+02,0.2414750D+02,
57 + 0.2575750D+02,0.2715250D+02,0.2836750D+02,0.2941100D+02,
58 + 0.3031000D+02,0.3109250D+02,0.3176300D+02,0.3232750D+02,
59 + 0.3281200D+02,0.3323200D+02,0.3358350D+02,0.3387750D+02,
60 + 0.3412650D+02,0.3434000D+02,0.3451900D+02,0.3466250D+02,
61 + 0.3477480D+02,0.3486355D+02,0.3493355D+02,0.3498775D+02,
62 + 0.3503010D+02,0.3506360D+02,0.3509020D+02,0.3511185D+02,
63 + 0.2740000D-01,0.5140000D-01,0.7395000D-01,0.9600000D-01,
64 + 0.1181500D+00,0.1409000D+00,0.1666000D+00,0.1984500D+00,
65 + 0.2393500D+00,0.2939500D+00,0.3679500D+00,0.4589500D+00,
66 + 0.5573000D+00,0.6593000D+00,0.7673000D+00,0.8848000D+00,
67 + 0.1017800D+01,0.1172300D+01,0.1348300D+01,0.1540800D+01,
68 + 0.1740300D+01,0.1937800D+01,0.2126300D+01,0.2299800D+01,
69 + 0.2453300D+01,0.2586300D+01,0.2702300D+01,0.2801900D+01,
70 + 0.2887550D+01,0.2962050D+01,0.3025900D+01,0.3079700D+01,
71 + 0.3125850D+01,0.3165850D+01,0.3199350D+01,0.3227400D+01,
72 + 0.3251150D+01,0.3271500D+01,0.3288600D+01,0.3302300D+01,
73 + 0.3312995D+01,0.3321445D+01,0.3328110D+01,0.3333270D+01,
74 + 0.3337305D+01,0.3340500D+01,0.3343035D+01,0.1334510D+01,
75 + 0.2435000D-02,0.4570000D-02,0.6575000D-02,0.8535000D-02,
76 + 0.1050500D-01,0.1253000D-01,0.1481500D-01,0.1764500D-01,
77 + 0.2128000D-01,0.2613500D-01,0.3272000D-01,0.4081000D-01,
78 + 0.4957000D-01,0.5866000D-01,0.6827000D-01,0.7875500D-01,
79 + 0.9065500D-01,0.1044050D+00,0.1200050D+00,0.1371050D+00,
80 + 0.1548550D+00,0.1724050D+00,0.1891550D+00,0.2045550D+00,
81 + 0.2182050D+00,0.2300550D+00,0.2403600D+00,0.2492200D+00,
82 + 0.2568350D+00,0.2634550D+00,0.2691300D+00,0.2739150D+00,
83 + 0.2780200D+00,0.2815750D+00,0.2845500D+00,0.2870400D+00,
84 + 0.2891500D+00,0.2909600D+00,0.2924750D+00,0.2936900D+00,
85 + 0.2946425D+00,0.2953940D+00,0.2959865D+00,0.2964455D+00,
86 + 0.2968045D+00,0.2970885D+00,0.2973140D+00,0.2974975D+00,
87 + 0.1740000D-03,0.3265000D-03,0.4695000D-03,0.6090000D-03,
88 + 0.7495000D-03,0.8940000D-03,0.1057000D-02,0.1259000D-02,
89 + 0.1518000D-02,0.1863500D-02,0.2332500D-02,0.2909000D-02,
90 + 0.3533000D-02,0.4180500D-02,0.4865000D-02,0.5610500D-02,
91 + 0.6455500D-02,0.7435000D-02,0.8550000D-02,0.9770000D-02,
92 + 0.1103500D-01,0.1229000D-01,0.1348500D-01,0.1458000D-01,
93 + 0.1555100D-01,0.1639550D-01,0.1713150D-01,0.1776300D-01,
94 + 0.1830600D-01,0.1877800D-01,0.1918200D-01,0.1952250D-01,
95 + 0.1981500D-01,0.2006850D-01,0.2028050D-01,0.2045801D-01,
96 + 0.2060851D-01,0.2073750D-01,0.2084566D-01,0.2093241D-01,
97 + 0.2100026D-01,0.2105381D-01,0.2109606D-01,0.2112876D-01,
98 + 0.2115431D-01,0.2117456D-01,0.2119066D-01,0.2120376D-01,
99 + 0.4885000D-05,0.9170000D-05,0.1319500D-04,0.1713000D-04,
100 + 0.2108000D-04,0.2513500D-04,0.2971500D-04,0.3539500D-04,
101 + 0.4268500D-04,0.5242500D-04,0.6562500D-04,0.8182500D-04,
102 + 0.9937500D-04,0.1175750D-03,0.1368250D-03,0.1578250D-03,
103 + 0.1816250D-03,0.2091750D-03,0.2404750D-03,0.2747750D-03,
104 + 0.3103250D-03,0.3454750D-03,0.3790250D-03,0.4098750D-03,
105 + 0.4372250D-03,0.4609750D-03,0.4816750D-03,0.4994750D-03,
106 + 0.5147750D-03,0.5280750D-03,0.5394750D-03,0.5490700D-03,
107 + 0.5572950D-03,0.5644250D-03,0.5703950D-03,0.5753899D-03,
108 + 0.5796200D-03,0.5832500D-03,0.5862950D-03,0.5887350D-03,
109 + 0.5906400D-03,0.5921450D-03,0.5933350D-03,0.5942565D-03,
110 + 0.5949755D-03,0.5955440D-03,0.5959955D-03,0.5963635D-03,
111 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
112 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
113 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
114 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
115 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
116 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
117 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
118 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
119 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
120 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
121 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
122 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
123 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
124 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
125 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
126 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
127 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
128 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
129 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
130 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
131 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
132 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
133 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
134 + 0.0000000D+00,0.0000000D+00,0.0000000D+00,0.0000000D+00,
135 + 0.9525000D-05,0.1785500D-04,0.2567500D-04,0.3332000D-04,
136 + 0.4100000D-04,0.4889000D-04,0.5779500D-04,0.6881000D-04,
137 + 0.8296000D-04,0.1018600D-03,0.1275100D-03,0.1590600D-03,
138 + 0.1932100D-03,0.2286100D-03,0.2660100D-03,0.3067600D-03,
139 + 0.3529600D-03,0.4065100D-03,0.4674100D-03,0.5341100D-03,
140 + 0.6032600D-03,0.6717100D-03,0.7370100D-03,0.7970100D-03,
141 + 0.8501600D-03,0.8963600D-03,0.9366100D-03,0.9711101D-03,
142 + 0.1000810D-02,0.1026660D-02,0.1048810D-02,0.1067460D-02,
143 + 0.1083460D-02,0.1097310D-02,0.1108910D-02,0.1118640D-02,
144 + 0.1126865D-02,0.1133915D-02,0.1139830D-02,0.1144570D-02,
145 + 0.1148275D-02,0.1151200D-02,0.1153510D-02,0.1155300D-02,
146 + 0.1156700D-02,0.1162205D-02,0.1170990D-02,0.1178145D-02,
147 + 0.9360000D-04,0.1757000D-03,0.2528000D-03,0.3281500D-03,
148 + 0.4038500D-03,0.4816500D-03,0.5694500D-03,0.6784000D-03,
149 + 0.8184000D-03,0.1004900D-02,0.1257900D-02,0.1568900D-02,
150 + 0.1905400D-02,0.2254400D-02,0.2623400D-02,0.3025400D-02,
151 + 0.3480900D-02,0.4008900D-02,0.4609400D-02,0.5266900D-02,
152 + 0.5948400D-02,0.6622900D-02,0.7266400D-02,0.7857900D-02,
153 + 0.8381900D-02,0.8836901D-02,0.9233401D-02,0.9573901D-02,
154 + 0.9866901D-02,0.1012140D-01,0.1033940D-01,0.1052340D-01,
155 + 0.1068140D-01,0.1081840D-01,0.1093290D-01,0.1102855D-01,
156 + 0.1110965D-01,0.1117920D-01,0.1123750D-01,0.1128425D-01,
157 + 0.1132085D-01,0.1134975D-01,0.1137255D-01,0.1139020D-01,
158 + 0.1140400D-01,0.1141492D-01,0.1142358D-01,0.1143063D-01,
159 + 0.2500000D-03,0.4690000D-03,0.6745000D-03,0.8755000D-03,
160 + 0.1077500D-02,0.1285000D-02,0.1519500D-02,0.1810000D-02,
161 + 0.2182500D-02,0.2679500D-02,0.3353500D-02,0.4182000D-02,
162 + 0.5079000D-02,0.6010000D-02,0.6994000D-02,0.8064000D-02,
163 + 0.9279000D-02,0.1068900D-01,0.1228900D-01,0.1403900D-01,
164 + 0.1585400D-01,0.1765400D-01,0.1937400D-01,0.2095400D-01,
165 + 0.2235400D-01,0.2356900D-01,0.2462600D-01,0.2553350D-01,
166 + 0.2631400D-01,0.2699250D-01,0.2757350D-01,0.2806300D-01,
167 + 0.2848350D-01,0.2884800D-01,0.2915300D-01,0.2940850D-01,
168 + 0.2962500D-01,0.2981050D-01,0.2996600D-01,0.3009050D-01,
169 + 0.3018780D-01,0.3026480D-01,0.3032550D-01,0.3037250D-01,
170 + 0.3040925D-01,0.3043835D-01,0.3046145D-01,0.3048025D-01,
171 + 0.3585000D-03,0.6725000D-03,0.9675000D-03,0.1256000D-02,
172 + 0.1545500D-02,0.1843000D-02,0.2179000D-02,0.2595500D-02,
173 + 0.3130000D-02,0.3843500D-02,0.4813500D-02,0.6003500D-02,
174 + 0.7288500D-02,0.8623499D-02,0.1003850D-01,0.1157850D-01,
175 + 0.1331850D-01,0.1533350D-01,0.1762850D-01,0.2014350D-01,
176 + 0.2274850D-01,0.2532850D-01,0.2779350D-01,0.3005850D-01,
177 + 0.3206350D-01,0.3380350D-01,0.3531851D-01,0.3661850D-01,
178 + 0.3773851D-01,0.3871351D-01,0.3954751D-01,0.4025051D-01,
179 + 0.4085401D-01,0.4137701D-01,0.4181501D-01,0.4218151D-01,
180 + 0.4249151D-01,0.4275751D-01,0.4298101D-01,0.4316001D-01,
181 + 0.4330001D-01,0.4341061D-01,0.4349771D-01,0.4356516D-01,
182 + 0.4361791D-01,0.4365961D-01,0.4369271D-01,0.4371971D-01,
183 + 0.1685000D-03,0.3160000D-03,0.4545000D-03,0.5900000D-03,
184 + 0.7260000D-03,0.8655000D-03,0.1023000D-02,0.1218500D-02,
185 + 0.1469500D-02,0.1804500D-02,0.2259000D-02,0.2817500D-02,
186 + 0.3422000D-02,0.4049500D-02,0.4712999D-02,0.5434999D-02,
187 + 0.6252999D-02,0.7203000D-02,0.8283000D-02,0.9462999D-02,
188 + 0.1068800D-01,0.1190300D-01,0.1306300D-01,0.1412800D-01,
189 + 0.1507000D-01,0.1588850D-01,0.1660150D-01,0.1721300D-01,
190 + 0.1773900D-01,0.1819650D-01,0.1858850D-01,0.1891850D-01,
191 + 0.1920150D-01,0.1944700D-01,0.1965250D-01,0.1982450D-01,
192 + 0.1997050D-01,0.2009550D-01,0.2020010D-01,0.2028410D-01,
193 + 0.2034985D-01,0.2040175D-01,0.2044265D-01,0.2047435D-01,
194 + 0.2049915D-01,0.2051875D-01,0.2053430D-01,0.2054695D-01 /
195
196 DATA AE_ABI /
197 + 0.3645000E-01,0.5211000E-01,0.5913000E-01,0.6203000E-01,
198 + 0.6304000E-01,0.6340450E-01,0.6353275E-01,0.6358405E-01,
199 + 0.6361780E-01,0.6364885E-01,0.6367990E-01,0.6371365E-01,
200 + 0.6375955E-01,0.6383380E-01,0.6392965E-01,0.6403360E-01,
201 + 0.6414810E-01,0.6426660E-01,0.6438010E-01,0.6448960E-01,
202 + 0.6459510E-01,0.6468170E-01,0.6474110E-01,0.6478295E-01,
203 + 0.6481670E-01,0.6484775E-01,0.6487610E-01,0.6490240E-01,
204 + 0.3510000E-01,0.5018000E-01,0.5694000E-01,0.5973500E-01,
205 + 0.6071000E-01,0.6106100E-01,0.6118450E-01,0.6123390E-01,
206 + 0.6126640E-01,0.6129630E-01,0.6132620E-01,0.6135870E-01,
207 + 0.6140290E-01,0.6147440E-01,0.6156670E-01,0.6166680E-01,
208 + 0.6177730E-01,0.6189180E-01,0.6200130E-01,0.6210680E-01,
209 + 0.6220820E-01,0.6229140E-01,0.6234860E-01,0.6238890E-01,
210 + 0.6242140E-01,0.6245130E-01,0.6247860E-01,0.6250395E-01,
211 + 0.3375000E-01,0.4825000E-01,0.5475000E-01,0.5744000E-01,
212 + 0.5838000E-01,0.5871750E-01,0.5883625E-01,0.5888375E-01,
213 + 0.5891500E-01,0.5894375E-01,0.5897250E-01,0.5900375E-01,
214 + 0.5904625E-01,0.5911500E-01,0.5920375E-01,0.5930000E-01,
215 + 0.5940650E-01,0.5951700E-01,0.5962200E-01,0.5972300E-01,
216 + 0.5982050E-01,0.5990050E-01,0.5995550E-01,0.5999425E-01,
217 + 0.6002550E-01,0.6005425E-01,0.6008050E-01,0.6010490E-01,
218 + 0.3240000E-01,0.4632000E-01,0.5256000E-01,0.5514000E-01,
219 + 0.5604000E-01,0.5636400E-01,0.5647800E-01,0.5652360E-01,
220 + 0.5655360E-01,0.5658120E-01,0.5660880E-01,0.5663880E-01,
221 + 0.5667960E-01,0.5674560E-01,0.5683080E-01,0.5692320E-01,
222 + 0.5702520E-01,0.5713070E-01,0.5723140E-01,0.5732860E-01,
223 + 0.5742220E-01,0.5749900E-01,0.5755180E-01,0.5758900E-01,
224 + 0.5761900E-01,0.5764660E-01,0.5767180E-01,0.5769520E-01,
225 + 0.3240000E-01,0.4632000E-01,0.5256000E-01,0.5514000E-01,
226 + 0.5604000E-01,0.5636400E-01,0.5647800E-01,0.5652360E-01,
227 + 0.5655360E-01,0.5658120E-01,0.5660880E-01,0.5663880E-01,
228 + 0.5667960E-01,0.5674560E-01,0.5683080E-01,0.5692320E-01,
229 + 0.5702520E-01,0.5713070E-01,0.5723140E-01,0.5732860E-01,
230 + 0.5742220E-01,0.5749900E-01,0.5755180E-01,0.5758900E-01,
231 + 0.5761900E-01,0.5764660E-01,0.5767180E-01,0.5769520E-01,
232 + 0.3105000E-01,0.4439000E-01,0.5037000E-01,0.5284500E-01,
233 + 0.5371000E-01,0.5402050E-01,0.5412975E-01,0.5417345E-01,
234 + 0.5420220E-01,0.5422865E-01,0.5425510E-01,0.5428385E-01,
235 + 0.5432295E-01,0.5438620E-01,0.5446785E-01,0.5455640E-01,
236 + 0.5465390E-01,0.5475485E-01,0.5485145E-01,0.5494460E-01,
237 + 0.5503430E-01,0.5510790E-01,0.5515850E-01,0.5519415E-01,
238 + 0.5522290E-01,0.5524935E-01,0.5527350E-01,0.5529590E-01,
239 + 0.2700000E-01,0.3860000E-01,0.4380000E-01,0.4595000E-01,
240 + 0.4670000E-01,0.4697000E-01,0.4706500E-01,0.4710300E-01,
241 + 0.4712800E-01,0.4715100E-01,0.4717400E-01,0.4719900E-01,
242 + 0.4723300E-01,0.4728800E-01,0.4735900E-01,0.4743600E-01,
243 + 0.4752100E-01,0.4760900E-01,0.4769300E-01,0.4777400E-01,
244 + 0.4785200E-01,0.4791600E-01,0.4796000E-01,0.4799100E-01,
245 + 0.4801600E-01,0.4803900E-01,0.4806000E-01,0.4807950E-01,
246 + 0.2430000E-01,0.3474000E-01,0.3942000E-01,0.4135500E-01,
247 + 0.4203000E-01,0.4227300E-01,0.4235850E-01,0.4239270E-01,
248 + 0.4241520E-01,0.4243590E-01,0.4245660E-01,0.4247910E-01,
249 + 0.4250970E-01,0.4255920E-01,0.4262310E-01,0.4269240E-01,
250 + 0.4276890E-01,0.4284810E-01,0.4292370E-01,0.4299660E-01,
251 + 0.4306680E-01,0.4312440E-01,0.4316400E-01,0.4319190E-01,
252 + 0.4321440E-01,0.4323510E-01,0.4325400E-01,0.4327155E-01,
253 + 0.2255000E-01,0.3225500E-01,0.3659500E-01,0.3838900E-01,
254 + 0.3901500E-01,0.3924050E-01,0.3931985E-01,0.3935155E-01,
255 + 0.3937240E-01,0.3939160E-01,0.3941080E-01,0.3943165E-01,
256 + 0.3946005E-01,0.3950600E-01,0.3956530E-01,0.3962960E-01,
257 + 0.3970055E-01,0.3977400E-01,0.3984415E-01,0.3991180E-01,
258 + 0.3997695E-01,0.4003040E-01,0.4006715E-01,0.4009305E-01,
259 + 0.4011390E-01,0.4013310E-01,0.4015065E-01,0.4016695E-01,
260 + 0.2130000E-01,0.3044500E-01,0.3455500E-01,0.3625450E-01,
261 + 0.3684700E-01,0.3706050E-01,0.3713575E-01,0.3716575E-01,
262 + 0.3718550E-01,0.3720370E-01,0.3722190E-01,0.3724165E-01,
263 + 0.3726850E-01,0.3731195E-01,0.3736805E-01,0.3742890E-01,
264 + 0.3749605E-01,0.3756555E-01,0.3763190E-01,0.3769590E-01,
265 + 0.3775750E-01,0.3780805E-01,0.3784280E-01,0.3786725E-01,
266 + 0.3788700E-01,0.3790520E-01,0.3792180E-01,0.3793720E-01,
267 + 0.2025000E-01,0.2895000E-01,0.3285000E-01,0.3446250E-01,
268 + 0.3502500E-01,0.3522750E-01,0.3529875E-01,0.3532725E-01,
269 + 0.3534600E-01,0.3536325E-01,0.3538050E-01,0.3539925E-01,
270 + 0.3542475E-01,0.3546600E-01,0.3551925E-01,0.3557700E-01,
271 + 0.3564075E-01,0.3570675E-01,0.3576975E-01,0.3583050E-01,
272 + 0.3588900E-01,0.3593700E-01,0.3597000E-01,0.3599325E-01,
273 + 0.3601200E-01,0.3602925E-01,0.3604500E-01,0.3605960E-01,
274 + 0.1920000E-01,0.2745500E-01,0.3114500E-01,0.3267050E-01,
275 + 0.3320300E-01,0.3339470E-01,0.3346215E-01,0.3348915E-01,
276 + 0.3350690E-01,0.3352320E-01,0.3353950E-01,0.3355725E-01,
277 + 0.3358140E-01,0.3362045E-01,0.3367085E-01,0.3372550E-01,
278 + 0.3378585E-01,0.3384835E-01,0.3390800E-01,0.3396550E-01,
279 + 0.3402090E-01,0.3406635E-01,0.3409760E-01,0.3411965E-01,
280 + 0.3413740E-01,0.3415370E-01,0.3416860E-01,0.3418245E-01 /
281
282 DATA BATM / 1222.6562D0,1144.9069D0,1305.5948D0,540.1778D0,0.D0 /
283 DATA CATM / 994186.38D0,878153.55D0,636143.04D0,772170.16D0,1.D-9/
284
285 DATA LAHG / 2200.0D2,4.0D5,1.0D6,4.0D6,1.0D7 /
286 DATA Lm /3.70D5,3.85D5,4.0D5,4.17D5,4.17D5,4.35D5,5.0D5,5.56D5,
287 + 5.99D5,6.33D5, 6.67D5,7.04D5 /
288 DATA LONG /
289 + 280.0,
290 + 300.0,
291 + 320.0,
292 + 340.0,
293 + 360.0,
294 + 380.0,
295 + 400.0,
296 + 450.0,
297 + 500.0,
298 + 550.0,
299 + 600.0,
300 + 650.0 /
301 DATA PI / 3.141592654D0 /
302
303c-- Take same Earth radius as in CORSIKA
304 parameter (rt = 6371315.D2)
305
306c-- Scale-height for Exponential density profile
307 parameter (hscale = 7.4d5)
308
309***********************************************************************
310*
311* SCATTERING PARAMETERS FOR RAYLEIGH:
312*
313* MEAN FREE PATH FOR SCATTERING RAYLEIGH (g/cm^2)
314 PARAMETER (XR=2.970D3)
315***********************************************************************
316
317c-- Observation level at La Palma
318 parameter (obslev = 2200.d2)
319
320 T_Ray = 1.0
321 T_Mie = 1.0
322 T_Oz = 1.0
323
324c-- Height calculated using an obslev = H(LaPalma) <> 0
325c H = -RT + SQRT(RT**2 + (height/COS(theta))**2 +
326c + (2.0D0*RT*height))
327 h = -rt + sqrt((rt+obslev)**2 +
328 + ((height-obslev)/cos(theta))**2 +
329 + (2.0d0*(rt+obslev)*(height-obslev)))
330
331 ROW = AINT(((H+1.)/1.0E5))
332
333***********************************************************************
334*
335* LARGE ZENITH ANGLE FACTOR (AIR MASS):
336
337c fs : air mass factor = path_lenght(za) / path_lenght(vertical)
338c at point of emission of cherenkov photon
339c => pure geometric correction on
340c absorption lenght measured at vertical height [km^-1]
341c--
342c ma = (1.D0/H)*(SQRT((RT*COS(theta))**2+(2*RT*H)+H**2)-
343c + (RT*COS(theta)))
344c--
345c mb = (1.d0/(h-obslev)) *
346c & ( -(rt+obslev)*cos(theta)
347c & +sqrt( (rt+h)**2 - ((rt+obslev)*sin(theta))**2) )
348c--
349
350c-- Air mass "m" calcualted using a one-exponential density
351c-- profile for the atmosphere, rho = rho_0 exp(-h/Hs)
352c-- with Hs = hscale = 7.4 km
353
354 Rcos2 = rt * cos(theta)**2
355 Rsin2 = rt * sin(theta)**2
356
357 x1 = sqrt((2. * obslev + Rcos2) / (2. * hscale))
358 x2 = sqrt((2. * h + Rcos2) / (2. * hscale))
359 x3 = sqrt((2. * obslev + rt) / (2. * hscale))
360 x4 = sqrt((2. * h + rt) / (2. * hscale))
361
362 e1 = derfc(x1)
363 e2 = derfc(x2)
364 e3 = derfc(x3)
365 e4 = derfc(x4)
366
367 m = exp(-Rsin2 / (2. * hscale)) * ((e1 - e2) / (e3 - e4))
368
369**********************************************************************
370*
371* RAYLEIGH SCATTERING
372
373 RHOTOT = 0.0
374 do 11 i=1,5
375 RHOF(i) = 0
376 11 continue
377
378 DO 91 I=2,5
379 IF ( H .LT. LAHG(I) ) THEN
380 RHOTOT = RHOTOT +
381 + BATM(I-1)*(EXP(-LAHG(I-1)/CATM(I-1) ) -
382 + EXP(-H/CATM(I-1)))
383 GOTO 92
384 ELSE
385 RHOTOT = RHOTOT +
386 + BATM(I-1)*(EXP(-LAHG(I-1)/CATM(I-1) ) -
387 + EXP(-LAHG(I)/CATM(I-1)))
388 ENDIF
389 91 CONTINUE
390
391
392 92 RHO_FI = m*RHOTOT
393
394 T_Ray = EXP(-(RHO_FI/XR)*(400.D0/wavelength)**4)
395
396
397***************************************************************************
398* *
399* MIE ABSORPTION: WE USE FOR HEIGHTS LOWER THAN 10 Km AN EXACT FORMULA*
400* FOR THE AEROSOL DENSITY, AND WE USE A TABLE FOR HIGHTS HIGHER THAN *
401* 10 Km. WE TAKE THE TABLE FROM THE ALTERMAN ARTICLE. *
402***************************************************************************
403
404
405***************************************************************************
406* *
407* MIE SCATTERING PARAMETERS *
408 Hm = 1.2D5
409***************************************************************************
410
411 IF (ROW.GT.27) THEN
412 ROW=28
413 ENDIF
414
415
416 CON_MI = 2
417
418 1001 IF (ABS(LONG(CON_MI)-wavelength).LT.0.01) THEN
419 TOT_AE =AE_ABI(ROW,CON_MI)
420 ELSEIF (LONG(CON_MI).GT.wavelength) THEN
421 A = (AE_ABI(ROW,CON_MI)-AE_ABI(ROW,CON_MI-1))/
422 + (LONG(CON_MI) - LONG(CON_MI-1))
423 B = AE_ABI(ROW,CON_MI) - (A*LONG(CON_MI))
424 TOT_AE = A*wavelength + B
425 ELSE
426 CON_MI = CON_MI +1
427 GOTO 1001
428 ENDIF
429
430 T_Mie = EXP(-(m*TOT_AE))
431
432
433
434***********************************************************************
435* *
436* OZONE ABSORPTION: We used the Elterman table. *
437* *
438***********************************************************************
439 IF (ROW.GT.47) THEN
440 ROW = 47
441 ENDIF
442
443 CON_OZ = 2
444 2001 IF (LONG(CON_OZ).GT.wavelength) THEN
445 A = (OZ_ABI(ROW,CON_OZ)-OZ_ABI(ROW,CON_OZ-1))/
446 + (LONG(CON_OZ) - LONG(CON_OZ-1))
447 B = OZ_ABI(ROW,CON_OZ) - (A*LONG(CON_OZ))
448 TOT_OZ = (A* wavelength)+B
449 ELSEIF (ABS(LONG(CON_OZ)-wavelength).LT.0.01) THEN
450 TOT_OZ = OZ_ABI(ROW,CON_OZ)
451 ELSE
452 CON_OZ = CON_OZ + 1
453 GOTO 2001
454 ENDIF
455
456 T_Oz = EXP(-(m*TOT_OZ))
457
458************************************************************************
459* *
460* TOTAL TRANSMISSION OF THE ATMOSPHERE: (RAYLEIGH + MIE + OZONE) *
461************************************************************************
462
463 tr_atmos = T_Ray*T_Mie*T_Oz
464
465 RETURN
466
467 END
468
469c @endcode
470
471c EOF
Note: See TracBrowser for help on using the repository browser.