1 | c**************************************************************************** |
---|
2 | c |
---|
3 | c Merging of different parameters definitions for new NLTE 15um param |
---|
4 | c |
---|
5 | c jul 2012 fgg+malv |
---|
6 | c**************************************************************************** |
---|
7 | c *** Old mz1d.par *** |
---|
8 | ! Grids parameters : |
---|
9 | |
---|
10 | integer nztabul ! # points in tabulation of Tesc & VC (ISO) |
---|
11 | parameter ( nztabul=79 ) |
---|
12 | |
---|
13 | ! NLTE parameters: |
---|
14 | |
---|
15 | integer nltot ! incluye el actual # alt in NLTE module |
---|
16 | parameter ( nltot=20 ) ! y el # alturas del Tstar110 |
---|
17 | |
---|
18 | integer nl ! actual # alt in NLTE module & C.Matrix |
---|
19 | parameter ( nl=12 ) |
---|
20 | integer nl2 ! = nl-2, needed for matrix inversion (mmh2) |
---|
21 | parameter ( nl2=nl-2 ) |
---|
22 | |
---|
23 | integer nzy |
---|
24 | parameter ( nzy = (nl-1)*4 + 1 ) ! Fine grid for C.Matrix |
---|
25 | |
---|
26 | integer nl_cts ! actual # alt para Tstar110 |
---|
27 | parameter ( nl_cts = 2 + nltot-nl ) |
---|
28 | integer nzy_cts ! fine grid for transmit calculation |
---|
29 | parameter ( nzy_cts = (nl_cts-1)*4 + 1 ) |
---|
30 | |
---|
31 | |
---|
32 | ! Other NLTE parameters: |
---|
33 | integer nisot ! number of isotopes considered |
---|
34 | integer nb ! number of bands included |
---|
35 | parameter ( nisot=4, nb=41 ) |
---|
36 | |
---|
37 | integer nhist ! # of temps in histogr. |
---|
38 | parameter ( nhist = 36 ) ! (get it from histograms!) |
---|
39 | |
---|
40 | integer nbox_max |
---|
41 | parameter ( nbox_max = 4 ) ! max.# boxes in histogram |
---|
42 | |
---|
43 | |
---|
44 | c *** Old tcr_15um.h *** |
---|
45 | |
---|
46 | integer itt_cza ! Selection of NLTE scheme |
---|
47 | parameter ( itt_cza = 13 ) |
---|
48 | |
---|
49 | real Ptop_atm, Pbottom_atm ! Upper and lower limits of |
---|
50 | ! NLTE model |
---|
51 | parameter ( Ptop_atm = 6.e-10 , Pbottom_atm = 8.e-6 ) |
---|
52 | |
---|
53 | |
---|
54 | real*8 rf19,rf20,rf21a,rf21b,rf21c,rf33bc |
---|
55 | parameter ( rf19 = 1.d0, rf20 = 1.d0, rf21a = 1.d0) |
---|
56 | parameter ( rf21b = 1.d0, rf21c = 1.d0, rf33bc = 1.d0 ) |
---|
57 | |
---|
58 | |
---|
59 | c *** Old bloque_dlvr11.f *** |
---|
60 | |
---|
61 | real nu(nisot,8) |
---|
62 | c data |
---|
63 | data nu(1,1),nu(1,2) /667.3801, 1335.1317/ |
---|
64 | data nu(2,1)/662.3734/ |
---|
65 | data nu(3,1)/648.4784/ |
---|
66 | data nu(4,1)/664.7289/ |
---|
67 | |
---|
68 | real nu12_0200,nu12_1000 |
---|
69 | parameter (nu12_0200 = 1285.4087) |
---|
70 | parameter (nu12_1000 = 1388.1847) |
---|
71 | |
---|
72 | integer indexisot(nisot) |
---|
73 | data indexisot/26,28,36,27/ |
---|
74 | |
---|
75 | ! ctes en el sistema cgs |
---|
76 | real*8 vlight, ee, hplanck, gamma |
---|
77 | parameter (vlight = 2.9979245e10) |
---|
78 | parameter (ee = 1.43876866) |
---|
79 | parameter (hplanck = 6.6260755e-27) |
---|
80 | parameter (gamma = 1.191043934e-5) |
---|
81 | |
---|
82 | |
---|
83 | ! datos de marte |
---|
84 | real imr(nisot) |
---|
85 | data imr / 0.987, 0.00408, 0.0112, 0.000742 / |
---|
86 | |
---|
87 | |
---|
88 | |
---|
89 | |
---|