source: CONFIG/publications/ICOLMDZORINCA_CO2_Transport_GMD_2023/INCA/src/INCA_MOD/input_data_tables_mod.F90 @ 6610

Last change on this file since 6610 was 6610, checked in by acosce, 10 months ago

INCA used for ICOLMDZORINCA_CO2_Transport_GMD_2023
File size: 15.8 KB
Line 
1!$Id: input_data_tables_mod.F90 104 2008-12-23 10:28:51Z acosce $
2!! =========================================================================
3!! INCA - INteraction with Chemistry and Aerosols
4!!
5!! Copyright Laboratoire des Sciences du Climat et de l'Environnement (LSCE)
6!!           Unite mixte CEA-CNRS-UVSQ
7!!
8!! Contributors to this INCA subroutine:
9!!
10!! Peter Hess, NCAR
11!! Jean-Francois Lamarque, NCAR
12!!
13!!
14!! Anne Cozic, LSCE, anne.cozic@cea.fr
15!! Yann Meurdesoif, LSCE, yann.meurdesoif@cea.fr
16!!
17!! This software is a computer program whose purpose is to simulate the
18!! atmospheric gas phase and aerosol composition. The model is designed to be
19!! used within a transport model or a general circulation model. This version
20!! of INCA was designed to be coupled to the LMDz GCM. LMDz-INCA accounts
21!! for emissions, transport (resolved and sub-grid scale), photochemical
22!! transformations, and scavenging (dry deposition and washout) of chemical
23!! species and aerosols interactively in the GCM. Several versions of the INCA
24!! model are currently used depending on the envisaged applications with the
25!! chemistry-climate model.
26!!
27!! This software is governed by the CeCILL  license under French law and
28!! abiding by the rules of distribution of free software.  You can  use,
29!! modify and/ or redistribute the software under the terms of the CeCILL
30!! license as circulated by CEA, CNRS and INRIA at the following URL
31!! "http://www.cecill.info".
32!!
33!! As a counterpart to the access to the source code and  rights to copy,
34!! modify and redistribute granted by the license, users are provided only
35!! with a limited warranty  and the software's author,  the holder of the
36!! economic rights,  and the successive licensors  have only  limited
37!! liability.
38!!
39!! In this respect, the user's attention is drawn to the risks associated
40!! with loading,  using,  modifying and/or developing or reproducing the
41!! software by the user in light of its specific status of free software,
42!! that may mean  that it is complicated to manipulate,  and  that  also
43!! therefore means  that it is reserved for developers  and  experienced
44!! professionals having in-depth computer knowledge. Users are therefore
45!! encouraged to load and test the software's suitability as regards their
46!! requirements in conditions enabling the security of their systems and/or
47!! data to be ensured and,  more generally, to use and operate it in the
48!! same conditions as regards security.
49!!
50!! The fact that you are presently reading this means that you have had
51!! knowledge of the CeCILL license and that you accept its terms.
52!! =========================================================================
53#include <inca_define.h>
54
55
56
57MODULE INPUT_DATA_TABLES
58
59  USE drydep_parameters
60  USE species_names
61
62  !
63  ! Table 1 from Wesely, Atmos. Environment, 1989, p1293
64  ! Table 2 from Sheih, microfiche PB86-218104 and Walcek, Atmos. Environment, 1986, p949
65  ! Table 3-5 compiled by P. Hess
66  !
67  ! index #1 : season
68  !           1 -> midsummer with lush vegetation
69  !           2 -> autumn with unharvested cropland
70  !           3 -> late autumn after frost, no snow
71  !           4 -> winter, snow on ground, and subfreezing
72  !           5 -> transitional spring with partially green short annuals
73  !
74  ! index #2 : landuse type
75  !           1 -> urban land
76  !           2 -> agricultural land
77  !           3 -> range land
78  !           4 -> deciduous forest
79  !           5 -> coniferous forest
80  !           6 -> mixed forest including wetland
81  !           7 -> water, both salt and fresh
82  !           8 -> barren land, mostly desert
83  !           9 -> nonforested wetland
84  !           10 -> mixed agricultural and range land
85  !           11 -> rocky open areas with low growing shrubs
86  !
87  ! JFL August 2000
88  !
89  REAL, DIMENSION(5,11), SAVE :: ri,rlu,rac,rgss,rgso,rcls,rclo
90!$OMP THREADPRIVATE(ri,rlu,rac,rgss,rgso,rcls,rclo)
91  !
92  data ri  (1,1:11) /1.e36,  60., 120.,  70., 130., 100.,1.e36,1.e36,  80., 100., 150./
93  data rlu (1,1:11) /1.e36,2000.,2000.,2000.,2000.,2000.,1.e36,1.e36,2500.,2000.,4000./
94  data rac (1,1:11) / 100., 200., 100.,2000.,2000.,2000.,   0.,   0., 300., 150., 200./
95  data rgss(1,1:11) / 400., 150., 350., 500., 500., 100.,   0.,1000.,   0., 220., 400./
96  data rgso(1,1:11) / 300., 150., 200., 200., 200., 300.,1000., 400.,1000., 180., 200./
97  data rcls(1,1:11) /1.e36,2000.,2000.,2000.,2000.,2000.,1.e36,1.e36,2500.,2000.,4000./
98  data rclo(1,1:11) /1.e36,1000.,1000.,1000.,1000.,1000.,1.e36,1.e36,1000.,1000.,1000./
99  !
100  data ri  (2,1:11) /1.e36,1.e36,1.e36,1.e36, 250., 500.,1.e36,1.e36,1.e36,1.e36,1.e36/
101  data rlu (2,1:11) /1.e36,9000.,9000.,9000.,4000.,8000.,1.e36,1.e36,9000.,9000.,9000./
102  data rac (2,1:11) / 100., 150., 100.,1500.,2000.,1700.,   0.,   0., 200., 120., 140./
103  data rgss(2,1:11) / 400., 200., 350., 500., 500., 100.,   0.,1000.,   0., 300., 400./
104  data rgso(2,1:11) / 300., 150., 200., 200., 200., 300.,1000., 400., 800., 180., 200./
105  data rcls(2,1:11) /1.e36,9000.,9000.,9000.,2000.,4000.,1.e36,1.e36,9000.,9000.,9000./
106  data rclo(2,1:11) /1.e36, 400., 400., 400.,1000., 600.,1.e36,1.e36, 400., 400., 400./
107  !
108  data ri  (3,1:11) /1.e36,1.e36,1.e36,1.e36, 250., 500.,1.e36,1.e36,1.e36,1.e36,1.e36/
109  data rlu (3,1:11) /1.e36,1.e36,9000.,9000.,4000.,8000.,1.e36,1.e36,9000.,9000.,9000./
110  data rac (3,1:11) / 100.,  10., 100.,1000.,2000.,1500.,   0.,   0., 100.,  50., 120./
111  data rgss(3,1:11) / 400., 150., 350., 500., 500., 200.,   0.,1000.,   0., 200., 400./
112  data rgso(3,1:11) / 300., 150., 200., 200., 200., 300.,1000., 400.,1000., 180., 200./
113  data rcls(3,1:11) /1.e36,1.e36,9000.,9000.,3000.,6000.,1.e36,1.e36,9000.,9000.,9000./
114  data rclo(3,1:11) /1.e36,1000., 400., 400.,1000., 600.,1.e36,1.e36, 800., 600., 600./
115  !
116  data ri  (4,1:11) /1.e36,1.e36,1.e36,1.e36, 400., 800.,1.e36,1.e36,1.e36,1.e36,1.e36/
117  data rlu (4,1:11) /1.e36,1.e36,1.e36,1.e36,6000.,9000.,1.e36,1.e36,9000.,9000.,9000./
118  data rac (4,1:11) / 100.,  10.,  10.,1000.,2000.,1500.,   0.,   0.,  50.,  10.,  50./
119  data rgss(4,1:11) / 100., 100., 100., 100., 100., 100.,   0.,1000., 100., 100.,  50./
120  data rgso(4,1:11) / 600.,3500.,3500.,3500.,3500.,3500.,1000., 400.,3500.,3500.,3500./
121  data rcls(4,1:11) /1.e36,1.e36,1.e36,9000., 200., 400.,1.e36,1.e36,9000.,1.e36,9000./
122  data rclo(4,1:11) /1.e36,1000.,1000., 400.,1500., 600.,1.e36,1.e36, 800.,1000., 800./
123  !
124  data ri  (5,1:11) /1.e36, 120., 240., 140., 250., 190.,1.e36,1.e36, 160., 200., 300./
125  data rlu (5,1:11) /1.e36,4000.,4000.,4000.,2000.,3000.,1.e36,1.e36,4000.,4000.,8000./
126  data rac (5,1:11) / 100.,  50.,  80.,1200.,2000.,1500.,   0.,   0., 200.,  60., 120./
127  data rgss(5,1:11) / 500., 150., 350., 500., 500., 200.,   0.,1000.,   0., 250., 400./
128  data rgso(5,1:11) / 300., 150., 200., 200., 200., 300.,1000., 400.,1000., 180., 200./
129  data rcls(5,1:11) /1.e36,4000.,4000.,4000.,2000.,3000.,1.e36,1.e36,4000.,4000.,8000./
130  data rclo(5,1:11) /1.e36,1000., 500., 500.,1500., 700.,1.e36,1.e36, 600., 800., 800./
131  !
132  ! roughness length
133  !
134  REAL, DIMENSION(5,11), SAVE :: z0
135!$OMP THREADPRIVATE(z0)
136  !
137  data z0  (1,1:11) /1.000,0.250,0.050,1.000,1.000,1.000,0.0006,0.002,0.150,0.100,0.100/
138  data z0  (2,1:11) /1.000,0.100,0.050,1.000,1.000,1.000,0.0006,0.002,0.100,0.080,0.080/
139  data z0  (3,1:11) /1.000,0.005,0.050,1.000,1.000,1.000,0.0006,0.002,0.100,0.020,0.060/
140  data z0  (4,1:11) /1.000,0.001,0.001,1.000,1.000,1.000,0.0006,0.002,0.001,0.001,0.040/
141  data z0  (5,1:11) /1.000,0.030,0.020,1.000,1.000,1.000,0.0006,0.002,0.010,0.030,0.060/
142  !
143  ! chemical data
144  !
145
146  INTEGER, PARAMETER :: n_species_table = 67
147  INTEGER, DIMENSION(ndep), SAVE  :: spec_map
148!$OMP THREADPRIVATE(spec_map)
149
150  REAL, DIMENSION(n_species_table*6), SAVE :: dheff
151  REAL, DIMENSION(n_species_table), SAVE   :: dfoxd
152  CHARACTER*20, DIMENSION(n_species_table), SAVE :: species_name_table
153!$OMP THREADPRIVATE(dheff, dfoxd)
154!$OMP THREADPRIVATE(species_name_table)
155  !
156
157  !--------------------------------------------------------------
158  ! Modif Anne 15/11/2005 : unification des tables entre
159  ! toutes les versions (AER CH4 CH4_AER NMHC NMHC_AER et GES)
160  !--------------------------------------------------------------
161  data species_name_table / 'O3'           &
162       ,'H2O2'          &
163       ,'HO'            &               
164       ,'OH'            &
165       ,'HO2'           &
166       ,'CO'            &
167       ,'CH4'           &
168       ,'CH3O2'         &
169       ,'CH3OOH'        &
170       ,'CH2O'          &
171       ,'CH2O2'         &
172       ,'NO'            &
173       ,'NO2'           &
174       ,'HNO2'          &
175       ,'HNO3'          &
176       ,'CO2'           &
177       ,'NH3'           &
178       ,'N2O5'          &
179       ,'NO3'           &
180       ,'CH3OH'         &
181       ,'HNO4'          &
182       ,'O1D'           &
183       ,'C2H6'          &
184       ,'CH3CH2OO'      &
185       ,'CH3CHOOCH2OH'  &
186       ,'CH3COCHOOCH2OH'&
187       ,'CH3CH2CH2CH3'  &
188       ,'CH3CHO'        &
189       ,'C2H5OOH'       &
190       ,'C3H6'          &
191       ,'CH3CHOOHCH2OH' &
192       ,'C2H4'          &
193       ,'PAN'           &  !'CH3COOONO2'
194       ,'CH3COOOH'      &
195       ,'C9H13CH3'      &
196       ,'CHOCHO'        &
197       ,'CH3COCHO'      &
198       ,'CH2OHCHO'      &
199       ,'CH3COOO'       &
200       ,'C3H8'          &
201       ,'CH3CHOOCH3'    &
202       ,'CH3COCH3'      &
203       ,'CH3CHOOHCH3'   &
204       ,'CH3COCH2OO'    &
205       ,'CH3COCH2O'     &
206       ,'CH3COCH2OOH'   &
207       ,'Radon'         &
208       ,'MCF'           &
209       ,'H2'            &
210       ,'C2H2'          &
211       ,'ISOP'          &
212       ,'APIN'          &
213       ,'C2H5OH'        &
214       ,'CH3COOH'       &
215       ,'MACR'          &
216       ,'MVK'           &
217       ,'MEK'           &
218       ,'ONITR'         &
219       ,'ONITU'         &
220       ,'SO2'          &
221       ,'DMS'           &
222       ,'H2S'           &
223       ,'DMSO'          &
224       ,'APp1g'         &
225       ,'APp2g'         &
226       ,'ARp1g'         &
227       ,'ARp2g'         &
228       /
229
230  !
231
232  !t
233  !    data in  [mol/l/atm]
234  !
235  !                 He       dH/R   Ks1     dH/R   Ks2     dH/R
236  !t
237  data dheff /1.10e-02, 2300.,0.     ,    0.,0.     ,    0.  &
238       ,8.33e+04, 7379.,2.2e-12,-3730.,0.     ,    0.   &
239       ,3.00e+01,    0.,0.     ,    0.,0.     ,    0.   & 
240       ,3.00e+01, 4500.,0.     ,    0.,0.     ,    0.   &
241       ,2.00e+03, 6600.,3.5e-05,    0.,0.     ,    0.   &
242       ,9.90e-04, 1300.,0.     ,    0.,0.     ,    0.   &
243       ,1.40e-03, 1600.,0.     ,    0.,0.     ,    0.   &
244       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
245       ,3.11e+02, 5241.,0.     ,    0.,0.     ,    0.   &
246       ,3.20e+03, 6800.,0.     ,    0.,0.     ,    0.   &
247       ,5.53e+03, 5700.,1.8e-04,-1510.,0.     ,    0.   &
248       ,1.90e-03, 1480.,0.     ,    0.,0.     ,    0.   &
249       ,1.20e-02, 2500.,0.     ,    0.,0.     ,    0.   &
250       ,5.00e+01, 4900.,0.     ,    0.,0.     ,    0.   &
251       ,0.      ,    0.,2.4e+06, 8700.,0.     ,    0.   &
252       ,3.40e-02, 2420.,4.5e-07,-1000.,3.6e-11,-1760.   &
253       ,5.90e+01, 4200.,1.7e-05, -450.,1.0e-14,-6716.   &
254       ,2.14e+00, 3362.,0.     ,    0.,0.     ,    0.   &
255       ,2.00e+00, 2000.,0.     ,    0.,0.     ,    0.   &
256       ,2.20e+02, 5200.,0.     ,    0.,0.     ,    0.   &
257       ,1.20e+04, 6900.,0.     ,    0.,0.     ,    0.   &
258       ,1.00e-16,    0.,0.     ,    0.,0.     ,    0.   &
259       ,1.90e-03, 2300.,0.     ,    0.,0.     ,    0.   &
260       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
261       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
262       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
263       ,1.70e-03,    0.,0.     ,    0.,0.     ,    0.   &
264       ,1.40e+01, 5600.,0.     ,    0.,0.     ,    0.   &
265       ,3.36e+02, 5995.,0.     ,    0.,0.     ,    0.   &
266       ,7.40e-03, 3400.,0.     ,    0.,0.     ,    0.   &
267       ,2.20e+02, 5653.,0.     ,    0.,0.     ,    0.   &
268       ,4.70e-03, 1800.,0.     ,    0.,0.     ,    0.   &
269       ,2.80e+00, 6500.,0.     ,    0.,0.     ,    0.   &
270       ,8.37e+02, 5308.,1.8e-04,-1510.,0.     ,    0.   &
271       ,1.70e-03,    0.,0.     ,    0.,0.     ,    0.   &
272       ,3.00e+05,    0.,0.     ,    0.,0.     ,    0.   &
273       ,3.71e+03, 7541.,0.     ,    0.,0.     ,    0.   &
274       ,4.14e+04, 4630.,0.     ,    0.,0.     ,    0.   &
275       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
276       ,1.45e-03, 2700.,0.     ,    0.,0.     ,    0.   &
277       ,3.00e+06,    0.,0.     ,    0.,0.     ,    0.   &
278       ,3.00e+01, 4600.,0.     ,    0.,0.     ,    0.   &
279       ,3.36e+02, 5995.,0.     ,    0.,0.     ,    0.   &
280       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
281       ,7.47e+00, 5241.,0.     ,    0.,0.     ,    0.   &
282       ,3.36e+02, 5995.,0.     ,    0.,0.     ,    0.   &
283       ,9.30e-03, 2600.,0.     ,    0.,0.     ,    0.   &
284       ,5.90e-02, 3900.,0.     ,    0.,0.     ,    0.   &
285       ,7.80e-04,  500.,0.     ,    0.,0.     ,    0.   &
286       ,4.10e-02, 1800.,0.     ,    0.,0.     ,    0.   &
287       ,2.80e-02,    0.,0.     ,    0.,0.     ,    0.   &
288       ,4.90e-02,    0.,0.     ,    0.,0.     ,    0.   &
289       ,1.90e+02, 6600.,0.     ,    0.,0.     ,    0.   &
290       ,4.10e+03, 6300.,0.     ,    0.,0.     ,    0.   &
291       ,6.50e+00,    0.,0.     ,    0.,0.     ,    0.   &
292       ,2.10e+01, 7800.,0.     ,    0.,0.     ,    0.   &
293       ,2.00e+01, 5000.,0.     ,    0.,0.     ,    0.   &
294       ,7.90e-01, 5400.,0.     ,    0.,0.     ,    0.   &
295       ,7.90e-01, 5400.,0.     ,    0.,0.     ,    0.   &
296       ,1.20e+00, 3200.,1.7e-02, 2090.,6.0e-08, 1120.   & 
297       ,4.80e-01, 3100.,0.     ,    0.,0.     ,    0.   & 
298       ,1.00e-03, 2300.,5.7e-08,    0.,1.3e-13,    0.   & 
299       ,5.00e+04,    0.,0.     ,    0.,0.     ,    0.   &   
300       ,1.00e+04,    0.,0.     ,    0.,0.     ,    0.   &
301       ,1.00e+03,    0.,0.     ,    0.,0.     ,    0.   &
302       ,1.00e+04,    0.,0.     ,    0.,0.     ,    0.   &
303       ,1.00e+03,    0.,0.     ,    0.,0.     ,    0.   &
304       /
305  !
306
307  data dfoxd / 1.0     &
308       ,1.0      &
309       ,1.0      &
310       ,1.       & 
311       ,0.1      &
312       ,1.e-36   &
313       ,1.e-36   &
314       ,0.1      &
315       ,0.1      &
316       ,1.e-36   &
317       ,0.1      &
318       ,1.e-36   &
319       ,0.1      &
320       ,1.e-36   &
321       ,1.e-36   &
322       ,1.e-36   &
323       ,1.e-36   &
324       ,0.1      &
325       ,0.1      &
326       ,1.e-36   &
327       ,0.1      &
328       ,1.E-36   &
329       ,1.e-36   &
330       ,0.1      &
331       ,0.1      &
332       ,0.1      &
333       ,1.e-36   &
334       ,1.e-36   &
335       ,0.1      &
336       ,1.e-36   &
337       ,0.1      &
338       ,1.e-36   &
339       ,1.e-36   &
340       ,0.1      &
341       ,1.e-36   &
342       ,1.e-36   &
343       ,1.e-36   &
344       ,1.e-36   &
345       ,0.1      &
346       ,1.e-36   &
347       ,0.1      &
348       ,0.1      &
349       ,0.1      &
350       ,0.1      &
351       ,0.1      &
352       ,0.1      &
353       ,1.e-36   &
354       ,1.e-36   &
355       ,1.e-36   &
356       ,1.e-36   &
357       ,1.0      &
358       ,0.1      &
359       ,1.e-36   &
360       ,1.e-36   &
361       ,1.e-36   &
362       ,0.1      &
363       ,0.1      &
364       ,1.0      &
365       ,0.1      &
366       ,1.e-36   &     
367       ,1.e-36   &     
368       ,1.e-36   &
369       ,.1       &             
370       ,1.e-36   &
371       ,1.e-36   &
372       ,1.e-36   &
373       ,1.e-36   &
374       /
375
376  !
377
378
379
380#ifdef NMHC
381  data spec_map / &
382     id_O3, id_HO2, id_OH, id_NO, id_NO2, id_NO3, id_HNO2, id_HNO3,        &
383     id_HNO4, id_N2O5, id_H2O2, id_CH3OOH, id_CH2O, id_CH3OH, id_C2H5OH,   &
384     id_CH3CHO, id_CH3COCHO, id_CH3COCH3, id_MEK, id_MVK, id_CH3COOH,      &
385     id_CH3COOOH, id_C2H5OOH, id_PAN, id_ONITR, id_ONITU, id_ISOP,         &
386     id_C2H6, id_C3H6, id_C2H4, id_C3H8, id_C2H2, id_APIN, id_MACR         &
387#ifdef AER
388     ,id_SO2, id_DMS, id_H2S, id_DMSO, id_NH3 &
389!SOA
390     ,id_APp1g &
391     ,id_APp2g &
392     ,id_ARp1g &
393     ,id_ARp2g &
394#endif
395     /
396#endif
397
398#if !defined(DUSS) && defined(AERONLY)
399  data spec_map / id_SO2, id_NO2, id_NO3, id_HNO3, id_N2O5, id_DMS, id_H2S, id_DMSO, id_NH3 /
400#endif
401
402#ifdef GES
403  data spec_map /id_CO /
404#endif
405
406end module INPUT_DATA_TABLES
Note: See TracBrowser for help on using the repository browser.