Context Navigation

einsteincoefs.h @ 472

Last change on this file since 472 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

Rev	Line
[414]	1	c**********************************************************************
	2	c*
	3	c* Einstein coefficients
	4	c*
	5	c* fgg mar 2010 (adapted from MALV get_a subroutine)
	6	c*
	7	c**********************************************************************
	8	c
	9	! isotope 626 for this isotope the values are from mlp
	10	c
	11	real*8 a1_103, a1_102, a1_101, a1_023, a1_022, a1_021
	12	c
	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)
	19	c
	20	c
	21	real*8 a1_003, a1_002, a1_301, a1_221, a1_141, a1_061
	22	c
	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)
	29	c
	30	c
	31	real*8 a1_201, a1_121, a1_041, a1_220, a1_130, a1_050
	32	c
	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)
	39	c
	40	c
	41	real*8 a1_120, a1_55, a1_44, a1_300, a1_a43, a1_210, a1_200
	42	c
	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)
	50	c
	51	c
	52	real*8 a1_140, a1_060, a1_0510, a1_040, a1_0420, a1_1220
	53	c
	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)
	60	c
	61	c
	62	real*8 a1_1111_1110,a1_1111_0110, a1_0311_0310,a1_0311_0110
	63	c
	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)
	68	c
	69	c
	70	real*8 a1_0111_010, a1_0201_0200, a1_0201_000, a1_0221_0220
	71	c
	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)
	76	c
	77	c
	78	real*8 a1_1001_1000, a1_1001_000, a1_001_000,a1_010_000
	79	c
	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)
	84	c
	85	c
	86	real*8 a1_0421_1220, a1_0421_0420, a1_0421_000
	87	c
	88	parameter (a1_0421_000 = 0.66d0)
	89	parameter (a1_0421_0420 = 360.5d0)
	90	parameter (a1_0421_1220 = 0.129d0)
	91	c
	92	c
	93	real*8 a1_1221_0111, a1_1221_1220, a1_1221_0420
	94	c
	95	parameter (a1_1221_0420 = 0.01227d0)
	96	parameter (a1_1221_1220 = 355.2d0)
	97	parameter (a1_1221_0111 = 0.0219d0)
	98	c
	99	c
	100	real*8 a1_0510_000, a1_0510_0420, a1_0510_0400
	101	c
	102	parameter (a1_0510_000 = 4.d-4)
	103	parameter (a1_0510_0420 = 1.086d0)
	104	parameter (a1_0510_0400 = 2.406d0)
	105	c
	106	c
	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)
	111	c
	112	c
	113	real*8 a1_1310_000 , a1_1310_0420 , a1_1310_0400
	114	c
	115	parameter (a1_1310_000 = 4.45d-4)
	116	parameter (a1_1310_0420 = 0.7535d0)
	117	parameter (a1_1310_0400 = 0.21d0)
	118	c
	119	c
	120	real*8 a1_1310_1200 , a1_1310_1220 , a1_1310_2000
	121	c
	122	parameter (a1_1310_1200 = 2.264d0)
	123	parameter (a1_1310_1220 = 0.3609d0)
	124	parameter (a1_1310_2000 = 0.074d0)
	125	c
	126	c
	127	real*8 a1_1310_0111,a1_1310_1000,a1_1310_0220,a1_1310_0200
	128	c
	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)
	133	c
	134	c
	135	real*8 a1_2110_000 , a1_2110_0420 , a1_2110_0400
	136	c
	137	parameter (a1_2110_000 = 3.265d-4)
	138	parameter (a1_2110_0420 = 6.48d-3)
	139	parameter (a1_2110_0400 = 1.54d-3)
	140	c
	141	c
	142	real*8 a1_2110_1200,a1_2110_1220,a1_2110_2000,a1_2110_0111
	143	c
	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)
	148	c
	149	c
	150	real*8 a1_2110_1000 , a1_2110_0220 , a1_2110_0200
	151	c
	152	parameter (a1_2110_1000 = 3.445d-2)
	153	parameter (a1_2110_0220 = 1.976d-2)
	154	parameter (a1_2110_0200 = 6.539d-3)
	155	c
	156	c
	157	real*8 a1_0403_0402, a1_0403_0202
	158	c
	159	parameter (a1_0403_0402 = 728.7d0)
	160	parameter (a1_0403_0202 = 28.9d0)
	161	c
	162	real*8 a1_1203_1202, a1_1203_1002, a1_1203_0202
	163	c
	164	c
	165	parameter (a1_1203_1202 = 716.3d0)
	166	parameter (a1_1203_1002 = 27.25d0)
	167	parameter (a1_1203_0202 = 43.68d0)
	168	c
	169	c
	170	real*8 a1_2003_2002, a1_2003_1002, a1_2003_0202
	171	c
	172	parameter (a1_2003_2002 = 995.7d0)
	173	parameter (a1_2003_1002 = 83.27d0)
	174	parameter (a1_2003_0202 = 1.15d0)
	175	c
	176	c
	177	real*8 a1_0402_0401, a1_0402_0201
	178	c
	179	parameter (a1_0402_0401 = 700.0d0)
	180	parameter (a1_0402_0201 = 28.97d0)
	181	c
	182	c
	183	real*8 a1_1202_1201, a1_1202_1001, a1_1202_0201
	184	c
	185	parameter (a1_1202_1201 = 687.0d0)
	186	parameter (a1_1202_1001 = 26.07d0)
	187	parameter (a1_1202_0201 = 39.57d0)
	188	c
	189	c
	190	real*8 a1_2002_2001, a1_2002_1001
	191	c
	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.

Context Navigation

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

23/11/11 == FGG + MALV

Download in other formats: