source: trunk/LMDZ.MARS/libf/phymars/einsteincoefs.h @ 426

Last change on this file since 426 was 414, checked in by aslmd, 14 years ago

LMDZ.MARS : NEW NLTE MODEL FROM GRANADA AMIGOS

23/11/11 == FGG + MALV

New parameterization of the NLTE 15 micron cooling. The old parameterization is kept as an option, including or not variable atomic oxygen concentration. A new flag is introduced in callphys.def, nltemodel, to select which parameterization wants to be used (new one, old one with variable [O], or old one with fixed [O], see below). Includes many new subroutines and commons in phymars. Some existing routines are also modified:

-physiq.F. Call to the new subroutine NLTE_leedat in first call. Call to nltecool modified to allow for variable atomic oxygen. Depending on the value of nltemodel, the new subroutine NLTEdlvr09_TCOOL is called instead of nltecool.

-inifis.F. Reading of nltemodel is added.

-callkeys.h Declaration of nltemodel is added.

The following lines should be added to callphys.def (ideally after setting callnlte):

# NLTE 15um scheme to use.
# 0-> Old scheme, static oxygen
# 1-> Old scheme, dynamic oxygen
# 2-> New scheme
nltemodel = 2

A new directory, NLTEDAT, has to be included in datagcm.

Improvements into NLTE NIR heating parameterization to take into account variability with O/CO2 ratio and SZA. A new subroutine, NIR_leedat.F, and a new common, NIRdata.h, are included. A new flag, nircorr, is added in callphys.def, to include or not these corrections. The following files are modified:

-nirco2abs.F: nq and pq are added in the arguments. The corrections factors are interpolated to the GCM grid and included in the clculation. A new subroutine, interpnir, is added at the end of the file.

-physiq.F: Call to NIR_leedat added at first call. Modified call to nirco2abs

-inifis: Reading new flag nircorr.

-callkeys.h: Declaration of nircorr.

The following lines have to be added to callphys.def (ideally after callnirco2):

# NIR NLTE correction for variable SZA and O/CO2?
# matters only if callnirco2=T
# 0-> no correction
# 1-> include correction
nircorr=1

A new data file, NIRcorrection_feb2011.dat, has to be included in datagcm.

