[1310] | 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 | |
---|