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