! #define TRACEBACK write (gol,'("in ",a," (",a,", line",i5,")")') rname, __FILE__, __LINE__; call goErr #define IF_NOTOK_RETURN(action) if (status/=0) then; TRACEBACK; action; return; end if #define IF_ERROR_RETURN(action) if (status> 0) then; TRACEBACK; action; return; end if ! #include "tm5.inc" ! module reaction_data implicit none public integer, parameter :: nreac=112 !number of thermal reactions integer, parameter :: nreacw=4 !number of aqueous phase reactions integer, parameter :: nthigh=320 integer, parameter :: ntlow=165 integer, parameter :: ntemp=nthigh-ntlow integer, parameter :: nrat=nreac+50 character(len=10), dimension(nreac),parameter :: ratnam=(/& 'RNOO3 ','RHO2NO ','RMO2NO ','RNO2OH ','ROHHNO3 ',& 'RNO2O3 ','RNONO3 ','RNO2NO3 ','RN2O5 ','RHNO4OH ',& 'RNO2HO2 ','RHNO4M ','RODM ','RH2OOD ','RO3HO2 ',& 'RCOOH ','RO3OH ','RHPOH ','RFRMOH ','RCH4OH ',& 'ROHMPER ','ROHROOH ','RMO2HO2a ','RMO2HO2b ','RMO2MO2 ',& 'RHO2OH ','RHO2HO2 ','RN2O5L ','RH2OH ','RC41 ',& 'RC43 ','RC44 ','RC46 ','RC47 ','RC48 ',& 'RC49 ','RC50 ','RC52 ','RC53 ','RC54 ',& 'RC57 ','RC58 ','RC59 ','RC61 ','RC62 ',& 'RC73 ','RC76 ','RC77 ','RC78 ','RC79 ',& 'RC80 ','RC81 ','RC82 ','RC83 ','RC84 ',& 'RC85 ','Rdmsoha ','Rdmsohb ','Rdmsno3 ','Rso2oh ',& 'RCnh3so4 ','RCnh3oh ','RCnh2oh ','RCnh2no ','RCnh2no2 ',& 'RCnh2ho2 ','RCnh2o2 ','RCnh2o3 ','ROHCH3OH ','ROHHCOOH ',& 'RNO3HO2 ','RNO3MO2 ','RNO3C2O3 ','RNO3XO2 ','RMCOOH ',& 'ROHC2H6 ','ROHETHOH ','ROHC3H8 ','ROHC3H6 ','RO3C3H6 ',& 'RNO3C3H6 ','ROHTERP ','RO3TERP ','RNO3TERP ','ROHISPD ',& 'RO3ISPD ','RNO3ISPD ','DRN222 ','RCO3PO3 ','RCO3PO2 ',& 'ROHACET ','RACO2HO2 ','RACO2MO2 ','RACO2NO ','RACO2XO2 ',& 'RXO2XO2N ','RXO2N ','RC3H7O2NO ','C3H7O2HO2 ','RHYPO2NO ',& 'RHYPO2HO2 ','RN2O5AER ','RNO3_AER ','RHO2_AER ','RHO2_L ',& 'RNH2O2NO ','RNH2O2O3 ','RNH2O2HO2 ','RHONO ','ROHHONO ',& 'RMO2NO2 ','RMO2NO2M '/) character(len=8), dimension(nreacw),parameter :: rwnam = & (/'kso2o3 ','kso2hp ','knh3so4 ','void '/) ! ! reaction rates ! integer, parameter :: KNOO3=1 integer, parameter :: KHO2NO=2 integer, parameter :: KMO2NO=3 integer, parameter :: KNO2OH=4 integer, parameter :: KOHHNO3=5 integer, parameter :: KNO2O3=6 integer, parameter :: KNONO3=7 integer, parameter :: KNO2NO3=8 integer, parameter :: KN2O5=9 integer, parameter :: KHNO4OH=10 integer, parameter :: KNO2HO2=11 integer, parameter :: KHNO4M=12 integer, parameter :: KODM=13 integer, parameter :: KH2OOD=14 integer, parameter :: KO3HO2=15 integer, parameter :: KCOOH=16 integer, parameter :: KO3OH=17 integer, parameter :: KHPOH=18 integer, parameter :: KFRMOH=19 integer, parameter :: KCH4OH=20 integer, parameter :: KOHMPER=21 integer, parameter :: KOHROOH=22 integer, parameter :: KMO2HO2A=23 integer, parameter :: KMO2HO2B=24 integer, parameter :: KMO2MO2=25 integer, parameter :: KHO2OH=26 integer, parameter :: KHO2HO2=27 integer, parameter :: KN2O5L=28 integer, parameter :: KH2OH=29 integer, parameter :: KC41=30 integer, parameter :: KC43=31 integer, parameter :: KC44=32 integer, parameter :: KC46=33 integer, parameter :: KC47=34 integer, parameter :: KC48=35 integer, parameter :: KC49=36 integer, parameter :: KC50=37 integer, parameter :: KC52=38 integer, parameter :: KC53=39 integer, parameter :: KC54=40 integer, parameter :: KC57=41 integer, parameter :: KC58=42 integer, parameter :: KC59=43 integer, parameter :: KC61=44 integer, parameter :: KC62=45 integer, parameter :: KC73=46 integer, parameter :: KC76=47 integer, parameter :: KC77=48 integer, parameter :: KC78=49 integer, parameter :: KC79=50 integer, parameter :: KC80=51 integer, parameter :: KC81=52 integer, parameter :: KC82=53 integer, parameter :: KC83=54 integer, parameter :: KC84=55 integer, parameter :: KC85=56 integer, parameter :: Kdmsoha=57 integer, parameter :: Kdmsohb=58 integer, parameter :: Kdmsno3=59 integer, parameter :: Kso2oh=60 integer, parameter :: Knh3so4=61 integer, parameter :: Knh3oh=62 integer, parameter :: knh2oh=63 integer, parameter :: Knh2no=64 integer, parameter :: Knh2no2=65 integer, parameter :: Knh2ho2=66 integer, parameter :: Knh2o2=67 integer, parameter :: Knh2o3=68 ! additional oxidation of BVOC integer, parameter :: kohch3oh=69 integer, parameter :: kohhcooh=70 ! ! additional NO3 reactions ! integer, parameter :: kno3ho2=71 integer, parameter :: kno3mo2=72 integer, parameter :: kno3c2o3=73 integer, parameter :: kno3xo2=74 ! ! integer, parameter :: kohmcooh=75 integer, parameter :: kohc2h6=76 integer, parameter :: kohethoh=77 ! ! C3 ! integer, parameter :: kohc3h8=78 integer, parameter :: kohc3h6=79 integer, parameter :: ko3c3h6=80 integer, parameter :: kno3c3h6=81 ! ! TERP ! integer, parameter :: kohterp=82 integer, parameter :: ko3terp=83 integer, parameter :: kno3terp=84 integer, parameter :: kohispd=85 integer, parameter :: ko3ispd=86 integer, parameter :: kno3ispd=87 ! integer, parameter :: krn222=88 ! for JO2 integer, parameter :: ko3po3=89 integer, parameter :: ko3po2=90 ! Acetone integer, parameter :: kohacet=91 integer, parameter :: kaco2ho2=92 integer, parameter :: kaco2mo2=93 integer, parameter :: kaco2no=94 integer, parameter :: kaco2xo2=95 ! JEW : ACO2 + XO2 -> (same as for XO2N) ! ! From CB05 ! integer, parameter :: kxo2xo2n=96 integer, parameter :: kxo2n=97 ! ! extra peroxy reactions to form acetone integer, parameter :: kc3h7o2no=98 integer, parameter :: kc3h7o2ho2=99 integer, parameter :: khypo2no=100 integer, parameter :: khypo2ho2=101 ! ! extra heterogeneous recations ! integer, parameter :: kn2o5_aer=102 integer, parameter :: kno3_aer=103 integer, parameter :: kho2_aer=104 integer, parameter :: kho2_l=105 integer, parameter :: knh2o2no=106 integer, parameter :: knh2o2o3=107 integer, parameter :: knh2o2ho2=108 ! ! HONO ! integer, parameter :: khono=109 integer, parameter :: kohhono=110 ! ! methyl peroxy nitrate ! integer, parameter :: kmo2no2=111 integer, parameter :: kch3o2no2m=112 ! ! additional parameters needed for 3-body and other reactions ! integer, parameter :: KNO2OHA=nreac+1 integer, parameter :: KNO2OHB=nreac+2 integer, parameter :: KOHHNO3A=nreac+3 integer, parameter :: KOHHNO3B=nreac+4 integer, parameter :: KOHHNO3C=nreac+5 integer, parameter :: KNO2NO3A=nreac+6 integer, parameter :: KNO2NO3B=nreac+7 integer, parameter :: KN2O5A=nreac+8 integer, parameter :: KN2O5B=nreac+9 integer, parameter :: KNO2HO2A=nreac+10 integer, parameter :: KNO2HO2B=nreac+11 integer, parameter :: KHNO4A=nreac+12 integer, parameter :: KHNO4B=nreac+13 integer, parameter :: KHONOA=nreac+14 integer, parameter :: KHONOB=nreac+15 integer, parameter :: KHO2HO2A=nreac+16 integer, parameter :: KHO2HO2B=nreac+17 integer, parameter :: KHO2HO2C=nreac+18 integer, parameter :: KC47A=nreac+19 integer, parameter :: KC47B=nreac+20 integer, parameter :: KC48A=nreac+21 integer, parameter :: KC48B=nreac+22 integer, parameter :: KC61A=nreac+23 integer, parameter :: KC61B=nreac+24 integer, parameter :: kso2oha=nreac+25 integer, parameter :: kso2ohb=nreac+26 integer, parameter :: kdmsohc=nreac+27 !JEW new parameters needed for the new CO + OH formulation integer, parameter :: kcooha=nreac+28 integer, parameter :: kcoohb=nreac+29 integer, parameter :: kcoohc=nreac+30 integer, parameter :: kcoohd=nreac+31 !JEW new parameters needed for the update to OH + organics integer, parameter :: kohmcho=nreac+32 integer, parameter :: kohmch2cho=nreac+33 integer, parameter :: kno3mcho=nreac+34 integer, parameter :: kno3mch2cho=nreac+35 ! JEW ; OLE reactions now calcuated using MVK and MACR rates to ! account for link of OLE to ISOP oxidation . For NO3 ! these are fed directly into to rr in chem_rates.F90 integer, parameter :: kohole=nreac+36 integer, parameter :: ko3ole=nreac+37 integer, parameter :: kno3ole=nreac+38 integer, parameter :: kohaceta=nreac+39 integer, parameter :: kohacetb=nreac+40 integer, parameter :: kohc3h6a=nreac+41 integer, parameter :: kohc3h6b=nreac+42 integer, parameter :: ko3mvk=nreac+43 integer, parameter :: ko3macr=nreac+44 integer, parameter :: kohmvk=nreac+45 integer, parameter :: kohmacr=nreac+46 integer, parameter :: kmo2no2a=nreac+47 integer, parameter :: kmo2no2b=nreac+48 integer, parameter :: kch3o2no2a=nreac+49 integer, parameter :: kch3o2no2b=nreac+50 integer, parameter :: kso2hp=2 integer, parameter :: kso2o3=1 ! rates_lut : for reaction rate look up table !WP! real,dimension(nrat,ntemp) :: rates_lut ! for N2O5 removal on aerosol parameterisation: ! nr_interval : number of interval for intergration over aerosol size integer,parameter :: nr_interval=13 ! np_interval : number of interval for integration over pressure integer,parameter :: np_interval=11 integer,parameter :: ncomponent=2 ! aerdens : density of aerosols(g/cm3) (water 1.0 ) real,parameter :: aerdens=1.8 ! rra : lookup table radius.. real,dimension(nr_interval):: rra = & (/1.0,1.25,1.5,1.75,2.0,2.25,2.5,2.75,3.,4.,5.,10.,15./) ! hetrem : heterogeneous removal coefficient [s-1] ! for 1 ppbv NH4HSO4 at given aerosol radius and pressure real,dimension(nr_interval,np_interval,ncomponent) :: hetrem real,dimension(180,12) :: h2_lat data h2_lat(1:180,1) /545.6847, 545.6823, 545.6799, 545.6778, 545.6754, 545.6731, 545.6709, 545.6686, 545.6663, 545.664, 545.6618, 545.6594, 545.6572, 545.655, & 545.6526, 545.6503, 545.6481, 545.6458, 545.6435, 545.6413, 545.639, 545.6367, 545.6344, 545.6321, 545.6298, 545.6032, 545.5523, 545.5013, & 545.4503, 545.3994, 545.3484, 545.2974, 545.2463, 545.1953, 545.1444, 545.0934, 545.0425, 544.9916, 544.9406, 544.8895, 544.8386, 544.7876, & 544.7366, 544.6857, 544.6347, 544.5837, 544.5327, 544.4818, 544.4308, 544.3798, 544.3288, 544.2778, 544.2268, 544.1759, 544.1249, 544.0739, & 544.0229, 543.972, 543.921, 543.8699, 543.819, 543.768, 543.717, 543.6661, 543.6151, 543.5641, 543.5131, 543.4622, 543.4113, 543.3602, & 543.3093, 543.2583, 543.2073, 543.1562, 543.1053, 543.0543, 542.7801, 542.2826, 541.7851, 541.2877, 540.7902, 540.2926, 539.7951, 539.2977, & 538.8001, 538.3026, 537.8052, 537.3076, 536.8102, 536.3127, 535.8152, 535.3176, 534.8202, 534.3227, 533.8251, 533.3276, 532.8302, 532.3327, & 531.8352, 531.3376, 530.8402, 530.3428, 529.8451, 529.3477, 528.8503, 528.3527, 527.8552, 527.3578, 526.8602, 526.3627, 525.7585, 525.0474, & 524.3364, 523.6252, 522.9142, 522.2032, 521.4922, 520.7811, 520.0701, 519.3591, 518.648, 517.937, 517.2259, 516.5148, 515.8038, 515.0928, & 514.3817, 513.6707, 512.9596, 512.2486, 511.4813, 510.6577, 509.834, 509.0103, 508.1867, 507.3629, 506.5393, 505.7157, 504.892, 504.0684, & 503.2448, 502.4212, 501.5975, 500.7738, 499.9502, 499.1265, 498.303, 497.4792, 496.6556, 495.8319, 495.0083, 494.1847, 493.361, 492.5374, & 491.7137, 490.8901, 490.0663, 489.2428, 488.4191, 487.5955, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, & 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837, 487.1837/ data h2_lat(1:180,2) /546.3832, 546.4077, 546.4323, 546.4571, 546.4817, 546.5063, 546.531, 546.5557, 546.5803, 546.605, 546.6297, 546.6543, 546.679, 546.7036, & 546.7283, 546.7529, 546.7776, 546.8022, 546.8269, 546.8516, 546.8762, 546.9009, 546.9255, 546.9503, 546.9749, 546.9255, 546.8022, 546.6788, & 546.5555, 546.4321, 546.3088, 546.1854, 546.062, 545.9386, 545.8153, 545.692, 545.5686, 545.4454, 545.3219, 545.1984, 545.0752, 544.9518, & 544.8284, 544.7051, 544.5818, 544.4584, 544.335, 544.2117, 544.0883, 543.965, 543.8416, 543.7181, 543.5948, 543.4715, 543.3481, 543.2247, & 543.1014, 542.9781, 542.8547, 542.7313, 542.608, 542.4846, 542.3612, 542.2379, 542.1144, 541.9911, 541.8677, 541.7444, 541.6212, 541.4977, & 541.3744, 541.251, 541.1276, 541.0042, 540.8809, 540.7575, 540.5496, 540.2568, 539.9642, 539.6716, 539.3789, 539.0861, 538.7935, 538.5009, & 538.2081, 537.9154, 537.6228, 537.3301, 537.0374, 536.7448, 536.4521, 536.1594, 535.8667, 535.574, 535.2813, 534.9886, 534.696, 534.4033, & 534.1107, 533.8179, 533.5253, 533.2327, 532.9398, 532.6472, 532.3546, 532.0619, 531.7692, 531.4766, 531.1838, 530.8912, 530.3141, 529.4526, & 528.591, 527.7293, 526.8678, 526.0063, 525.1448, 524.2833, 523.4216, 522.5601, 521.6986, 520.837, 519.9755, 519.114, 518.2524, 517.3909, & 516.5293, 515.6678, 514.8062, 513.9447, 513.2912, 512.8458, 512.4003, 511.9549, 511.5095, 511.064, 510.6187, 510.1732, 509.7278, 509.2823, & 508.837, 508.3916, 507.9461, 507.5007, 507.0553, 506.6098, 506.1645, 505.719, 505.2736, 504.8281, 504.3828, 503.9373, 503.4919, 503.0464, & 502.6011, 502.1556, 501.7101, 501.2648, 500.8193, 500.374, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, & 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512, 500.1512/ data h2_lat(1:180,3) /544.8172, 544.8079, 544.7986, 544.7896, 544.7803, 544.7711, 544.7619, 544.7527, 544.7435, 544.7343, 544.7252, 544.7159, 544.7068, 544.6976, & 544.6884, 544.6791, 544.67, 544.6608, 544.6516, 544.6424, 544.6332, 544.624, 544.6147, 544.6057, 544.5964, 544.5524, 544.4738, 544.3951, & 544.3164, 544.2377, 544.1591, 544.0804, 544.0016, 543.9229, 543.8442, 543.7656, 543.6869, 543.6084, 543.5296, 543.4508, 543.3723, 543.2935, & 543.2148, 543.1362, 543.0575, 542.9788, 542.9001, 542.8214, 542.7427, 542.6641, 542.5853, 542.5066, 542.4279, 542.3493, 542.2706, 542.1919, & 542.1133, 542.0346, 541.9559, 541.8771, 541.7985, 541.7198, 541.6411, 541.5624, 541.4837, 541.4051, 541.3264, 541.2477, 541.1691, 541.0903, & 541.0117, 540.933, 540.8542, 540.7755, 540.6968, 540.6182, 540.4709, 540.2548, 540.0388, 539.8229, 539.6068, 539.3907, 539.1746, 538.9587, & 538.7426, 538.5266, 538.3106, 538.0945, 537.8785, 537.6625, 537.4464, 537.2303, 537.0143, 536.7983, 536.5823, 536.3661, 536.1502, 535.9341, & 535.7181, 535.5021, 535.2861, 535.0701, 534.8539, 534.6379, 534.422, 534.2059, 533.9899, 533.7739, 533.5577, 533.3418, 532.848, 532.0765, & 531.3049, 530.5332, 529.7617, 528.9902, 528.2187, 527.4471, 526.6755, 525.904, 525.1324, 524.3609, 523.5894, 522.8177, 522.0462, 521.2747, & 520.5031, 519.7316, 518.96, 518.1885, 517.6494, 517.3428, 517.0361, 516.7293, 516.4227, 516.116, 515.8094, 515.5027, 515.196, 514.8893, & 514.5828, 514.2761, 513.9694, 513.6627, 513.3561, 513.0494, 512.7428, 512.4361, 512.1294, 511.8227, 511.5161, 511.2094, 510.9027, 510.5961, & 510.2895, 509.9828, 509.676, 509.3695, 509.0627, 508.7562, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, & 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028, 508.6028/ data h2_lat(1:180,4) /541.3416, 541.2589, 541.1764, 541.094, 541.0115, 540.9289, 540.8464, 540.764, 540.6815, 540.5989, 540.5165, 540.434, 540.3515, 540.269, & 540.1865, 540.1039, 540.0216, 539.939, 539.8565, 539.7741, 539.6915, 539.609, 539.5266, 539.4441, 539.3615, 539.3212, 539.3232, 539.325, & 539.3269, 539.3288, 539.3307, 539.3326, 539.3344, 539.3362, 539.3381, 539.34, 539.3419, 539.3439, 539.3456, 539.3474, 539.3494, 539.3513, & 539.3531, 539.3551, 539.3569, 539.3588, 539.3607, 539.3625, 539.3644, 539.3663, 539.3681, 539.3699, 539.3718, 539.3738, 539.3756, 539.3775, & 539.3794, 539.3813, 539.3832, 539.3849, 539.3869, 539.3887, 539.3906, 539.3925, 539.3943, 539.3962, 539.3981, 539.4, 539.402, 539.4037, & 539.4056, 539.4075, 539.4094, 539.4111, 539.4131, 539.415, 539.356, 539.2363, 539.1165, 538.9967, 538.8769, 538.757, 538.6372, 538.5175, & 538.3976, 538.2778, 538.1581, 538.0382, 537.9184, 537.7986, 537.6788, 537.559, 537.4392, 537.3194, 537.1995, 537.0797, 536.96, 536.8401, & 536.7203, 536.6005, 536.4808, 536.361, 536.241, 536.1213, 536.0016, 535.8817, 535.7619, 535.6422, 535.5222, 535.4025, 535.0604, 534.4958, & 533.9312, 533.3666, 532.802, 532.2375, 531.673, 531.1085, 530.5438, 529.9793, 529.4148, 528.8502, 528.2856, 527.7211, 527.1566, 526.592, & 526.0274, 525.463, 524.8983, 524.3338, 523.8645, 523.4905, 523.1164, 522.7423, 522.3683, 521.9942, 521.6202, 521.2462, 520.8721, 520.4981,& 520.1241, 519.7501, 519.376, 519.002, 518.6279, 518.2538, 517.8799, 517.5058, 517.1318, 516.7576, 516.3837, 516.0096, 515.6356, 515.2616, & 514.8876, 514.5135, 514.1394, 513.7655, 513.3914, 513.0173, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, & 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303, 512.8303/ data h2_lat(1:180,5) /534.4771, 534.4755, 534.4741, 534.4728, 534.4714, 534.47, 534.4686, 534.4672, 534.4658, 534.4644, 534.4631, 534.4616, 534.4603, 534.4589, & 534.4575, 534.4561, 534.4548, 534.4533, 534.4519, 534.4506, 534.4491, 534.4478, 534.4464, 534.445, 534.4435, 534.4635, 534.5049, 534.5462, & 534.5875, 534.6289, 534.6702, 534.7116, 534.7528, 534.7941, 534.8355, 534.8769, 534.9182, 534.9597, 535.0009, 535.0421, 535.0836, 535.1249, & 535.1663, 535.2076, 535.249, 535.2902, 535.3316, 535.3729, 535.4142, 535.4556, 535.4969, 535.5382, 535.5795, 535.621, 535.6622, 535.7036, & 535.7449, 535.7863, 535.8276, 535.8688, 535.9103, 535.9515, 535.9929, 536.0342, 536.0756, 536.1169, 536.1582, 536.1996, 536.241, 536.2822, & 536.3236, 536.365, 536.4062, 536.4475, 536.4889, 536.5302, 536.5631, 536.5875, 536.6119, 536.6363, 536.6607, 536.6849, 536.7094, 536.7338, & 536.7581, 536.7825, 536.8069, 536.8312, 536.8556, 536.8799, 536.9044, 536.9286, 536.9531, 536.9775, 537.0018, 537.0261, 537.0505, 537.0749, & 537.0993, 537.1236, 537.148, 537.1725, 537.1967, 537.2211, 537.2455, 537.2698, 537.2943, 537.3187, 537.3429, 537.3674, 537.1082, 536.5656, & 536.0229, 535.4801, 534.9374, 534.3949, 533.8522, 533.3096, 532.7668, 532.2242, 531.6816, 531.1389, 530.5963, 530.0536, 529.5109, 528.9683, & 528.4256, 527.883, 527.3403, 526.7977, 526.3599, 526.0271, 525.6943, 525.3614, 525.0286, 524.6957, 524.363, 524.0302, 523.6973, 523.3645, & 523.0318, 522.699, 522.3661, 522.0333, 521.7005, 521.3676, 521.035, 520.702, 520.3693, 520.0363, 519.7036, 519.3708, 519.0379, 518.7051, & 518.3724, 518.0396, 517.7066, 517.374, 517.0411, 516.7083, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, & 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418, 516.5418/ data h2_lat(1:180,6) /528.9905, 529.0618, 529.1333, 529.2048, 529.2762, 529.3477, 529.4192, 529.4907, 529.5621, 529.6336, 529.7051, 529.7766, 529.848, 529.9195, & 529.991, 530.0624, 530.134, 530.2054, 530.2769, 530.3484, 530.4197, 530.4912, 530.5627, 530.6342, 530.7056, 530.7879, 530.881, 530.9739, & 531.067, 531.16, 531.2531, 531.3461, 531.4391, 531.532, 531.6252, 531.7182, 531.8112, 531.9044, 531.9973, 532.0903, 532.1834, 532.2764, & 532.3694, 532.4626, 532.5555, 532.6485, 532.7416, 532.8346, 532.9277, 533.0207, 533.1137, 533.2067, 533.2997, 533.3928, 533.4858, 533.5789, & 533.6719, 533.765, 533.858, 533.9509, 534.044, 534.137, 534.23, 534.3231, 534.4161, 534.5092, 534.6022, 534.6952, 534.7884, 534.8812, 534.9743,& 535.0674, 535.1603, 535.2533, 535.3464, 535.4395, 535.5298, 535.6174, 535.705, 535.7927, 535.8802, 535.9678, 536.0554, 536.1431, 536.2306, & 536.3182, 536.4058, 536.4934, 536.581, 536.6686, 536.7562, 536.8438, 536.9314, 537.019, 537.1066, 537.1941, 537.2818, 537.3694, 537.457, & 537.5446, 537.6322, 537.7198, 537.8073, 537.895, 537.9826, 538.0702, 538.1578, 538.2455, 538.3329, 538.4206, 538.2084, 537.6962, 537.184, & 536.6717, 536.1595, 535.6475, 535.1353, 534.6232, 534.111, 533.5988, 533.0867, 532.5745, 532.0623, 531.5502, 531.038, 530.5258, 530.0137, & 529.5015, 528.9893, 528.4772, 527.9321, 527.3542, 526.7761, 526.1981, 525.6202, 525.0421, 524.4642, 523.8862, 523.3082, 522.7302, 522.1523, & 521.5743, 520.9963, 520.4182, 519.8403, 519.2622, 518.6844, 518.1063, 517.5283, 516.9502, 516.3724, 515.7943, 515.2163, 514.6384, 514.0604, & 513.4824, 512.9043, 512.3265, 511.7484, 511.1704, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, & 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814, 510.8814/ data h2_lat(1:180,7) /528.1784, 528.1115, 528.0446, 527.978, 527.9111, 527.8442, 527.7775, 527.7108, 527.6439, 527.5771, 527.5104, 527.4435, 527.3768, 527.3101, & 527.2432, 527.1764, 527.1097, 527.0428, 526.9761, 526.9094, 526.8425, 526.7757, 526.7089, 526.6422, 526.5753, 526.6195, 526.7747, 526.9299, & 527.085, 527.2402, 527.3954, 527.5505, 527.7056, 527.8607, 528.0159, 528.1711, 528.3263, 528.4816, 528.6366, 528.7917, 528.947, 529.1021, & 529.2572, 529.4125, 529.5676, 529.7227, 529.8779, 530.033, 530.1882, 530.3434, 530.4985, 530.6536, 530.8088, 530.964, 531.1191, 531.2743, & 531.4294, 531.5846, 531.7398, 531.8948, 532.05, 532.2052, 532.3604, 532.5155, 532.6707, 532.8259, 532.981, 533.1362, 533.2914, 533.4464, & 533.6017, 533.7568, 533.9119, 534.067, 534.2222, 534.3774, 534.4937, 534.5709, 534.6481, 534.7255, 534.8027, 534.8798, 534.9571, 535.0344, & 535.1116, 535.1888, 535.2661, 535.3433, 535.4206, 535.4978, 535.5751, 535.6523, 535.7296, 535.8068, 535.884, 535.9612, 536.0385, 536.1158, & 536.193, 536.2702, 536.3475, 536.4248, 536.5019, 536.5792, 536.6565, 536.7337, 536.811, 536.8883, 536.9654, 537.0427, 536.8989, 536.5341, & 536.1691, 535.8041, 535.4393, 535.0745, 534.7097, 534.3448, 533.9799, 533.6151, 533.2501, 532.8853, 532.5204, 532.1555, 531.7906, 531.4258, & 531.0609, 530.6961, 530.3311, 529.9663, 529.149, 527.8792, 526.6093, 525.3394, 524.0696, 522.7997, 521.53, 520.2601, 518.9903, 517.7205, & 516.4507, 515.1809, 513.9109, 512.6411, 511.3713, 510.1015, 508.8318, 507.5618, 506.292, 505.0221, 503.7524, 502.4825, 501.2126, 499.9428, & 498.6731, 497.4032, 496.1333, 494.8636, 493.5937, 492.3239, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, & 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889, 491.6889/ data h2_lat(1:180,8) /528.6539, 528.5361, 528.4184, 528.3008, 528.183, 528.0653, 527.9477, 527.83, 527.7123, 527.5946, 527.477, 527.3592, 527.2416, 527.1239, & 527.0062, 526.8884, 526.7709, 526.6531, 526.5355, 526.4178, 526.3, 526.1824, 526.0647, 525.947, 525.8293, 525.834, 525.9611, 526.088, & 526.215, 526.342, 526.4691, 526.5961, 526.723, 526.85, 526.9771, 527.1041, 527.2311, 527.3583, 527.4852, 527.6121, 527.7392, 527.8661, & 527.9932, 528.1202, 528.2473, 528.3742, 528.5012, 528.6282, 528.7552, 528.8823, 529.0093, 529.1362, 529.2632, 529.3903, 529.5173, 529.6443, & 529.7714, 529.8984, 530.0254, 530.1523, 530.2794, 530.4063, 530.5334, 530.6604, 530.7874, 530.9144, 531.0414, 531.1684, 531.2955, 531.4224, & 531.5495, 531.6765, 531.8035, 531.9304, 532.0574, 532.1845, 532.2175, 532.1564, 532.0953, 532.0342, 531.9731, 531.9119, 531.8508, 531.7898, & 531.7286, 531.6674, 531.6064, 531.5452, 531.4841, 531.423, 531.3619, 531.3007, 531.2397, 531.1786, 531.1174, 531.0562, 530.9952, 530.9341, & 530.873, 530.8118, 530.7507, 530.6897, 530.6284, 530.5674, 530.5063, 530.4451, 530.384, 530.3231, 530.2618, 530.2007, 529.99, 529.6296, & 529.2692, 528.9088, 528.5483, 528.1881, 527.8278, 527.4673, 527.1069, 526.7466, 526.3862, 526.0258, 525.6655, 525.3051, 524.9447, 524.5844, & 524.2239, 523.8636, 523.5032, 523.1429, 522.1282, 520.4591, 518.79, 517.1209, 515.4518, 513.7827, 512.1137, 510.4446, 508.7755, 507.1065, & 505.4375, 503.7684, 502.0993, 500.4302, 498.7611, 497.092, 495.423, 493.7539, 492.0848, 490.4157, 488.7467, 487.0776, 485.4085, 483.7394, & 482.0704, 480.4013, 478.7322, 477.0632, 475.394, 473.725, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, & 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905, 472.8905/ data h2_lat(1:180,9) /527.6614, 527.6564, 527.6516, 527.6469, 527.6419, 527.637, 527.6323, 527.6274, 527.6226, 527.6178, 527.613, 527.6081, 527.6033, 527.5985, & 527.5937, 527.5888, 527.5841, 527.5792, 527.5743, 527.5696, 527.5647, 527.5599, 527.5551, 527.5503, 527.5454, 527.606, 527.7322, 527.8583, & 527.9844, 528.1105, 528.2367, 528.3629, 528.4888, 528.6149, 528.7411, 528.8672, 528.9933, 529.1196, 529.2456, 529.3716, 529.4979, 529.624, & 529.7501, 529.8763, 530.0024, 530.1284, 530.2546, 530.3807, 530.5068, 530.6329, 530.759, 530.8851, 531.0112, 531.1375, 531.2635, 531.3896, & 531.5157, 531.6419, 531.768, 531.894, 532.0202, 532.1463, 532.2725, 532.3986, 532.5247, 532.6508, 532.7769, 532.903, 533.0293, 533.1552, & 533.2814, 533.4076, 533.5336, 533.6597, 533.7858, 533.912, 533.8359, 533.5576, 533.2794, 533.0011, 532.7228, 532.4444, 532.1661, 531.8879, & 531.6096, 531.3312, 531.053, 530.7747, 530.4964, 530.2181, 529.9398, 529.6615, 529.3832, 529.105, 528.8266, 528.5483, 528.27, 527.9918, & 527.7134, 527.4351, 527.1569, 526.8787, 526.6002, 526.322, 526.0437, 525.7654, 525.4871, 525.2089, 524.9305, 524.6522, 524.1245, 523.3474, & 522.5701, 521.7928, 521.0157, 520.2386, 519.4614, 518.6842, 517.907, 517.1299, 516.3527, 515.5755, 514.7983, 514.0211, 513.244, 512.4667, & 511.6896, 510.9124, 510.1352, 509.3581, 508.2245, 506.7346, 505.2447, 503.7547, 502.2649, 500.7749, 499.285, 497.7951, 496.3052, 494.8153, & 493.3254, 491.8355, 490.3456, 488.8556, 487.3657, 485.8758, 484.386, 482.896, 481.4061, 479.916, 478.4262, 476.9363, 475.4463, 473.9565, & 472.4666, 470.9767, 469.4866, 467.9969, 466.5069, 465.017, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, & 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272, 464.272/ data h2_lat(1:180,10) / 528.5276, 528.5571, 528.5867, 528.6165, 528.646, 528.6756, 528.7053, 528.7349, 528.7645, 528.7942, 528.8239, 528.8534, 528.8831, 528.9127, & 528.9424, 528.9719, 529.0017, 529.0312, 529.0609, 529.0906, 529.1202, 529.1498, 529.1794, 529.2091, 529.2386, 529.3627, 529.5812, 529.7996, & 530.0179, 530.2364, 530.4548, 530.6732, 530.8915, 531.1099, 531.3283, 531.5468, 531.7651, 531.9837, 532.202, 532.4203, 532.6388, 532.8572, & 533.0756, 533.2941, 533.5125, 533.7308, 533.9492, 534.1676, 534.386, 534.6044, 534.8228, 535.0411, 535.2596, 535.478, 535.6964, 535.9148, & 536.1332, 536.3517, 536.5701, 536.7884, 537.0068, 537.2252, 537.4436, 537.662, 537.8804, 538.0989, 538.3172, 538.5356, 538.7542, 538.9724, & 539.1909, 539.4093, 539.6277, 539.8459, 540.0645, 540.2828, 540.1044, 539.5292, 538.9539, 538.3786, 537.8033, 537.2279, 536.6526, 536.0775, & 535.502, 534.9268, 534.3515, 533.7762, 533.2009, 532.6255, 532.0503, 531.4749, 530.8997, 530.3244, 529.749, 529.1737, 528.5984, 528.0232, & 527.4479, 526.8725, 526.2973, 525.722, 525.1466, 524.5714, 523.9961, 523.4208, 522.8455, 522.2703, 521.6948, 521.1196, 520.3181, 519.2905, & 518.2628, 517.235, 516.2073, 515.1797, 514.152, 513.1243, 512.0966, 511.069, 510.0413, 509.0136, 507.9859, 506.9582, 505.9305, 504.9028, & 503.8752, 502.8475, 501.8198, 500.7921, 499.6252, 498.3189, 497.0126, 495.7062, 494.4, 493.0936, 491.7874, 490.4811, 489.1748, 487.8685, & 486.5623, 485.256, 483.9496, 482.6433, 481.3371, 480.0308, 478.7246, 477.4182, 476.1119, 474.8055, 473.4994, 472.1931, 470.8867, 469.5804, & 468.2742, 466.9679, 465.6615, 464.3554, 463.049, 461.7428, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, & 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896, 461.0896/ data h2_lat(1:180,11) /533.6163, 533.6492, 533.6822, 533.7154, 533.7484, 533.7814, 533.8145, 533.8475, 533.8806, 533.9136, 533.9468, 533.9797, 534.0129, 534.0459, & 534.079, 534.112, 534.1451, 534.1782, 534.2112, 534.2443, 534.2773, 534.3104, 534.3434, 534.3765, 534.4095, 534.5023, 534.655, 534.8076, & 534.9602, 535.1129, 535.2654, 535.4181, 535.5706, 535.7231, 535.8759, 536.0284, 536.1811, 536.3338, 536.4863, 536.6388, 536.7916, 536.9441, & 537.0967, 537.2495, 537.402, 537.5546, 537.7072, 537.8598, 538.0124, 538.165, 538.3176, 538.4702, 538.6228, 538.7755, 538.928, 539.0806, & 539.2333, 539.3859, 539.5385, 539.691, 539.8437, 539.9963, 540.1489, 540.3015, 540.454, 540.6068, 540.7593, 540.9119, 541.0647, 541.2171, & 541.3698, 541.5224, 541.675, 541.8275, 541.9802, 542.1328, 541.882, 541.2277, 540.5734, 539.9191, 539.2648, 538.6104, 537.9561, 537.3019, & 536.6475, 535.9932, 535.339, 534.6846, 534.0303, 533.376, 532.7217, 532.0674, 531.4131, 530.7589, 530.1045, 529.4501, 528.7959, 528.1416, & 527.4873, 526.8329, 526.1787, 525.5245, 524.87, 524.2158, 523.5615, 522.9072, 522.2529, 521.5987, 520.9442, 520.29, 519.4531, 518.4333, & 517.4136, 516.3939, 515.3741, 514.3546, 513.3349, 512.3152, 511.2954, 510.2758, 509.2561, 508.2365, 507.2168, 506.197, 505.1773, 504.1577, & 503.138, 502.1183, 501.0986, 500.0789, 499.0187, 497.9179, 496.8171, 495.7162, 494.6155, 493.5146, 492.4139, 491.313, 490.2122, 489.1115, & 488.0107, 486.91, 485.8091, 484.7083, 483.6075, 482.5067, 481.406, 480.3051, 479.2043, 478.1034, 477.0027, 475.9019, 474.8011, 473.7003, & 472.5995, 471.4987, 470.3978, 469.2972, 468.1963, 467.0956, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, & 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451, 466.5451/ data h2_lat(1:180,12) /539.8811, 539.8895, 539.8982, 539.9069, 539.9155, 539.9241, 539.9327, 539.9414, 539.95, 539.9587, 539.9673, 539.976, 539.9846, 539.9933, & 540.0019, 540.0105, 540.0192, 540.0278, 540.0364, 540.0451, 540.0537, 540.0623, 540.071, 540.0797, 540.0883, 540.1172, 540.1665, 540.2157, & 540.265, 540.3142, 540.3635, 540.4127, 540.4619, 540.5111, 540.5604, 540.6097, 540.6589, 540.7083, 540.7574, 540.8066, 540.856, 540.9052, & 540.9544, 541.0037, 541.0529, 541.1021, 541.1514, 541.2006, 541.2499, 541.2992, 541.3484, 541.3976, 541.4468, 541.4962, 541.5453, 541.5946, & 541.6439, 541.6931, 541.7424, 541.7915, 541.8408, 541.89, 541.9393, 541.9886, 542.0378, 542.0871, 542.1363, 542.1855, 542.2349, 542.284, & 542.3333, 542.3826, 542.4318, 542.4809, 542.5302, 542.5795, 542.3266, 541.7715, 541.2165, 540.6615, 540.1064, 539.5512, 538.9962, 538.4412, & 537.886, 537.331, 536.7759, 536.2208, 535.6658, 535.1107, 534.5557, 534.0005, 533.4455, 532.8904, 532.3353, 531.7802, 531.2252, 530.6701, & 530.1151, 529.5599, 529.0049, 528.4499, 527.8947, 527.3397, 526.7847, 526.2295, 525.6744, 525.1195, 524.5643, 524.0093, 523.2617, 522.3215, & 521.3813, 520.4409, 519.5007, 518.5607, 517.6205, 516.6803, 515.74, 514.7999, 513.8597, 512.9195, 511.9793, 511.0391, 510.0989, 509.1587, & 508.2185, 507.2784, 506.3381, 505.398, 504.4366, 503.4541, 502.4716, 501.489, 500.5066, 499.524, 498.5415, 497.559, 496.5765, 495.594, & 494.6116, 493.6291, 492.6465, 491.664, 490.6815, 489.699, 488.7166, 487.734, 486.7515, 485.7689, 484.7865, 483.804, 482.8214, 481.8389, & 480.8565, 479.8739, 478.8914, 477.909, 476.9264, 475.944, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, & 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527, 475.4527/ end module reaction_data