Small changes to the molecular diffusion scheme to fix the number of species considered, to avoid problems when compiling with more than 15 tracers (for example, when CH4 is included). Modified subroutines: aeronomars/moldiff.F and aeronomars/moldiffcoeff.F
File size: 5.1 KB
Line 
1c**********************************************************************
2c*
3c* Einstein coefficients
4c*
5c* fgg mar 2010 (adapted from MALV get_a subroutine)
6c*
7c**********************************************************************
8c                                               
9! isotope  626          for this isotope the values are from mlp   
10c
11      real*8 a1_103, a1_102, a1_101, a1_023, a1_022, a1_021
12c
13  parameter (a1_103       = 1260.d0)
14  parameter (a1_102       = 820.d0)
15  parameter (a1_101       = 410.d0)
16  parameter (a1_023       = 1220.d0)
17  parameter (a1_022       = 814.d0)
18  parameter (a1_021       = 406.d0)
19c
20c
21      real*8    a1_003, a1_002, a1_301, a1_221, a1_141, a1_061
22c
23  parameter (a1_003       = 1200.d0)
24  parameter (a1_002       = 780.d0)
25  parameter (a1_301       = 400.d0)
26  parameter (a1_221       = 400.d0)
27  parameter (a1_141       = 400.d0)
28  parameter (a1_061       = 400.d0)
29c
30c
31      real*8 a1_201, a1_121, a1_041, a1_220, a1_130, a1_050
32c
33  parameter (a1_201       = 424.d0)
34  parameter (a1_121       = 418.d0)
35  parameter (a1_041       = 418.d0)
36  parameter (a1_220       = 1.43d0)
37  parameter (a1_130       = 2.3d0)
38  parameter (a1_050       = 7.15d0)
39c
40c
41      real*8 a1_120, a1_55, a1_44, a1_300, a1_a43, a1_210, a1_200
42c
43  parameter (a1_120       = 2.08d0)
44  parameter (a1_55        = (5.67d0+a1_120)/2.d0)
45  parameter (a1_44        = 3.68d0)
46  parameter (a1_300       = 4.3d0)
47  parameter (a1_a43       = 2.44d0)
48  parameter (a1_210       = a1_a43+a1_44)
49  parameter (a1_200       = 4.53d0)
50c
51c
52      real*8 a1_140, a1_060, a1_0510, a1_040,  a1_0420, a1_1220
53c
54  parameter (a1_140       = 2.27d0)
55  parameter (a1_060       = 2.32d0)
56  parameter (a1_0510      = a1_0510_0420+a1_0510_0400)
57  parameter (a1_040       = 2.35d0)
58  parameter (a1_0420      = 3.3d0)
59  parameter (a1_1220      = 4.71d0)
60c
61c
62      real*8 a1_1111_1110,a1_1111_0110, a1_0311_0310,a1_0311_0110
63c
64  parameter (a1_1111_1110 = 367.0d0)
65  parameter (a1_1111_0110 = 17.0d0)
66  parameter (a1_0311_0310 = 370.8d0)
67  parameter (a1_0311_0110 = 9.65d0)
68c
69c
70      real*8 a1_0111_010, a1_0201_0200, a1_0201_000, a1_0221_0220
71c
72  parameter (a1_0111_010  = 388.4d0)
73  parameter (a1_0201_000  = 10.39d0)
74  parameter (a1_0201_0200 = 383.756d0)
75  parameter (a1_0221_0220 = 386.56d0)
76c
77c
78      real*8 a1_1001_1000, a1_1001_000, a1_001_000,a1_010_000
79c
80  parameter (a1_1001_1000 = 381.222d0)
81  parameter (a1_1001_000  = 16.68d0)
82  parameter (a1_001_000   = 389.34d0)
83  parameter (a1_010_000   = 1.3546d0)
84c
85c
86      real*8 a1_0421_1220, a1_0421_0420, a1_0421_000
87c
88  parameter (a1_0421_000  = 0.66d0)
89  parameter (a1_0421_0420 = 360.5d0)
90  parameter (a1_0421_1220 = 0.129d0)
91c
92c
93      real*8 a1_1221_0111, a1_1221_1220, a1_1221_0420
94c
95  parameter (a1_1221_0420 = 0.01227d0)
96  parameter (a1_1221_1220 = 355.2d0)
97  parameter (a1_1221_0111 = 0.0219d0)
98c
99c
100      real*8 a1_0510_000, a1_0510_0420, a1_0510_0400
101c
102  parameter (a1_0510_000  = 4.d-4)
103  parameter (a1_0510_0420 = 1.086d0)
104  parameter (a1_0510_0400 = 2.406d0)
105c
106c
107      real*8 a1_0510_1000, a1_0510_0220, a1_0510_0200
108  parameter (a1_0510_1000 = 3.44d-5)
109  parameter (a1_0510_0220 = 6.339d-4)
110  parameter (a1_0510_0200 = 1.5d-3)
111c
112c
113      real*8 a1_1310_000 , a1_1310_0420 , a1_1310_0400
114c
115  parameter (a1_1310_000  = 4.45d-4)
116  parameter (a1_1310_0420 = 0.7535d0)
117  parameter (a1_1310_0400 = 0.21d0)
118c
119c
120      real*8 a1_1310_1200 , a1_1310_1220 , a1_1310_2000
121c
122  parameter (a1_1310_1200 = 2.264d0)
123  parameter (a1_1310_1220 = 0.3609d0)
124  parameter (a1_1310_2000 = 0.074d0)
125c
126c
127      real*8 a1_1310_0111,a1_1310_1000,a1_1310_0220,a1_1310_0200
128c
129  parameter (a1_1310_0111 = 1.64d-4)
130  parameter (a1_1310_1000 = 1.433d-3)
131  parameter (a1_1310_0220 = 2.4d-4)
132  parameter (a1_1310_0200 = 9.719d-3)
133c
134c
135      real*8 a1_2110_000 , a1_2110_0420 , a1_2110_0400
136c
137  parameter (a1_2110_000  = 3.265d-4)
138  parameter (a1_2110_0420 = 6.48d-3)
139  parameter (a1_2110_0400 = 1.54d-3)
140c
141c
142      real*8 a1_2110_1200,a1_2110_1220,a1_2110_2000,a1_2110_0111
143c
144  parameter (a1_2110_1200 = 0.0738d0)
145  parameter (a1_2110_1220 = 1.957d0)
146  parameter (a1_2110_2000 = 3.289d0)
147  parameter (a1_2110_0111 = 2.68d-4)
148c
149c
150      real*8 a1_2110_1000 , a1_2110_0220 , a1_2110_0200
151c
152  parameter (a1_2110_1000 = 3.445d-2)
153  parameter (a1_2110_0220 = 1.976d-2)
154  parameter (a1_2110_0200 = 6.539d-3)
155c
156c
157      real*8 a1_0403_0402, a1_0403_0202
158c
159  parameter (a1_0403_0402 = 728.7d0)
160  parameter (a1_0403_0202 = 28.9d0)
161c
162      real*8 a1_1203_1202, a1_1203_1002, a1_1203_0202
163c
164c
165  parameter (a1_1203_1202 = 716.3d0)
166  parameter (a1_1203_1002 = 27.25d0)
167  parameter (a1_1203_0202 = 43.68d0)
168c
169c
170      real*8 a1_2003_2002, a1_2003_1002, a1_2003_0202
171c
172  parameter (a1_2003_2002 = 995.7d0)
173  parameter (a1_2003_1002 = 83.27d0)
174  parameter (a1_2003_0202 = 1.15d0)
175c
176c
177      real*8 a1_0402_0401, a1_0402_0201
178c
179  parameter (a1_0402_0401 = 700.0d0)
180  parameter (a1_0402_0201 = 28.97d0)
181c
182c
183      real*8 a1_1202_1201, a1_1202_1001, a1_1202_0201
184c
185  parameter (a1_1202_1201 = 687.0d0) 
186  parameter (a1_1202_1001 =  26.07d0)
187  parameter (a1_1202_0201 =  39.57d0)
188c
189c
190      real*8 a1_2002_2001, a1_2002_1001
191c
192  parameter (a1_2002_2001 = 686.1d0)
193  parameter (a1_2002_1001 = 54.67d0)
194
Note: See TracBrowser for help on using the repository browser.