1 MODULE module_data_racm
3 ! ****************************************************************
7 ! Generated by KPP - symbolic chemistry Kinetics PreProcessor
8 ! KPP is developed at CGRER labs University of Iowa by
9 ! Valeriu Damian & Adrian Sandu
11 ! File : racm_kpp_seulex.h
12 ! Time : Thu Apr 15 16:55:52 2004
14 ! ****************************************************************
19 PARAMETER (epsilc_racm=1.E-12)
20 ! NPHOT - The number of photolytic reactions
24 ! NSPEC - The number of species involved
28 ! NVAR - The number of Variable species
31 ! NVARACT - The number of Active species
33 PARAMETER (nvaract=69)
34 ! NRAD - The number of Radical species
37 ! NFIX - The number of Fixed species
40 ! NREACT - The number of reactions
43 PARAMETER (nreact=237)
44 ! NVARST - Starting of variables in conc. vect.
47 ! NRADST - Starting of radicals in conc. vect.
51 ! NFIXST - Starting of fixed in conc. vect.
54 ! NONZERO_V - Number of nonzero variable elements
56 PARAMETER (nonzero_v=925)
57 ! LU_NONZERO_V - Number of nonzero variable LU elements
59 PARAMETER (lu_nonzero_v=1051)
60 ! NONZERO_R - Number of nonzero radical elements
63 PARAMETER (nonzero_r=0)
64 ! CNVAR - (NVAR+1) The number of elements in compressed row format
66 PARAMETER (cnvar=nvar+1)
67 ! CNRAD - (NRAD+1) The number of elements in compressed row format
69 PARAMETER (cnrad=nrad+1)
70 ! NLOOKAT - number of species to look at
72 PARAMETER (nlookat=77)
73 ! NMASS - number of atoms to check mass balance
76 ! NX - X grid dimension
79 ! NY - Y grid dimension
83 ! NZ - Z grid dimension
86 ! MAX_XYZ - maximum grid dimension
89 ! NS - number of species to transport
94 PARAMETER (pi=3.14159265358979)
96 ! Declaration for global variables
100 ! VAR - Concentrations of variable species (global)
102 ! RAD - Concentrations of radical species (global)
104 ! FIX - Concentrations of fixed species (global)
106 !RS EQUIVALENCE (c(1+nvarst),var(1))
107 !RS EQUIVALENCE (c(1+nradst),rad(1))
108 !RS EQUIVALENCE (c(1+nfixst),fix(1))
110 !RS commonBs which have to be checked for WRF
112 !RS To be eliminated later
114 ! TIME - current integration time
116 !RS COMMON /gdata/time
117 ! SUN - light intensity
119 !RS COMMON /gdata/sun
120 ! C_DEFAULT - Default concentration for all species
121 !RS REAL*8 c_default(nspec)
122 !RS COMMON /gdata/c_default
123 ! C - Concentration for all species
126 ! TSTART - integration start time
128 !RS COMMON /gdata/tstart
129 ! TEND - integration end time
132 !RS COMMON /gdata/tend
133 ! SMASS - names of atoms for mass balance
134 !RS CHARACTER*12 smass(1)
135 !RS COMMON /chargdata/smass
138 !RS COMMON /gdata/temp
141 !RS COMMON /gdata/pres
142 ! DT - integration step
146 ! LOOKAT - indexes of species to look at
147 !RS INTEGER lookat(77)
148 !RS COMMON /intgdata/lookat
150 ! SLOOKAT - names of species to look at
151 !RS CHARACTER*12 slookat(77)
152 !RS COMMON /chargdata/slookat
153 ! NMONITOR - number of species to monitor
155 !RS PARAMETER (nmonitor=3)
156 ! MONITOR - indexes of species to monitor
157 !RS INTEGER monitor(3)
158 !RS COMMON /intgdata/monitor
159 ! SMONITOR - names of species to monitor
161 !RS CHARACTER*12 smonitor(3)
162 !RS COMMON /chargdata/smonitor
167 !RS COMMON blocks which ARE WRF compatible
168 ! RTOLS - (scalar) relative tolerance
170 !RS COMMON /gdata/rtols
171 ! ATOL - Absolute tolerance
172 !RS REAL*8 atol(nspec)
173 !RS COMMON /gdata/atol
174 ! RTOL - Relative tolerance
175 !RS REAL*8 rtol(nspec)
176 !RS COMMON /gdata/rtol
177 ! STEPMIN - minimum allowed intergation step
179 !RS COMMON /gdata/stepmin
180 ! STEPMAX - maximum allowed integration step
182 !RS COMMON /gdata/stepmax
183 ! CFACTOR - Conversion factor
185 !RS COMMON /gdata/cfactor
186 ! TRANS - indexes of species to transport
188 !RS COMMON /intgdata/trans
189 ! STRANS - names of species to transport
190 !RS CHARACTER*12 strans(1)
191 !RS COMMON /chargdata/strans
192 ! SEQN - equation names
193 CHARACTER*55 seqn(237)
194 !RS COMMON /chargdata/seqn
199 ! Indeces declaration for variable species
220 PARAMETER (i_n2o5=10)
232 PARAMETER (i_hno4=16)
236 PARAMETER (i_hono=18)
238 PARAMETER (i_h2o2=19)
242 PARAMETER (i_addt=21)
244 PARAMETER (i_addx=22)
246 PARAMETER (i_hket=23)
250 PARAMETER (i_addc=25)
254 PARAMETER (i_hno3=27)
266 PARAMETER (i_dien=33)
270 PARAMETER (i_tpan=35)
272 PARAMETER (i_etep=36)
277 PARAMETER (i_oltp=38)
279 PARAMETER (i_olip=39)
281 PARAMETER (i_mgly=40)
283 PARAMETER (i_cslp=41)
287 PARAMETER (i_limp=43)
289 PARAMETER (i_hc5p=44)
291 PARAMETER (i_hc8p=45)
293 PARAMETER (i_tolp=46)
295 PARAMETER (i_xylp=47)
297 PARAMETER (i_hcho=48)
299 PARAMETER (i_apip=49)
301 PARAMETER (i_isop=50)
303 PARAMETER (i_macr=51)
305 PARAMETER (i_hc3p=52)
311 PARAMETER (i_tco3=55)
319 PARAMETER (i_olnn=59)
321 PARAMETER (i_olnd=60)
323 PARAMETER (i_ethp=61)
329 PARAMETER (i_ketp=64)
333 PARAMETER (i_aco3=66)
337 PARAMETER (i_onit=68)
350 ! Indexes declaration for radical species
352 PARAMETER (ir_addt=1)
354 PARAMETER (ir_addx=2)
356 PARAMETER (ir_addc=3)
358 PARAMETER (ir_etep=4)
360 PARAMETER (ir_oltp=5)
362 PARAMETER (ir_olip=6)
364 PARAMETER (ir_cslp=7)
366 PARAMETER (ir_limp=8)
368 PARAMETER (ir_hc5p=9)
370 PARAMETER (ir_hc8p=10)
372 PARAMETER (ir_tolp=11)
374 PARAMETER (ir_xylp=12)
376 PARAMETER (ir_apip=13)
378 PARAMETER (ir_isop=14)
381 PARAMETER (ir_hc3p=15)
383 PARAMETER (ir_ethp=16)
385 PARAMETER (ir_o3p=17)
387 PARAMETER (ir_tco3=18)
390 PARAMETER (ir_mo2=19)
392 PARAMETER (ir_o1d=20)
394 PARAMETER (ir_olnn=21)
396 PARAMETER (ir_olnd=22)
398 PARAMETER (ir_rpho=23)
400 PARAMETER (ir_xo2=24)
402 PARAMETER (ir_ketp=25)
404 ! Indexes declaration for fixed species
415 ! User defined variables
419 !RS COMMON /gdata/stepstart
421 ! End user defined variables
423 ! ****************************************************************
425 ! Sparse Data Header File
427 ! Generated by KPP - symbolic chemistry Kinetics PreProcessor
428 ! KPP is developed at CGRER labs University of Iowa by
429 ! Valeriu Damian & Adrian Sandu
432 ! File : racm_kpp_seulex_s.h
433 ! Time : Thu Apr 15 16:55:52 2004
434 ! Working directory : /home/haas/Chemie/kpp_1.2/kpp-1.2/examples/racm
435 ! Equation file : racm_kpp_seulex.k
436 ! Output root filename : racm_kpp_seulex
438 ! ****************************************************************
444 !RS! IROW_V - row indexes for the Jacobian of variables
445 !RS INTEGER irow_v(nonzero_v)
446 !RS COMMON /sdata/irow_v
447 !RS! ICOL_V - column indexes for the Jacobian of variables
448 !RS INTEGER icol_v(nonzero_v)
449 !RS COMMON /sdata/icol_v
450 !RS! CROW_V - compressed row indexes for the Jacobian of variables
451 !RS INTEGER crow_v(cnvar)
452 !RS COMMON /sdata/crow_v
453 !RS! DIAG_V - diagonal indexes for the Jacobian of variables
454 !RS INTEGER diag_v(cnvar)
455 !RS COMMON /sdata/diag_v
456 !RS! LU_IROW_V - row indexes for the LU Jacobian of variables
457 !RS INTEGER lu_irow_v(lu_nonzero_v)
458 !RS COMMON /sdata/lu_irow_v
459 !RS! LU_ICOL_V - column indexes for the LU Jacobian of variables
460 INTEGER lu_icol_v(1051)
462 !RS COMMON /sdata/lu_icol_v
463 !RS! LU_CROW_V - compressed row indexes for the LU Jacobian of variables
464 INTEGER lu_crow_v(74)
465 !RS COMMON /sdata/lu_crow_v
466 !RS! LU_DIAG_V - diagonal indexes for the LU Jacobian of variables
467 INTEGER lu_diag_v(74)
468 !RS COMMON /sdata/lu_diag_v
469 !RS! IROW_R - row indexes for the Jacobian of radicals
470 INTEGER irow_r(nonzero_r+1)
471 !RS COMMON /sdata/irow_r
472 !RS! ICOL_R - column indexes for the Jacobian of radicals
473 INTEGER icol_r(nonzero_r+1)
474 !RS COMMON /sdata/icol_r
475 !RS! CROW_R - compressed row indexes for the Jacobian of radicals
476 INTEGER crow_r(cnrad)
477 !RS COMMON /sdata/crow_r
478 !RS! DIAG_R - diagonal indexes for the Jacobian of radicals
479 INTEGER diag_r(cnrad)
480 !RS COMMON /sdata/diag_r
488 !RS DATA lookat/1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, &
489 !RS 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, &
490 !RS 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, &
491 !RS 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, &
492 !RS 69, 70, 71, 72, 73, 74, 75, 76, 77/
495 !RS DATA slookat/'SULF', 'CO2', 'ORA1', 'ORA2', 'SO2', 'O1D', 'HC5', 'TOL', &
496 !RS 'XYL', 'N2O5', 'HC8', 'HC3', 'ETH', 'CH4', 'UDD', 'HNO4', 'OP1', &
497 !RS 'HONO', 'H2O2', 'PHO', 'ADDT', 'ADDX', 'HKET', 'ETE', 'ADDC', 'PAA', &
498 !RS 'HNO3', 'CO', 'API', 'LIM', 'PAN', 'CSL', 'DIEN', 'GLY', 'TPAN', &
499 !RS 'ETEP', 'ISO', 'OLTP', 'OLIP', 'MGLY', 'CSLP', 'KET', 'LIMP', 'HC5P', &
500 !RS 'HC8P', 'TOLP', 'XYLP', 'HCHO', 'APIP', 'ISOP', 'MACR', 'HC3P', 'ALD', &
501 !RS 'DCB', 'TCO3', 'XO2', 'OLT', 'OLI', 'OLNN', 'OLND', 'ETHP', 'O3P', &
502 !RS 'O3', 'KETP', 'MO2', 'ACO3', 'HO', 'ONIT', 'HO2', 'NO3', 'OP2', 'NO', &
503 !RS 'NO2', 'H2O', 'N2', 'O2', 'H2'/
506 !RS DATA monitor/63, 72, 73/
509 !RS DATA smonitor/'O3', 'NO', 'NO2'/
515 DATA (seqn(i),i=1,24)/' NO2 --> O3P + NO ' &
516 , ' O3 --> O1D + O2 ', &
517 ' O3 --> O3P + O2 ', &
518 ' HONO --> HO + NO ', &
519 ' HNO3 --> HO + NO2 ', &
520 ' HNO4 --> 0HO + 1HO2 + 0NO3 + 1NO2', &
521 ' NO3 --> NO + O2 ', &
522 ' NO3 --> O3P + NO2 ', &
524 ' HCHO --> CO + H2 ', &
525 ' HCHO --> CO + 2HO2 ', &
526 ' ALD --> CO + MO2 + HO2', &
527 ' OP1 --> HCHO + HO + HO2', &
528 ' OP2 --> ALD + HO + HO2', &
529 ' PAA --> MO2 + HO ', &
530 ' KET --> ETHP + ACO3', &
531 ' GLY --> 2CO + 0HCHO + 1H2', &
532 ' GLY --> 2CO + 0HCHO + 1HO2 + 0H2', &
533 ' MGLY --> CO + ACO3 + HO2', &
534 ' DCB --> TCO3 + HO2 ', &
535 ' ONIT --> 1KET + 0ALD + HO2 + NO2', &
536 ' MACR --> CO + HCHO + ACO3 + HO2', &
537 ' HKET --> HCHO + ACO3 + HO2', &
540 DATA (seqn(i),i=25,48)/ &
541 ' O3P + O3 --> 2O2 ', &
542 ' O1D + N2 --> O3P + N2 ', &
543 ' O1D + O2 --> O3P + O2 ', &
544 ' O1D + H2O --> 2HO ', &
545 ' O3 + HO --> HO2 + O2 ', &
546 ' O3 + HO2 --> HO + 2O2 ', &
547 ' HO + HO2 --> H2O + O2 ', &
548 ' H2O2 + HO --> HO2 + H2O ', &
549 ' 2HO2 --> H2O2 + O2 ', &
550 ' 2HO2 + H2O --> H2O2 + H2O + O2', &
551 ' O3P + NO --> NO2 ', &
552 ' O3P + NO2 --> NO + O2 ', &
553 ' O3P + NO2 --> NO3 ', &
554 ' HO + NO --> HONO ', &
555 ' HO + NO2 --> HNO3 ', &
556 ' HO + NO3 --> HO2 + NO2 ', &
557 ' HO2 + NO --> HO + NO2 ', &
558 ' HO2 + NO2 --> HNO4 ', &
559 ' HNO4 --> HO2 + NO2 ', &
560 ' HO2 + NO3 --> 0HNO3 + 1HO + 1NO2 + O2', &
561 ' HONO + HO --> NO2 + H2O ', &
562 ' HNO3 + HO --> NO3 + H2O ', &
563 ' HNO4 + HO --> NO2 + H2O + O2', &
564 ' O3 + NO --> NO2 + O2 '/
566 DATA (seqn(i),i=49,72)/ &
567 ' O3 + NO2 --> NO3 + O2 ', &
568 ' 2NO + O2 --> 2NO2 ', &
569 ' NO3 + NO --> 2NO2 ', &
570 ' NO3 + NO2 --> NO + NO2 + O2', &
571 ' NO3 + NO2 --> N2O5 ', &
572 ' N2O5 --> NO3 + NO2 ', &
573 ' 2NO3 --> 2NO2 + O2 ', &
574 ' HO + H2 --> HO2 + H2O ', &
575 ' SO2 + HO --> SULF + HO2 ', &
576 ' CO + HO --> CO2 + HO2 ', &
577 ' ISO + O3P --> 0CO + 0HCHO + 0DCB + 0XO2 + 1OLT + 0HO', &
578 ' MACR + O3P --> ALD ', &
579 ' CH4 + HO --> MO2 + H2O ', &
580 ' ETH + HO --> ETHP + H2O ', &
581 ' HC3 + HO --> 0ORA1 + 0CO + 0GLY + 0HCHO + 1HC3P + 0', &
582 ' HC5 + HO --> 0KET + 1HC5P + 0HO2 + H2O', &
583 ' HC8 + HO --> 0HKET + 1HC8P + 0ALD + 0HO2 + H2O', &
584 ' ETE + HO --> ETEP ', &
585 ' OLT + HO --> OLTP ', &
586 ' OLI + HO --> OLIP ', &
587 ' DIEN + HO --> ISOP ', &
588 ' ISO + HO --> ISOP ', &
589 ' API + HO --> APIP ', &
590 ' LIM + HO --> LIMP '/
592 DATA (seqn(i),i=73,96)/ &
593 ' TOL + HO --> 1ADDT + 0XO2 + 0HO2', &
594 ' XYL + HO --> 1ADDX + 0XO2 + 0HO2', &
595 ' CSL + HO --> 0PHO + 1ADDC + 0XO2 + 0HO2', &
596 ' HCHO + HO --> CO + HO2 + H2O', &
597 ' ALD + HO --> ACO3 + H2O ', &
598 ' KET + HO --> KETP + H2O ', &
599 ' HKET + HO --> MGLY + HO2 + H2O', &
600 ' GLY + HO --> 2CO + HO2 + H2O', &
601 ' MGLY + HO --> CO + ACO3 + H2O', &
602 ' MACR + HO --> 0HKET + 0CO + 0MGLY + 0HCHO + 1TCO3 + ', &
603 ' DCB + HO --> 0UDD + 0GLY + 0MGLY + 0TCO3 + 0XO2 + 0', &
604 ' UDD + HO --> 0KET + 1ALD + HO2', &
605 ' OP1 + HO --> 0HCHO + 1MO2 + 0HO', &
606 ' HO + OP2 --> 0KET + 0HC3P + 0ALD + 0XO2 + 0HO', &
607 ' PAA + HO --> 0HCHO + 0XO2 + 1ACO3 + 0HO2', &
608 ' PAN + HO --> HCHO + XO2 + NO3 + H2O', &
609 ' TPAN + HO --> 1HKET + 0PAN + 0HCHO + XO2 + 0HO2 + 1N', &
610 ' HO + ONIT --> HC3P + NO2 + H2O', &
611 ' HCHO + NO3 --> HNO3 + CO + HO2', &
612 ' ALD + NO3 --> HNO3 + ACO3', &
613 ' GLY + NO3 --> HNO3 + 2CO + HO2', &
614 ' MGLY + NO3 --> HNO3 + CO + ACO3', &
615 ' MACR + NO3 --> 0HNO3 + 1CO + 0TCO3 + 1OLNN', &
616 ' DCB + NO3 --> 0HNO3 + 0GLY + 0MGLY + 0KET + 0ALD + 0'/
619 DATA (seqn(i),i=97,120)/ &
620 ' CSL + NO3 --> PHO + HNO3 ', &
621 ' ETE + NO3 --> 1OLNN + 0OLND', &
622 ' OLT + NO3 --> 0OLNN + 1OLND', &
623 ' OLI + NO3 --> 0OLNN + 1OLND', &
624 ' DIEN + NO3 --> 1MACR + 1OLNN + 0OLND', &
625 ' ISO + NO3 --> 1MACR + 1OLNN + 0OLND', &
626 ' API + NO3 --> 0OLNN + 1OLND', &
627 ' LIM + NO3 --> 0OLNN + 1OLND', &
628 ' TPAN + NO3 --> 0PAN + 0HCHO + XO2 + 1ONIT + 1NO3 + 0N', &
629 ' ETE + O3 --> 0ORA1 + 0CO + HCHO + 0HO + 0HO2 + 0H2', &
630 ' OLT + O3 --> 0ORA1 + 0ORA2 + 0ETH + 0CH4 + 0H2O2 + ', &
631 ' OLI + O3 --> 0ORA2 + 0ETH + 0CH4 + 0H2O2 + 0CO + 0K', &
632 ' DIEN + O3 --> 0ORA1 + 0H2O2 + 0CO + 1HCHO + 0MACR + ', &
633 ' ISO + O3 --> 0ORA1 + 0H2O2 + 0CO + 1HCHO + 0MACR + ', &
634 ' API + O3 --> 0H2O2 + 0CO + 1KET + 1ALD + 0ETHP + 0K', &
635 ' LIM + O3 --> 0ORA1 + 0ORA2 + 0H2O2 + 0CO + 0HCHO + ', &
636 ' MACR + O3 --> 0ORA1 + 0ORA2 + 1CO + 1MGLY + 0HCHO + ', &
637 ' DCB + O3 --> 0ORA1 + 0ORA2 + 0PAA + 1CO + 0GLY + 1M', &
638 ' TPAN + O3 --> 0ORA1 + 0CO + 0PAN + 1HCHO + 1ACO3 + 0', &
639 ' PHO + NO2 --> 0CSL + ONIT', &
640 ' PHO + HO2 --> CSL ', &
641 ' ADDT + NO2 --> HONO + CSL ', &
642 ' ADDT + O2 --> 0CSL + 1TOLP + 0HO2', &
643 ' ADDT + O3 --> CSL + HO '/
645 DATA (seqn(i),i=121,144)/ &
646 ' ADDX + NO2 --> HONO + CSL ', &
647 ' ADDX + O2 --> 0CSL + 1XYLP + 0HO2', &
648 ' ADDX + O3 --> CSL + HO ', &
649 ' ADDC + NO2 --> HONO + CSL ', &
650 ' ADDC + O2 --> 0CSL + 1CSLP + 0HO2', &
651 ' ADDC + O3 --> CSL + HO ', &
652 ' ACO3 + NO2 --> PAN ', &
653 ' PAN --> ACO3 + NO2 ', &
654 ' TCO3 + NO2 --> TPAN ', &
655 ' TPAN --> TCO3 + NO2 ', &
656 ' MO2 + NO --> HCHO + HO2 + NO2', &
657 ' ETHP + NO --> ALD + HO2 + NO2', &
658 ' HC3P + NO --> 0GLY + 1KET + 0HCHO + 0ALD + 0XO2 + 0E', &
659 ' HC5P + NO --> 1KET + 0HCHO + 0ALD + 0XO2 + 0ETHP + 0', &
660 ' HC8P + NO --> 1KET + 0ALD + 0XO2 + 0ETHP + 0ONIT + 1', &
661 ' ETEP + NO --> 2HCHO + 0ALD + HO2 + NO2', &
662 ' OLTP + NO --> 0KET + HCHO + 1ALD + HO2 + NO2', &
663 ' OLIP + NO --> 0KET + 2ALD + HO2 + NO2', &
664 ' ISOP + NO --> 1HCHO + 0MACR + 0OLT + 0ONIT + 1HO2 + ', &
665 ' APIP + NO --> 1KET + 1ALD + 0ONIT + 1HO2 + 1NO2', &
666 ' LIMP + NO --> 0HCHO + 0MACR + 0OLI + 0ONIT + 1HO2 + ', &
667 ' TOLP + NO --> 1GLY + 1MGLY + 0DCB + 0ONIT + 1HO2 + 1', &
668 ' XYLP + NO --> 0GLY + 1MGLY + 1DCB + 0ONIT + 1HO2 + 1', &
669 ' CSLP + NO --> GLY + MGLY + HO2 + NO2'/
671 DATA (seqn(i),i=145,168)/ &
672 ' ACO3 + NO --> MO2 + NO2 ', &
673 ' TCO3 + NO --> HCHO + ACO3 + NO2', &
674 ' KETP + NO --> 1MGLY + 0ALD + 0XO2 + 0ACO3 + 1HO2 + N', &
675 ' OLNN + NO --> ONIT + HO2 + NO2', &
676 ' OLND + NO --> 0KET + 0HCHO + 1ALD + 2NO2', &
677 ' MO2 + HO2 --> OP1 ', &
678 ' ETHP + HO2 --> OP2 ', &
679 ' HC3P + HO2 --> OP2 ', &
680 ' HC5P + HO2 --> OP2 ', &
681 ' HC8P + HO2 --> OP2 ', &
682 ' ETEP + HO2 --> OP2 ', &
683 ' OLTP + HO2 --> OP2 ', &
684 ' OLIP + HO2 --> OP2 ', &
685 ' ISOP + HO2 --> OP2 ', &
686 ' APIP + HO2 --> OP2 ', &
687 ' LIMP + HO2 --> OP2 ', &
688 ' TOLP + HO2 --> OP2 ', &
689 ' XYLP + HO2 --> OP2 ', &
690 ' CSLP + HO2 --> OP2 ', &
691 ' ACO3 + HO2 --> PAA ', &
692 ' ACO3 + HO2 --> ORA2 + O3 ', &
693 ' TCO3 + HO2 --> OP2 ', &
694 ' TCO3 + HO2 --> ORA2 + O3 ', &
695 ' KETP + HO2 --> OP2 '/
697 DATA (seqn(i),i=169,192)/ &
698 ' OLNN + HO2 --> ONIT ', &
699 ' OLND + HO2 --> ONIT ', &
700 ' 2MO2 --> 1HCHO + 1HO2', &
701 ' ETHP + MO2 --> 1HCHO + 1ALD + HO2', &
702 ' HC3P + MO2 --> 0GLY + 0MGLY + 0KET + 1HCHO + 1ALD + 0', &
703 ' HC5P + MO2 --> 0KET + 1HCHO + 1ALD + 0XO2 + 0ETHP + 0', &
704 ' HC8P + MO2 --> 0KET + 1HCHO + 0ALD + 0XO2 + 0ETHP + 1', &
705 ' ETEP + MO2 --> 2HCHO + 0ALD + HO2', &
706 ' OLTP + MO2 --> 0KET + 1HCHO + 1ALD + HO2', &
707 ' OLIP + MO2 --> 0KET + 1HCHO + 1ALD + HO2', &
708 ' ISOP + MO2 --> 1HCHO + 1MACR + 0OLT + 0OLI + HO2', &
709 ' APIP + MO2 --> KET + HCHO + ALD + 2HO2', &
710 ' LIMP + MO2 --> 1HCHO + 1MACR + 0OLI + 2HO2', &
711 ' TOLP + MO2 --> 1GLY + 0MGLY + HCHO + DCB + HO2', &
712 ' XYLP + MO2 --> 0GLY + 1MGLY + HCHO + DCB + HO2', &
713 ' CSLP + MO2 --> GLY + MGLY + HCHO + 2HO2', &
714 ' MO2 + ACO3 --> HCHO + MO2 + HO2', &
715 ' MO2 + ACO3 --> ORA2 + HCHO', &
716 ' TCO3 + MO2 --> 2HCHO + ACO3 + HO2', &
717 ' TCO3 + MO2 --> ORA2 + HCHO', &
718 ' KETP + MO2 --> 0HKET + 0MGLY + 1HCHO + 0ALD + 0XO2 + ', &
719 ' OLNN + MO2 --> 1HCHO + ONIT + HO2', &
720 ' OLND + MO2 --> 0KET + 1HCHO + 1ALD + 0ONIT + 0HO2 + 0', &
721 ' ETHP + ACO3 --> 0ORA2 + ALD + 0MO2 + 0HO2'/
723 DATA (seqn(i),i=193,216)/ &
724 'HC3P + ACO3 --> 0ORA2 + 0GLY + 0MGLY + 0KET + 0HCHO + ', &
725 'HC5P + ACO3 --> 0ORA2 + 0KET + 0HCHO + 1ALD + 0XO2 + 0', &
726 'HC8P + ACO3 --> 0ORA2 + 1KET + 0ALD + 0XO2 + 0ETHP + 1', &
727 'ETEP + ACO3 --> 0ORA2 + 1HCHO + 1ALD + 0MO2 + 0HO2', &
728 'OLTP + ACO3 --> 0ORA2 + 0KET + 1HCHO + 1ALD + 1MO2 + 1', &
729 'OLIP + ACO3 --> 0ORA2 + 1KET + 1ALD + 1MO2 + 1HO2', &
730 'ISOP + ACO3 --> 0ORA2 + 0HCHO + 1MACR + 0OLT + 1MO2 + ', &
731 ' APIP + ACO3 --> KET + ALD + MO2 + HO2', &
732 'LIMP + ACO3 --> 0HCHO + 1MACR + 0OLI + MO2 + HO2', &
733 ' TOLP + ACO3 --> 1GLY + 0MGLY + DCB + MO2 + HO2', &
734 ' XYLP + ACO3 --> 0GLY + 1MGLY + DCB + MO2 + HO2', &
735 ' CSLP + ACO3 --> GLY + MGLY + MO2 + HO2', &
736 ' 2ACO3 --> 2MO2 ', &
737 ' TCO3 + ACO3 --> HCHO + MO2 + ACO3', &
738 'KETP + ACO3 --> 0ORA2 + 1MGLY + 0KET + 0ALD + 0XO2 + 0', &
739 ' OLNN + ACO3 --> 0ORA2 + 0MO2 + ONIT + 0HO2', &
740 'OLND + ACO3 --> 0ORA2 + 0KET + 0HCHO + 1ALD + 1MO2 + 0', &
741 ' 2OLNN --> 2ONIT + HO2', &
742 'OLNN + OLND --> 0KET + 0HCHO + 1ALD + 2ONIT + 0HO2 + 0', &
743 ' 2OLND --> 0KET + 1HCHO + 1ALD + ONIT + NO2', &
744 ' MO2 + NO3 --> HCHO + HO2 + NO2', &
745 ' ETHP + NO3 --> ALD + HO2 + NO2', &
746 ' HC3P + NO3 --> 0GLY + 1KET + 0HCHO + 0ALD + 0XO2 + 0E', &
747 ' HC5P + NO3 --> 1KET + 0HCHO + 0ALD + 0XO2 + 0ETHP + 0'/
749 DATA (seqn(i),i=217,237)/ &
750 ' HC8P + NO3 --> 1KET + 0ALD + 1XO2 + 0ETHP + 1HO2 + NO', &
751 ' ETEP + NO3 --> 2HCHO + 0ALD + HO2 + NO2', &
752 ' OLTP + NO3 --> 0KET + HCHO + 1ALD + HO2 + NO2', &
753 ' OLIP + NO3 --> 0KET + 2ALD + HO2 + NO2', &
754 ' ISOP + NO3 --> 1HCHO + 1MACR + 0OLT + HO2 + NO2', &
755 ' APIP + NO3 --> KET + ALD + HO2 + NO2', &
756 ' LIMP + NO3 --> 0HCHO + 1MACR + 0OLI + HO2 + NO2', &
757 ' TOLP + NO3 --> 1GLY + 1MGLY + 0DCB + HO2 + NO2', &
758 ' XYLP + NO3 --> 1GLY + 1MGLY + DCB + HO2 + NO2', &
759 ' CSLP + NO3 --> GLY + MGLY + HO2 + NO2', &
760 ' ACO3 + NO3 --> MO2 + NO2 ', &
761 ' TCO3 + NO3 --> HCHO + ACO3 + NO2', &
762 ' KETP + NO3 --> 1MGLY + 0ALD + 0XO2 + 0ACO3 + 1HO2 + N', &
763 ' OLNN + NO3 --> ONIT + HO2 + NO2', &
764 ' OLND + NO3 --> 0KET + 0HCHO + 1ALD + 2NO2', &
765 ' XO2 + HO2 --> OP2 ', &
766 ' XO2 + MO2 --> HCHO + HO2 ', &
767 ' XO2 + ACO3 --> MO2 ', &
769 ' XO2 + NO --> NO2 ', &
770 ' XO2 + NO3 --> NO2 '/
775 DATA (lu_icol_v(i),i=1,252)/1, 5, 67, 2, 28, 67, 3, 12, 24, 30, 33, 35, &
776 37, 51, 54, 57, 63, 67, 4, 30, 36, 38, 39, 44, 45, 50, 51, 52, 54, 55, &
777 57, 58, 59, 60, 61, 63, 64, 65, 66, 69, 5, 67, 6, 63, 7, 67, 8, 67, 9, &
778 67, 10, 70, 73, 11, 67, 12, 67, 13, 57, 58, 63, 67, 14, 57, 58, 63, &
779 67, 15, 54, 67, 16, 67, 69, 73, 17, 65, 67, 69, 18, 21, 22, 25, 67, &
780 72, 73, 19, 29, 30, 33, 37, 57, 58, 63, 67, 69, 20, 32, 67, 69, 70, &
781 73, 8, 21, 63, 67, 73, 9, 22, 63, 67, 73, 11, 23, 35, 51, 64, 65, 67, &
782 24, 63, 67, 70, 25, 32, 63, 67, 73, 26, 54, 63, 66, 67, 69, 27, 32, &
783 34, 40, 48, 51, 53, 54, 67, 69, 70, 73, 12, 24, 28, 29, 30, 33, 34, &
784 35, 37, 40, 48, 51, 53, 54, 57, 58, 62, 63, 67, 70, 29, 63, 67, 70, &
785 30, 63, 67, 70, 31, 35, 63, 66, 67, 70, 73, 20, 21, 22, 25, 32, 63, &
786 67, 69, 70, 73, 33, 63, 67, 70, 12, 34, 41, 46, 47, 52, 54, 63, 65, &
787 66, 67, 70, 72, 35, 55, 63, 67, 70, 73, 24, 36, 63, 65, 66, 67, 69, &
788 70, 72, 37, 62, 63, 67, 70, 38, 57, 65, 66, 67, 69, 70, 72, 39, 58, &
789 65, 66, 67, 69, 70, 72, 23, 35, 40, 41, 46, 47, 51, 52, 54/
791 DATA (lu_icol_v(i),i=253,504)/55, 63, 64, 65, 66, 67, 70, 72, 73, 25, &
792 32, 41, 63, 65, 66, 67, 69, 70, 72, 73, 7, 15, 29, 38, 39, 42, 44, 45, &
793 49, 52, 54, 57, 58, 59, 60, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, &
794 30, 43, 63, 65, 66, 67, 69, 70, 72, 7, 44, 65, 66, 67, 69, 70, 72, 11, &
795 45, 65, 66, 67, 69, 70, 72, 21, 46, 63, 65, 66, 67, 69, 70, 72, 73, &
796 22, 47, 63, 65, 66, 67, 69, 70, 72, 73, 12, 17, 23, 24, 26, 30, 31, &
797 33, 34, 35, 36, 37, 38, 39, 41, 43, 44, 45, 46, 47, 48, 49, 50, 51, &
798 52, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 69, 70, &
799 72, 73, 29, 49, 63, 65, 66, 67, 69, 70, 72, 33, 37, 50, 62, 63, 65, &
800 66, 67, 69, 70, 72, 30, 33, 37, 43, 50, 51, 62, 63, 65, 66, 67, 69, &
801 70, 72, 12, 52, 65, 66, 67, 68, 69, 70, 71, 72, 11, 12, 15, 29, 36, &
802 38, 39, 44, 45, 49, 51, 52, 53, 54, 57, 58, 59, 60, 61, 62, 63, 64, &
803 65, 66, 67, 68, 69, 70, 71, 72, 37, 46, 47, 54, 62, 63, 65, 66, 67, &
804 69, 70, 72, 73, 35, 51, 54, 55, 62, 63, 65, 66, 67, 69, 70, 72, 73, 8, &
805 9, 26, 31, 32, 33, 35, 37, 44, 45, 51, 52, 54, 55, 56, 62, 63, 64, 65/
807 DATA (lu_icol_v(i),i=505,756)/66, 67, 68, 69, 70, 71, 72, 73, 30, 33, &
808 37, 50, 57, 62, 63, 65, 66, 67, 69, 70, 72, 43, 50, 58, 62, 63, 65, &
809 66, 67, 69, 70, 72, 24, 29, 30, 33, 37, 51, 57, 58, 59, 60, 62, 63, &
810 65, 66, 67, 69, 70, 72, 24, 29, 30, 33, 37, 57, 58, 59, 60, 62, 63, &
811 65, 66, 67, 69, 70, 72, 13, 29, 30, 42, 44, 45, 49, 52, 54, 57, 58, &
812 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 6, 33, 37, &
813 51, 62, 63, 65, 66, 67, 69, 70, 72, 73, 21, 22, 24, 25, 29, 30, 32, &
814 33, 35, 37, 51, 54, 55, 57, 58, 62, 63, 65, 66, 67, 69, 70, 72, 73, &
815 29, 30, 33, 37, 42, 44, 45, 49, 52, 54, 57, 58, 59, 60, 62, 63, 64, &
816 65, 66, 67, 68, 69, 70, 71, 72, 73, 14, 17, 26, 33, 36, 37, 38, 39, &
817 41, 43, 44, 45, 46, 47, 49, 50, 52, 53, 54, 55, 56, 57, 58, 59, 60, &
818 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 23, 26, 31, 33, &
819 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, &
820 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, &
821 70, 71, 72, 73, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17, 18, 19, 21, &
824 DATA (lu_icol_v(i),i=757,1008)/23, 24, 25, 26, 27, 28, 29, 30, 31, 32, &
825 33, 34, 35, 37, 40, 41, 42, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, &
826 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, &
827 71, 72, 73, 20, 32, 35, 43, 44, 45, 46, 47, 49, 50, 52, 55, 59, 60, &
828 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 5, 7, 8, 9, 11, 12, 15, &
829 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 32, 33, 34, 35, &
830 36, 37, 38, 39, 40, 41, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, &
831 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, &
832 71, 72, 73, 10, 16, 24, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, &
833 39, 40, 41, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, &
834 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, &
835 36, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 55, 56, 57, 58, &
836 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 18, 21, 22, 25, &
837 32, 36, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 52, 55, 56, 57, 58, &
838 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 10, 16, &
841 DATA (lu_icol_v(i),i=1009,1051)/22, 25, 27, 31, 32, 34, 35, 36, 38, 39, &
842 40, 41, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, &
843 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73/
846 DATA lu_crow_v/1, 4, 7, 19, 41, 43, 45, 47, 49, 51, 54, 56, 58, 63, 68, &
847 71, 75, 79, 86, 96, 102, 107, 112, 119, 123, 128, 134, 146, 166, 170, &
848 174, 181, 191, 195, 208, 214, 223, 228, 236, 244, 262, 273, 298, 307, &
849 315, 323, 333, 343, 386, 395, 406, 420, 430, 460, 473, 486, 513, 526, &
850 537, 555, 572, 598, 611, 635, 661, 699, 741, 804, 829, 890, 938, 968, &
854 DATA lu_diag_v/1, 4, 7, 19, 41, 43, 45, 47, 49, 51, 54, 56, 58, 63, 68, &
855 71, 75, 79, 86, 96, 103, 108, 113, 119, 123, 128, 134, 148, 166, 170, &
856 174, 185, 191, 196, 208, 215, 223, 228, 236, 246, 264, 278, 299, 308, &
857 316, 324, 334, 363, 387, 397, 411, 421, 442, 463, 476, 500, 517, 528, &
858 545, 563, 585, 602, 627, 651, 690, 733, 797, 823, 885, 934, 965, 1002, &
862 END MODULE module_data_racm