-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathshielding.f90
843 lines (821 loc) · 77.8 KB
/
shielding.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
!******************************************************************************
! MODULE: shielding
!******************************************************************************
!
! DESCRIPTION:
!> @brief Module that contain everything about shielding of H2 and CO molecules
!! H2 is only shielded by himself. CO is shielded both by H2 and himself.
!! Self shielding is computed in the ode_solver module in the subroutine
!! set_dependant_rates
!!
!
!******************************************************************************
module shielding
use numerical_types
implicit none
integer, parameter :: NL1=105
integer, parameter :: NL2=52
integer, parameter :: NL3=43
! H2 Shielding factors
! NH2 (cm-2)
real(double_precision), parameter, dimension(nl1) :: N1H2 = (/0.000d+00, 3.690d+11, &
3.715d+12, 3.948d+13, 1.233d+14, 2.536d+14, 4.342d+14, 6.653d+14, 6.689d+14, &
9.075d+14, 1.234d+15, 1.631d+15, 2.105d+15, 2.363d+15, 2.899d+15, 3.207d+15, 3.848d+15, 4.636d+15, &
5.547d+15, 6.604d+15, 7.855d+15, 9.368d+15, 1.122d+16, 1.352d+16, 1.643d+16, 2.017d+16, 2.515d+16, &
3.190d+16, 4.128d+16, 5.439d+16, 7.315d+16, 1.009d+17, 1.432d+17, 2.092d+17, 3.123d+17, 4.738d+17, &
5.388d+17, 8.935d+17, 1.381d+18, 2.164d+18, 3.330d+18, 5.024d+18, 7.404d+18, 9.029d+18, 1.316d+19, &
1.813d+19, 2.453d+19, 3.248d+19, 4.216d+19, 5.370d+19, 6.722d+19, 8.277d+19, 9.894d+19, 1.186d+20, &
1.404d+20, 1.644d+20, 1.908d+20, 2.197d+20, 2.510d+20, 2.849d+20, 3.214d+20, 3.604d+20, 4.019d+20, &
4.456d+20, 4.915d+20, 5.393d+20, 5.886d+20, 6.392d+20, 6.909d+20, 8.505d+20, 7.433d+20, 7.965d+20, &
9.056d+20, 9.627d+20, 1.011d+21, 1.068d+21, 1.125d+21, 1.185d+21, 1.250d+21, 1.327d+21, 1.428d+21, &
1.578d+21, 1.851d+21, 2.128d+21, 2.298d+21, 2.389d+21, 2.459d+21, 2.519d+21, 2.571d+21, 2.618d+21, &
2.707d+21, 2.790d+21, 2.887d+21, 3.001d+21, 3.139d+21, 3.303d+21, 3.497d+21, 3.722d+21, 3.983d+21, &
4.283d+21, 4.644d+21, 5.127d+21, 5.945d+21, 8.205d+21, 1.015d+22/)
! teta H2
real(double_precision), parameter, dimension(nl1) :: T1H2 = (/1.000d+00, 9.983d-01, 9.853d-01, &
8.761d-01, 7.199d-01, 5.728d-01, 4.455d-01, 3.431d-01, 3.418d-01, &
2.732d-01, 2.110d-01, 1.619d-01, 1.236d-01, 1.084d-01, 8.447d-02, 7.410d-02, 5.774d-02, 4.416d-02, &
3.390d-02, 2.625d-02, 2.048d-02, 1.606d-02, 1.264d-02, 9.987d-03, 7.937d-03, 6.343d-03, 5.088d-03, &
4.089d-03, 3.283d-03, 2.640d-03, 2.130d-03, 1.725d-03, 1.397d-03, 1.129d-03, 9.097d-04, 7.340d-04, &
6.883d-04, 5.377d-04, 4.352d-04, 3.475d-04, 2.771d-04, 2.205d-04, 1.753d-04, 1.549d-04, 1.210d-04, &
9.666d-05, 7.705d-05, 6.148d-05, 4.904d-05, 3.909d-05, 3.112d-05, 2.473d-05, 1.997d-05, 1.578d-05, &
1.244d-05, 9.769d-06, 7.634d-06, 5.932d-06, 4.581d-06, 3.515d-06, 2.679d-06, 2.029d-06, 1.527d-06, &
1.144d-06, 8.523d-07, 6.332d-07, 4.693d-07, 3.475d-07, 2.574d-07, 1.047d-07, 1.907d-07, 1.413d-07, &
7.739d-08, 5.677d-08, 4.386d-08, 3.227d-08, 2.385d-08, 1.750d-08, 1.248d-08, 8.389d-09, 5.026d-09, &
2.382d-09, 6.259d-10, 1.653d-10, 7.399d-11, 4.824d-11, 3.474d-11, 2.633d-11, 2.069d-11, 1.663d-11, &
1.099d-11, 7.506d-12, 4.825d-12, 2.864d-12, 1.534d-12, 7.324d-13, 3.087d-13, 1.135d-13, 3.591d-14, &
9.689d-15, 2.045d-15, 2.618d-16, 8.918d-18, 3.041d-21, 1.739d-23/)
! CO shidlding factors
! NCO (cm-2)
real(double_precision), parameter, dimension(nl2) :: N2CO = (/ 0.000d+00, 1.000d+12, &
1.650d+12, 2.995d+12, 5.979d+12, 1.313d+13, 3.172d+13, 8.429d+13, 2.464d+14, &
7.923d+14, 1.670d+15, 2.595d+15, 4.435d+15, 6.008d+15, 8.952d+15, 1.334d+16, 1.661d+16, 2.274d+16, &
3.115d+16, 4.266d+16, 5.843d+16, 8.002d+16, 1.096d+17, 1.501d+17, 2.055d+17, 2.815d+17, 4.241d+17, &
6.389d+17, 9.625d+17, 1.450d+18, 2.184d+18, 3.291d+18, 4.124d+18, 5.685d+18, 7.838d+18, 1.080d+19, &
1.285d+19, 1.681d+19, 2.199d+19, 2.538d+19, 3.222d+19, 4.091d+19, 5.193d+19, 5.893d+19, 7.356d+19, &
8.269d+19, 9.246d+19, 1.031d+20, 1.148d+20, 1.277d+20, 1.419d+20, 1.578d+20/)
! Teta CO
real(double_precision), parameter, dimension(nl2) :: T2CO = (/ 1.000d+00, 9.990d-01, &
9.981d-01, 9.961d-01, 9.912d-01, 9.815d-01, 9.601d-01, 9.113d-01, 8.094d-01, &
6.284d-01, 4.808d-01, 3.889d-01, 2.827d-01, 2.293d-01, 1.695d-01, 1.224d-01, 1.017d-01, 7.764d-02, &
5.931d-02, 4.546d-02, 3.506d-02, 2.728d-02, 2.143d-02, 1.700d-02, 1.360d-02, 1.094d-02, 8.273d-03, &
6.283d-03, 4.773d-03, 3.611d-03, 2.704d-03, 1.986d-03, 1.657d-03, 1.258d-03, 9.332d-04, 6.745d-04, &
5.596d-04, 4.123d-04, 2.982d-04, 2.490d-04, 1.827d-04, 1.324d-04, 9.473d-05, 7.891d-05, 5.668d-05, &
4.732d-05, 3.967d-05, 3.327d-05, 2.788d-05, 2.331d-05, 1.944d-05, 1.619d-05/)
!NH2 (cm-2)
real(double_precision), parameter, dimension(nl3) :: N2H2 = (/0.000d+00, 2.666d+13, &
3.801d+14, 6.634d+15, 8.829d+16, 9.268d+17, 1.007d+18, 2.021d+18, 3.036d+18, &
4.051d+18, 5.066d+18, 6.082d+18, 7.097d+18, 8.112d+18, 9.341d+18, 1.014d+19, 2.030d+19, 3.045d+19, &
4.061d+19, 5.076d+19, 6.092d+19, 7.107d+19, 8.123d+19, 9.353d+19, 1.015d+20, 2.031d+20, 3.047d+20, &
4.062d+20, 5.078d+20, 6.094d+20, 7.109d+20, 8.125d+20, 9.355d+20, 1.016d+21, 2.031d+21, 3.047d+21, &
4.063d+21, 5.078d+21, 6.094d+21, 7.110d+21, 8.125d+21, 9.355d+21, 1.016d+22/)
! Teta H2
real(double_precision), parameter, dimension(nl3) :: T2H2 = (/1.000d+00, 9.999d-01, &
9.893d-01, 9.678d-01, 9.465d-01, 9.137d-01, 9.121d-01, 8.966d-01, 8.862d-01, &
8.781d-01, 8.716d-01, 8.660d-01, 8.612d-01, 8.569d-01, 8.524d-01, 8.497d-01, 8.262d-01, 8.118d-01, &
8.010d-01, 7.921d-01, 7.841d-01, 7.769d-01, 7.702d-01, 7.626d-01, 7.579d-01, 7.094d-01, 6.712d-01, &
6.378d-01, 6.074d-01, 5.791d-01, 5.524d-01, 5.271d-01, 4.977d-01, 4.793d-01, 2.837d-01, 1.526d-01, &
7.774d-02, 3.952d-02, 2.093d-02, 1.199d-02, 7.666d-03, 5.333d-03, 4.666d-03/)
! AV
real(double_precision), parameter, dimension(nl3) :: AV2 = (/0.000d+00, 1.000d-07, &
1.000d-06, 1.000d-05, 1.000d-04, 1.000d-03, 1.086d-03, 2.171d-03, 3.257d-03, &
4.343d-03, 5.429d-03, 6.514d-03, 7.600d-03, 8.686d-03, 1.000d-02, 1.086d-02, 2.171d-02, 3.257d-02, &
4.343d-02, 5.429d-02, 6.514d-02, 7.600d-02, 8.686d-02, 1.000d-01, 1.086d-01, 2.171d-01, 3.257d-01, &
4.343d-01, 5.429d-01, 6.514d-01, 7.600d-01, 8.686d-01, 1.000d+00, 1.086d+00, 2.171d+00, 3.257d+00, &
4.343d+00, 5.429d+00, 6.514d+00, 7.600d+00, 8.686d+00, 1.000d+01, 1.086d+01/)
! teta AV
real(double_precision), parameter, dimension(nl3) :: T2AV = (/1.000d+00, 1.000d+00, &
1.000d+00, 1.000d+00, 9.991d-01, 9.887d-01, 9.877d-01, 9.755d-01, 9.638d-01, &
9.524d-01, 9.413d-01, 9.306d-01, 9.202d-01, 9.101d-01, 8.983d-01, 8.908d-01, 8.081d-01, 7.422d-01, &
6.875d-01, 6.404d-01, 5.990d-01, 5.620d-01, 5.284d-01, 4.916d-01, 4.696d-01, 2.757d-01, 1.705d-01, &
1.084d-01, 7.013d-02, 4.593d-02, 3.037d-02, 2.023d-02, 1.247d-02, 9.130d-03, 2.133d-04, 6.269d-06, &
2.108d-07, 7.633d-09, 2.873d-10, 1.102d-11, 4.273d-13, 8.412d-15, 6.510d-16/)
! -------- Shielding of CO from Visser et al. (2009) ---------
real(double_precision), parameter, dimension(47) :: NCO_CO = (/ &
1.000E+00, 1.000E+10, 1.585E+10, 2.512E+10, 3.981E+10, 6.310E+10, 1.000E+11, 1.585E+11, 2.512E+11, &
3.981E+11, 6.310E+11, 1.000E+12, 1.585E+12, 2.512E+12, 3.981E+12, 6.310E+12, 1.000E+13, 1.585E+13, &
2.512E+13, 3.981E+13, 6.310E+13, 1.000E+14, 1.585E+14, 2.512E+14, 3.981E+14, 6.310E+14, 1.000E+15, &
1.585E+15, 2.512E+15, 3.981E+15, 6.310E+15, 1.000E+16, 1.585E+16, 2.512E+16, 3.981E+16, 6.310E+16, &
1.000E+17, 1.585E+17, 2.512E+17, 3.981E+17, 6.310E+17, 1.000E+18, 1.585E+18, 2.512E+18, 3.981E+18, &
6.310E+18, 1.000E+19/)
real(double_precision), parameter, dimension(42) :: NH2_CO = (/ &
1.000E+00, 1.000E+15, 1.585E+15, 2.512E+15, 3.981E+15, 6.310E+15, 1.000E+16, 1.585E+16, 2.512E+16, &
3.981E+16, 6.310E+16, 1.000E+17, 1.585E+17, 2.512E+17, 3.981E+17, 6.310E+17, 1.000E+18, 1.585E+18, &
2.512E+18, 3.981E+18, 6.310E+18, 1.000E+19, 1.585E+19, 2.512E+19, 3.981E+19, 6.310E+19, 1.000E+20, &
1.585E+20, 2.512E+20, 3.981E+20, 6.310E+20, 1.000E+21, 1.585E+21, 2.512E+21, 3.981E+21, 6.310E+21, &
1.000E+22, 1.585E+22, 2.512E+22, 3.981E+22, 6.310E+22, 1.000E+23/)
real(double_precision), parameter, dimension(42,47) :: SHIELDING_FUNCTION_CO = reshape((/ &
1.000E+00, 8.985E-01, 8.966E-01, 8.949E-01, 8.932E-01, 8.915E-01, 8.899E-01, 8.855E-01, 8.834E-01, 8.814E-01, &
8.792E-01, 8.766E-01, 8.735E-01, 8.697E-01, 8.652E-01, 8.600E-01, 8.543E-01, 8.482E-01, 8.416E-01, 8.345E-01, &
8.265E-01, 8.176E-01, 8.073E-01, 7.949E-01, 7.784E-01, 7.553E-01, 7.223E-01, 6.758E-01, 6.127E-01, 5.310E-01, &
4.328E-01, 3.260E-01, 2.241E-01, 1.394E-01, 7.604E-02, 3.382E-02, 1.108E-02, 2.364E-03, 2.983E-04, 2.297E-05, &
1.546E-06, 3.938E-07, &
9.997E-01, 8.983E-01, 8.963E-01, 8.946E-01, 8.929E-01, 8.912E-01, 8.896E-01, 8.852E-01, 8.831E-01, 8.811E-01, &
8.789E-01, 8.764E-01, 8.732E-01, 8.694E-01, 8.649E-01, 8.597E-01, 8.540E-01, 8.479E-01, 8.414E-01, 8.342E-01, &
8.262E-01, 8.173E-01, 8.070E-01, 7.946E-01, 7.782E-01, 7.550E-01, 7.220E-01, 6.756E-01, 6.125E-01, 5.309E-01, &
4.327E-01, 3.260E-01, 2.241E-01, 1.394E-01, 7.604E-02, 3.382E-02, 1.108E-02, 2.364E-03, 2.982E-04, 2.297E-05, &
1.546E-06, 3.938E-07, &
9.995E-01, 8.981E-01, 8.962E-01, 8.944E-01, 8.927E-01, 8.911E-01, 8.895E-01, 8.850E-01, 8.829E-01, 8.809E-01, &
8.788E-01, 8.762E-01, 8.730E-01, 8.692E-01, 8.647E-01, 8.595E-01, 8.539E-01, 8.477E-01, 8.412E-01, 8.340E-01, &
8.260E-01, 8.171E-01, 8.069E-01, 7.944E-01, 7.780E-01, 7.549E-01, 7.219E-01, 6.754E-01, 6.124E-01, 5.308E-01, &
4.326E-01, 3.259E-01, 2.241E-01, 1.394E-01, 7.603E-02, 3.382E-02, 1.108E-02, 2.364E-03, 2.982E-04, 2.297E-05, &
1.546E-06, 3.938E-07, &
9.992E-01, 8.978E-01, 8.959E-01, 8.941E-01, 8.924E-01, 8.908E-01, 8.892E-01, 8.848E-01, 8.826E-01, 8.807E-01, &
8.785E-01, 8.759E-01, 8.728E-01, 8.689E-01, 8.644E-01, 8.593E-01, 8.536E-01, 8.475E-01, 8.409E-01, 8.337E-01, &
8.258E-01, 8.169E-01, 8.066E-01, 7.942E-01, 7.778E-01, 7.546E-01, 7.216E-01, 6.752E-01, 6.122E-01, 5.306E-01, &
4.325E-01, 3.258E-01, 2.241E-01, 1.393E-01, 7.603E-02, 3.381E-02, 1.108E-02, 2.364E-03, 2.982E-04, 2.296E-05, &
1.546E-06, 3.938E-07, &
9.988E-01, 8.974E-01, 8.955E-01, 8.937E-01, 8.920E-01, 8.904E-01, 8.888E-01, 8.843E-01, 8.822E-01, 8.802E-01, &
8.781E-01, 8.755E-01, 8.723E-01, 8.685E-01, 8.640E-01, 8.588E-01, 8.532E-01, 8.470E-01, 8.405E-01, 8.333E-01, &
8.254E-01, 8.164E-01, 8.062E-01, 7.938E-01, 7.774E-01, 7.543E-01, 7.213E-01, 6.749E-01, 6.119E-01, 5.304E-01, &
4.323E-01, 3.258E-01, 2.240E-01, 1.393E-01, 7.603E-02, 3.381E-02, 1.108E-02, 2.364E-03, 2.982E-04, 2.296E-05, &
1.546E-06, 3.938E-07, &
9.981E-01, 8.967E-01, 8.948E-01, 8.930E-01, 8.913E-01, 8.897E-01, 8.881E-01, 8.837E-01, 8.815E-01, 8.796E-01, &
8.774E-01, 8.748E-01, 8.716E-01, 8.678E-01, 8.633E-01, 8.582E-01, 8.525E-01, 8.464E-01, 8.398E-01, 8.326E-01, &
8.247E-01, 8.158E-01, 8.055E-01, 7.931E-01, 7.767E-01, 7.536E-01, 7.207E-01, 6.743E-01, 6.114E-01, 5.300E-01, &
4.321E-01, 3.256E-01, 2.240E-01, 1.393E-01, 7.602E-02, 3.381E-02, 1.108E-02, 2.364E-03, 2.982E-04, 2.296E-05, &
1.546E-06, 3.938E-07, &
9.970E-01, 8.956E-01, 8.937E-01, 8.920E-01, 8.903E-01, 8.886E-01, 8.870E-01, 8.826E-01, 8.805E-01, 8.785E-01, &
8.763E-01, 8.737E-01, 8.706E-01, 8.668E-01, 8.623E-01, 8.571E-01, 8.514E-01, 8.453E-01, 8.387E-01, 8.316E-01, &
8.236E-01, 8.147E-01, 8.045E-01, 7.921E-01, 7.757E-01, 7.527E-01, 7.198E-01, 6.735E-01, 6.107E-01, 5.295E-01, &
4.317E-01, 3.253E-01, 2.238E-01, 1.393E-01, 7.601E-02, 3.381E-02, 1.107E-02, 2.363E-03, 2.982E-04, 2.296E-05, &
1.545E-06, 3.938E-07, &
9.953E-01, 8.939E-01, 8.920E-01, 8.903E-01, 8.886E-01, 8.869E-01, 8.853E-01, 8.809E-01, 8.788E-01, 8.768E-01, &
8.746E-01, 8.720E-01, 8.689E-01, 8.651E-01, 8.606E-01, 8.554E-01, 8.497E-01, 8.436E-01, 8.371E-01, 8.299E-01, &
8.220E-01, 8.130E-01, 8.028E-01, 7.904E-01, 7.741E-01, 7.511E-01, 7.183E-01, 6.722E-01, 6.095E-01, 5.285E-01, &
4.310E-01, 3.250E-01, 2.237E-01, 1.392E-01, 7.599E-02, 3.380E-02, 1.107E-02, 2.363E-03, 2.982E-04, 2.296E-05, &
1.545E-06, 3.938E-07, &
9.926E-01, 8.913E-01, 8.894E-01, 8.876E-01, 8.859E-01, 8.843E-01, 8.827E-01, 8.782E-01, 8.761E-01, 8.741E-01, &
8.720E-01, 8.694E-01, 8.662E-01, 8.624E-01, 8.579E-01, 8.528E-01, 8.471E-01, 8.410E-01, 8.344E-01, 8.273E-01, &
8.193E-01, 8.104E-01, 8.002E-01, 7.879E-01, 7.716E-01, 7.487E-01, 7.160E-01, 6.701E-01, 6.077E-01, 5.271E-01, &
4.300E-01, 3.243E-01, 2.234E-01, 1.391E-01, 7.596E-02, 3.379E-02, 1.107E-02, 2.363E-03, 2.981E-04, 2.296E-05, &
1.545E-06, 3.938E-07, &
9.883E-01, 8.871E-01, 8.852E-01, 8.834E-01, 8.818E-01, 8.801E-01, 8.785E-01, 8.741E-01, 8.720E-01, 8.700E-01, &
8.678E-01, 8.652E-01, 8.621E-01, 8.583E-01, 8.538E-01, 8.486E-01, 8.429E-01, 8.368E-01, 8.303E-01, 8.231E-01, &
8.152E-01, 8.063E-01, 7.962E-01, 7.839E-01, 7.677E-01, 7.450E-01, 7.125E-01, 6.668E-01, 6.049E-01, 5.248E-01, &
4.283E-01, 3.234E-01, 2.229E-01, 1.390E-01, 7.591E-02, 3.378E-02, 1.107E-02, 2.362E-03, 2.980E-04, 2.296E-05, &
1.545E-06, 3.938E-07, &
9.817E-01, 8.807E-01, 8.788E-01, 8.770E-01, 8.753E-01, 8.737E-01, 8.721E-01, 8.676E-01, 8.655E-01, 8.635E-01, &
8.614E-01, 8.588E-01, 8.556E-01, 8.518E-01, 8.473E-01, 8.422E-01, 8.365E-01, 8.304E-01, 8.238E-01, 8.167E-01, &
8.088E-01, 8.000E-01, 7.899E-01, 7.777E-01, 7.617E-01, 7.391E-01, 7.069E-01, 6.618E-01, 6.005E-01, 5.212E-01, &
4.258E-01, 3.219E-01, 2.223E-01, 1.387E-01, 7.584E-02, 3.376E-02, 1.106E-02, 2.360E-03, 2.979E-04, 2.295E-05, &
1.545E-06, 3.938E-07, &
9.716E-01, 8.707E-01, 8.688E-01, 8.671E-01, 8.654E-01, 8.637E-01, 8.621E-01, 8.577E-01, 8.556E-01, 8.536E-01, &
8.515E-01, 8.489E-01, 8.457E-01, 8.419E-01, 8.374E-01, 8.323E-01, 8.266E-01, 8.205E-01, 8.140E-01, 8.069E-01, &
7.990E-01, 7.903E-01, 7.802E-01, 7.682E-01, 7.524E-01, 7.301E-01, 6.984E-01, 6.540E-01, 5.937E-01, 5.158E-01, &
4.220E-01, 3.196E-01, 2.212E-01, 1.384E-01, 7.573E-02, 3.372E-02, 1.105E-02, 2.358E-03, 2.977E-04, 2.295E-05, &
1.545E-06, 3.938E-07, &
9.563E-01, 8.558E-01, 8.539E-01, 8.521E-01, 8.504E-01, 8.488E-01, 8.472E-01, 8.428E-01, 8.406E-01, 8.387E-01, &
8.365E-01, 8.339E-01, 8.308E-01, 8.270E-01, 8.225E-01, 8.173E-01, 8.117E-01, 8.056E-01, 7.991E-01, 7.920E-01, &
7.842E-01, 7.756E-01, 7.657E-01, 7.538E-01, 7.383E-01, 7.166E-01, 6.856E-01, 6.423E-01, 5.835E-01, 5.076E-01, &
4.161E-01, 3.161E-01, 2.196E-01, 1.378E-01, 7.555E-02, 3.366E-02, 1.103E-02, 2.355E-03, 2.974E-04, 2.294E-05, &
1.545E-06, 3.938E-07, &
9.338E-01, 8.338E-01, 8.319E-01, 8.302E-01, 8.285E-01, 8.269E-01, 8.253E-01, 8.209E-01, 8.187E-01, 8.168E-01, &
8.146E-01, 8.120E-01, 8.089E-01, 8.051E-01, 8.006E-01, 7.955E-01, 7.898E-01, 7.838E-01, 7.773E-01, 7.703E-01, &
7.626E-01, 7.540E-01, 7.443E-01, 7.328E-01, 7.177E-01, 6.967E-01, 6.668E-01, 6.250E-01, 5.685E-01, 4.955E-01, &
4.075E-01, 3.110E-01, 2.172E-01, 1.369E-01, 7.528E-02, 3.357E-02, 1.101E-02, 2.350E-03, 2.970E-04, 2.292E-05, &
1.545E-06, 3.938E-07, &
9.021E-01, 8.030E-01, 8.011E-01, 7.993E-01, 7.976E-01, 7.960E-01, 7.944E-01, 7.900E-01, 7.879E-01, 7.859E-01, &
7.838E-01, 7.812E-01, 7.781E-01, 7.743E-01, 7.698E-01, 7.647E-01, 7.591E-01, 7.531E-01, 7.466E-01, 7.397E-01, &
7.321E-01, 7.237E-01, 7.143E-01, 7.032E-01, 6.888E-01, 6.688E-01, 6.404E-01, 6.008E-01, 5.473E-01, 4.784E-01, &
3.952E-01, 3.036E-01, 2.137E-01, 1.356E-01, 7.486E-02, 3.343E-02, 1.097E-02, 2.342E-03, 2.962E-04, 2.290E-05, &
1.544E-06, 3.938E-07, &
8.599E-01, 7.621E-01, 7.602E-01, 7.585E-01, 7.568E-01, 7.551E-01, 7.536E-01, 7.492E-01, 7.471E-01, 7.451E-01, &
7.429E-01, 7.404E-01, 7.372E-01, 7.335E-01, 7.290E-01, 7.239E-01, 7.183E-01, 7.124E-01, 7.060E-01, 6.992E-01, &
6.918E-01, 6.836E-01, 6.746E-01, 6.640E-01, 6.504E-01, 6.317E-01, 6.053E-01, 5.686E-01, 5.192E-01, 4.557E-01, &
3.788E-01, 2.937E-01, 2.089E-01, 1.338E-01, 7.422E-02, 3.322E-02, 1.090E-02, 2.330E-03, 2.951E-04, 2.286E-05, &
1.544E-06, 3.938E-07, &
8.080E-01, 7.122E-01, 7.103E-01, 7.086E-01, 7.069E-01, 7.053E-01, 7.037E-01, 6.993E-01, 6.972E-01, 6.953E-01, &
6.931E-01, 6.906E-01, 6.874E-01, 6.837E-01, 6.793E-01, 6.742E-01, 6.686E-01, 6.627E-01, 6.565E-01, 6.498E-01, &
6.425E-01, 6.347E-01, 6.261E-01, 6.162E-01, 6.037E-01, 5.865E-01, 5.624E-01, 5.292E-01, 4.847E-01, 4.276E-01, &
3.584E-01, 2.810E-01, 2.025E-01, 1.312E-01, 7.327E-02, 3.289E-02, 1.081E-02, 2.310E-03, 2.933E-04, 2.280E-05, &
1.542E-06, 3.938E-07, &
7.498E-01, 6.569E-01, 6.551E-01, 6.533E-01, 6.517E-01, 6.501E-01, 6.485E-01, 6.442E-01, 6.421E-01, 6.401E-01, &
6.380E-01, 6.355E-01, 6.324E-01, 6.286E-01, 6.242E-01, 6.192E-01, 6.137E-01, 6.079E-01, 6.017E-01, 5.952E-01, &
5.881E-01, 5.806E-01, 5.725E-01, 5.633E-01, 5.519E-01, 5.364E-01, 5.149E-01, 4.854E-01, 4.461E-01, 3.959E-01, &
3.350E-01, 2.661E-01, 1.947E-01, 1.277E-01, 7.189E-02, 3.239E-02, 1.066E-02, 2.281E-03, 2.905E-04, 2.270E-05, &
1.541E-06, 3.937E-07, &
6.900E-01, 6.014E-01, 5.996E-01, 5.979E-01, 5.962E-01, 5.947E-01, 5.931E-01, 5.888E-01, 5.867E-01, 5.848E-01, &
5.827E-01, 5.802E-01, 5.771E-01, 5.734E-01, 5.690E-01, 5.640E-01, 5.586E-01, 5.529E-01, 5.468E-01, 5.405E-01, &
5.337E-01, 5.265E-01, 5.189E-01, 5.104E-01, 5.000E-01, 4.861E-01, 4.670E-01, 4.410E-01, 4.067E-01, 3.632E-01, &
3.103E-01, 2.497E-01, 1.854E-01, 1.232E-01, 6.995E-02, 3.165E-02, 1.043E-02, 2.236E-03, 2.863E-04, 2.255E-05, &
1.538E-06, 3.937E-07, &
6.323E-01, 5.494E-01, 5.476E-01, 5.460E-01, 5.444E-01, 5.428E-01, 5.413E-01, 5.371E-01, 5.350E-01, 5.331E-01, &
5.310E-01, 5.285E-01, 5.255E-01, 5.218E-01, 5.175E-01, 5.126E-01, 5.072E-01, 5.016E-01, 4.957E-01, 4.895E-01, &
4.830E-01, 4.761E-01, 4.689E-01, 4.611E-01, 4.516E-01, 4.391E-01, 4.221E-01, 3.991E-01, 3.691E-01, 3.313E-01, &
2.854E-01, 2.324E-01, 1.748E-01, 1.176E-01, 6.729E-02, 3.060E-02, 1.011E-02, 2.171E-03, 2.801E-04, 2.233E-05, &
1.534E-06, 3.937E-07, &
5.777E-01, 5.021E-01, 5.004E-01, 4.988E-01, 4.972E-01, 4.957E-01, 4.942E-01, 4.901E-01, 4.881E-01, 4.862E-01, &
4.842E-01, 4.817E-01, 4.787E-01, 4.751E-01, 4.709E-01, 4.660E-01, 4.608E-01, 4.553E-01, 4.496E-01, 4.436E-01, &
4.373E-01, 4.308E-01, 4.240E-01, 4.168E-01, 4.081E-01, 3.966E-01, 3.812E-01, 3.607E-01, 3.342E-01, 3.010E-01, &
2.609E-01, 2.143E-01, 1.629E-01, 1.106E-01, 6.380E-02, 2.915E-02, 9.646E-03, 2.078E-03, 2.713E-04, 2.200E-05, &
1.527E-06, 3.937E-07, &
5.250E-01, 4.578E-01, 4.562E-01, 4.546E-01, 4.531E-01, 4.517E-01, 4.503E-01, 4.463E-01, 4.443E-01, 4.425E-01, &
4.405E-01, 4.381E-01, 4.352E-01, 4.317E-01, 4.276E-01, 4.229E-01, 4.178E-01, 4.125E-01, 4.069E-01, 4.012E-01, &
3.952E-01, 3.891E-01, 3.828E-01, 3.760E-01, 3.680E-01, 3.575E-01, 3.434E-01, 3.249E-01, 3.012E-01, 2.718E-01, &
2.364E-01, 1.953E-01, 1.496E-01, 1.024E-01, 5.941E-02, 2.725E-02, 9.033E-03, 1.955E-03, 2.595E-04, 2.154E-05, &
1.516E-06, 3.936E-07, &
4.720E-01, 4.137E-01, 4.122E-01, 4.107E-01, 4.094E-01, 4.080E-01, 4.067E-01, 4.028E-01, 4.010E-01, 3.993E-01, &
3.974E-01, 3.951E-01, 3.923E-01, 3.889E-01, 3.849E-01, 3.804E-01, 3.755E-01, 3.704E-01, 3.651E-01, 3.597E-01, &
3.542E-01, 3.485E-01, 3.427E-01, 3.364E-01, 3.290E-01, 3.194E-01, 3.066E-01, 2.898E-01, 2.686E-01, 2.426E-01, &
2.115E-01, 1.753E-01, 1.349E-01, 9.286E-02, 5.415E-02, 2.490E-02, 8.269E-03, 1.802E-03, 2.446E-04, 2.089E-05, &
1.500E-06, 3.935E-07, &
4.177E-01, 3.681E-01, 3.667E-01, 3.655E-01, 3.642E-01, 3.630E-01, 3.618E-01, 3.582E-01, 3.565E-01, 3.549E-01, &
3.531E-01, 3.509E-01, 3.483E-01, 3.451E-01, 3.414E-01, 3.371E-01, 3.325E-01, 3.277E-01, 3.228E-01, 3.178E-01, &
3.127E-01, 3.076E-01, 3.023E-01, 2.966E-01, 2.898E-01, 2.811E-01, 2.695E-01, 2.546E-01, 2.358E-01, 2.129E-01, &
1.857E-01, 1.543E-01, 1.192E-01, 8.230E-02, 4.814E-02, 2.219E-02, 7.385E-03, 1.625E-03, 2.268E-04, 2.002E-05, &
1.477E-06, 3.933E-07, &
3.614E-01, 3.201E-01, 3.190E-01, 3.179E-01, 3.169E-01, 3.159E-01, 3.148E-01, 3.114E-01, 3.099E-01, 3.085E-01, &
3.068E-01, 3.049E-01, 3.025E-01, 2.996E-01, 2.962E-01, 2.923E-01, 2.880E-01, 2.836E-01, 2.792E-01, 2.747E-01, &
2.701E-01, 2.656E-01, 2.609E-01, 2.559E-01, 2.498E-01, 2.420E-01, 2.317E-01, 2.186E-01, 2.022E-01, 1.825E-01, &
1.592E-01, 1.325E-01, 1.025E-01, 7.094E-02, 4.160E-02, 1.922E-02, 6.427E-03, 1.434E-03, 2.066E-04, 1.886E-05, &
1.442E-06, 3.931E-07, &
3.028E-01, 2.694E-01, 2.685E-01, 2.677E-01, 2.669E-01, 2.661E-01, 2.653E-01, 2.622E-01, 2.609E-01, 2.597E-01, &
2.583E-01, 2.566E-01, 2.546E-01, 2.520E-01, 2.490E-01, 2.456E-01, 2.418E-01, 2.379E-01, 2.340E-01, 2.300E-01, &
2.261E-01, 2.222E-01, 2.182E-01, 2.138E-01, 2.085E-01, 2.017E-01, 1.929E-01, 1.816E-01, 1.678E-01, 1.513E-01, &
1.320E-01, 1.098E-01, 8.509E-02, 5.903E-02, 3.472E-02, 1.612E-02, 5.434E-03, 1.236E-03, 1.841E-04, 1.736E-05, &
1.393E-06, 3.928E-07, &
2.434E-01, 2.177E-01, 2.171E-01, 2.165E-01, 2.159E-01, 2.154E-01, 2.148E-01, 2.120E-01, 2.111E-01, 2.101E-01, &
2.090E-01, 2.076E-01, 2.060E-01, 2.039E-01, 2.014E-01, 1.985E-01, 1.953E-01, 1.920E-01, 1.886E-01, 1.853E-01, &
1.820E-01, 1.787E-01, 1.753E-01, 1.716E-01, 1.671E-01, 1.614E-01, 1.540E-01, 1.448E-01, 1.335E-01, 1.202E-01, &
1.048E-01, 8.726E-02, 6.767E-02, 4.704E-02, 2.779E-02, 1.299E-02, 4.441E-03, 1.034E-03, 1.592E-04, 1.548E-05, &
1.329E-06, 3.923E-07, &
1.871E-01, 1.683E-01, 1.679E-01, 1.676E-01, 1.672E-01, 1.669E-01, 1.665E-01, 1.642E-01, 1.635E-01, 1.629E-01, &
1.620E-01, 1.611E-01, 1.598E-01, 1.582E-01, 1.563E-01, 1.541E-01, 1.516E-01, 1.489E-01, 1.461E-01, 1.433E-01, &
1.406E-01, 1.378E-01, 1.350E-01, 1.320E-01, 1.283E-01, 1.236E-01, 1.177E-01, 1.104E-01, 1.016E-01, 9.136E-02, &
7.956E-02, 6.618E-02, 5.133E-02, 3.576E-02, 2.123E-02, 1.003E-02, 3.492E-03, 8.355E-04, 1.325E-04, 1.326E-05, &
1.251E-06, 3.918E-07, &
1.387E-01, 1.256E-01, 1.254E-01, 1.252E-01, 1.250E-01, 1.248E-01, 1.246E-01, 1.227E-01, 1.223E-01, 1.219E-01, &
1.214E-01, 1.207E-01, 1.199E-01, 1.188E-01, 1.175E-01, 1.159E-01, 1.140E-01, 1.120E-01, 1.098E-01, 1.075E-01, &
1.053E-01, 1.031E-01, 1.008E-01, 9.827E-02, 9.528E-02, 9.158E-02, 8.696E-02, 8.136E-02, 7.474E-02, 6.707E-02, &
5.830E-02, 4.841E-02, 3.752E-02, 2.616E-02, 1.561E-02, 7.466E-03, 2.653E-03, 6.521E-04, 1.059E-04, 1.093E-05, &
1.168E-06, 3.913E-07, &
1.012E-01, 9.224E-02, 9.214E-02, 9.204E-02, 9.193E-02, 9.182E-02, 9.170E-02, 9.024E-02, 9.003E-02, 8.978E-02, &
8.948E-02, 8.911E-02, 8.861E-02, 8.795E-02, 8.712E-02, 8.608E-02, 8.481E-02, 8.334E-02, 8.169E-02, 7.995E-02, &
7.817E-02, 7.637E-02, 7.452E-02, 7.248E-02, 7.009E-02, 6.719E-02, 6.364E-02, 5.941E-02, 5.447E-02, 4.880E-02, &
4.233E-02, 3.506E-02, 2.712E-02, 1.892E-02, 1.134E-02, 5.480E-03, 1.982E-03, 4.978E-04, 8.236E-05, 8.801E-06, &
1.092E-06, 3.908E-07, &
7.401E-02, 6.786E-02, 6.781E-02, 6.776E-02, 6.770E-02, 6.764E-02, 6.758E-02, 6.640E-02, 6.629E-02, 6.616E-02, &
6.601E-02, 6.581E-02, 6.554E-02, 6.517E-02, 6.469E-02, 6.406E-02, 6.326E-02, 6.227E-02, 6.110E-02, 5.981E-02, &
5.843E-02, 5.701E-02, 5.552E-02, 5.388E-02, 5.199E-02, 4.972E-02, 4.701E-02, 4.382E-02, 4.013E-02, 3.591E-02, &
3.109E-02, 2.569E-02, 1.982E-02, 1.381E-02, 8.295E-03, 4.038E-03, 1.479E-03, 3.767E-04, 6.326E-05, 7.086E-06, &
1.032E-06, 3.904E-07, &
5.467E-02, 5.034E-02, 5.031E-02, 5.028E-02, 5.025E-02, 5.022E-02, 5.018E-02, 4.917E-02, 4.911E-02, 4.905E-02, &
4.897E-02, 4.887E-02, 4.874E-02, 4.854E-02, 4.828E-02, 4.791E-02, 4.743E-02, 4.680E-02, 4.602E-02, 4.510E-02, &
4.407E-02, 4.297E-02, 4.179E-02, 4.048E-02, 3.899E-02, 3.723E-02, 3.515E-02, 3.274E-02, 2.996E-02, 2.678E-02, &
2.315E-02, 1.907E-02, 1.466E-02, 1.020E-02, 6.124E-03, 2.990E-03, 1.102E-03, 2.830E-04, 4.825E-05, 5.751E-06, &
9.867E-07, 3.901E-07, &
4.075E-02, 3.767E-02, 3.765E-02, 3.763E-02, 3.761E-02, 3.759E-02, 3.757E-02, 3.667E-02, 3.664E-02, 3.661E-02, &
3.657E-02, 3.652E-02, 3.644E-02, 3.634E-02, 3.619E-02, 3.598E-02, 3.569E-02, 3.530E-02, 3.479E-02, 3.417E-02, &
3.343E-02, 3.260E-02, 3.167E-02, 3.064E-02, 2.946E-02, 2.808E-02, 2.649E-02, 2.465E-02, 2.255E-02, 2.014E-02, &
1.737E-02, 1.426E-02, 1.092E-02, 7.568E-03, 4.536E-03, 2.213E-03, 8.166E-04, 2.105E-04, 3.638E-05, 4.678E-06, &
9.497E-07, 3.898E-07, &
3.063E-02, 2.844E-02, 2.842E-02, 2.841E-02, 2.839E-02, 2.838E-02, 2.837E-02, 2.759E-02, 2.757E-02, 2.755E-02, &
2.753E-02, 2.749E-02, 2.745E-02, 2.738E-02, 2.729E-02, 2.717E-02, 2.699E-02, 2.675E-02, 2.643E-02, 2.602E-02, &
2.550E-02, 2.488E-02, 2.417E-02, 2.335E-02, 2.242E-02, 2.134E-02, 2.010E-02, 1.869E-02, 1.708E-02, 1.523E-02, &
1.310E-02, 1.073E-02, 8.184E-03, 5.648E-03, 3.374E-03, 1.643E-03, 6.052E-04, 1.558E-04, 2.723E-05, 3.829E-06, &
9.197E-07, 3.896E-07, &
2.323E-02, 2.167E-02, 2.166E-02, 2.165E-02, 2.164E-02, 2.162E-02, 2.161E-02, 2.096E-02, 2.095E-02, 2.094E-02, &
2.092E-02, 2.090E-02, 2.087E-02, 2.083E-02, 2.077E-02, 2.069E-02, 2.058E-02, 2.043E-02, 2.023E-02, 1.995E-02, &
1.960E-02, 1.915E-02, 1.861E-02, 1.796E-02, 1.722E-02, 1.637E-02, 1.540E-02, 1.430E-02, 1.304E-02, 1.161E-02, &
9.965E-03, 8.134E-03, 6.184E-03, 4.253E-03, 2.534E-03, 1.232E-03, 4.528E-04, 1.162E-04, 2.054E-05, 3.200E-06, &
8.973E-07, 3.894E-07, &
1.775E-02, 1.664E-02, 1.663E-02, 1.662E-02, 1.661E-02, 1.661E-02, 1.660E-02, 1.608E-02, 1.607E-02, 1.606E-02, &
1.605E-02, 1.603E-02, 1.601E-02, 1.598E-02, 1.594E-02, 1.589E-02, 1.582E-02, 1.572E-02, 1.559E-02, 1.541E-02, &
1.516E-02, 1.484E-02, 1.443E-02, 1.393E-02, 1.334E-02, 1.267E-02, 1.190E-02, 1.103E-02, 1.004E-02, 8.915E-03, &
7.636E-03, 6.216E-03, 4.713E-03, 3.233E-03, 1.924E-03, 9.352E-04, 3.439E-04, 8.813E-05, 1.583E-05, 2.755E-06, &
8.813E-07, 3.893E-07, &
1.362E-02, 1.282E-02, 1.282E-02, 1.281E-02, 1.281E-02, 1.280E-02, 1.280E-02, 1.239E-02, 1.238E-02, 1.237E-02, &
1.237E-02, 1.235E-02, 1.234E-02, 1.232E-02, 1.229E-02, 1.226E-02, 1.221E-02, 1.214E-02, 1.205E-02, 1.193E-02, &
1.176E-02, 1.152E-02, 1.122E-02, 1.084E-02, 1.038E-02, 9.843E-03, 9.231E-03, 8.540E-03, 7.759E-03, 6.873E-03, &
5.872E-03, 4.768E-03, 3.606E-03, 2.470E-03, 1.470E-03, 7.164E-04, 2.644E-04, 6.808E-05, 1.253E-05, 2.446E-06, &
8.701E-07, 3.893E-07, &
1.044E-02, 9.868E-03, 9.865E-03, 9.861E-03, 9.858E-03, 9.854E-03, 9.851E-03, 9.538E-03, 9.533E-03, 9.529E-03, &
9.523E-03, 9.515E-03, 9.504E-03, 9.490E-03, 9.471E-03, 9.445E-03, 9.411E-03, 9.363E-03, 9.298E-03, 9.209E-03, &
9.086E-03, 8.919E-03, 8.697E-03, 8.411E-03, 8.056E-03, 7.635E-03, 7.150E-03, 6.601E-03, 5.983E-03, 5.286E-03, &
4.504E-03, 3.647E-03, 2.751E-03, 1.882E-03, 1.122E-03, 5.497E-04, 2.048E-04, 5.344E-05, 1.023E-05, 2.236E-06, &
8.624E-07, 3.892E-07, &
7.963E-03, 7.556E-03, 7.554E-03, 7.551E-03, 7.549E-03, 7.546E-03, 7.544E-03, 7.308E-03, 7.305E-03, 7.302E-03, &
7.297E-03, 7.291E-03, 7.284E-03, 7.273E-03, 7.259E-03, 7.239E-03, 7.213E-03, 7.177E-03, 7.129E-03, 7.063E-03, &
6.973E-03, 6.851E-03, 6.689E-03, 6.477E-03, 6.209E-03, 5.885E-03, 5.506E-03, 5.074E-03, 4.587E-03, 4.039E-03, &
3.430E-03, 2.768E-03, 2.083E-03, 1.423E-03, 8.510E-04, 4.205E-04, 1.589E-04, 4.250E-05, 8.603E-06, 2.094E-06, &
8.574E-07, 3.891E-07, &
6.037E-03, 5.743E-03, 5.741E-03, 5.739E-03, 5.737E-03, 5.735E-03, 5.733E-03, 5.558E-03, 5.555E-03, 5.553E-03, &
5.549E-03, 5.545E-03, 5.539E-03, 5.530E-03, 5.519E-03, 5.504E-03, 5.483E-03, 5.456E-03, 5.419E-03, 5.369E-03, &
5.302E-03, 5.212E-03, 5.093E-03, 4.938E-03, 4.740E-03, 4.496E-03, 4.207E-03, 3.871E-03, 3.490E-03, 3.062E-03, &
2.589E-03, 2.081E-03, 1.561E-03, 1.067E-03, 6.409E-04, 3.201E-04, 1.233E-04, 3.411E-05, 7.410E-06, 1.996E-06, &
8.539E-07, 3.891E-07, &
4.541E-03, 4.325E-03, 4.323E-03, 4.322E-03, 4.320E-03, 4.319E-03, 4.317E-03, 4.189E-03, 4.187E-03, 4.185E-03, &
4.182E-03, 4.178E-03, 4.173E-03, 4.167E-03, 4.158E-03, 4.146E-03, 4.130E-03, 4.108E-03, 4.079E-03, 4.041E-03, &
3.991E-03, 3.925E-03, 3.838E-03, 3.725E-03, 3.581E-03, 3.402E-03, 3.186E-03, 2.931E-03, 2.637E-03, 2.305E-03, &
1.940E-03, 1.552E-03, 1.161E-03, 7.942E-04, 4.796E-04, 2.425E-04, 9.539E-05, 2.751E-05, 6.479E-06, 1.922E-06, &
8.513E-07, 3.890E-07, &
3.378E-03, 3.223E-03, 3.222E-03, 3.220E-03, 3.219E-03, 3.218E-03, 3.216E-03, 3.123E-03, 3.121E-03, 3.120E-03, &
3.117E-03, 3.114E-03, 3.111E-03, 3.105E-03, 3.098E-03, 3.089E-03, 3.076E-03, 3.059E-03, 3.037E-03, 3.008E-03, &
2.971E-03, 2.922E-03, 2.859E-03, 2.778E-03, 2.675E-03, 2.546E-03, 2.388E-03, 2.198E-03, 1.975E-03, 1.721E-03, &
1.443E-03, 1.150E-03, 8.580E-04, 5.872E-04, 3.565E-04, 1.823E-04, 7.329E-05, 2.211E-05, 5.689E-06, 1.857E-06, &
8.490E-07, 3.890E-07, &
2.470E-03, 2.363E-03, 2.362E-03, 2.361E-03, 2.360E-03, 2.359E-03, 2.358E-03, 2.290E-03, 2.289E-03, 2.288E-03, &
2.286E-03, 2.284E-03, 2.281E-03, 2.277E-03, 2.271E-03, 2.264E-03, 2.254E-03, 2.241E-03, 2.225E-03, 2.204E-03, &
2.176E-03, 2.141E-03, 2.096E-03, 2.038E-03, 1.965E-03, 1.874E-03, 1.761E-03, 1.623E-03, 1.459E-03, 1.270E-03, &
1.061E-03, 8.429E-04, 6.273E-04, 4.290E-04, 2.612E-04, 1.348E-04, 5.529E-05, 1.752E-05, 4.977E-06, 1.796E-06, &
8.467E-07, 3.889E-07, &
1.759E-03, 1.690E-03, 1.689E-03, 1.689E-03, 1.688E-03, 1.687E-03, 1.686E-03, 1.638E-03, 1.637E-03, 1.636E-03, &
1.634E-03, 1.633E-03, 1.630E-03, 1.628E-03, 1.624E-03, 1.618E-03, 1.611E-03, 1.602E-03, 1.590E-03, 1.575E-03, &
1.556E-03, 1.531E-03, 1.500E-03, 1.460E-03, 1.409E-03, 1.346E-03, 1.266E-03, 1.169E-03, 1.052E-03, 9.161E-04, &
7.650E-04, 6.062E-04, 4.495E-04, 3.062E-04, 1.862E-04, 9.649E-05, 4.038E-05, 1.356E-05, 4.328E-06, 1.737E-06, &
8.444E-07, 3.887E-07, &
1.210E-03, 1.169E-03, 1.169E-03, 1.168E-03, 1.168E-03, 1.167E-03, 1.167E-03, 1.133E-03, 1.132E-03, 1.131E-03, &
1.130E-03, 1.129E-03, 1.128E-03, 1.126E-03, 1.123E-03, 1.119E-03, 1.115E-03, 1.108E-03, 1.100E-03, 1.090E-03, &
1.077E-03, 1.061E-03, 1.040E-03, 1.013E-03, 9.785E-04, 9.352E-04, 8.811E-04, 8.145E-04, 7.343E-04, 6.405E-04, &
5.354E-04, 4.238E-04, 3.130E-04, 2.119E-04, 1.281E-04, 6.645E-05, 2.845E-05, 1.029E-05, 3.767E-06, 1.684E-06, &
8.421E-07, 3.885E-07, &
8.046E-04, 7.815E-04, 7.811E-04, 7.808E-04, 7.805E-04, 7.802E-04, 7.799E-04, 7.572E-04, 7.567E-04, 7.562E-04, &
7.557E-04, 7.549E-04, 7.540E-04, 7.527E-04, 7.509E-04, 7.486E-04, 7.456E-04, 7.415E-04, 7.364E-04, 7.299E-04, &
7.215E-04, 7.108E-04, 6.969E-04, 6.792E-04, 6.566E-04, 6.280E-04, 5.923E-04, 5.484E-04, 4.956E-04, 4.336E-04, &
3.635E-04, 2.882E-04, 2.126E-04, 1.432E-04, 8.610E-05, 4.483E-05, 1.982E-05, 7.862E-06, 3.331E-06, 1.640E-06, &
8.397E-07, 3.881E-07, &
5.240E-04, 5.112E-04, 5.110E-04, 5.108E-04, 5.106E-04, 5.104E-04, 5.103E-04, 4.958E-04, 4.955E-04, 4.952E-04, &
4.948E-04, 4.943E-04, 4.937E-04, 4.929E-04, 4.919E-04, 4.904E-04, 4.885E-04, 4.859E-04, 4.827E-04, 4.785E-04, &
4.731E-04, 4.662E-04, 4.573E-04, 4.459E-04, 4.313E-04, 4.129E-04, 3.899E-04, 3.619E-04, 3.281E-04, 2.884E-04, &
2.432E-04, 1.941E-04, 1.440E-04, 9.745E-05, 5.885E-05, 3.108E-05, 1.437E-05, 6.317E-06, 3.040E-06, 1.608E-06, &
8.374E-07, 3.875E-07 /), (/ 42,47 /))
! -------- Shielding of N2 from Visser et al. (2009) ---------
real(double_precision), parameter, dimension(46) :: nn2_n2 = (/ &
1.000e+10, 1.585e+10, 2.512e+10, 3.981e+10, 6.310e+10, 1.000e+11, 1.585e+11, 2.512e+11, 3.981e+11, &
6.310e+11, 1.000e+12, 1.585e+12, 2.512e+12, 3.981e+12, 6.310e+12, 1.000e+13, 1.585e+13, 2.512e+13, &
3.981e+13, 6.310e+13, 1.000e+14, 1.585e+14, 2.512e+14, 3.981e+14, 6.310e+14, 1.000e+15, 1.585e+15, &
2.512e+15, 3.981e+15, 6.310e+15, 1.000e+16, 1.585e+16, 2.512e+16, 3.981e+16, 6.310e+16, 1.000e+17, &
1.585e+17, 2.512e+17, 3.981e+17, 6.310e+17, 1.000e+18, 1.585e+18, 2.512e+18, 3.981e+18, 6.310e+18, &
1.000e+19/)
real(double_precision), parameter, dimension(46) :: nh2_n2 = (/ &
1.000e+14, 1.585e+14, 2.512e+14, 3.981e+14, 6.310e+14, 1.000e+15, 1.585e+15, 2.512e+15, 3.981e+15, &
6.310e+15, 1.000e+16, 1.585e+16, 2.512e+16, 3.981e+16, 6.310e+16, 1.000e+17, 1.585e+17, 2.512e+17, &
3.981e+17, 6.310e+17, 1.000e+18, 1.585e+18, 2.512e+18, 3.981e+18, 6.310e+18, 1.000e+19, 1.585e+19, &
2.512e+19, 3.981e+19, 6.310e+19, 1.000e+20, 1.585e+20, 2.512e+20, 3.981e+20, 6.310e+20, 1.000e+21, &
1.585e+21, 2.512e+21, 3.981e+21, 6.310e+21, 1.000e+22, 1.585e+22, 2.512e+22, 3.981e+22, 6.310e+22, &
1.000e+23/)
!real(double_precision), parameter, dimension(46,46) :: shielding_function_n2_nh1e14 = reshape((/ &
!9.981E-01, 9.974E-01, 9.966E-01, 9.956E-01, 9.947E-01, 9.937E-01, 9.929E-01, 9.921E-01, 9.914E-01, 9.909E-01, &
!9.905E-01, 9.900E-01, 9.894E-01, 9.885E-01, 9.872E-01, 9.854E-01, 9.829E-01, 9.797E-01, 9.759E-01, 9.712E-01, &
!9.653E-01, 9.573E-01, 9.462E-01, 9.308E-01, 9.094E-01, 8.799E-01, 8.404E-01, 7.898E-01, 7.284E-01, 6.582E-01, &
!5.825E-01, 5.069E-01, 4.398E-01, 3.862E-01, 3.403E-01, 2.903E-01, 2.293E-01, 1.601E-01, 9.273E-02, 4.133E-02, &
!1.339E-02, 3.307E-03, 6.703E-04, 8.300E-05, 3.461E-06, 2.305E-08, &
!9.980E-01, 9.973E-01, 9.965E-01, 9.955E-01, 9.945E-01, 9.936E-01, 9.928E-01, 9.920E-01, 9.913E-01, 9.908E-01, &
!9.904E-01, 9.899E-01, 9.893E-01, 9.884E-01, 9.871E-01, 9.852E-01, 9.828E-01, 9.796E-01, 9.758E-01, 9.711E-01, &
!9.652E-01, 9.571E-01, 9.460E-01, 9.306E-01, 9.093E-01, 8.798E-01, 8.403E-01, 7.896E-01, 7.283E-01, 6.581E-01, &
!5.824E-01, 5.068E-01, 4.397E-01, 3.861E-01, 3.403E-01, 2.902E-01, 2.293E-01, 1.600E-01, 9.271E-02, 4.133E-02, &
!1.339E-02, 3.306E-03, 6.703E-04, 8.300E-05, 3.461E-06, 2.305E-08, &
!9.978E-01, 9.971E-01, 9.963E-01, 9.953E-01, 9.943E-01, 9.934E-01, 9.926E-01, 9.918E-01, 9.911E-01, 9.906E-01, &
!9.902E-01, 9.897E-01, 9.891E-01, 9.882E-01, 9.869E-01, 9.850E-01, 9.826E-01, 9.794E-01, 9.756E-01, 9.709E-01, &
!9.650E-01, 9.569E-01, 9.458E-01, 9.305E-01, 9.091E-01, 8.796E-01, 8.401E-01, 7.895E-01, 7.281E-01, 6.579E-01, &
!5.822E-01, 5.067E-01, 4.396E-01, 3.860E-01, 3.402E-01, 2.901E-01, 2.292E-01, 1.600E-01, 9.268E-02, 4.131E-02, &
!1.339E-02, 3.305E-03, 6.702E-04, 8.300E-05, 3.461E-06, 2.305E-08, &
!9.975E-01, 9.968E-01, 9.959E-01, 9.950E-01, 9.940E-01, 9.931E-01, 9.922E-01, 9.914E-01, 9.908E-01, 9.903E-01, &
!9.899E-01, 9.894E-01, 9.888E-01, 9.878E-01, 9.866E-01, 9.847E-01, 9.823E-01, 9.791E-01, 9.753E-01, 9.706E-01, &
!9.647E-01, 9.566E-01, 9.455E-01, 9.302E-01, 9.088E-01, 8.793E-01, 8.398E-01, 7.892E-01, 7.278E-01, 6.577E-01, &
!5.820E-01, 5.065E-01, 4.394E-01, 3.858E-01, 3.400E-01, 2.900E-01, 2.291E-01, 1.599E-01, 9.263E-02, 4.129E-02, &
!1.338E-02, 3.304E-03, 6.700E-04, 8.299E-05, 3.461E-06, 2.305E-08, &
!9.970E-01, 9.963E-01, 9.954E-01, 9.945E-01, 9.935E-01, 9.926E-01, 9.917E-01, 9.909E-01, 9.903E-01, 9.898E-01, &
!9.894E-01, 9.889E-01, 9.883E-01, 9.874E-01, 9.861E-01, 9.842E-01, 9.818E-01, 9.786E-01, 9.748E-01, 9.701E-01, &
!9.642E-01, 9.562E-01, 9.451E-01, 9.297E-01, 9.083E-01, 8.789E-01, 8.394E-01, 7.888E-01, 7.274E-01, 6.573E-01, &
!5.816E-01, 5.061E-01, 4.391E-01, 3.856E-01, 3.398E-01, 2.898E-01, 2.290E-01, 1.598E-01, 9.255E-02, 4.125E-02, &
!1.337E-02, 3.302E-03, 6.698E-04, 8.298E-05, 3.460E-06, 2.305E-08, &
!9.962E-01, 9.955E-01, 9.947E-01, 9.937E-01, 9.927E-01, 9.918E-01, 9.910E-01, 9.902E-01, 9.895E-01, 9.890E-01, &
!9.886E-01, 9.881E-01, 9.875E-01, 9.866E-01, 9.853E-01, 9.835E-01, 9.810E-01, 9.779E-01, 9.740E-01, 9.693E-01, &
!9.634E-01, 9.554E-01, 9.443E-01, 9.289E-01, 9.076E-01, 8.781E-01, 8.387E-01, 7.881E-01, 7.268E-01, 6.567E-01, &
!5.811E-01, 5.056E-01, 4.386E-01, 3.851E-01, 3.394E-01, 2.895E-01, 2.287E-01, 1.596E-01, 9.243E-02, 4.120E-02, &
!1.335E-02, 3.299E-03, 6.695E-04, 8.296E-05, 3.460E-06, 2.304E-08, &
!9.950E-01, 9.943E-01, 9.934E-01, 9.925E-01, 9.915E-01, 9.906E-01, 9.897E-01, 9.889E-01, 9.883E-01, 9.878E-01, &
!9.873E-01, 9.869E-01, 9.862E-01, 9.853E-01, 9.841E-01, 9.822E-01, 9.798E-01, 9.766E-01, 9.728E-01, 9.681E-01, &
!9.622E-01, 9.542E-01, 9.431E-01, 9.277E-01, 9.064E-01, 8.770E-01, 8.375E-01, 7.870E-01, 7.257E-01, 6.557E-01, &
!5.802E-01, 5.048E-01, 4.379E-01, 3.845E-01, 3.388E-01, 2.889E-01, 2.283E-01, 1.593E-01, 9.224E-02, 4.111E-02, &
!1.332E-02, 3.293E-03, 6.690E-04, 8.293E-05, 3.459E-06, 2.303E-08, &
!9.930E-01, 9.923E-01, 9.915E-01, 9.905E-01, 9.896E-01, 9.886E-01, 9.878E-01, 9.870E-01, 9.863E-01, 9.858E-01, &
!9.854E-01, 9.850E-01, 9.843E-01, 9.834E-01, 9.821E-01, 9.803E-01, 9.779E-01, 9.747E-01, 9.709E-01, 9.662E-01, &
!9.603E-01, 9.523E-01, 9.412E-01, 9.259E-01, 9.046E-01, 8.752E-01, 8.358E-01, 7.853E-01, 7.241E-01, 6.542E-01, &
!5.788E-01, 5.036E-01, 4.368E-01, 3.835E-01, 3.379E-01, 2.881E-01, 2.276E-01, 1.588E-01, 9.195E-02, 4.097E-02, &
!1.328E-02, 3.285E-03, 6.682E-04, 8.288E-05, 3.457E-06, 2.302E-08, &
!9.900E-01, 9.893E-01, 9.885E-01, 9.875E-01, 9.865E-01, 9.856E-01, 9.848E-01, 9.840E-01, 9.833E-01, 9.828E-01, &
!9.824E-01, 9.819E-01, 9.813E-01, 9.804E-01, 9.791E-01, 9.773E-01, 9.749E-01, 9.718E-01, 9.679E-01, 9.633E-01, &
!9.574E-01, 9.494E-01, 9.383E-01, 9.230E-01, 9.017E-01, 8.724E-01, 8.330E-01, 7.826E-01, 7.216E-01, 6.518E-01, &
!5.766E-01, 5.016E-01, 4.350E-01, 3.819E-01, 3.364E-01, 2.869E-01, 2.266E-01, 1.580E-01, 9.150E-02, 4.076E-02, &
!1.321E-02, 3.273E-03, 6.671E-04, 8.281E-05, 3.455E-06, 2.301E-08, &
!9.853E-01, 9.846E-01, 9.838E-01, 9.828E-01, 9.819E-01, 9.810E-01, 9.801E-01, 9.793E-01, 9.787E-01, 9.782E-01, &
!9.778E-01, 9.773E-01, 9.767E-01, 9.758E-01, 9.745E-01, 9.727E-01, 9.703E-01, 9.672E-01, 9.634E-01, 9.587E-01, &
!9.528E-01, 9.449E-01, 9.338E-01, 9.185E-01, 8.973E-01, 8.680E-01, 8.287E-01, 7.785E-01, 7.176E-01, 6.482E-01, &
!5.733E-01, 4.986E-01, 4.323E-01, 3.794E-01, 3.342E-01, 2.849E-01, 2.250E-01, 1.569E-01, 9.081E-02, 4.044E-02, &
!1.311E-02, 3.255E-03, 6.653E-04, 8.269E-05, 3.450E-06, 2.298E-08, &
!9.782E-01, 9.775E-01, 9.767E-01, 9.757E-01, 9.748E-01, 9.738E-01, 9.730E-01, 9.722E-01, 9.716E-01, 9.711E-01, &
!9.707E-01, 9.702E-01, 9.696E-01, 9.687E-01, 9.674E-01, 9.657E-01, 9.633E-01, 9.602E-01, 9.564E-01, 9.518E-01, &
!9.459E-01, 9.380E-01, 9.269E-01, 9.117E-01, 8.905E-01, 8.614E-01, 8.223E-01, 7.722E-01, 7.117E-01, 6.426E-01, &
!5.683E-01, 4.941E-01, 4.283E-01, 3.758E-01, 3.309E-01, 2.820E-01, 2.226E-01, 1.552E-01, 8.979E-02, 3.998E-02, &
!1.297E-02, 3.229E-03, 6.626E-04, 8.252E-05, 3.444E-06, 2.294E-08, &
!9.674E-01, 9.668E-01, 9.659E-01, 9.650E-01, 9.640E-01, 9.631E-01, 9.623E-01, 9.615E-01, 9.609E-01, 9.604E-01, &
!9.600E-01, 9.595E-01, 9.589E-01, 9.580E-01, 9.568E-01, 9.551E-01, 9.527E-01, 9.496E-01, 9.459E-01, 9.413E-01, &
!9.355E-01, 9.276E-01, 9.166E-01, 9.014E-01, 8.804E-01, 8.514E-01, 8.125E-01, 7.628E-01, 7.028E-01, 6.344E-01, &
!5.607E-01, 4.874E-01, 4.223E-01, 3.704E-01, 3.260E-01, 2.778E-01, 2.192E-01, 1.527E-01, 8.831E-02, 3.930E-02, &
!1.276E-02, 3.192E-03, 6.588E-04, 8.224E-05, 3.434E-06, 2.287E-08, &
!9.516E-01, 9.509E-01, 9.501E-01, 9.492E-01, 9.482E-01, 9.473E-01, 9.465E-01, 9.458E-01, 9.451E-01, 9.446E-01, &
!9.442E-01, 9.438E-01, 9.432E-01, 9.423E-01, 9.411E-01, 9.394E-01, 9.371E-01, 9.341E-01, 9.304E-01, 9.259E-01, &
!9.201E-01, 9.123E-01, 9.014E-01, 8.864E-01, 8.655E-01, 8.367E-01, 7.982E-01, 7.491E-01, 6.898E-01, 6.223E-01, &
!5.498E-01, 4.776E-01, 4.136E-01, 3.626E-01, 3.190E-01, 2.717E-01, 2.143E-01, 1.492E-01, 8.623E-02, 3.837E-02, &
!1.248E-02, 3.142E-03, 6.535E-04, 8.183E-05, 3.418E-06, 2.277E-08, &
!9.287E-01, 9.281E-01, 9.273E-01, 9.264E-01, 9.254E-01, 9.246E-01, 9.238E-01, 9.230E-01, 9.224E-01, 9.219E-01, &
!9.215E-01, 9.211E-01, 9.205E-01, 9.197E-01, 9.185E-01, 9.169E-01, 9.146E-01, 9.117E-01, 9.082E-01, 9.037E-01, &
!8.980E-01, 8.903E-01, 8.795E-01, 8.647E-01, 8.440E-01, 8.156E-01, 7.777E-01, 7.293E-01, 6.711E-01, 6.051E-01, &
!5.342E-01, 4.639E-01, 4.015E-01, 3.517E-01, 3.093E-01, 2.632E-01, 2.075E-01, 1.444E-01, 8.339E-02, 3.712E-02, &
!1.212E-02, 3.078E-03, 6.463E-04, 8.120E-05, 3.393E-06, 2.260E-08, &
!8.967E-01, 8.961E-01, 8.953E-01, 8.944E-01, 8.935E-01, 8.927E-01, 8.919E-01, 8.912E-01, 8.905E-01, 8.901E-01, &
!8.897E-01, 8.893E-01, 8.888E-01, 8.880E-01, 8.869E-01, 8.853E-01, 8.832E-01, 8.804E-01, 8.770E-01, 8.727E-01, &
!8.671E-01, 8.595E-01, 8.489E-01, 8.343E-01, 8.140E-01, 7.862E-01, 7.491E-01, 7.019E-01, 6.452E-01, 5.811E-01, &
!5.127E-01, 4.449E-01, 3.848E-01, 3.368E-01, 2.959E-01, 2.517E-01, 1.982E-01, 1.379E-01, 7.959E-02, 3.546E-02, &
!1.164E-02, 2.996E-03, 6.363E-04, 8.023E-05, 3.355E-06, 2.235E-08, &
!8.534E-01, 8.528E-01, 8.521E-01, 8.512E-01, 8.504E-01, 8.496E-01, 8.488E-01, 8.481E-01, 8.475E-01, 8.471E-01, &
!8.467E-01, 8.464E-01, 8.459E-01, 8.451E-01, 8.441E-01, 8.427E-01, 8.407E-01, 8.381E-01, 8.349E-01, 8.308E-01, &
!8.254E-01, 8.180E-01, 8.077E-01, 7.934E-01, 7.737E-01, 7.466E-01, 7.106E-01, 6.649E-01, 6.104E-01, 5.492E-01, &
!4.841E-01, 4.197E-01, 3.626E-01, 3.171E-01, 2.783E-01, 2.365E-01, 1.861E-01, 1.293E-01, 7.465E-02, 3.333E-02, &
!1.104E-02, 2.889E-03, 6.224E-04, 7.878E-05, 3.295E-06, 2.195E-08, &
!7.978E-01, 7.973E-01, 7.965E-01, 7.958E-01, 7.950E-01, 7.942E-01, 7.935E-01, 7.929E-01, 7.923E-01, 7.919E-01, &
!7.916E-01, 7.912E-01, 7.908E-01, 7.901E-01, 7.892E-01, 7.880E-01, 7.862E-01, 7.839E-01, 7.808E-01, 7.770E-01, &
!7.718E-01, 7.647E-01, 7.548E-01, 7.410E-01, 7.220E-01, 6.959E-01, 6.614E-01, 6.178E-01, 5.662E-01, 5.085E-01, &
!4.477E-01, 3.878E-01, 3.348E-01, 2.924E-01, 2.563E-01, 2.175E-01, 1.710E-01, 1.187E-01, 6.849E-02, 3.068E-02, &
!1.028E-02, 2.751E-03, 6.026E-04, 7.659E-05, 3.205E-06, 2.136E-08, &
!7.307E-01, 7.302E-01, 7.295E-01, 7.288E-01, 7.281E-01, 7.274E-01, 7.268E-01, 7.262E-01, 7.257E-01, 7.254E-01, &
!7.251E-01, 7.248E-01, 7.244E-01, 7.238E-01, 7.230E-01, 7.219E-01, 7.204E-01, 7.184E-01, 7.157E-01, 7.121E-01, &
!7.072E-01, 7.005E-01, 6.910E-01, 6.779E-01, 6.598E-01, 6.350E-01, 6.024E-01, 5.615E-01, 5.134E-01, 4.603E-01, &
!4.047E-01, 3.503E-01, 3.021E-01, 2.635E-01, 2.306E-01, 1.954E-01, 1.533E-01, 1.063E-01, 6.131E-02, 2.756E-02, &
!9.363E-03, 2.573E-03, 5.745E-04, 7.335E-05, 3.071E-06, 2.047E-08, &
!6.549E-01, 6.545E-01, 6.539E-01, 6.533E-01, 6.527E-01, 6.521E-01, 6.516E-01, 6.511E-01, 6.506E-01, 6.503E-01, &
!6.500E-01, 6.498E-01, 6.494E-01, 6.490E-01, 6.483E-01, 6.474E-01, 6.462E-01, 6.444E-01, 6.420E-01, 6.388E-01, &
!6.343E-01, 6.280E-01, 6.192E-01, 6.069E-01, 5.898E-01, 5.667E-01, 5.364E-01, 4.987E-01, 4.548E-01, 4.069E-01, &
!3.574E-01, 3.092E-01, 2.666E-01, 2.323E-01, 2.029E-01, 1.716E-01, 1.343E-01, 9.292E-02, 5.355E-02, 2.415E-02, &
!8.327E-03, 2.354E-03, 5.360E-04, 6.875E-05, 2.879E-06, 1.920E-08, &
!5.746E-01, 5.742E-01, 5.737E-01, 5.732E-01, 5.727E-01, 5.722E-01, 5.718E-01, 5.714E-01, 5.710E-01, 5.707E-01, &
!5.705E-01, 5.703E-01, 5.700E-01, 5.696E-01, 5.691E-01, 5.683E-01, 5.673E-01, 5.658E-01, 5.637E-01, 5.608E-01, &
!5.568E-01, 5.510E-01, 5.429E-01, 5.315E-01, 5.159E-01, 4.947E-01, 4.670E-01, 4.330E-01, 3.939E-01, 3.517E-01, &
!3.088E-01, 2.674E-01, 2.307E-01, 2.008E-01, 1.751E-01, 1.478E-01, 1.154E-01, 7.960E-02, 4.578E-02, 2.068E-02, &
!7.220E-03, 2.094E-03, 4.855E-04, 6.253E-05, 2.619E-06, 1.747E-08, &
!4.935E-01, 4.931E-01, 4.928E-01, 4.924E-01, 4.920E-01, 4.916E-01, 4.913E-01, 4.910E-01, 4.907E-01, 4.904E-01, &
!4.902E-01, 4.900E-01, 4.898E-01, 4.894E-01, 4.890E-01, 4.884E-01, 4.875E-01, 4.862E-01, 4.844E-01, 4.818E-01, &
!4.782E-01, 4.731E-01, 4.658E-01, 4.555E-01, 4.414E-01, 4.224E-01, 3.979E-01, 3.679E-01, 3.339E-01, 2.977E-01, &
!2.614E-01, 2.268E-01, 1.960E-01, 1.707E-01, 1.486E-01, 1.251E-01, 9.743E-02, 6.698E-02, 3.840E-02, 1.733E-02, &
!6.097E-03, 1.802E-03, 4.236E-04, 5.472E-05, 2.292E-06, 1.529E-08, &
!4.141E-01, 4.138E-01, 4.135E-01, 4.132E-01, 4.129E-01, 4.126E-01, 4.124E-01, 4.121E-01, 4.119E-01, 4.117E-01, &
!4.115E-01, 4.113E-01, 4.111E-01, 4.108E-01, 4.104E-01, 4.098E-01, 4.091E-01, 4.080E-01, 4.064E-01, 4.042E-01, &
!4.010E-01, 3.965E-01, 3.901E-01, 3.811E-01, 3.688E-01, 3.523E-01, 3.311E-01, 3.056E-01, 2.768E-01, 2.466E-01, &
!2.166E-01, 1.883E-01, 1.631E-01, 1.421E-01, 1.236E-01, 1.038E-01, 8.066E-02, 5.527E-02, 3.156E-02, 1.420E-02, &
!5.007E-03, 1.495E-03, 3.543E-04, 4.581E-05, 1.917E-06, 1.279E-08, &
!3.379E-01, 3.376E-01, 3.374E-01, 3.371E-01, 3.369E-01, 3.366E-01, 3.364E-01, 3.362E-01, 3.360E-01, 3.358E-01, &
!3.356E-01, 3.355E-01, 3.353E-01, 3.350E-01, 3.347E-01, 3.342E-01, 3.335E-01, 3.326E-01, 3.313E-01, 3.293E-01, &
!3.266E-01, 3.228E-01, 3.173E-01, 3.097E-01, 2.994E-01, 2.857E-01, 2.682E-01, 2.473E-01, 2.238E-01, 1.993E-01, &
!1.751E-01, 1.523E-01, 1.321E-01, 1.152E-01, 1.001E-01, 8.390E-02, 6.502E-02, 4.442E-02, 2.528E-02, 1.134E-02, &
!3.995E-03, 1.197E-03, 2.843E-04, 3.673E-05, 1.535E-06, 1.023E-08, &
!2.666E-01, 2.664E-01, 2.662E-01, 2.659E-01, 2.657E-01, 2.655E-01, 2.653E-01, 2.651E-01, 2.649E-01, 2.648E-01, &
!2.646E-01, 2.645E-01, 2.643E-01, 2.641E-01, 2.638E-01, 2.633E-01, 2.628E-01, 2.620E-01, 2.608E-01, 2.592E-01, &
!2.570E-01, 2.538E-01, 2.493E-01, 2.433E-01, 2.350E-01, 2.243E-01, 2.107E-01, 1.945E-01, 1.763E-01, 1.570E-01, &
!1.379E-01, 1.199E-01, 1.039E-01, 9.053E-02, 7.856E-02, 6.576E-02, 5.087E-02, 3.468E-02, 1.970E-02, 8.825E-03, &
!3.112E-03, 9.333E-04, 2.214E-04, 2.850E-05, 1.188E-06, 7.907E-09, &
!2.030E-01, 2.028E-01, 2.026E-01, 2.024E-01, 2.022E-01, 2.020E-01, 2.019E-01, 2.017E-01, 2.016E-01, 2.014E-01, &
!2.013E-01, 2.011E-01, 2.010E-01, 2.007E-01, 2.005E-01, 2.001E-01, 1.997E-01, 1.990E-01, 1.981E-01, 1.968E-01, &
!1.949E-01, 1.924E-01, 1.890E-01, 1.843E-01, 1.782E-01, 1.703E-01, 1.605E-01, 1.488E-01, 1.354E-01, 1.210E-01, &
!1.063E-01, 9.228E-02, 7.987E-02, 6.945E-02, 6.016E-02, 5.028E-02, 3.883E-02, 2.644E-02, 1.502E-02, 6.748E-03, &
!2.392E-03, 7.200E-04, 1.703E-04, 2.181E-05, 9.059E-07, 6.013E-09, &
!1.505E-01, 1.503E-01, 1.501E-01, 1.499E-01, 1.497E-01, 1.496E-01, 1.494E-01, 1.493E-01, 1.491E-01, 1.490E-01, &
!1.489E-01, 1.488E-01, 1.486E-01, 1.484E-01, 1.482E-01, 1.479E-01, 1.475E-01, 1.469E-01, 1.462E-01, 1.451E-01, &
!1.437E-01, 1.418E-01, 1.392E-01, 1.358E-01, 1.315E-01, 1.261E-01, 1.194E-01, 1.114E-01, 1.022E-01, 9.186E-02, &
!8.091E-02, 7.020E-02, 6.066E-02, 5.265E-02, 4.551E-02, 3.796E-02, 2.927E-02, 1.992E-02, 1.133E-02, 5.120E-03, &
!1.834E-03, 5.553E-04, 1.310E-04, 1.669E-05, 6.906E-07, 4.570E-09, &
!1.105E-01, 1.104E-01, 1.102E-01, 1.100E-01, 1.098E-01, 1.097E-01, 1.096E-01, 1.094E-01, 1.093E-01, 1.092E-01, &
!1.091E-01, 1.089E-01, 1.088E-01, 1.087E-01, 1.085E-01, 1.082E-01, 1.079E-01, 1.074E-01, 1.068E-01, 1.059E-01, &
!1.048E-01, 1.033E-01, 1.014E-01, 9.901E-02, 9.602E-02, 9.237E-02, 8.798E-02, 8.275E-02, 7.654E-02, 6.929E-02, &
!6.130E-02, 5.326E-02, 4.598E-02, 3.984E-02, 3.436E-02, 2.859E-02, 2.200E-02, 1.495E-02, 8.519E-03, 3.874E-03, &
!1.402E-03, 4.273E-04, 1.005E-04, 1.275E-05, 5.258E-07, 3.470E-09, &
!8.157E-02, 8.141E-02, 8.124E-02, 8.108E-02, 8.093E-02, 8.080E-02, 8.068E-02, 8.057E-02, 8.046E-02, 8.035E-02, &
!8.025E-02, 8.014E-02, 8.003E-02, 7.989E-02, 7.972E-02, 7.950E-02, 7.922E-02, 7.884E-02, 7.833E-02, 7.766E-02, &
!7.678E-02, 7.564E-02, 7.421E-02, 7.247E-02, 7.037E-02, 6.790E-02, 6.498E-02, 6.149E-02, 5.725E-02, 5.215E-02, &
!4.635E-02, 4.037E-02, 3.486E-02, 3.015E-02, 2.593E-02, 2.151E-02, 1.650E-02, 1.120E-02, 6.383E-03, 2.918E-03, &
!1.066E-03, 3.263E-04, 7.648E-05, 9.666E-06, 3.976E-07, 2.619E-09, &
!6.054E-02, 6.040E-02, 6.025E-02, 6.010E-02, 5.997E-02, 5.986E-02, 5.976E-02, 5.966E-02, 5.956E-02, 5.947E-02, &
!5.938E-02, 5.929E-02, 5.919E-02, 5.907E-02, 5.893E-02, 5.875E-02, 5.852E-02, 5.821E-02, 5.780E-02, 5.727E-02, &
!5.659E-02, 5.573E-02, 5.467E-02, 5.341E-02, 5.193E-02, 5.022E-02, 4.822E-02, 4.580E-02, 4.281E-02, 3.916E-02, &
!3.494E-02, 3.052E-02, 2.637E-02, 2.277E-02, 1.951E-02, 1.611E-02, 1.231E-02, 8.332E-03, 4.752E-03, 2.184E-03, &
!8.049E-04, 2.472E-04, 5.772E-05, 7.268E-06, 2.984E-07, 1.963E-09, &
!4.499E-02, 4.488E-02, 4.475E-02, 4.463E-02, 4.453E-02, 4.443E-02, 4.435E-02, 4.426E-02, 4.419E-02, 4.411E-02, &
!4.404E-02, 4.397E-02, 4.389E-02, 4.379E-02, 4.368E-02, 4.353E-02, 4.335E-02, 4.311E-02, 4.280E-02, 4.240E-02, &
!4.189E-02, 4.126E-02, 4.049E-02, 3.960E-02, 3.856E-02, 3.737E-02, 3.595E-02, 3.422E-02, 3.206E-02, 2.941E-02, &
!2.633E-02, 2.306E-02, 1.994E-02, 1.718E-02, 1.466E-02, 1.204E-02, 9.154E-03, 6.170E-03, 3.518E-03, 1.626E-03, &
!6.049E-04, 1.863E-04, 4.327E-05, 5.426E-06, 2.224E-07, 1.462E-09, &
!3.353E-02, 3.344E-02, 3.335E-02, 3.325E-02, 3.317E-02, 3.310E-02, 3.304E-02, 3.298E-02, 3.292E-02, 3.286E-02, &
!3.281E-02, 3.276E-02, 3.269E-02, 3.262E-02, 3.254E-02, 3.243E-02, 3.230E-02, 3.212E-02, 3.190E-02, 3.161E-02, &
!3.126E-02, 3.081E-02, 3.028E-02, 2.966E-02, 2.894E-02, 2.808E-02, 2.706E-02, 2.578E-02, 2.418E-02, 2.223E-02, &
!1.996E-02, 1.753E-02, 1.518E-02, 1.305E-02, 1.108E-02, 9.040E-03, 6.829E-03, 4.579E-03, 2.607E-03, 1.211E-03, &
!4.538E-04, 1.397E-04, 3.220E-05, 4.016E-06, 1.643E-07, 1.080E-09, &
!2.509E-02, 2.503E-02, 2.497E-02, 2.491E-02, 2.486E-02, 2.481E-02, 2.476E-02, 2.472E-02, 2.469E-02, 2.465E-02, &
!2.461E-02, 2.458E-02, 2.454E-02, 2.449E-02, 2.443E-02, 2.436E-02, 2.427E-02, 2.415E-02, 2.400E-02, 2.381E-02, &
!2.358E-02, 2.328E-02, 2.292E-02, 2.250E-02, 2.199E-02, 2.137E-02, 2.060E-02, 1.964E-02, 1.844E-02, 1.699E-02, &
!1.530E-02, 1.350E-02, 1.171E-02, 1.005E-02, 8.484E-03, 6.872E-03, 5.151E-03, 3.430E-03, 1.945E-03, 9.043E-04, &
!3.397E-04, 1.040E-04, 2.367E-05, 2.930E-06, 1.195E-07, 7.855E-10, &
!1.885E-02, 1.881E-02, 1.878E-02, 1.875E-02, 1.872E-02, 1.869E-02, 1.866E-02, 1.864E-02, 1.862E-02, 1.860E-02, &
!1.858E-02, 1.856E-02, 1.853E-02, 1.850E-02, 1.847E-02, 1.842E-02, 1.837E-02, 1.829E-02, 1.820E-02, 1.808E-02, &
!1.793E-02, 1.774E-02, 1.751E-02, 1.722E-02, 1.686E-02, 1.640E-02, 1.583E-02, 1.512E-02, 1.422E-02, 1.314E-02, &
!1.189E-02, 1.053E-02, 9.164E-03, 7.857E-03, 6.594E-03, 5.298E-03, 3.935E-03, 2.597E-03, 1.462E-03, 6.769E-04, &
!2.529E-04, 7.628E-05, 1.706E-05, 2.091E-06, 8.502E-08, 5.594E-10, &
!1.420E-02, 1.419E-02, 1.417E-02, 1.416E-02, 1.414E-02, 1.413E-02, 1.411E-02, 1.410E-02, 1.409E-02, 1.408E-02, &
!1.407E-02, 1.406E-02, 1.404E-02, 1.403E-02, 1.401E-02, 1.398E-02, 1.394E-02, 1.389E-02, 1.384E-02, 1.376E-02, &
!1.367E-02, 1.355E-02, 1.340E-02, 1.321E-02, 1.296E-02, 1.264E-02, 1.223E-02, 1.171E-02, 1.105E-02, 1.026E-02, &
!9.327E-03, 8.308E-03, 7.258E-03, 6.219E-03, 5.192E-03, 4.137E-03, 3.043E-03, 1.988E-03, 1.108E-03, 5.076E-04, &
!1.869E-04, 5.504E-05, 1.201E-05, 1.452E-06, 5.885E-08, 3.882E-10, &
!1.072E-02, 1.071E-02, 1.070E-02, 1.070E-02, 1.069E-02, 1.068E-02, 1.067E-02, 1.066E-02, 1.066E-02, 1.065E-02, &
!1.064E-02, 1.064E-02, 1.063E-02, 1.062E-02, 1.060E-02, 1.058E-02, 1.056E-02, 1.053E-02, 1.049E-02, 1.044E-02, &
!1.039E-02, 1.031E-02, 1.021E-02, 1.009E-02, 9.924E-03, 9.710E-03, 9.428E-03, 9.062E-03, 8.596E-03, 8.021E-03, &
!7.343E-03, 6.585E-03, 5.779E-03, 4.955E-03, 4.118E-03, 3.256E-03, 2.374E-03, 1.535E-03, 8.454E-04, 3.819E-04, &
!1.376E-04, 3.920E-05, 8.278E-06, 9.845E-07, 3.972E-08, 2.630E-10, &
!8.086E-03, 8.081E-03, 8.075E-03, 8.069E-03, 8.063E-03, 8.057E-03, 8.051E-03, 8.045E-03, 8.040E-03, 8.035E-03, &
!8.030E-03, 8.025E-03, 8.018E-03, 8.010E-03, 7.999E-03, 7.985E-03, 7.967E-03, 7.944E-03, 7.916E-03, 7.883E-03, &
!7.842E-03, 7.790E-03, 7.724E-03, 7.640E-03, 7.532E-03, 7.389E-03, 7.200E-03, 6.952E-03, 6.632E-03, 6.231E-03, &
!5.749E-03, 5.194E-03, 4.584E-03, 3.937E-03, 3.262E-03, 2.564E-03, 1.855E-03, 1.190E-03, 6.495E-04, 2.894E-04, &
!1.019E-04, 2.796E-05, 5.683E-06, 6.621E-07, 2.657E-08, 1.769E-10, &
!6.092E-03, 6.088E-03, 6.083E-03, 6.078E-03, 6.074E-03, 6.069E-03, 6.065E-03, 6.060E-03, 6.056E-03, 6.052E-03, &
!6.048E-03, 6.043E-03, 6.038E-03, 6.031E-03, 6.022E-03, 6.011E-03, 5.996E-03, 5.979E-03, 5.957E-03, 5.931E-03, &
!5.899E-03, 5.859E-03, 5.810E-03, 5.749E-03, 5.672E-03, 5.573E-03, 5.444E-03, 5.276E-03, 5.059E-03, 4.785E-03, &
!4.450E-03, 4.053E-03, 3.599E-03, 3.099E-03, 2.565E-03, 2.008E-03, 1.447E-03, 9.247E-04, 5.023E-04, 2.218E-04, &
!7.672E-05, 2.039E-05, 3.998E-06, 4.562E-07, 1.820E-08, 1.217E-10, &
!4.583E-03, 4.580E-03, 4.576E-03, 4.573E-03, 4.569E-03, 4.566E-03, 4.562E-03, 4.559E-03, 4.556E-03, 4.553E-03, &
!4.549E-03, 4.545E-03, 4.540E-03, 4.535E-03, 4.527E-03, 4.518E-03, 4.507E-03, 4.493E-03, 4.475E-03, 4.454E-03, &
!4.428E-03, 4.397E-03, 4.358E-03, 4.311E-03, 4.254E-03, 4.182E-03, 4.091E-03, 3.974E-03, 3.826E-03, 3.640E-03, &
!3.410E-03, 3.130E-03, 2.798E-03, 2.418E-03, 2.002E-03, 1.564E-03, 1.125E-03, 7.189E-04, 3.906E-04, 1.722E-04, &
!5.899E-05, 1.534E-05, 2.922E-06, 3.274E-07, 1.298E-08, 8.685E-11, &
!3.450E-03, 3.448E-03, 3.445E-03, 3.442E-03, 3.440E-03, 3.437E-03, 3.435E-03, 3.432E-03, 3.430E-03, 3.427E-03, &
!3.424E-03, 3.421E-03, 3.417E-03, 3.412E-03, 3.406E-03, 3.399E-03, 3.390E-03, 3.379E-03, 3.365E-03, 3.348E-03, &
!3.328E-03, 3.303E-03, 3.272E-03, 3.236E-03, 3.192E-03, 3.139E-03, 3.074E-03, 2.993E-03, 2.891E-03, 2.763E-03, &
!2.604E-03, 2.406E-03, 2.165E-03, 1.880E-03, 1.558E-03, 1.216E-03, 8.738E-04, 5.588E-04, 3.043E-04, 1.344E-04, &
!4.585E-05, 1.173E-05, 2.181E-06, 2.402E-07, 9.446E-09, 6.310E-11, &
!2.601E-03, 2.599E-03, 2.597E-03, 2.595E-03, 2.593E-03, 2.591E-03, 2.589E-03, 2.587E-03, 2.585E-03, 2.583E-03, &
!2.581E-03, 2.578E-03, 2.575E-03, 2.571E-03, 2.567E-03, 2.562E-03, 2.555E-03, 2.546E-03, 2.535E-03, 2.522E-03, &
!2.506E-03, 2.487E-03, 2.464E-03, 2.436E-03, 2.404E-03, 2.366E-03, 2.319E-03, 2.262E-03, 2.192E-03, 2.103E-03, &
!1.992E-03, 1.851E-03, 1.674E-03, 1.460E-03, 1.213E-03, 9.467E-04, 6.791E-04, 4.340E-04, 2.367E-04, 1.047E-04, &
!3.566E-05, 9.012E-06, 1.639E-06, 1.774E-07, 6.904E-09, 4.594E-11, &
!1.959E-03, 1.958E-03, 1.956E-03, 1.955E-03, 1.953E-03, 1.952E-03, 1.950E-03, 1.949E-03, 1.948E-03, 1.946E-03, &
!1.944E-03, 1.942E-03, 1.940E-03, 1.937E-03, 1.934E-03, 1.930E-03, 1.925E-03, 1.919E-03, 1.911E-03, 1.901E-03, &
!1.890E-03, 1.875E-03, 1.858E-03, 1.839E-03, 1.816E-03, 1.789E-03, 1.757E-03, 1.717E-03, 1.669E-03, 1.607E-03, &
!1.528E-03, 1.426E-03, 1.296E-03, 1.136E-03, 9.467E-04, 7.392E-04, 5.295E-04, 3.376E-04, 1.838E-04, 8.127E-05, &
!2.763E-05, 6.922E-06, 1.238E-06, 1.320E-07, 5.088E-09, 3.371E-11, &
!1.477E-03, 1.475E-03, 1.474E-03, 1.473E-03, 1.472E-03, 1.471E-03, 1.470E-03, 1.469E-03, 1.468E-03, 1.467E-03, &
!1.466E-03, 1.465E-03, 1.463E-03, 1.461E-03, 1.459E-03, 1.457E-03, 1.453E-03, 1.449E-03, 1.444E-03, 1.437E-03, &
!1.429E-03, 1.419E-03, 1.407E-03, 1.393E-03, 1.378E-03, 1.359E-03, 1.338E-03, 1.311E-03, 1.277E-03, 1.234E-03, &
!1.178E-03, 1.104E-03, 1.008E-03, 8.880E-04, 7.436E-04, 5.820E-04, 4.165E-04, 2.645E-04, 1.432E-04, 6.296E-05, &
!2.132E-05, 5.313E-06, 9.402E-07, 9.914E-08, 3.790E-09, 2.507E-11, &
!1.115E-03, 1.115E-03, 1.114E-03, 1.113E-03, 1.112E-03, 1.112E-03, 1.111E-03, 1.111E-03, 1.110E-03, 1.109E-03, &
!1.109E-03, 1.108E-03, 1.107E-03, 1.106E-03, 1.104E-03, 1.103E-03, 1.100E-03, 1.098E-03, 1.094E-03, 1.090E-03, &
!1.084E-03, 1.078E-03, 1.070E-03, 1.061E-03, 1.050E-03, 1.038E-03, 1.023E-03, 1.005E-03, 9.822E-04, 9.520E-04, &
!9.119E-04, 8.585E-04, 7.881E-04, 6.979E-04, 5.876E-04, 4.618E-04, 3.307E-04, 2.091E-04, 1.122E-04, 4.878E-05, &
!1.640E-05, 4.068E-06, 7.137E-07, 7.444E-08, 2.826E-09, 1.872E-11, &
!8.388E-04, 8.383E-04, 8.379E-04, 8.375E-04, 8.370E-04, 8.366E-04, 8.362E-04, 8.359E-04, 8.355E-04, 8.351E-04, &
!8.347E-04, 8.342E-04, 8.337E-04, 8.330E-04, 8.323E-04, 8.313E-04, 8.300E-04, 8.283E-04, 8.261E-04, 8.232E-04, &
!8.196E-04, 8.153E-04, 8.100E-04, 8.039E-04, 7.968E-04, 7.885E-04, 7.784E-04, 7.660E-04, 7.500E-04, 7.288E-04, &
!7.005E-04, 6.622E-04, 6.109E-04, 5.441E-04, 4.610E-04, 3.643E-04, 2.616E-04, 1.650E-04, 8.772E-05, 3.769E-05, &
!1.257E-05, 3.104E-06, 5.386E-07, 5.537E-08, 2.087E-09, 1.390E-11, &
!6.198E-04, 6.196E-04, 6.193E-04, 6.190E-04, 6.187E-04, 6.185E-04, 6.182E-04, 6.180E-04, 6.177E-04, 6.175E-04, &
!6.173E-04, 6.170E-04, 6.167E-04, 6.164E-04, 6.159E-04, 6.154E-04, 6.146E-04, 6.136E-04, 6.122E-04, 6.103E-04, &
!6.080E-04, 6.050E-04, 6.015E-04, 5.973E-04, 5.923E-04, 5.865E-04, 5.794E-04, 5.708E-04, 5.597E-04, 5.450E-04, &
!5.251E-04, 4.982E-04, 4.617E-04, 4.135E-04, 3.527E-04, 2.806E-04, 2.025E-04, 1.278E-04, 6.764E-05, 2.883E-05, &
!9.559E-06, 2.353E-06, 4.025E-07, 4.060E-08, 1.518E-09, 1.021E-11, &
!4.439E-04, 4.437E-04, 4.436E-04, 4.434E-04, 4.432E-04, 4.430E-04, 4.429E-04, 4.427E-04, 4.426E-04, 4.424E-04, &
!4.423E-04, 4.421E-04, 4.420E-04, 4.418E-04, 4.415E-04, 4.412E-04, 4.407E-04, 4.401E-04, 4.392E-04, 4.381E-04, &
!4.365E-04, 4.346E-04, 4.322E-04, 4.293E-04, 4.259E-04, 4.218E-04, 4.169E-04, 4.109E-04, 4.032E-04, 3.932E-04, &
!3.796E-04, 3.611E-04, 3.360E-04, 3.025E-04, 2.597E-04, 2.081E-04, 1.513E-04, 9.600E-05, 5.090E-05, 2.168E-05, &
!7.188E-06, 1.765E-06, 2.967E-07, 2.926E-08, 1.086E-09, 7.384E-12 /), (/ 46,46 /))
real(double_precision), parameter, dimension(46,46) :: shielding_function_n2_nh1e20 = reshape((/ &
9.478E-01, 9.471E-01, 9.463E-01, 9.453E-01, 9.444E-01, 9.435E-01, 9.426E-01, 9.418E-01, 9.412E-01, 9.407E-01, &
9.402E-01, 9.398E-01, 9.392E-01, 9.383E-01, 9.370E-01, 9.352E-01, 9.328E-01, 9.298E-01, 9.260E-01, 9.216E-01, &
9.159E-01, 9.084E-01, 8.979E-01, 8.835E-01, 8.635E-01, 8.360E-01, 7.989E-01, 7.511E-01, 6.927E-01, 6.252E-01, &
5.518E-01, 4.782E-01, 4.130E-01, 3.618E-01, 3.192E-01, 2.732E-01, 2.172E-01, 1.528E-01, 8.954E-02, 4.043E-02, &
1.326E-02, 3.300E-03, 6.701E-04, 8.298E-05, 3.460E-06, 2.305E-08, &
9.477E-01, 9.470E-01, 9.462E-01, 9.452E-01, 9.442E-01, 9.433E-01, 9.425E-01, 9.417E-01, 9.410E-01, 9.405E-01, &
9.401E-01, 9.397E-01, 9.390E-01, 9.382E-01, 9.369E-01, 9.351E-01, 9.327E-01, 9.296E-01, 9.259E-01, 9.214E-01, &
9.158E-01, 9.082E-01, 8.978E-01, 8.834E-01, 8.634E-01, 8.359E-01, 7.988E-01, 7.510E-01, 6.926E-01, 6.251E-01, &
5.517E-01, 4.781E-01, 4.129E-01, 3.618E-01, 3.191E-01, 2.732E-01, 2.171E-01, 1.528E-01, 8.952E-02, 4.043E-02, &
1.326E-02, 3.299E-03, 6.700E-04, 8.298E-05, 3.460E-06, 2.304E-08, &
9.475E-01, 9.468E-01, 9.460E-01, 9.450E-01, 9.440E-01, 9.431E-01, 9.423E-01, 9.415E-01, 9.408E-01, 9.403E-01, &
9.399E-01, 9.395E-01, 9.389E-01, 9.380E-01, 9.367E-01, 9.349E-01, 9.325E-01, 9.294E-01, 9.257E-01, 9.212E-01, &
9.156E-01, 9.080E-01, 8.976E-01, 8.832E-01, 8.632E-01, 8.357E-01, 7.986E-01, 7.508E-01, 6.924E-01, 6.250E-01, &
5.516E-01, 4.780E-01, 4.128E-01, 3.617E-01, 3.190E-01, 2.731E-01, 2.170E-01, 1.528E-01, 8.949E-02, 4.041E-02, &
1.326E-02, 3.298E-03, 6.699E-04, 8.297E-05, 3.460E-06, 2.304E-08, &
9.472E-01, 9.465E-01, 9.456E-01, 9.447E-01, 9.437E-01, 9.428E-01, 9.420E-01, 9.412E-01, 9.405E-01, 9.400E-01, &
9.396E-01, 9.392E-01, 9.385E-01, 9.377E-01, 9.364E-01, 9.346E-01, 9.322E-01, 9.291E-01, 9.254E-01, 9.209E-01, &
9.153E-01, 9.077E-01, 8.973E-01, 8.829E-01, 8.629E-01, 8.354E-01, 7.983E-01, 7.505E-01, 6.922E-01, 6.247E-01, &
5.513E-01, 4.778E-01, 4.126E-01, 3.615E-01, 3.189E-01, 2.729E-01, 2.169E-01, 1.527E-01, 8.944E-02, 4.039E-02, &
1.325E-02, 3.297E-03, 6.698E-04, 8.297E-05, 3.460E-06, 2.304E-08, &
9.467E-01, 9.460E-01, 9.452E-01, 9.442E-01, 9.432E-01, 9.423E-01, 9.415E-01, 9.407E-01, 9.400E-01, 9.395E-01, &
9.391E-01, 9.387E-01, 9.380E-01, 9.372E-01, 9.359E-01, 9.341E-01, 9.317E-01, 9.286E-01, 9.249E-01, 9.205E-01, &
9.148E-01, 9.073E-01, 8.968E-01, 8.824E-01, 8.625E-01, 8.349E-01, 7.979E-01, 7.501E-01, 6.917E-01, 6.244E-01, &
5.510E-01, 4.774E-01, 4.124E-01, 3.612E-01, 3.186E-01, 2.727E-01, 2.168E-01, 1.526E-01, 8.937E-02, 4.035E-02, &
1.324E-02, 3.295E-03, 6.696E-04, 8.296E-05, 3.459E-06, 2.304E-08, &
9.459E-01, 9.452E-01, 9.444E-01, 9.434E-01, 9.425E-01, 9.416E-01, 9.407E-01, 9.399E-01, 9.392E-01, 9.388E-01, &
9.383E-01, 9.379E-01, 9.373E-01, 9.364E-01, 9.351E-01, 9.333E-01, 9.309E-01, 9.279E-01, 9.242E-01, 9.197E-01, &
9.141E-01, 9.065E-01, 8.961E-01, 8.817E-01, 8.617E-01, 8.342E-01, 7.971E-01, 7.494E-01, 6.911E-01, 6.237E-01, &
5.504E-01, 4.769E-01, 4.119E-01, 3.608E-01, 3.183E-01, 2.724E-01, 2.165E-01, 1.524E-01, 8.925E-02, 4.030E-02, &
1.322E-02, 3.292E-03, 6.693E-04, 8.294E-05, 3.459E-06, 2.304E-08, &
9.447E-01, 9.440E-01, 9.431E-01, 9.422E-01, 9.412E-01, 9.403E-01, 9.395E-01, 9.387E-01, 9.380E-01, 9.375E-01, &
9.371E-01, 9.367E-01, 9.361E-01, 9.352E-01, 9.339E-01, 9.321E-01, 9.297E-01, 9.267E-01, 9.230E-01, 9.185E-01, &
9.129E-01, 9.053E-01, 8.949E-01, 8.805E-01, 8.606E-01, 8.331E-01, 7.960E-01, 7.483E-01, 6.900E-01, 6.228E-01, &
5.495E-01, 4.761E-01, 4.112E-01, 3.602E-01, 3.177E-01, 2.719E-01, 2.161E-01, 1.520E-01, 8.906E-02, 4.021E-02, &
1.319E-02, 3.286E-03, 6.688E-04, 8.291E-05, 3.458E-06, 2.303E-08, &
9.427E-01, 9.421E-01, 9.412E-01, 9.403E-01, 9.393E-01, 9.384E-01, 9.376E-01, 9.368E-01, 9.361E-01, 9.356E-01, &
9.352E-01, 9.348E-01, 9.341E-01, 9.333E-01, 9.320E-01, 9.302E-01, 9.278E-01, 9.248E-01, 9.211E-01, 9.166E-01, &
9.110E-01, 9.035E-01, 8.930E-01, 8.787E-01, 8.587E-01, 8.313E-01, 7.943E-01, 7.466E-01, 6.884E-01, 6.213E-01, &
5.482E-01, 4.749E-01, 4.101E-01, 3.591E-01, 3.167E-01, 2.711E-01, 2.154E-01, 1.516E-01, 8.877E-02, 4.007E-02, &
1.315E-02, 3.278E-03, 6.680E-04, 8.286E-05, 3.456E-06, 2.302E-08, &
9.397E-01, 9.391E-01, 9.382E-01, 9.373E-01, 9.363E-01, 9.354E-01, 9.346E-01, 9.338E-01, 9.331E-01, 9.326E-01, &
9.322E-01, 9.318E-01, 9.312E-01, 9.303E-01, 9.290E-01, 9.273E-01, 9.249E-01, 9.218E-01, 9.181E-01, 9.137E-01, &
9.081E-01, 9.006E-01, 8.902E-01, 8.758E-01, 8.559E-01, 8.285E-01, 7.915E-01, 7.440E-01, 6.859E-01, 6.189E-01, &
5.460E-01, 4.729E-01, 4.083E-01, 3.576E-01, 3.153E-01, 2.698E-01, 2.144E-01, 1.508E-01, 8.832E-02, 3.986E-02, &
1.308E-02, 3.266E-03, 6.668E-04, 8.279E-05, 3.454E-06, 2.300E-08, &
9.351E-01, 9.344E-01, 9.336E-01, 9.327E-01, 9.317E-01, 9.308E-01, 9.300E-01, 9.292E-01, 9.285E-01, 9.280E-01, &
9.276E-01, 9.272E-01, 9.266E-01, 9.257E-01, 9.244E-01, 9.227E-01, 9.203E-01, 9.173E-01, 9.136E-01, 9.092E-01, &
9.036E-01, 8.961E-01, 8.857E-01, 8.714E-01, 8.515E-01, 8.241E-01, 7.873E-01, 7.399E-01, 6.820E-01, 6.153E-01, &
5.427E-01, 4.700E-01, 4.057E-01, 3.551E-01, 3.131E-01, 2.679E-01, 2.128E-01, 1.497E-01, 8.763E-02, 3.955E-02, &
1.298E-02, 3.248E-03, 6.650E-04, 8.267E-05, 3.450E-06, 2.297E-08, &
9.280E-01, 9.273E-01, 9.265E-01, 9.256E-01, 9.246E-01, 9.237E-01, 9.229E-01, 9.221E-01, 9.215E-01, 9.210E-01, &
9.206E-01, 9.201E-01, 9.195E-01, 9.187E-01, 9.174E-01, 9.157E-01, 9.133E-01, 9.103E-01, 9.067E-01, 9.023E-01, &
8.967E-01, 8.892E-01, 8.789E-01, 8.646E-01, 8.448E-01, 8.175E-01, 7.809E-01, 7.337E-01, 6.761E-01, 6.098E-01, &
5.377E-01, 4.655E-01, 4.016E-01, 3.515E-01, 3.098E-01, 2.650E-01, 2.105E-01, 1.480E-01, 8.662E-02, 3.908E-02, &
1.284E-02, 3.222E-03, 6.624E-04, 8.249E-05, 3.443E-06, 2.293E-08, &
9.174E-01, 9.167E-01, 9.159E-01, 9.150E-01, 9.140E-01, 9.131E-01, 9.123E-01, 9.115E-01, 9.109E-01, 9.104E-01, &
9.100E-01, 9.095E-01, 9.090E-01, 9.081E-01, 9.069E-01, 9.052E-01, 9.029E-01, 8.999E-01, 8.963E-01, 8.919E-01, &
8.864E-01, 8.789E-01, 8.686E-01, 8.544E-01, 8.348E-01, 8.077E-01, 7.712E-01, 7.244E-01, 6.673E-01, 6.016E-01, &
5.303E-01, 4.588E-01, 3.957E-01, 3.462E-01, 3.050E-01, 2.608E-01, 2.071E-01, 1.456E-01, 8.515E-02, 3.841E-02, &
1.264E-02, 3.185E-03, 6.586E-04, 8.222E-05, 3.433E-06, 2.287E-08, &
9.017E-01, 9.010E-01, 9.002E-01, 8.993E-01, 8.984E-01, 8.975E-01, 8.967E-01, 8.959E-01, 8.953E-01, 8.948E-01, &
8.944E-01, 8.940E-01, 8.934E-01, 8.925E-01, 8.914E-01, 8.897E-01, 8.874E-01, 8.845E-01, 8.810E-01, 8.767E-01, &
8.712E-01, 8.638E-01, 8.536E-01, 8.395E-01, 8.200E-01, 7.931E-01, 7.571E-01, 7.108E-01, 6.544E-01, 5.897E-01, &
5.194E-01, 4.492E-01, 3.872E-01, 3.385E-01, 2.981E-01, 2.548E-01, 2.022E-01, 1.421E-01, 8.308E-02, 3.748E-02, &
1.236E-02, 3.135E-03, 6.533E-04, 8.181E-05, 3.417E-06, 2.276E-08, &
8.791E-01, 8.784E-01, 8.776E-01, 8.767E-01, 8.758E-01, 8.749E-01, 8.741E-01, 8.734E-01, 8.728E-01, 8.723E-01, &
8.719E-01, 8.715E-01, 8.709E-01, 8.701E-01, 8.690E-01, 8.674E-01, 8.652E-01, 8.624E-01, 8.589E-01, 8.547E-01, &
8.493E-01, 8.420E-01, 8.319E-01, 8.180E-01, 7.988E-01, 7.723E-01, 7.368E-01, 6.912E-01, 6.360E-01, 5.726E-01, &
5.041E-01, 4.356E-01, 3.751E-01, 3.277E-01, 2.884E-01, 2.464E-01, 1.955E-01, 1.373E-01, 8.025E-02, 3.623E-02, &
1.199E-02, 3.071E-03, 6.460E-04, 8.117E-05, 3.392E-06, 2.260E-08, &
8.474E-01, 8.468E-01, 8.460E-01, 8.451E-01, 8.442E-01, 8.434E-01, 8.426E-01, 8.419E-01, 8.413E-01, 8.408E-01, &
8.405E-01, 8.401E-01, 8.395E-01, 8.388E-01, 8.377E-01, 8.362E-01, 8.341E-01, 8.314E-01, 8.281E-01, 8.240E-01, &
8.187E-01, 8.116E-01, 8.017E-01, 7.880E-01, 7.691E-01, 7.432E-01, 7.084E-01, 6.640E-01, 6.103E-01, 5.490E-01, &
4.828E-01, 4.169E-01, 3.587E-01, 3.130E-01, 2.753E-01, 2.350E-01, 1.864E-01, 1.308E-01, 7.648E-02, 3.459E-02, &
1.152E-02, 2.989E-03, 6.361E-04, 8.021E-05, 3.354E-06, 2.234E-08, &
8.047E-01, 8.041E-01, 8.033E-01, 8.025E-01, 8.017E-01, 8.009E-01, 8.001E-01, 7.994E-01, 7.989E-01, 7.984E-01, &
7.981E-01, 7.977E-01, 7.972E-01, 7.965E-01, 7.955E-01, 7.941E-01, 7.922E-01, 7.897E-01, 7.866E-01, 7.827E-01, &
7.776E-01, 7.706E-01, 7.610E-01, 7.477E-01, 7.293E-01, 7.041E-01, 6.704E-01, 6.276E-01, 5.760E-01, 5.174E-01, &
4.546E-01, 3.921E-01, 3.369E-01, 2.937E-01, 2.580E-01, 2.201E-01, 1.744E-01, 1.224E-01, 7.159E-02, 3.247E-02, &
1.091E-02, 2.883E-03, 6.221E-04, 7.876E-05, 3.294E-06, 2.195E-08, &
7.499E-01, 7.494E-01, 7.487E-01, 7.479E-01, 7.471E-01, 7.464E-01, 7.457E-01, 7.450E-01, 7.445E-01, 7.441E-01, &
7.438E-01, 7.434E-01, 7.430E-01, 7.424E-01, 7.415E-01, 7.403E-01, 7.385E-01, 7.363E-01, 7.334E-01, 7.297E-01, &
7.249E-01, 7.182E-01, 7.089E-01, 6.961E-01, 6.784E-01, 6.542E-01, 6.220E-01, 5.812E-01, 5.324E-01, 4.775E-01, &
4.189E-01, 3.608E-01, 3.097E-01, 2.695E-01, 2.364E-01, 2.015E-01, 1.595E-01, 1.119E-01, 6.551E-02, 2.984E-02, &
1.016E-02, 2.744E-03, 6.023E-04, 7.656E-05, 3.204E-06, 2.135E-08, &
6.841E-01, 6.836E-01, 6.829E-01, 6.822E-01, 6.815E-01, 6.809E-01, 6.802E-01, 6.797E-01, 6.792E-01, 6.788E-01, &
6.785E-01, 6.782E-01, 6.779E-01, 6.773E-01, 6.766E-01, 6.755E-01, 6.740E-01, 6.721E-01, 6.695E-01, 6.661E-01, &
6.615E-01, 6.552E-01, 6.464E-01, 6.342E-01, 6.174E-01, 5.945E-01, 5.641E-01, 5.259E-01, 4.807E-01, 4.301E-01, &
3.768E-01, 3.242E-01, 2.778E-01, 2.414E-01, 2.114E-01, 1.799E-01, 1.423E-01, 9.974E-02, 5.843E-02, 2.675E-02, &
9.247E-03, 2.567E-03, 5.743E-04, 7.333E-05, 3.070E-06, 2.046E-08, &
6.102E-01, 6.097E-01, 6.092E-01, 6.086E-01, 6.080E-01, 6.074E-01, 6.069E-01, 6.064E-01, 6.060E-01, 6.056E-01, &
6.054E-01, 6.051E-01, 6.048E-01, 6.044E-01, 6.037E-01, 6.029E-01, 6.016E-01, 6.000E-01, 5.977E-01, 5.947E-01, &
5.905E-01, 5.846E-01, 5.764E-01, 5.650E-01, 5.492E-01, 5.278E-01, 4.997E-01, 4.647E-01, 4.235E-01, 3.782E-01, &
3.308E-01, 2.844E-01, 2.436E-01, 2.113E-01, 1.847E-01, 1.569E-01, 1.239E-01, 8.674E-02, 5.083E-02, 2.339E-02, &
8.217E-03, 2.348E-03, 5.358E-04, 6.873E-05, 2.878E-06, 1.919E-08, &
5.325E-01, 5.321E-01, 5.316E-01, 5.311E-01, 5.306E-01, 5.302E-01, 5.297E-01, 5.293E-01, 5.290E-01, 5.287E-01, &
5.285E-01, 5.283E-01, 5.280E-01, 5.276E-01, 5.271E-01, 5.264E-01, 5.254E-01, 5.240E-01, 5.220E-01, 5.193E-01, &
5.155E-01, 5.102E-01, 5.027E-01, 4.922E-01, 4.777E-01, 4.582E-01, 4.327E-01, 4.012E-01, 3.647E-01, 3.250E-01, &
2.841E-01, 2.444E-01, 2.093E-01, 1.814E-01, 1.583E-01, 1.342E-01, 1.057E-01, 7.390E-02, 4.327E-02, 1.997E-02, &
7.119E-03, 2.088E-03, 4.853E-04, 6.251E-05, 2.618E-06, 1.747E-08, &
4.549E-01, 4.546E-01, 4.542E-01, 4.538E-01, 4.534E-01, 4.531E-01, 4.527E-01, 4.524E-01, 4.521E-01, 4.519E-01, &
4.517E-01, 4.515E-01, 4.513E-01, 4.510E-01, 4.505E-01, 4.499E-01, 4.491E-01, 4.479E-01, 4.462E-01, 4.438E-01, &
4.405E-01, 4.357E-01, 4.289E-01, 4.195E-01, 4.066E-01, 3.891E-01, 3.666E-01, 3.389E-01, 3.073E-01, 2.735E-01, &
2.391E-01, 2.060E-01, 1.768E-01, 1.533E-01, 1.335E-01, 1.129E-01, 8.879E-02, 6.189E-02, 3.616E-02, 1.670E-02, &
6.007E-03, 1.797E-03, 4.235E-04, 5.470E-05, 2.291E-06, 1.529E-08, &
3.799E-01, 3.797E-01, 3.794E-01, 3.790E-01, 3.788E-01, 3.785E-01, 3.782E-01, 3.780E-01, 3.777E-01, 3.775E-01, &
3.774E-01, 3.772E-01, 3.770E-01, 3.767E-01, 3.763E-01, 3.758E-01, 3.751E-01, 3.741E-01, 3.726E-01, 3.705E-01, &
3.676E-01, 3.634E-01, 3.575E-01, 3.493E-01, 3.380E-01, 3.229E-01, 3.035E-01, 2.801E-01, 2.535E-01, 2.254E-01, &
1.972E-01, 1.703E-01, 1.465E-01, 1.271E-01, 1.106E-01, 9.337E-02, 7.323E-02, 5.090E-02, 2.964E-02, 1.366E-02, &
4.930E-03, 1.490E-03, 3.541E-04, 4.579E-05, 1.917E-06, 1.279E-08, &
3.088E-01, 3.086E-01, 3.083E-01, 3.080E-01, 3.078E-01, 3.076E-01, 3.074E-01, 3.072E-01, 3.070E-01, 3.068E-01, &
3.066E-01, 3.065E-01, 3.063E-01, 3.060E-01, 3.057E-01, 3.053E-01, 3.046E-01, 3.037E-01, 3.025E-01, 3.007E-01, &
2.982E-01, 2.946E-01, 2.896E-01, 2.827E-01, 2.732E-01, 2.607E-01, 2.448E-01, 2.256E-01, 2.041E-01, 1.814E-01, &
1.587E-01, 1.373E-01, 1.183E-01, 1.027E-01, 8.927E-02, 7.525E-02, 5.889E-02, 4.083E-02, 2.371E-02, 1.090E-02, &
3.933E-03, 1.193E-03, 2.842E-04, 3.672E-05, 1.534E-06, 1.023E-08, &
2.429E-01, 2.427E-01, 2.425E-01, 2.422E-01, 2.420E-01, 2.418E-01, 2.417E-01, 2.415E-01, 2.413E-01, 2.411E-01, &
2.410E-01, 2.408E-01, 2.407E-01, 2.404E-01, 2.402E-01, 2.398E-01, 2.392E-01, 2.385E-01, 2.374E-01, 2.360E-01, &
2.339E-01, 2.309E-01, 2.268E-01, 2.212E-01, 2.137E-01, 2.039E-01, 1.916E-01, 1.768E-01, 1.601E-01, 1.424E-01, &
1.246E-01, 1.077E-01, 9.280E-02, 8.052E-02, 6.992E-02, 5.885E-02, 4.598E-02, 3.182E-02, 1.845E-02, 8.479E-03, &
3.063E-03, 9.306E-04, 2.213E-04, 2.849E-05, 1.188E-06, 7.905E-09, &
1.845E-01, 1.843E-01, 1.841E-01, 1.839E-01, 1.837E-01, 1.836E-01, 1.834E-01, 1.833E-01, 1.831E-01, 1.830E-01, &
1.828E-01, 1.827E-01, 1.825E-01, 1.823E-01, 1.821E-01, 1.818E-01, 1.813E-01, 1.807E-01, 1.798E-01, 1.786E-01, &
1.769E-01, 1.746E-01, 1.714E-01, 1.671E-01, 1.615E-01, 1.543E-01, 1.454E-01, 1.347E-01, 1.225E-01, 1.093E-01, &
9.571E-02, 8.266E-02, 7.112E-02, 6.162E-02, 5.343E-02, 4.490E-02, 3.503E-02, 2.423E-02, 1.406E-02, 6.481E-03, &
2.355E-03, 7.179E-04, 1.702E-04, 2.181E-05, 9.057E-07, 6.012E-09, &
1.365E-01, 1.363E-01, 1.361E-01, 1.359E-01, 1.358E-01, 1.356E-01, 1.355E-01, 1.353E-01, 1.352E-01, 1.351E-01, &
1.350E-01, 1.348E-01, 1.347E-01, 1.345E-01, 1.343E-01, 1.340E-01, 1.336E-01, 1.331E-01, 1.324E-01, 1.314E-01, &
1.301E-01, 1.283E-01, 1.259E-01, 1.228E-01, 1.188E-01, 1.138E-01, 1.078E-01, 1.005E-01, 9.212E-02, 8.266E-02, &
7.258E-02, 6.269E-02, 5.388E-02, 4.661E-02, 4.034E-02, 3.384E-02, 2.636E-02, 1.822E-02, 1.059E-02, 4.916E-03, &
1.805E-03, 5.537E-04, 1.309E-04, 1.669E-05, 6.904E-07, 4.569E-09, &
1.001E-01, 9.996E-02, 9.978E-02, 9.961E-02, 9.945E-02, 9.931E-02, 9.918E-02, 9.905E-02, 9.893E-02, 9.882E-02, &
9.871E-02, 9.860E-02, 9.847E-02, 9.832E-02, 9.813E-02, 9.789E-02, 9.757E-02, 9.713E-02, 9.655E-02, 9.576E-02, &
9.471E-02, 9.333E-02, 9.156E-02, 8.931E-02, 8.655E-02, 8.320E-02, 7.919E-02, 7.443E-02, 6.877E-02, 6.216E-02, &
5.484E-02, 4.746E-02, 4.079E-02, 3.525E-02, 3.044E-02, 2.547E-02, 1.980E-02, 1.367E-02, 7.963E-03, 3.720E-03, &
1.380E-03, 4.260E-04, 1.004E-04, 1.275E-05, 5.257E-07, 3.469E-09, &
7.385E-02, 7.370E-02, 7.354E-02, 7.338E-02, 7.324E-02, 7.312E-02, 7.300E-02, 7.289E-02, 7.278E-02, 7.268E-02, &
7.259E-02, 7.248E-02, 7.237E-02, 7.224E-02, 7.208E-02, 7.187E-02, 7.160E-02, 7.124E-02, 7.076E-02, 7.013E-02, &
6.930E-02, 6.824E-02, 6.690E-02, 6.528E-02, 6.334E-02, 6.106E-02, 5.838E-02, 5.519E-02, 5.134E-02, 4.670E-02, &
4.142E-02, 3.596E-02, 3.096E-02, 2.673E-02, 2.302E-02, 1.920E-02, 1.487E-02, 1.025E-02, 5.972E-03, 2.804E-03, &
1.049E-03, 3.253E-04, 7.644E-05, 9.662E-06, 3.975E-07, 2.618E-09, &
5.482E-02, 5.468E-02, 5.454E-02, 5.440E-02, 5.428E-02, 5.417E-02, 5.407E-02, 5.397E-02, 5.388E-02, 5.380E-02, &
5.371E-02, 5.363E-02, 5.353E-02, 5.342E-02, 5.328E-02, 5.311E-02, 5.289E-02, 5.259E-02, 5.221E-02, 5.171E-02, &
5.106E-02, 5.025E-02, 4.926E-02, 4.808E-02, 4.671E-02, 4.513E-02, 4.328E-02, 4.107E-02, 3.836E-02, 3.505E-02, &
3.123E-02, 2.722E-02, 2.348E-02, 2.027E-02, 1.740E-02, 1.444E-02, 1.114E-02, 7.648E-03, 4.455E-03, 2.101E-03, &
7.928E-04, 2.465E-04, 5.768E-05, 7.265E-06, 2.983E-07, 1.963E-09, &
4.078E-02, 4.067E-02, 4.055E-02, 4.043E-02, 4.033E-02, 4.024E-02, 4.016E-02, 4.008E-02, 4.001E-02, 3.994E-02, &
3.987E-02, 3.980E-02, 3.972E-02, 3.963E-02, 3.952E-02, 3.938E-02, 3.921E-02, 3.898E-02, 3.868E-02, 3.830E-02, &
3.782E-02, 3.722E-02, 3.650E-02, 3.566E-02, 3.470E-02, 3.359E-02, 3.228E-02, 3.070E-02, 2.873E-02, 2.633E-02, &
2.356E-02, 2.062E-02, 1.784E-02, 1.538E-02, 1.316E-02, 1.086E-02, 8.325E-03, 5.689E-03, 3.309E-03, 1.568E-03, &
5.963E-04, 1.857E-04, 4.324E-05, 5.424E-06, 2.223E-07, 1.462E-09, &
3.045E-02, 3.037E-02, 3.028E-02, 3.019E-02, 3.011E-02, 3.004E-02, 2.998E-02, 2.992E-02, 2.987E-02, 2.982E-02, &
2.976E-02, 2.971E-02, 2.965E-02, 2.958E-02, 2.950E-02, 2.940E-02, 2.927E-02, 2.911E-02, 2.889E-02, 2.862E-02, &
2.828E-02, 2.786E-02, 2.736E-02, 2.677E-02, 2.609E-02, 2.530E-02, 2.435E-02, 2.318E-02, 2.172E-02, 1.996E-02, &
1.792E-02, 1.576E-02, 1.367E-02, 1.179E-02, 1.004E-02, 8.228E-03, 6.262E-03, 4.251E-03, 2.465E-03, 1.171E-03, &
4.478E-04, 1.393E-04, 3.218E-05, 4.014E-06, 1.642E-07, 1.079E-09, &
2.288E-02, 2.282E-02, 2.276E-02, 2.271E-02, 2.265E-02, 2.261E-02, 2.257E-02, 2.253E-02, 2.249E-02, 2.246E-02, &
2.242E-02, 2.239E-02, 2.235E-02, 2.230E-02, 2.225E-02, 2.218E-02, 2.209E-02, 2.198E-02, 2.183E-02, 2.165E-02, &
2.142E-02, 2.114E-02, 2.080E-02, 2.040E-02, 1.992E-02, 1.934E-02, 1.863E-02, 1.775E-02, 1.666E-02, 1.534E-02, &
1.383E-02, 1.223E-02, 1.065E-02, 9.183E-03, 7.783E-03, 6.332E-03, 4.776E-03, 3.214E-03, 1.852E-03, 8.782E-04, &
3.357E-04, 1.036E-04, 2.364E-05, 2.928E-06, 1.194E-07, 7.851E-10, &
1.728E-02, 1.725E-02, 1.722E-02, 1.718E-02, 1.715E-02, 1.713E-02, 1.710E-02, 1.708E-02, 1.706E-02, 1.704E-02, &
1.702E-02, 1.700E-02, 1.698E-02, 1.695E-02, 1.691E-02, 1.687E-02, 1.681E-02, 1.674E-02, 1.665E-02, 1.654E-02, &
1.639E-02, 1.621E-02, 1.599E-02, 1.572E-02, 1.538E-02, 1.495E-02, 1.442E-02, 1.376E-02, 1.294E-02, 1.196E-02, &
1.084E-02, 9.634E-03, 8.429E-03, 7.270E-03, 6.132E-03, 4.948E-03, 3.694E-03, 2.459E-03, 1.403E-03, 6.599E-04, &
2.501E-04, 7.602E-05, 1.704E-05, 2.089E-06, 8.496E-08, 5.590E-10, &
1.310E-02, 1.309E-02, 1.307E-02, 1.305E-02, 1.304E-02, 1.303E-02, 1.301E-02, 1.300E-02, 1.299E-02, 1.298E-02, &
1.297E-02, 1.296E-02, 1.294E-02, 1.293E-02, 1.291E-02, 1.288E-02, 1.284E-02, 1.280E-02, 1.274E-02, 1.267E-02, &
1.258E-02, 1.247E-02, 1.232E-02, 1.214E-02, 1.191E-02, 1.161E-02, 1.123E-02, 1.074E-02, 1.014E-02, 9.413E-03, &
8.579E-03, 7.674E-03, 6.744E-03, 5.819E-03, 4.886E-03, 3.909E-03, 2.887E-03, 1.898E-03, 1.068E-03, 4.960E-04, &
1.849E-04, 5.482E-05, 1.198E-05, 1.451E-06, 5.880E-08, 3.878E-10, &
9.942E-03, 9.935E-03, 9.926E-03, 9.918E-03, 9.909E-03, 9.900E-03, 9.892E-03, 9.885E-03, 9.878E-03, 9.872E-03, &
9.865E-03, 9.858E-03, 9.849E-03, 9.839E-03, 9.825E-03, 9.807E-03, 9.783E-03, 9.752E-03, 9.715E-03, 9.670E-03, &
9.613E-03, 9.541E-03, 9.448E-03, 9.328E-03, 9.174E-03, 8.973E-03, 8.710E-03, 8.370E-03, 7.941E-03, 7.417E-03, &
6.807E-03, 6.130E-03, 5.413E-03, 4.674E-03, 3.908E-03, 3.102E-03, 2.269E-03, 1.474E-03, 8.181E-04, 3.735E-04, &
1.360E-04, 3.901E-05, 8.257E-06, 9.827E-07, 3.967E-08, 2.627E-10, &
7.538E-03, 7.533E-03, 7.527E-03, 7.521E-03, 7.515E-03, 7.509E-03, 7.503E-03, 7.498E-03, 7.493E-03, 7.488E-03, &
7.483E-03, 7.477E-03, 7.470E-03, 7.462E-03, 7.452E-03, 7.438E-03, 7.420E-03, 7.397E-03, 7.370E-03, 7.338E-03, &
7.298E-03, 7.248E-03, 7.185E-03, 7.104E-03, 7.002E-03, 6.867E-03, 6.691E-03, 6.460E-03, 6.165E-03, 5.799E-03, &
5.362E-03, 4.864E-03, 4.318E-03, 3.733E-03, 3.111E-03, 2.455E-03, 1.781E-03, 1.147E-03, 6.295E-04, 2.829E-04, &
1.005E-04, 2.778E-05, 5.663E-06, 6.604E-07, 2.652E-08, 1.766E-10, &
5.710E-03, 5.705E-03, 5.701E-03, 5.696E-03, 5.691E-03, 5.687E-03, 5.682E-03, 5.678E-03, 5.674E-03, 5.670E-03, &
5.666E-03, 5.661E-03, 5.656E-03, 5.649E-03, 5.640E-03, 5.629E-03, 5.615E-03, 5.597E-03, 5.576E-03, 5.550E-03, &
5.519E-03, 5.481E-03, 5.433E-03, 5.375E-03, 5.302E-03, 5.208E-03, 5.088E-03, 4.931E-03, 4.730E-03, 4.479E-03, &
4.174E-03, 3.815E-03, 3.405E-03, 2.949E-03, 2.454E-03, 1.928E-03, 1.393E-03, 8.926E-04, 4.871E-04, 2.166E-04, &
7.556E-05, 2.023E-05, 3.979E-06, 4.546E-07, 1.815E-08, 1.214E-10, &
4.321E-03, 4.318E-03, 4.315E-03, 4.311E-03, 4.308E-03, 4.304E-03, 4.301E-03, 4.298E-03, 4.294E-03, 4.291E-03, &
4.288E-03, 4.283E-03, 4.279E-03, 4.273E-03, 4.266E-03, 4.257E-03, 4.246E-03, 4.232E-03, 4.215E-03, 4.194E-03, &
4.169E-03, 4.138E-03, 4.101E-03, 4.056E-03, 4.001E-03, 3.933E-03, 3.847E-03, 3.738E-03, 3.600E-03, 3.428E-03, &
3.216E-03, 2.960E-03, 2.657E-03, 2.307E-03, 1.919E-03, 1.504E-03, 1.084E-03, 6.945E-04, 3.788E-04, 1.679E-04, &
5.798E-05, 1.518E-05, 2.903E-06, 3.258E-07, 1.293E-08, 8.659E-11, &
3.276E-03, 3.273E-03, 3.270E-03, 3.268E-03, 3.265E-03, 3.263E-03, 3.260E-03, 3.258E-03, 3.255E-03, 3.252E-03, &
3.249E-03, 3.246E-03, 3.242E-03, 3.238E-03, 3.232E-03, 3.225E-03, 3.216E-03, 3.205E-03, 3.191E-03, 3.175E-03, &
3.155E-03, 3.130E-03, 3.101E-03, 3.066E-03, 3.024E-03, 2.974E-03, 2.912E-03, 2.835E-03, 2.740E-03, 2.620E-03, &
2.472E-03, 2.289E-03, 2.065E-03, 1.799E-03, 1.496E-03, 1.171E-03, 8.431E-04, 5.402E-04, 2.951E-04, 1.309E-04, &
4.498E-05, 1.159E-05, 2.163E-06, 2.387E-07, 9.399E-09, 6.285E-11, &
2.487E-03, 2.485E-03, 2.483E-03, 2.481E-03, 2.479E-03, 2.477E-03, 2.475E-03, 2.473E-03, 2.471E-03, 2.469E-03, &
2.467E-03, 2.464E-03, 2.461E-03, 2.457E-03, 2.453E-03, 2.448E-03, 2.441E-03, 2.432E-03, 2.422E-03, 2.409E-03, &
2.393E-03, 2.375E-03, 2.352E-03, 2.326E-03, 2.295E-03, 2.258E-03, 2.214E-03, 2.160E-03, 2.093E-03, 2.009E-03, &
1.904E-03, 1.771E-03, 1.604E-03, 1.402E-03, 1.168E-03, 9.129E-04, 6.560E-04, 4.199E-04, 2.295E-04, 1.019E-04, &
3.490E-05, 8.876E-06, 1.621E-06, 1.758E-07, 6.858E-09, 4.569E-11, &
1.885E-03, 1.884E-03, 1.882E-03, 1.881E-03, 1.879E-03, 1.878E-03, 1.877E-03, 1.875E-03, 1.874E-03, 1.872E-03, &
1.870E-03, 1.868E-03, 1.866E-03, 1.864E-03, 1.860E-03, 1.856E-03, 1.852E-03, 1.845E-03, 1.838E-03, 1.828E-03, &
1.817E-03, 1.803E-03, 1.786E-03, 1.767E-03, 1.745E-03, 1.719E-03, 1.689E-03, 1.651E-03, 1.605E-03, 1.546E-03, &
1.470E-03, 1.372E-03, 1.248E-03, 1.094E-03, 9.135E-04, 7.142E-04, 5.121E-04, 3.268E-04, 1.782E-04, 7.900E-05, &
2.698E-05, 6.797E-06, 1.220E-06, 1.305E-07, 5.043E-09, 3.347E-11, &
1.428E-03, 1.427E-03, 1.426E-03, 1.425E-03, 1.424E-03, 1.423E-03, 1.422E-03, 1.421E-03, 1.420E-03, 1.419E-03, &
1.418E-03, 1.417E-03, 1.415E-03, 1.413E-03, 1.411E-03, 1.408E-03, 1.405E-03, 1.401E-03, 1.396E-03, 1.389E-03, &
1.381E-03, 1.371E-03, 1.360E-03, 1.347E-03, 1.332E-03, 1.314E-03, 1.293E-03, 1.268E-03, 1.235E-03, 1.194E-03, &
1.139E-03, 1.068E-03, 9.754E-04, 8.591E-04, 7.196E-04, 5.636E-04, 4.036E-04, 2.564E-04, 1.389E-04, 6.114E-05, &
2.077E-05, 5.199E-06, 9.234E-07, 9.765E-08, 3.746E-09, 2.483E-11, &
1.083E-03, 1.083E-03, 1.082E-03, 1.081E-03, 1.081E-03, 1.080E-03, 1.079E-03, 1.079E-03, 1.078E-03, 1.077E-03, &
1.077E-03, 1.076E-03, 1.075E-03, 1.074E-03, 1.073E-03, 1.071E-03, 1.069E-03, 1.066E-03, 1.063E-03, 1.058E-03, &
1.053E-03, 1.046E-03, 1.039E-03, 1.030E-03, 1.020E-03, 1.008E-03, 9.939E-04, 9.764E-04, 9.540E-04, 9.246E-04, &
8.856E-04, 8.337E-04, 7.653E-04, 6.777E-04, 5.707E-04, 4.486E-04, 3.213E-04, 2.031E-04, 1.089E-04, 4.734E-05, &
1.593E-05, 3.965E-06, 6.976E-07, 7.299E-08, 2.783E-09, 1.849E-11, &
8.178E-04, 8.173E-04, 8.169E-04, 8.164E-04, 8.160E-04, 8.156E-04, 8.152E-04, 8.148E-04, 8.144E-04, 8.140E-04, &
8.136E-04, 8.132E-04, 8.127E-04, 8.120E-04, 8.113E-04, 8.103E-04, 8.090E-04, 8.074E-04, 8.052E-04, 8.024E-04, &
7.989E-04, 7.946E-04, 7.895E-04, 7.835E-04, 7.766E-04, 7.684E-04, 7.587E-04, 7.466E-04, 7.310E-04, 7.104E-04, &
6.827E-04, 6.453E-04, 5.954E-04, 5.303E-04, 4.494E-04, 3.552E-04, 2.550E-04, 1.607E-04, 8.531E-05, 3.658E-05, &
1.219E-05, 3.013E-06, 5.235E-07, 5.398E-08, 2.045E-09, 1.367E-11, &
6.061E-04, 6.058E-04, 6.055E-04, 6.052E-04, 6.050E-04, 6.047E-04, 6.045E-04, 6.042E-04, 6.040E-04, 6.038E-04, &
6.035E-04, 6.033E-04, 6.030E-04, 6.026E-04, 6.022E-04, 6.016E-04, 6.009E-04, 5.999E-04, 5.985E-04, 5.967E-04, &
5.943E-04, 5.915E-04, 5.880E-04, 5.838E-04, 5.790E-04, 5.733E-04, 5.664E-04, 5.580E-04, 5.471E-04, 5.328E-04, &
5.135E-04, 4.871E-04, 4.515E-04, 4.045E-04, 3.450E-04, 2.745E-04, 1.980E-04, 1.248E-04, 6.591E-05, 2.800E-05, &
9.258E-06, 2.275E-06, 3.887E-07, 3.930E-08, 1.479E-09, 9.984E-12, &
4.350E-04, 4.349E-04, 4.347E-04, 4.345E-04, 4.343E-04, 4.342E-04, 4.340E-04, 4.339E-04, 4.337E-04, 4.336E-04, &
4.334E-04, 4.333E-04, 4.331E-04, 4.329E-04, 4.327E-04, 4.323E-04, 4.319E-04, 4.312E-04, 4.304E-04, 4.292E-04, &
4.277E-04, 4.258E-04, 4.235E-04, 4.206E-04, 4.173E-04, 4.133E-04, 4.085E-04, 4.027E-04, 3.952E-04, 3.854E-04, &
3.722E-04, 3.541E-04, 3.296E-04, 2.968E-04, 2.549E-04, 2.043E-04, 1.484E-04, 9.402E-05, 4.972E-05, 2.108E-05, &
6.958E-06, 1.700E-06, 2.845E-07, 2.808E-08, 1.049E-09, 7.170E-12 /), (/ 46,46 /))
contains
end module shielding