| 1 | !********************************************************************** |
|---|
| 2 | |
|---|
| 3 | ! param.cmn |
|---|
| 4 | |
|---|
| 5 | ! fgg mar 02 first version |
|---|
| 6 | !********************************************************************** |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | common/uv/ crscabsi2,c1_16,c17_24,c25_29,c30_31,c32, & |
|---|
| 10 | & c33,c34,c35,c36,jfotsout,co2crsc195,co2crsc295, & |
|---|
| 11 | & t0,fluxtop,freccen,jabsifotsintpar,e107,date_e107,e107_tab, & |
|---|
| 12 | & coefit0,coefit1,coefit2,coefit3,coefit4 |
|---|
| 13 | |
|---|
| 14 | |
|---|
| 15 | |
|---|
| 16 | real crscabsi2(nabs,16) !cross section |
|---|
| 17 | real c1_16(nz2,16) !Col. int. 1 (o2+o+h2+n) (cm^-2) |
|---|
| 18 | real c17_24(nz2) !Col. int. 17-24 (co2+o2+n2+no+co+no2) (cm^-2) |
|---|
| 19 | real c25_29(nz2) !Col. int. 25-29 (co2+o2+h2o+h2o2+no+co+no2) (cm^-2) |
|---|
| 20 | real c30_31(nz2) !Col. int. 30-31 (co2+o2+h2o+h2o2+no+no2) |
|---|
| 21 | real c32(nz2) !col. int. 32 (co2+h2o2+o2+no+no2) (cm^-2) |
|---|
| 22 | real c33(nz2) !col. int. 33 (h2o2+o2+no2) (cm^-2) |
|---|
| 23 | real c34(nz2) !col. int. 34 (h2o2+o2+o3+no2) (cm^-2) |
|---|
| 24 | real c35(nz2) !col. int. 35 (h2o2+o3+no2) (cm^-2) |
|---|
| 25 | real c36(nz2) !col. int. 36 (o3+no2) (cm^-2) |
|---|
| 26 | real jfotsout(ninter,nabs,nlayermx) |
|---|
| 27 | real co2crsc195(9) |
|---|
| 28 | real co2crsc295(9) |
|---|
| 29 | real t0(nz2) |
|---|
| 30 | real fluxtop(ninter) |
|---|
| 31 | real freccen(ninter) !representative wavelenght |
|---|
| 32 | real jabsifotsintpar(nz2,nabs,ninter) |
|---|
| 33 | real e107,date_e107(669),e107_tab(669) |
|---|
| 34 | real coefit0(ninter,nabs),coefit1(ninter,nabs) |
|---|
| 35 | real coefit2(ninter,nabs) |
|---|
| 36 | real coefit3(ninter,nabs),coefit4(ninter,nabs) |
|---|
| 37 | |
|---|
| 38 | |
|---|
| 39 | common/fotoquimica/ ch2, ch3, ch4, ch5, ch7,ch9,ch10,ch11, & |
|---|
| 40 | & ch13,ch14,ch15,ch18,ch19,ch20,ch21,ch22, & |
|---|
| 41 | & ch23,ch24,ch30,ch31,ch32,ch33,ch34,ch35, & |
|---|
| 42 | & ch36,ch37,ch38,ch39,ch40,ch41,ch42,ch43, & |
|---|
| 43 | & ch45,ch46,ch47,ch48,ch49,ch50,ch55,ch56, & |
|---|
| 44 | & ch57,ch58,ch59,ch62,ch63,ch64,ch65, & |
|---|
| 45 | & ch66,ch67,ch68,ch69,ch70,ch71,ch72, & |
|---|
| 46 | & ch73,ch74,ch75,ch76,ch85,ch86,ch87,rcoef, & |
|---|
| 47 | & jdistot,jdistot_b,jion |
|---|
| 48 | |
|---|
| 49 | real*8 ch2, ch3, ch4, ch5, ch7,ch9,ch10,ch11,ch13,ch14,ch15,ch18 |
|---|
| 50 | real*8 ch19,ch20,ch21,ch22,ch23,ch24,ch30,ch31,ch32,ch33,ch34 |
|---|
| 51 | real*8 ch35,ch36,ch37,ch38,ch39,ch40,ch41,ch42,ch43,ch45 |
|---|
| 52 | real*8 ch46,ch47,ch48,ch49,ch50,ch55,ch56,ch57,ch58,ch59,ch62 |
|---|
| 53 | real*8 ch63,ch64,ch65,ch66,ch67,ch68,ch69,ch70,ch71 |
|---|
| 54 | real*8 ch72,ch73,ch74,ch75,ch76,ch85,ch86,ch87 |
|---|
| 55 | real*8 rcoef(61,3) |
|---|
| 56 | !reaction rates |
|---|
| 57 | real jdistot(nabs,nlayermx) !photodissociation rates |
|---|
| 58 | real jdistot_b(nabs,nlayermx) |
|---|
| 59 | |
|---|
| 60 | real jion(nabs,nlayermx,4) |
|---|
| 61 | |
|---|
| 62 | |
|---|
| 63 | |
|---|
| 64 | |
|---|
| 65 | common/prod/ Pco2,Po2,Po3p,Pco,Ph,Poh,Pho2,Ph2,Ph2o,Po1d,Ph2o2, & |
|---|
| 66 | & Po3,Pn,Pno,Pno2,Pn2,Pn2d, & |
|---|
| 67 | & Pco2plus,Poplus,Po2plus,Pelect,Pcoplus,Pcplus, & |
|---|
| 68 | & Pnplus, Pnoplus, Pn2plus, Phplus, Phco2plus, & |
|---|
| 69 | & Pco2tot,Po2tot,Po3ptot,Pcotot,Phtot,Pohtot, & |
|---|
| 70 | & Pho2tot,Ph2tot,Ph2otot,Po1dtot,Ph2o2tot, & |
|---|
| 71 | & Po3tot,Pntot,Pnotot,Pno2tot,Pn2tot,Pn2dtot, & |
|---|
| 72 | & Pcoplustot,Pcplustot, & |
|---|
| 73 | & Pco2plustot,Poplustot,Po2plustot,Pelecttot, & |
|---|
| 74 | & Pnplustot,Pnoplustot,Pn2plustot,Phplustot, & |
|---|
| 75 | & Phco2plustot |
|---|
| 76 | |
|---|
| 77 | |
|---|
| 78 | real*8 Pco2(nlayermx,nreact),Po2(nlayermx,nreact) |
|---|
| 79 | real*8 Po3p(nlayermx,nreact) |
|---|
| 80 | real*8 Pco(nlayermx,nreact),Ph(nlayermx,nreact) |
|---|
| 81 | real*8 Poh(nlayermx,nreact) |
|---|
| 82 | real*8 Pho2(nlayermx,nreact),Ph2(nlayermx,nreact) |
|---|
| 83 | real*8 Ph2o(nlayermx,nreact) |
|---|
| 84 | real*8 Po1d(nlayermx,nreact),Ph2o2(nlayermx,nreact) |
|---|
| 85 | real*8 Po3(nlayermx,nreact),Pn(nlayermx,nreact) |
|---|
| 86 | real*8 Pno(nlayermx,nreact) |
|---|
| 87 | real*8 Pno2(nlayermx,nreact),Pn2(nlayermx,nreact) |
|---|
| 88 | real*8 Pn2d(nlayermx,nreact) |
|---|
| 89 | real*8 Pco2plus(nlayermx,nreact),Poplus(nlayermx,nreact) |
|---|
| 90 | real*8 Po2plus(nlayermx,nreact), Pelect(nlayermx,nreact) |
|---|
| 91 | real*8 Pcoplus(nlayermx,nreact),Pcplus(nlayermx,nreact) |
|---|
| 92 | real*8 Pnplus(nlayermx,nreact),Pnoplus(nlayermx,nreact) |
|---|
| 93 | real*8 Pn2plus(nlayermx,nreact),Phplus(nlayermx,nreact) |
|---|
| 94 | real*8 Phco2plus(nlayermx,nreact) |
|---|
| 95 | real*8 Pco2tot(nlayermx),Po2tot(nlayermx) |
|---|
| 96 | real*8 Po3ptot(nlayermx),Pcotot(nlayermx) |
|---|
| 97 | real*8 Phtot(nlayermx) |
|---|
| 98 | real*8 Pohtot(nlayermx),Pho2tot(nlayermx) |
|---|
| 99 | real*8 Ph2tot(nlayermx),Ph2otot(nlayermx) |
|---|
| 100 | real*8 Po1dtot(nlayermx),Ph2o2tot(nlayermx) |
|---|
| 101 | real*8 Po3tot(nlayermx),Pntot(nlayermx),Pnotot(nlayermx) |
|---|
| 102 | real*8 Pno2tot(nlayermx),Pn2tot(nlayermx),Pn2dtot(nlayermx) |
|---|
| 103 | real*8 Pco2plustot(nlayermx), Poplustot(nlayermx) |
|---|
| 104 | real*8 Po2plustot(nlayermx),Pelecttot(nlayermx) |
|---|
| 105 | real*8 Pcoplustot(nlayermx), Pcplustot(nlayermx) |
|---|
| 106 | real*8 Pnplustot(nlayermx), Pnoplustot(nlayermx) |
|---|
| 107 | real*8 Pn2plustot(nlayermx),Phplustot(nlayermx) |
|---|
| 108 | real*8 Phco2plustot(nlayermx) |
|---|
| 109 | |
|---|
| 110 | |
|---|
| 111 | common/loss/ Lco2,Lo2,Lo3p,Lco,Lh,Loh,Lho2,Lh2,Lh2o,Lo1d,Lh2o2, & |
|---|
| 112 | & Lo3,Ln,Lno,Lno2,Ln2,Ln2d,Lelect, & |
|---|
| 113 | & Lco2plus,Lo2plus,Loplus, Lcoplus, Lcplus, & |
|---|
| 114 | & Lnplus,Lnoplus,Ln2plus, Lhplus, Lhco2plus, & |
|---|
| 115 | & Lco2tot,Lo2tot,Lo3ptot,Lcotot,Lhtot,Lohtot, & |
|---|
| 116 | & Lho2tot,Lh2tot,Lh2otot,Lo1dtot,Lh2o2tot, & |
|---|
| 117 | & Lo3tot,Lntot,Lnotot,Lno2tot,Ln2tot,Ln2dtot, & |
|---|
| 118 | & Lcoplustot,Lcplustot,Lco2plustot,Loplustot, & |
|---|
| 119 | & Lo2plustot,Lelecttot,Lnplustot,Lnoplustot, & |
|---|
| 120 | & Ln2plustot,Lhplustot,Lhco2plustot |
|---|
| 121 | |
|---|
| 122 | |
|---|
| 123 | real*8 Lco2(nlayermx,nreact),Lo2(nlayermx,nreact) |
|---|
| 124 | real*8 Lo3p(nlayermx,nreact) |
|---|
| 125 | real*8 Lco(nlayermx,nreact),Lh(nlayermx,nreact) |
|---|
| 126 | real*8 Loh(nlayermx,nreact) |
|---|
| 127 | real*8 Lho2(nlayermx,nreact),Lh2(nlayermx,nreact) |
|---|
| 128 | real*8 Lh2o(nlayermx,nreact) |
|---|
| 129 | real*8 Lo1d(nlayermx,nreact),Lh2o2(nlayermx,nreact) |
|---|
| 130 | real*8 Lo3(nlayermx,nreact),Ln(nlayermx,nreact) |
|---|
| 131 | real*8 Lno(nlayermx,nreact) |
|---|
| 132 | real*8 Lno2(nlayermx,nreact),Ln2(nlayermx,nreact) |
|---|
| 133 | real*8 Ln2d(nlayermx,nreact) |
|---|
| 134 | real*8 Lco2plus(nlayermx,nreact), Loplus(nlayermx,nreact) |
|---|
| 135 | real*8 Lo2plus(nlayermx,nreact), Lelect(nlayermx,nreact) |
|---|
| 136 | real*8 Lcoplus(nlayermx,nreact), Lcplus(nlayermx,nreact) |
|---|
| 137 | real*8 Lnplus(nlayermx,nreact), Lnoplus(nlayermx,nreact) |
|---|
| 138 | real*8 Ln2plus(nlayermx,nreact), Lhplus(nlayermx,nreact) |
|---|
| 139 | real*8 Lhco2plus(nlayermx,nreact) |
|---|
| 140 | real*8 Lco2tot(nlayermx),Lo2tot(nlayermx) |
|---|
| 141 | real*8 Lo3ptot(nlayermx),Lcotot(nlayermx) |
|---|
| 142 | real*8 Lhtot(nlayermx),Lohtot(nlayermx) |
|---|
| 143 | real*8 Lho2tot(nlayermx),Lh2tot(nlayermx) |
|---|
| 144 | real*8 Lh2otot(nlayermx),Lo1dtot(nlayermx) |
|---|
| 145 | real*8 Lh2o2tot(nlayermx) |
|---|
| 146 | real*8 Lo3tot(nlayermx),Lntot(nlayermx),Lnotot(nlayermx) |
|---|
| 147 | real*8 Lno2tot(nlayermx),Ln2tot(nlayermx) |
|---|
| 148 | real*8 Ln2dtot(nlayermx) |
|---|
| 149 | real*8 Lco2plustot(nlayermx), Loplustot(nlayermx) |
|---|
| 150 | real*8 Lo2plustot(nlayermx), Lelecttot(nlayermx) |
|---|
| 151 | real*8 Lcoplustot(nlayermx), Lcplustot(nlayermx) |
|---|
| 152 | real*8 Lnplustot(nlayermx), Lnoplustot(nlayermx) |
|---|
| 153 | real*8 Ln2plustot(nlayermx), Lhplustot(nlayermx) |
|---|
| 154 | real*8 Lhco2plustot(nlayermx) |
|---|
| 155 | |
|---|
| 156 | |
|---|
| 157 | |
|---|
| 158 | common/vida/ tminco2,tmino2,tmino3p,tminco,tminh,tminoh, & |
|---|
| 159 | & tminho2,tminh2,tminh2o,tmino1d,tminh2o2,tmino3,tminn, & |
|---|
| 160 | & tminno,tminno2,tminn2,tminn2d, & |
|---|
| 161 | & tminco2plus,tminoplus,tmino2plus,tmincoplus, & |
|---|
| 162 | & tmincplus,tminnplus,tminnoplus, & |
|---|
| 163 | & tminn2plus,tminhplus,tminhco2plus |
|---|
| 164 | |
|---|
| 165 | |
|---|
| 166 | |
|---|
| 167 | real*8 tminco2(nlayermx),tmino2(nlayermx) |
|---|
| 168 | real*8 tmino3p(nlayermx),tminco(nlayermx) |
|---|
| 169 | real*8 tminh(nlayermx),tminoh(nlayermx) |
|---|
| 170 | real*8 tminho2(nlayermx),tminh2(nlayermx) |
|---|
| 171 | real*8 tminh2o(nlayermx),tmino1d(nlayermx),tminh2o2(nlayermx) |
|---|
| 172 | real*8 tmino3(nlayermx),tminn(nlayermx),tminno(nlayermx) |
|---|
| 173 | real*8 tminno2(nlayermx),tminn2(nlayermx),tminn2d(nlayermx) |
|---|
| 174 | real*8 tminco2plus(nlayermx),tminoplus(nlayermx) |
|---|
| 175 | real*8 tmino2plus(nlayermx) |
|---|
| 176 | real*8 tmincoplus(nlayermx) |
|---|
| 177 | real*8 tmincplus(nlayermx) |
|---|
| 178 | real*8 tminnplus(nlayermx),tminnoplus(nlayermx) |
|---|
| 179 | real*8 tminn2plus(nlayermx) |
|---|
| 180 | real*8 tminhplus(nlayermx) |
|---|
| 181 | real*8 tminhco2plus(nlayermx) |
|---|
| 182 | |
|---|
| 183 | |
|---|
| 184 | common/flujoctes/ fluxtophr,ct1,p1,ct2,p2 |
|---|
| 185 | |
|---|
| 186 | real fluxtophr(ninter) |
|---|
| 187 | real ct1(ninter),p1(ninter),ct2(ninter),p2(ninter) |
|---|
| 188 | |
|---|
| 189 | common/phdisef/ efdisco2, efdiso2, efdish2, efdish2o, & |
|---|
| 190 | & efdish2o2,efdiso3,efdiso,efdisn,efdish,efdisno,efdisn2, & |
|---|
| 191 | & efdisno2,efdisco,efionco2,efionn2,efionco, & |
|---|
| 192 | & efiono3p,efionn,efionno,efionh |
|---|
| 193 | |
|---|
| 194 | |
|---|
| 195 | real efdisco2(ninter), efdiso2(ninter), efdish2o(ninter) |
|---|
| 196 | real efdish2o2(ninter), efdish2(ninter), efdiso3(ninter) |
|---|
| 197 | real efdiso(ninter), efdisn(ninter), efdish(ninter) |
|---|
| 198 | real efdisno(ninter), efdisn2(ninter), efdisno2(ninter) |
|---|
| 199 | real efdisco(ninter) |
|---|
| 200 | real efionco2(ninter,4) |
|---|
| 201 | real efionn2(ninter,2) |
|---|
| 202 | real efionco(ninter,3) |
|---|
| 203 | real efiono3p(ninter),efionn(ninter) |
|---|
| 204 | real efionno(ninter),efionh(ninter) |
|---|
| 205 | |
|---|
| 206 | |
|---|
| 207 | |
|---|
| 208 | |
|---|
| 209 | |
|---|
| 210 | |
|---|
