Changeset 635 for trunk/LMDZ.MARS/libf/aeronomars/param_read.F
- Timestamp:
- Apr 26, 2012, 3:22:19 PM (13 years ago)
- File:
-
- 1 edited
-
trunk/LMDZ.MARS/libf/aeronomars/param_read.F (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/aeronomars/param_read.F
r38 r635 5 5 6 6 c common variables and constants 7 8 #include "param.h" 9 #include "param_v3.h" 10 #include "datafile.h" 11 12 7 include "dimensions.h" 8 include "dimphys.h" 9 include 'param.h' 10 include 'param_v4.h' 11 include 'datafile.h' 12 13 13 14 c local variables 14 15 15 16 integer i,j,k,inter !indexes 16 integer ierr 17 integer ierr,lnblnk 17 18 real nada 18 logical,save :: firstcall=.true. 19 20 19 20 21 21 c*************************+PROGRAM STARTS************************** 22 22 23 23 24 c data 25 if (firstcall) then 26 anchint(:)=(/4.9,5.3,6.1,12.8,2.0,6.9,11.8,0.4,1.9,4.9,0.9,4.0, 27 $6.9,4.9,4.9,10.2,5.6,1.9,0.9,0.9,0.9,0.9,3.7 28 $,13.1,3.9,6.0,4.8,23.9,8.4,8.5,26.8,7.5,127.6/) 29 30 crscabsi2(1,1:ninter2)=(/0.0,4.42e-18,7.51e-18 31 $,1.498e-17,2.344e-17,2.388e-17,2.927E-17,3.161e-17,3.161e-17, 32 $3.4e-17,3.4e-17,2.588e-17,2.596e-17,2.248e-17,3.183e-17, 33 $1.284e-17/) 34 35 crscabsi2(2,1:ninter2)=(/2.736e-19,1.055e-18,4.789e-18 36 $,1.063e-17,1.559e-17,1.698e-17,2.164e-17,2.241e-17,2.427e-17, 37 $2.444e-17,2.502e-17,2.516e-17,2.907e-17,3.649e-17,3.066e-17, 38 $2.093e-17/) 39 40 crscabsi2(3,1:ninter2)=(/2.603e-19,6.985e-19,3.97e-18 41 $,6.621e-18,8.433e-18,9.011e-18,9.862e-18,1.033e-17,1.012e-17, 42 $1.033e-17,1.033e-17,1.034e-17,9.0e-18,3.564e-18,3.572e-18 43 $,3.47E-18/) 44 45 freccen(:)=(/3.4,7.5,14.5,23.0,30.3,34.1,49.6,50.5,52.5,56.0, 24 c data for the UV heating tabulation 25 26 data (crscabsi2(1,j),j=1,16) /5.61031E-19,1.59677E-18,4.7072E-18, 27 $ 1.48254e-17,2.07445e-17,2.573e-17,2.901e-17,3.083e-17, 28 $ 3.217e-17,3.539e-17,3.658e-17,3.63e-17,3.41239e-17, 29 $ 2.71019e-17,4.93677e-17,1.64e-17/ 30 31 data (crscabsi2(2,j),j=1,16) /0.27250E-18,0.11650E-17,0.39250E-17, 32 $ 0.10630E-16,0.15590E-16,0.17180E-16,0.19270E-16,0.22860E-16, 33 $ 0.24270E-16,0.24440E-16,0.25020E-16,0.26600E-16,0.25400E-16, 34 $ 0.35800E-16,0.25590E-16,0.16740E-16/ 35 36 data (crscabsi2(3,j),j=1,16) /0.2776E-18,0.9792E-18,0.3313E-17, 37 $ 0.6621E-17,0.8481E-17,0.9146E-17,0.9414E-17,0.1039E-16, 38 $ 0.1012E-16,0.1033E-16,0.1033E-16,0.1033E-16,0.8268E-17, 39 $ 0.6563E-17,0.3506E-17,0.3470E-17/ 40 41 data (crscabsi2(5,j),j=1,16) /.5E-20,.1077607E-19,.5670491E-19, 42 $ .3322716E-18,.1054509E-17,.1700005E-17,.3171188E-17, 43 $ .4734241E-17,.5108741E-17,.6022236E-17,.6741537E-17, 44 $ .7277079E-17,.9070787E-17,.9708916E-17,.4026281E-17,0.0/ 45 46 data (crscabsi2(8,j),j=1,16) /0.0, 7.44175e-19, 2.23167e-18, 47 $ 8.46200e-18,1.18275e-17,1.54900e-17,2.32475e-17,2.41373e-17, 48 $ 2.55482e-17,2.38431e-17,2.28600e-17,2.35067e-17,2.56000e-17, 49 $ 2.64636e-17,2.86260e-17,3.26561e-17/ 50 51 data(crscabsi2(9,j),j=1,16) /3.48182e-20,3.37038e-19,1.03077e-18, 52 $ 4.01364e-18,6.45e-18,7.8e-18,1.0e-17,1.13500e-17,1.15500e-17, 53 $ 1.18000e-17,1.17500e-17,1.16000e-17,1.28667e-17,1.18500e-17, 54 $ 1.11000e-17,9.50000e-18/ 55 56 data(crscabsi2(10,j),j=1,16) /0.0,9.39833e-19,2.87714e-18, 57 $ 9.66900e-18,1.37063e-17,1.61820e-17,2.30450e-17,2.63373e-17, 58 $ 2.63773e-17,2.67677e-17,2.64100e-17,2.53000e-17,2.18100e-17, 59 $ 2.04941e-17,2.28160e-17,2.93550e-17/ 60 61 data(crscabsi2(11,j),j=1,16) /0.0,9.58555e-19,2.52767e-18, 62 $ 8.29700e-18,1.21850e-17,1.40500e-17,1.97025e-17,2.12018e-17, 63 $ 2.14673e-17,2.20331e-17,2.21500e-17,2.21600e-17,2.33200e-17, 64 $ 2.67800e-17,2.56400e-17,3.58561e-17/ 65 66 data(crscabsi2(12,j),j=1,16) /0.0,1.0e-20,2.5e-20,1.30400e-19, 67 $ 2.93800e-19,4.36000e-19,8.49400e-19,1.29400e-18,1.40500e-18, 68 $ 1.67600e-18,1.93400e-18,2.12200e-18,2.75800e-18,3.48400e-18, 69 $ 4.17200e-18,5.26000e-18/ 70 71 data(crscabsi2(13,j),j=1,16) /0.0,1.60e-18,4.99111e-18,1.48496e-17 72 $ ,2.17395e-17,2.55857e-17,2.87754e-17,3.65571e-17,3.85691e-17, 73 $ 4.16286e-17,4.15117e-17,4.05901e-17,3.64000e-17,2.99670e-17, 74 $ 2.46796e-17,2.51789e-17/ 75 76 data freccen /3.4,7.5,14.5,23.0,30.3,34.1,49.6,50.5,52.5,56.0, 46 77 $59.0,61.5,68.7,73.1,78.4,83.1,92.4,97.5,99.3,100.1,100.7,102.1, 47 78 $104.5,116.8,121.3,127.0,130.6,153.7,162.8,171.4 48 $,195.6,206.3,273.5/) 49 50 co2crsc195(:)=(/3.691e-19,4.44216e-20,3.86945e-19,5.94208e-19, 51 $2.93217e-19,7.58769e-20,8.60192e-21,4.20007e-24,2.29996e-26/) 52 53 co2crsc295(:)=(/3.691e-19,5.21572e-20,4.23488e-19,6.54728e-19, 54 $3.30227e-19,1.03183e-19,1.55722e-20,1.72317e-23,7.0e-25/) 55 56 firstcall=.false. 57 endif ! of if (firstcall) 79 $,195.6,206.3,222.0,236.0,289.0,600./ 80 81 data co2crsc195/2.05864e-17,5.90557e-20,3.1027e-19,6.70653e-19, 82 $4.55132e-19,8.87122e-20,1.32138e-20,7.22244e-23,2.88002e-26/ 83 84 data co2crsc295/2.05897e-17,6.71104e-20,3.45509e-19,7.45711e-19, 85 $4.82752e-19,1.11594e-19,1.98308e-20,1.3853e-22,2.1414e-25/ 58 86 59 87 c Reads tabulated functions 60 88 61 open(210,status='old',file=trim(datafile)//'/EUVDAT/coln.dat', 62 & iostat=ierr) 63 64 IF (ierr.NE.0) THEN 89 !Tabulated column amount 90 open(210, status = 'old', 91 c $file=datafile(1:lnblnk(datafile))//'/EUVDAT/coln.dat',iostat=ierr) 92 $file=datafile 93 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/coln.dat',iostat=ierr) 94 95 IF (ierr.NE.0) THEN 65 96 write(*,*)'cant find directory EUVDAT and content coln.dat' 66 write(*,*)'(in aeronomars/param_read .F)'67 write(*,*)'It should be in :', trim(datafile),'/'97 write(*,*)'(in aeronomars/param_read_iono.F)' 98 write(*,*)'It should be in :', datafile,'/' 68 99 write(*,*)'1) You can change this directory address in ' 69 100 write(*,*)' file phymars/datafile.h' … … 73 104 STOP 74 105 ENDIF 75 76 open(220,file=trim(datafile)//'/EUVDAT/j2_an.dat') 77 open(230,file=trim(datafile)//'/EUVDAT/j3_an.dat') 78 open(240,file=trim(datafile)//'/EUVDAT/j1_an.dat') 79 open(250,file=trim(datafile)//'/EUVDAT/j2_bn.dat') 80 open(260,file=trim(datafile)//'/EUVDAT/j2_cn.dat') 81 open(270,file=trim(datafile)//'/EUVDAT/j1_bn.dat') 82 open(280,file=trim(datafile)//'/EUVDAT/j1_cn.dat') 83 open(290,file=trim(datafile)//'/EUVDAT/j1_dn.dat') 84 open(150,file=trim(datafile)//'/EUVDAT/j4n.dat') 85 open(160,file=trim(datafile)//'/EUVDAT/j5n.dat') 86 open(170,file=trim(datafile)//'/EUVDAT/j6n.dat') 87 88 do i=210,290,10 89 read(i,*) 90 read(i,*) 91 end do 92 93 do i=150,170,10 94 read(i,*) 95 read(i,*) 96 end do 97 98 do i=nz2,1,-1 99 read(210,*) c23(i),(c123(j,i),j=2,ninter2),c12(i),c1(i),ch2o2(i) 100 end do 101 102 do i=nz2,1,-1 103 read(220,*) (jabsifotsint(j,2,i),j=1,ninter2) 104 end do 105 106 do i=nz2,1,-1 107 read(230,*) (jabsifotsint(j,3,i),j=1,ninter2) 108 end do 109 110 do i=nz2,1,-1 111 read(240,*) (jabsifotsint(j,1,i),j=2,ninter2) 112 end do 113 114 do i=nz2,1,-1 115 read(250,*) (jabsifotsint(j,2,i),j=17,24) 116 end do 117 118 do i=nz2,1,-1 119 read(260,*) (jabsifotsint(j,2,i),j=25,31) 120 end do 121 122 do i=nz2,1,-1 123 read(270,*) (jabsifotsint(j,1,i),j=17,24) 124 end do 125 126 do i=nz2,1,-1 127 read(280,*) (jabsifotsint(j,1,i),j=25,31) 128 end do 129 130 do i=nz2,1,-1 131 read(290,*) jabsifotsint(32,1,i) 132 end do 133 134 do i=nz2,1,-1 135 read(160,*) (jabsifotsint(j,5,i),j=1,15) 136 end do 137 138 do i=nz2,1,-1 139 read(150,*) (jabsifotsint(j,4,i),j=25,31) 140 end do 141 142 do i=nz2,1,-1 143 read(170,*) (jabsifotsint(j,6,i),j=25,33) 144 end do 145 146 do i=210,290,10 106 107 !Tabulated photoabsorption coefficients 108 open(220,file=datafile 109 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/j2_an.dat') 110 open(230,file=datafile 111 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/j3_an.dat') 112 open(240,file=datafile 113 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/j1_an.dat') 114 open(250,file=datafile 115 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/j2_bn.dat') 116 open(260,file=datafile 117 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/j2_cn.dat') 118 open(300,file=datafile 119 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j2_dn.dat') 120 open(270,file=datafile 121 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j1_bn.dat') 122 open(280,file=datafile 123 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j1_cn.dat') 124 open(290,file=datafile 125 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j1_dn.dat') 126 open(150,file=datafile 127 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j4n.dat') 128 open(160,file=datafile 129 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j5n.dat') 130 open(170,file=datafile 131 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j6n.dat') 132 open(180,file=datafile 133 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j7n.dat') 134 open(390,file=datafile 135 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j8_an.dat') 136 open(400,file=datafile 137 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j8_bn.dat') 138 open(410,file=datafile 139 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j9n.dat') 140 open(420,file=datafile 141 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j10_an.dat') 142 open(430,file=datafile 143 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j10_bn.dat') 144 open(440,file=datafile 145 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j10_cn.dat') 146 open(450,file=datafile 147 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j11_an.dat') 148 open(460,file=datafile 149 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j11_bn.dat') 150 open(470,file=datafile 151 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j11_cn.dat') 152 open(480,file=datafile 153 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j12n.dat') 154 open(490,file=datafile 155 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j13_an.dat') 156 open(500,file=datafile 157 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j13_bn.dat') 158 open(510,file=datafile 159 $ (1:lnblnk(datafile))//'/EUVDAT//param_v5/j13_cn.dat') 160 161 162 do i=210,300,10 163 read(i,*) 164 read(i,*) 165 end do 166 167 do i=150,180,10 168 read(i,*) 169 read(i,*) 170 end do 171 172 do i=390,510,10 173 read(i,*) 174 read(i,*) 175 enddo 176 177 do i=nz2,1,-1 178 read(210,*) (c1_16(i,j),j=1,16),c17_24(i),c25_29(i),c30_31(i), 179 $ c32(i),c33(i),c34(i),c35(i),c36(i) 180 end do 181 182 do i=nz2,1,-1 183 read(220,*) (jabsifotsintpar(j,2,i),j=1,16) 184 end do 185 186 do i=nz2,1,-1 187 read(230,*) (jabsifotsintpar(j,3,i),j=1,16) 188 end do 189 190 do i=nz2,1,-1 191 read(240,*) (jabsifotsintpar(j,1,i),j=1,16) 192 end do 193 194 do i=nz2,1,-1 195 read(250,*) (jabsifotsintpar(j,2,i),j=17,24) 196 end do 197 198 199 do i=nz2,1,-1 200 read(260,*) (jabsifotsintpar(j,2,i),j=25,31) 201 end do 202 203 do i=nz2,1,-1 204 read(270,*) (jabsifotsintpar(j,1,i),j=17,24) 205 end do 206 207 do i=nz2,1,-1 208 read(280,*) (jabsifotsintpar(j,1,i),j=25,31) 209 end do 210 211 do i=nz2,1,-1 212 read(290,*) jabsifotsintpar(32,1,i) 213 end do 214 215 do i=nz2,1,-1 216 read(300,*) (jabsifotsintpar(j,2,i),j=32,34) 217 end do 218 219 do i=nz2,1,-1 220 read(160,*) (jabsifotsintpar(j,5,i),j=1,15) 221 end do 222 223 do i=nz2,1,-1 224 read(150,*) (jabsifotsintpar(j,4,i),j=25,31) 225 end do 226 227 do i=nz2,1,-1 228 read(170,*) (jabsifotsintpar(j,6,i),j=25,35) 229 end do 230 231 do i=nz2,1,-1 232 read(180,*) (jabsifotsintpar(j,7,i),j=34,36) 233 end do 234 235 do i=nz2,1,-1 236 read(390,*) (jabsifotsintpar(j,8,i),j=2,16) 237 enddo 238 239 do i=nz2,1,-1 240 read(400,*) (jabsifotsintpar(j,8,i),j=17,24) 241 enddo 242 243 do i=nz2,1,-1 244 read(410,*) (jabsifotsintpar(j,9,i),j=1,16) 245 enddo 246 247 do i=nz2,1,-1 248 read(420,*) (jabsifotsintpar(j,10,i),j=2,16) 249 enddo 250 251 do i=nz2,1,-1 252 read(430,*) (jabsifotsintpar(j,10,i),j=17,24) 253 enddo 254 255 do i=nz2,1,-1 256 read(440,*) (jabsifotsintpar(j,10,i),j=25,32) 257 enddo 258 259 do i=nz2,1,-1 260 read(450,*) (jabsifotsintpar(j,11,i),j=2,16) 261 enddo 262 263 do i=nz2,1,-1 264 read(460,*) (jabsifotsintpar(j,11,i),j=17,24) 265 enddo 266 267 do i=nz2,1,-1 268 read(470,*) (jabsifotsintpar(j,11,i),j=25,29) 269 enddo 270 271 do i=nz2,1,-1 272 read(480,*) (jabsifotsintpar(j,12,i),j=2,16) 273 enddo 274 275 do i=nz2,1,-1 276 read(490,*) (jabsifotsintpar(j,13,i),j=2,16) 277 enddo 278 279 do i=nz2,1,-1 280 read(500,*) (jabsifotsintpar(j,13,i),j=17,24) 281 enddo 282 283 do i=nz2,1,-1 284 read(510,*) (jabsifotsintpar(j,13,i),j=25,36) 285 enddo 286 287 do i=210,300,10 147 288 close(i) 148 289 end do 149 290 150 do i=150,1 70,10291 do i=150,180,10 151 292 close(i) 152 293 end do 153 294 154 155 c reads t0 156 157 open(120,file=trim(datafile)//'/EUVDAT/t0.dat') 158 do i=1,201 159 read(120,*)t0(i) 160 end do 161 close(120) 162 163 open(100,file=trim(datafile)//'/EUVDAT/flujo.dat') 295 do i=390,510,10 296 close(i) 297 enddo 298 299 300 c set t0 301 302 do i=1,nz2 303 t0(i)=195. 304 end do 305 306 164 307 do i=1,ninter 165 read(100,*) inter,fluxtophr(i) 308 fluxtop(i)=1. 309 end do 310 311 !Parameters for the variation of the solar flux with 11 years cycle 312 open(100,file=datafile 313 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/varflujo.dat') 314 read(100,*) 315 do i=1,24 316 read(100,*) inter,ct1(i),p1(i),ct2(i),p2(i),nada 166 317 end do 167 318 close(100) 168 319 169 170 open(99,file=trim(datafile)//'/EUVDAT/varflujo.dat') 171 read(99,*) 172 do i=1,24 173 read(99,*) inter,ct1(i),p1(i),ct2(i),p2(i),nada 174 end do 175 close(99) 176 177 c eficiencias de disociacion (de Torr et al, 1979) 178 179 do inter=1,11 180 efdisco2(inter) = 0. 181 efdiso2(inter) = 0. 182 end do 183 320 c dissociation and ionization efficiencies 321 322 do inter=1,ninter 323 efdisco2(inter)=0. 324 efdiso2(inter)=0. 325 efdish2(inter)=0. 326 efdish2o(inter)=0. 327 efdish2o2(inter)=0. 328 efdiso3(inter)=0. 329 efdisco(inter)=0. 330 efdisn2(inter)=0. 331 efdisno(inter)=0. 332 efdisno2(inter)=0. 333 efionco2(inter,1)=0. 334 efionco2(inter,2)=0. 335 efionco2(inter,3)=0. 336 efionco2(inter,4)=0. 337 efiono3p(inter)=0. 338 efionn2(inter,1)=0. 339 efionn2(inter,2)=0. 340 efionco(inter,1)=0. 341 efionco(inter,2)=0. 342 efionco(inter,3)=0. 343 efionn(inter)=0. 344 efionh(inter)=0. 345 efionno(inter)=0. 346 enddo 347 348 349 c CO2, O2, NO 350 351 open(120,file=datafile 352 $ (1:lnblnk(datafile))//'/EUVDAT/param_v5/efdis_inter.dat') 353 read(120,*) 354 ! do i=1,21 355 ! read(120,*)inter,efdisco2(inter),efdiso2(inter),efdisno(inter) 356 do inter=8,28 357 read(120,*)i,efdisco2(inter),efdiso2(inter),efdisno(inter) 358 enddo 359 do inter=29,ninter 360 efdisco2(inter)=1. 361 efdiso2(inter)=1. 362 efdisno(inter)=1. 363 enddo 364 365 366 c N2 367 368 efdisn2(15)=0.1 369 do inter=16,ninter 370 efdisn2(inter)=1. 371 enddo 372 373 374 c CO 375 376 efdisco(16)=0.5 377 do inter=17,ninter 378 efdisco(inter)=1. 379 enddo 380 381 382 c O, N, H 383 384 do inter=1,ninter 385 efdiso(inter)=0. 386 efdisn(inter)=0. 387 efdish(inter)=0. 388 enddo 389 390 391 c H2O, H2O2, O3, NO2 392 393 do inter=25,31 394 efdish2o(inter)=1. 395 enddo 396 do inter=25,35 397 efdish2o2(inter)=1. 398 enddo 399 do inter=34,36 400 efdiso3(inter)=1. 401 enddo 402 do inter=27,36 403 efdisno2(inter)=1. 404 enddo 184 405 do inter=1,15 185 efdish2(inter) = 1. 186 end do 187 188 efdisco2(12) = 0.183 189 efdiso2(12) = 0.003 190 191 efdisco2(13) = 0.163 192 efdiso2(13) = 0.170 193 194 efdisco2(14) = 0.243 195 efdiso2(14) =0.180 196 197 efdisco2(15) = 0.323 198 efdiso2(15) =0.653 199 200 efdisco2(16) = 0.235 201 efdiso2(16) =0.616 202 203 do inter=17,32 204 efdisco2(inter) = 1.0 205 end do 206 207 efdiso2(17) = 0.399 208 do inter=18,20 209 efdiso2(inter) = 0.261 210 end do 211 212 do inter=21,23 213 efdiso2(inter) = 0.755 214 end do 215 216 217 do inter=24,31 218 efdiso2(inter) = 1. 219 end do 220 221 do inter=25,31 222 efdish2o(inter) = 1. 223 end do 224 225 do inter=25,33 226 efdish2o2(inter) = 1. 227 end do 406 efdish2(inter)=1. 407 enddo 408 409 !4 possible channels for CO2 ionization 410 do inter=14,16 411 efionco2(inter,1)=1.-efdisco2(inter) 412 enddo 413 efionco2(13,1)=0.805*(1.-efdisco2(13)) 414 efionco2(13,2)=0.195*(1.-efdisco2(13)) 415 do inter=11,12 416 efionco2(inter,3)=1.-efdisco2(inter) 417 enddo 418 efionco2(10,3)=0.9*(1.-efdisco2(10)) 419 efionco2(10,4)=0.1*(1.-efdisco2(10)) 420 do inter=2,9 421 efionco2(inter,4)=1.-efdisco2(inter) 422 enddo 423 424 !For O(3p) total ionization under 91.1 nm 425 do inter=1,16 426 efiono3p(inter)=1.d0 427 enddo 428 429 !2 channels for N2 ionization 430 do inter=9,15 431 efionn2(inter,1)=1.-efdisn2(inter) 432 enddo 433 do inter=2,8 434 efionn2(inter,2)=1.-efdisn2(inter) 435 enddo 436 437 !3 channels for CO ionization 438 do inter=11,16 439 efionco(inter,1)=1.-efdisco(inter) 440 enddo 441 efionco(10,1)=0.87*(1.-efdisco(10)) 442 efionco(10,2)=0.13*(1.-efdisco(10)) 443 do inter=8,9 444 efionco(inter,2)=1.-efdisco(inter) 445 enddo 446 efionco(7,2)=0.1*(1.-efdisco(7)) 447 efionco(7,3)=0.9*(1.-efdisco(7)) 448 do inter=2,6 449 efionco(inter,3)=1.-efdisco(inter) 450 enddo 451 452 !Total ionization under 85 nm for N 453 do inter=1,16 454 efionn(inter)=1. 455 enddo 456 457 !NO 458 do inter=2,28 459 efionno(inter)=1.-efdisno(inter) 460 enddo 461 462 !Total ionization under 90 nm for H 463 do inter=3,16 464 efionh(inter)=1. 465 enddo 228 466 229 467 … … 232 470 233 471 end 472
Note: See TracChangeset
for help on using the changeset viewer.
