source: trunk/LMDZ.MARS/libf/phymars/nlte_rates.h @ 482

Last change on this file since 482 was 414, checked in by aslmd, 13 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: 3.4 KB
Line 
1        common/rates_vt/ k7a(4),k7b(4), k7ap(4),k7bp(4), 
2     @          k3aa(4),k3ab(4),k3ac(4), k3aap(4),k3abp(4),k3acp(4),
3     @          k3ba(4),k3bb(4),k3bc(4), k3bap(4),k3bbp(4),k3bcp(4),
4     @          k19aa(4),k19ab(4),k19ac(4), k19aap(4),k19abp(4),k19acp(4),
5     @          k19ba(4),k19bb(4),k19bc(4), k19bap(4),k19bbp(4),k19bcp(4),
6     @          k19ca(4),k19cb(4),k19cc(4), k19cap(4),k19cbp(4),k19ccp(4),
7     @          k20a(4),k20b(4),k20c(4), k20ap(4),k20bp(4),k20cp(4),
8     @          k27a,k27b,k27c, k27ap,k27bp,k27cp
9
10        real*8 k7a,k7b, k7ap,k7bp
11        real*8 k3aa,k3ab,k3ac, k3aap,k3abp,k3acp
12        real*8 k3ba,k3bb,k3bc, k3bap,k3bbp,k3bcp
13        real*8 k19aa,k19ab,k19ac, k19aap,k19abp,k19acp
14        real*8 k19ba,k19bb,k19bc, k19bap,k19bbp,k19bcp
15        real*8 k19ca,k19cb,k19cc, k19cap,k19cbp,k19ccp
16        real*8 k20a,k20b,k20c, k20ap,k20bp,k20cp
17        real*8 k27a,k27b,k27c, k27ap,k27bp,k27cp
18
19        common/rates_vv/ k1(4),k1p(4), 
20     @          k2a,k2b, k2x,k2y,k2z, k2xp,k2yp,k2zp,
21     @          k6,k6p, k6a(2:4),k6b(2:4), k6ap(2:4),k6bp(2:4),
22     @          k21a,k21ap, k21a1(2:4),k21a2(2:4), k21a1p(2:4),k21a2p(2:4),
23     @          k21b(4),k21c(4), k21bp(4),k21cp(4),
24     @          k31,k32,
25     @          k33a1,k33a2,k33b1,k33b2,k33c, 
26     @          k33a1p(2:4),k33a2p(2:4),k33b1p(2:4),k33b2p(2:4),k33cp(2:4),
27     @          k28a,k28b,k28c, k28ap,k28bp,k28cp,
28     @          k26a,k26b,k26c,k26d, k26ap(4),k26bp(4),k26cp(4),k26dp(4), 
29     @          k41p_taylor, k41p_shved, k41p_starr_hannock, 
30     @          k41_1,k41p_1, k41_2,k41p_2, k42,k42p
31
32        real*8 k1,k1p
33        real*8 k2a,k2b, k2x,k2y,k2z, k2xp,k2yp,k2zp
34        real*8 k6,k6p, k6a,k6b, k6ap,k6bp
35        real*8 k21a,k21ap, k21a1,k21a2, k21a1p,k21a2p
36        real*8 k21b,k21c, k21bp,k21cp
37        real*8 k31,k32
38        real*8 k33a1,k33a2,k33b1,k33b2,k33c
39        real*8 k33a1p,k33a2p,k33b1p,k33b2p,k33cp
40        real*8 k28a,k28b,k28c, k28ap,k28bp,k28cp
41        real*8 k26a,k26b,k26c,k26d, k26ap,k26bp,k26cp,k26dp
42        real*8 k41p_taylor, k41p_shved, k41p_starr_hannock
43        real*8 k41_1,k41p_1, k41_2,k41p_2, k42,k42p
44
45
46        common/rates_k26isot/ k26a21,k26c21,k26d21, 
47     @     k26a22,k26c22,k26d22, k26a23,k26c23,k26d23, 
48     @     k26a24,k26c24,k26d24, 
49     @     k26a32,k26c32,k26d32, k26a33,k26c33,k26d33, 
50     @     k26a31,k26c31,k26d31, 
51     @     k26a34,k26c34,k26d34, k26a42,k26c42,k26d42, 
52     @     k26a41,k26c41,k26d41, 
53     @     k26a43,k26c43,k26d43, k26a44,k26c44,k26d44
54
55        real*8 k26a21,k26c21,k26d21, 
56     @     k26a22,k26c22,k26d22, k26a23,k26c23,k26d23, 
57     @     k26a24,k26c24,k26d24, 
58     @     k26a32,k26c32,k26d32, k26a33,k26c33,k26d33, 
59     @     k26a31,k26c31,k26d31, 
60     @     k26a34,k26c34,k26d34, k26a42,k26c42,k26d42, 
61     @     k26a41,k26c41,k26d41, 
62     @     k26a43,k26c43,k26d43, k26a44,k26c44,k26d44
63
64
65        common/rates_last/ k23k21c, k24k21c, k34k21c, 
66     @          k23k21cp, k24k21cp, k34k21cp, k43,k43p, k_vthcl
67
68        real*8 k23k21c,k24k21c,k34k21c, k23k21cp,k24k21cp,k34k21cp
69        real*8 k43,k43p, k_vthcl
70
71        common/rates_V09/ k41_3,k41p_3, k41_4,k41p_4, k41iso_1,k41iso_1p, 
72     @                    k41iso_2,k41iso_2p, k41iso_3,k41iso_3p, 
73     @                    k42b, k42c, k42bp, k42cp, k43iso,k43isop, 
74     @                    k44a,k44b,k44c,k44d, k44ap,k44bp,k44cp,k44dp, 
75     @                    k42iso,k42isop, k42isob,k42isobp
76        real*8  k41_3,k41p_3, k41_4,k41p_4, k41iso_1,k41iso_1p
77        real*8  k41iso_2,k41iso_2p, k41iso_3,k41iso_3p
78        real*8  k42b, k42c, k42bp, k42cp, k43iso,k43isop
79        real*8  k44a,k44b,k44c,k44d, k44ap,k44bp,k44cp,k44dp
80        real*8  k42iso,k42isop, k42isob,k42isobp
Note: See TracBrowser for help on using the repository browser.