source: trunk/LMDZ.MARS/libf/phymars/nlte_paramdef.h @ 1009

Last change on this file since 1009 was 757, checked in by emillour, 12 years ago

Mars GCM:

  • Improvement of the NLTE 15um scheme (for running with nltemodel = 2); now MUCH faster than previously (by a factor 5 or so):
  • Improvements included to the parameterization:
    • Cool-to-space calculation included above P(atm)=1e-10, with a soft merging to the full result (without the CTS approximation) below that level
    • exhaustive cleaning of the code, including FTNCHK and reordering of loops, subroutines and internal calls
    • simplification of the precomputed tables of CO2 bands' atmospheric transmittances
    • the two internal grids (the one used in the CTS region and the one below) have been , in order to reduce the CPU time consumption
    • reading of the spectroscopic histograms is made only once, at the beginning of the GCM, to avoid repetitive readings of ASCII files
    • F90 matrix operations (matmul,...) included.
  • Changes in routines:
    • removed nlte_leedat.F
    • updated nlte_calc.F, nlte_tcool.F, nlte_aux.F
    • updated nlte_commons.h, nlte_paramdef.h
    • added nlte_setup.F
  • Important: The input files (in the NLTEDAT directory) read as input by these routines have changed. the NLTEDAT directory should now on contain only the following files:

deltanu26.dat enelow27.dat hid26-3.dat parametp_Tstar_IAA1204.dat
deltanu27.dat enelow28.dat hid26-4.dat parametp_VC_IAA1204.dat
deltanu28.dat enelow36.dat hid27-1.dat
deltanu36.dat hid26-1.dat hid28-1.dat
enelow26.dat hid26-2.dat hid36-1.dat

FGG+MALV
File size: 2.6 KB
Line 
1c****************************************************************************
2c
3c       Merging of different parameters definitions for new NLTE 15um param
4c
5c       jul 2012    fgg+malv
6c****************************************************************************
7c *** 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
44c *** 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 = 3.e-10 , Pbottom_atm = 2.e-5 ) 
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
59c *** Old bloque_dlvr11.f ***
60
61        real nu(nisot,8)
62c 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
Note: See TracBrowser for help on using the repository browser.