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