MODULE mod_chem_spack_fexchem USE mod_chem_spack_rates, ONLY: rates ! subroutine IMPLICIT NONE PRIVATE PUBLIC :: fexchem ! subroutine CONTAINS SUBROUTINE fexchem(y,rk,chem,ngas,ijkbegin,ijkend,maxblock_size,nr) !------------------------------------------------------------------------ ! ! -- DESCRIPTION ! ! This routine computes the chemical production term. ! This routine is automatically generated by SPACK. ! Mechanism: ../Mechanism/RACM ! Species: ../Mechanism/ciRA77 ! !------------------------------------------------------------------------ ! ! -- INPUT VARIABLES ! ! Y: chemical concentrations. ! RK: kinetic rates. ! ZCSOURC: volumic emissions. ! ! -- INPUT/OUTPUT VARIABLES ! ! -- OUTPUT VARIABLES ! ! CHEM: array of chemical production terms. ! !------------------------------------------------------------------------ ! ! -- REMARKS ! !------------------------------------------------------------------------ ! ! -- MODIFICATIONS ! !------------------------------------------------------------------------ ! ! -- AUTHOR(S) ! ! SPACK. ! !------------------------------------------------------------------------ IMPLICIT NONE INTEGER , INTENT(IN) :: ngas INTEGER , INTENT(IN) :: ijkbegin INTEGER , INTENT(IN) :: ijkend INTEGER , INTENT(IN) :: maxblock_size INTEGER , INTENT(IN) :: nr DOUBLE PRECISION , INTENT(IN) :: rk(maxblock_size,nr) DOUBLE PRECISION , INTENT(IN) :: y(maxblock_size,ngas) DOUBLE PRECISION , INTENT(OUT) :: chem(maxblock_size,NGAS) DOUBLE PRECISION :: w(maxblock_size,nr) INTEGER :: ijk ! Compute reaction rates. call rates(rk,y,w,ngas,ijkbegin,ijkend,maxblock_size,nr) ! Chemical production terms. DO ijk=ijkbegin,ijkend chem(ijk, 1) = - w(ijk, 2) & - w(ijk, 3) & + w(ijk, 24) & - w(ijk, 25) & - w(ijk, 29) & - w(ijk, 30) & - w(ijk, 48) & - w(ijk, 49) & - w(ijk,106) & - w(ijk,107) & - w(ijk,108) & - w(ijk,109) & - w(ijk,110) & - w(ijk,111) & - w(ijk,112) & - w(ijk,113) & - w(ijk,114) & - w(ijk,115) & - w(ijk,120) & - w(ijk,123) & - w(ijk,126) & + w(ijk,165) & + w(ijk,167) chem(ijk, 2) = - w(ijk, 9) & - w(ijk, 32) & + w(ijk, 33) & + w(ijk, 34) & + 0.6000000000000000D-02 * w(ijk,107) & + 0.1100000000000000D-01 * w(ijk,108) & + 0.1000000000000000D-02 * w(ijk,109) & + 0.1000000000000000D-02 * w(ijk,110) & + 0.2000000000000000D-01 * w(ijk,111) & + 0.2000000000000000D-01 * w(ijk,112) chem(ijk, 3) = + w(ijk, 1) & + w(ijk, 4) & + w(ijk, 7) & - w(ijk, 35) & + w(ijk, 36) & - w(ijk, 38) & - w(ijk, 41) & - w(ijk, 48) & - 0.2000000000000000D+01 * w(ijk, 50) & - w(ijk, 51) & + w(ijk, 52) & - w(ijk,131) & - w(ijk,132) & - w(ijk,133) & - w(ijk,134) & - w(ijk,135) & - w(ijk,136) & - w(ijk,137) & - w(ijk,138) & - w(ijk,139) & - w(ijk,140) & - w(ijk,141) & - w(ijk,142) & - w(ijk,143) & - w(ijk,144) & - w(ijk,145) & - w(ijk,146) & - w(ijk,147) & - w(ijk,148) & - w(ijk,149) & - w(ijk,236) chem(ijk, 4) = - w(ijk, 1) & + w(ijk, 5) & + 0.6500000000000000D+00 * w(ijk, 6) & + w(ijk, 8) & + w(ijk, 21) & + w(ijk, 35) & - w(ijk, 36) & - w(ijk, 37) & - w(ijk, 39) & + w(ijk, 40) & + w(ijk, 41) & - w(ijk, 42) & + w(ijk, 43) & + 0.7000000000000000D+00 * w(ijk, 44) & + w(ijk, 45) & + w(ijk, 47) & + w(ijk, 48) & - w(ijk, 49) & + 0.2000000000000000D+01 * w(ijk, 50) & + 0.2000000000000000D+01 * w(ijk, 51) & - w(ijk, 53) & + w(ijk, 54) & + 0.2000000000000000D+01 * w(ijk, 55) & + w(ijk, 90) & + 0.5000000000000000D+00 * w(ijk, 96) & + 0.4000000000000000D+00 * w(ijk,105) & + 0.7000000000000000D+00 * w(ijk,115) & - w(ijk,116) & - w(ijk,118) & - w(ijk,121) & - w(ijk,124) & - w(ijk,127) & + w(ijk,128) & - w(ijk,129) & + w(ijk,130) & + w(ijk,131) & + w(ijk,132) & + 0.9409999999999999D+00 * w(ijk,133) & + 0.8760000000000000D+00 * w(ijk,134) & + 0.7390000000000000D+00 * w(ijk,135) & + w(ijk,136) & + w(ijk,137) & + w(ijk,138) & + 0.8470000000000000D+00 * w(ijk,139) & + 0.8000000000000000D+00 * w(ijk,140) & + 0.6500000000000000D+00 * w(ijk,141) & + 0.9500000000000000D+00 * w(ijk,142) & + 0.9500000000000000D+00 * w(ijk,143) & + w(ijk,144) & + w(ijk,145) & + w(ijk,146) & + w(ijk,147) & + w(ijk,148) & + 0.2000000000000000D+01 * w(ijk,149) & + 0.5000000000000000D+00 * w(ijk,191) & + 0.5160000000000000D+00 * w(ijk,209) & + 0.5000000000000000D+00 * w(ijk,211) & + w(ijk,212) & + w(ijk,213) & + w(ijk,214) & + w(ijk,215) & + w(ijk,216) & + w(ijk,217) & + w(ijk,218) & + w(ijk,219) & + w(ijk,220) & + w(ijk,221) & + w(ijk,222) & + w(ijk,223) & + w(ijk,224) & + w(ijk,225) & + w(ijk,226) & + w(ijk,227) & + w(ijk,228) & + w(ijk,229) & + w(ijk,230) & + 0.2000000000000000D+01 * w(ijk,231) & + w(ijk,236) & + w(ijk,237) chem(ijk, 5) = + 0.3500000000000000D+00 * w(ijk, 6) & - w(ijk, 7) & - w(ijk, 8) & + w(ijk, 37) & - w(ijk, 40) & - w(ijk, 44) & + w(ijk, 46) & + w(ijk, 49) & - w(ijk, 51) & - w(ijk, 52) & - w(ijk, 53) & + w(ijk, 54) & - 0.2000000000000000D+01 * w(ijk, 55) & + w(ijk, 88) & + 0.6000000000000000D+00 * w(ijk, 89) & - w(ijk, 91) & - w(ijk, 92) & - w(ijk, 93) & - w(ijk, 94) & - w(ijk, 95) & - w(ijk, 96) & - w(ijk, 97) & - w(ijk, 98) & - w(ijk, 99) & - w(ijk,100) & - w(ijk,101) & - w(ijk,102) & - w(ijk,103) & - w(ijk,104) & - 0.4000000000000000D+00 * w(ijk,105) & - w(ijk,213) & - w(ijk,214) & - w(ijk,215) & - w(ijk,216) & - w(ijk,217) & - w(ijk,218) & - w(ijk,219) & - w(ijk,220) & - w(ijk,221) & - w(ijk,222) & - w(ijk,223) & - w(ijk,224) & - w(ijk,225) & - w(ijk,226) & - w(ijk,227) & - w(ijk,228) & - w(ijk,229) & - w(ijk,230) & - w(ijk,231) & - w(ijk,237) chem(ijk, 6) = + w(ijk, 53) & - w(ijk, 54) chem(ijk, 7) = - w(ijk, 4) & + w(ijk, 38) & - w(ijk, 45) & + w(ijk,118) & + w(ijk,121) & + w(ijk,124) chem(ijk, 8) = - w(ijk, 5) & + w(ijk, 39) & + 0.3000000000000000D+00 * w(ijk, 44) & - w(ijk, 46) & + w(ijk, 91) & + w(ijk, 92) & + w(ijk, 93) & + w(ijk, 94) & + 0.2000000000000000D+00 * w(ijk, 95) & + 0.5000000000000000D+00 * w(ijk, 96) & + w(ijk, 97) chem(ijk, 9) = - w(ijk, 6) & + w(ijk, 42) & - w(ijk, 43) & - w(ijk, 47) chem(ijk, 10) = - w(ijk, 57) chem(ijk, 11) = + w(ijk, 57) chem(ijk, 12) = + w(ijk, 10) & + w(ijk, 11) & + w(ijk, 12) & + 0.1870000000000000D+01 * w(ijk, 17) & + 0.1550000000000000D+01 * w(ijk, 18) & + w(ijk, 19) & + w(ijk, 22) & - w(ijk, 58) & + 0.1000000000000000D-01 * w(ijk, 59) & + 0.3600000000000000D-01 * w(ijk, 63) & + w(ijk, 76) & + 0.2000000000000000D+01 * w(ijk, 80) & + w(ijk, 81) & + 0.4100000000000000D+00 * w(ijk, 82) & + w(ijk, 91) & + 0.2000000000000000D+01 * w(ijk, 93) & + w(ijk, 94) & + 0.8000000000000000D+00 * w(ijk, 95) & + 0.4300000000000000D+00 * w(ijk,106) & + 0.3700000000000000D+00 * w(ijk,107) & + 0.3000000000000000D+00 * w(ijk,108) & + 0.3600000000000000D+00 * w(ijk,109) & + 0.3600000000000000D+00 * w(ijk,110) & + 0.1400000000000000D+00 * w(ijk,111) & + 0.1400000000000000D+00 * w(ijk,112) & + 0.5400000000000000D+00 * w(ijk,113) & + 0.6600000000000000D+00 * w(ijk,114) & + 0.1300000000000000D+00 * w(ijk,115) chem(ijk, 13) = + w(ijk, 1) & + w(ijk, 3) & + w(ijk, 8) & - w(ijk, 24) & - w(ijk, 25) & + w(ijk, 26) & + w(ijk, 27) & - w(ijk, 35) & - w(ijk, 36) & - w(ijk, 37) & - w(ijk, 59) & - w(ijk, 60) & + 0.9000000000000000D-01 * w(ijk,109) & + 0.9000000000000000D-01 * w(ijk,110) chem(ijk, 14) = + w(ijk, 2) & - w(ijk, 26) & - w(ijk, 27) & - w(ijk, 28) chem(ijk, 15) = + w(ijk, 4) & + w(ijk, 5) & + 0.3500000000000000D+00 * w(ijk, 6) & + 0.2000000000000000D+01 * w(ijk, 9) & + w(ijk, 13) & + w(ijk, 14) & + w(ijk, 15) & + 0.2000000000000000D+01 * w(ijk, 28) & - w(ijk, 29) & + w(ijk, 30) & - w(ijk, 31) & - w(ijk, 32) & - w(ijk, 38) & - w(ijk, 39) & - w(ijk, 40) & + w(ijk, 41) & + 0.7000000000000000D+00 * w(ijk, 44) & - w(ijk, 45) & - w(ijk, 46) & - w(ijk, 47) & - w(ijk, 56) & - w(ijk, 57) & - w(ijk, 58) & + 0.2000000000000000D-01 * w(ijk, 59) & - w(ijk, 61) & - w(ijk, 62) & - 0.9640000000000000D+00 * w(ijk, 63) & - w(ijk, 64) & - w(ijk, 65) & - w(ijk, 66) & - w(ijk, 67) & - w(ijk, 68) & - w(ijk, 69) & - w(ijk, 70) & - w(ijk, 71) & - w(ijk, 72) & - w(ijk, 73) & - w(ijk, 74) & - w(ijk, 75) & - w(ijk, 76) & - w(ijk, 77) & - w(ijk, 78) & - w(ijk, 79) & - w(ijk, 80) & - w(ijk, 81) & - w(ijk, 82) & - w(ijk, 83) & - w(ijk, 84) & - 0.6500000000000000D+00 * w(ijk, 85) & - 0.5100000000000000D+00 * w(ijk, 86) & - w(ijk, 87) & - w(ijk, 88) & - w(ijk, 89) & - w(ijk, 90) & + 0.1200000000000000D+00 * w(ijk,106) & + 0.4000000000000000D+00 * w(ijk,107) & + 0.6300000000000000D+00 * w(ijk,108) & + 0.2800000000000000D+00 * w(ijk,109) & + 0.2800000000000000D+00 * w(ijk,110) & + 0.8500000000000000D+00 * w(ijk,111) & + 0.8500000000000000D+00 * w(ijk,112) & + 0.7000000000000001D-01 * w(ijk,113) & + 0.2100000000000000D+00 * w(ijk,114) & + 0.3600000000000000D-01 * w(ijk,115) & + w(ijk,120) & + w(ijk,123) & + w(ijk,126) chem(ijk, 16) = + 0.6500000000000000D+00 * w(ijk, 6) & + 0.2000000000000000D+01 * w(ijk, 11) & + w(ijk, 12) & + w(ijk, 13) & + w(ijk, 14) & + 0.8000000000000000D+00 * w(ijk, 18) & + w(ijk, 19) & + w(ijk, 20) & + w(ijk, 21) & + w(ijk, 22) & + w(ijk, 23) & + w(ijk, 29) & - w(ijk, 30) & - w(ijk, 31) & + w(ijk, 32) & - 0.2000000000000000D+01 * w(ijk, 33) & - 0.2000000000000000D+01 * w(ijk, 34) & + w(ijk, 40) & - w(ijk, 41) & - w(ijk, 42) & + w(ijk, 43) & - w(ijk, 44) & + w(ijk, 56) & + w(ijk, 57) & + w(ijk, 58) & + 0.2800000000000000D+00 * w(ijk, 59) & + 0.3810000000000000D+00 * w(ijk, 63) & + 0.2500000000000000D+00 * w(ijk, 64) & + 0.4900000000000000D-01 * w(ijk, 65) & + 0.1000000000000000D+00 * w(ijk, 73) & + 0.1000000000000000D+00 * w(ijk, 74) & + 0.5000000000000000D-01 * w(ijk, 75) & + w(ijk, 76) & + w(ijk, 79) & + w(ijk, 80) & + 0.4900000000000000D+00 * w(ijk, 82) & + 0.5000000000000000D+00 * w(ijk, 83) & + w(ijk, 84) & + 0.3500000000000000D+00 * w(ijk, 87) & + 0.4000000000000000D+00 * w(ijk, 89) & + w(ijk, 91) & + w(ijk, 93) & + 0.5000000000000000D+00 * w(ijk, 96) & + 0.2600000000000000D+00 * w(ijk,106) & + 0.2500000000000000D+00 * w(ijk,107) & + 0.2200000000000000D+00 * w(ijk,108) & + 0.3000000000000000D+00 * w(ijk,109) & + 0.3000000000000000D+00 * w(ijk,110) & + 0.1000000000000000D+00 * w(ijk,111) & + 0.1000000000000000D+00 * w(ijk,112) & + 0.2900000000000000D+00 * w(ijk,113) & + 0.2900000000000000D+00 * w(ijk,114) & + 0.8000000000000000D-01 * w(ijk,115) & - w(ijk,117) & + 0.2000000000000000D-01 * w(ijk,119) & + 0.2000000000000000D-01 * w(ijk,122) & + 0.2000000000000000D-01 * w(ijk,125) & + w(ijk,131) & + w(ijk,132) & + 0.7420000000000000D+00 * w(ijk,133) & + 0.5990000000000000D+00 * w(ijk,134) & + 0.6060000000000000D+00 * w(ijk,135) & + w(ijk,136) & + w(ijk,137) & + w(ijk,138) & + 0.8470000000000000D+00 * w(ijk,139) & + 0.8000000000000000D+00 * w(ijk,140) & + 0.6500000000000000D+00 * w(ijk,141) & + 0.9500000000000000D+00 * w(ijk,142) & + 0.9500000000000000D+00 * w(ijk,143) & + w(ijk,144) & + 0.7700000000000000D+00 * w(ijk,147) & + w(ijk,148) & - w(ijk,150) & - w(ijk,151) & - w(ijk,152) & - w(ijk,153) & - w(ijk,154) & - w(ijk,155) & - w(ijk,156) & - w(ijk,157) & - w(ijk,158) & - w(ijk,159) & - w(ijk,160) & - w(ijk,161) & - w(ijk,162) & - w(ijk,163) & - w(ijk,164) & - w(ijk,165) & - w(ijk,166) & - w(ijk,167) & - w(ijk,168) & - w(ijk,169) & - w(ijk,170) & + 0.6600000000000000D+00 * w(ijk,171) & + w(ijk,172) & + 0.9920000000000000D+00 * w(ijk,173) & + 0.9460000000000000D+00 * w(ijk,174) & + 0.9930000000000000D+00 * w(ijk,175) & + w(ijk,176) & + w(ijk,177) & + w(ijk,178) & + w(ijk,179) & + 0.2000000000000000D+01 * w(ijk,180) & + 0.2000000000000000D+01 * w(ijk,181) & + w(ijk,182) & + w(ijk,183) & + 0.2000000000000000D+01 * w(ijk,184) & + w(ijk,185) & + w(ijk,187) & + 0.8800000000000000D+00 * w(ijk,189) & + w(ijk,190) & + 0.5000000000000000D+00 * w(ijk,191) & + 0.5000000000000000D+00 * w(ijk,192) & + 0.4880000000000000D+00 * w(ijk,193) & + 0.4380000000000000D+00 * w(ijk,194) & + 0.4890000000000000D+00 * w(ijk,195) & + 0.5000000000000000D+00 * w(ijk,196) & + 0.5010000000000000D+00 * w(ijk,197) & + 0.5100000000000000D+00 * w(ijk,198) & + 0.5060000000000000D+00 * w(ijk,199) & + w(ijk,200) & + w(ijk,201) & + w(ijk,202) & + w(ijk,203) & + w(ijk,204) & + 0.3800000000000000D+00 * w(ijk,207) & + 0.5000000000000000D+00 * w(ijk,208) & + w(ijk,210) & + 0.5000000000000000D+00 * w(ijk,211) & + w(ijk,213) & + w(ijk,214) & + 0.7920000000000000D+00 * w(ijk,215) & + 0.6990000000000000D+00 * w(ijk,216) & + 0.8450000000000000D+00 * w(ijk,217) & + w(ijk,218) & + w(ijk,219) & + w(ijk,220) & + w(ijk,221) & + w(ijk,222) & + w(ijk,223) & + w(ijk,224) & + w(ijk,225) & + w(ijk,226) & + 0.7700000000000000D+00 * w(ijk,229) & + w(ijk,230) & - w(ijk,232) & + w(ijk,233) chem(ijk, 17) = - w(ijk, 61) & + 0.6000000000000000D-01 * w(ijk,107) & + 0.7000000000000001D-01 * w(ijk,108) chem(ijk, 18) = - w(ijk, 62) & + 0.3000000000000000D-01 * w(ijk,107) & + 0.6000000000000000D-01 * w(ijk,108) chem(ijk, 19) = - w(ijk, 63) chem(ijk, 20) = - w(ijk, 64) chem(ijk, 21) = - w(ijk, 65) chem(ijk, 22) = - w(ijk, 66) & - w(ijk, 98) & - w(ijk,106) chem(ijk, 23) = + 0.8600000000000000D+00 * w(ijk, 59) & - w(ijk, 67) & - w(ijk, 99) & - w(ijk,107) & + 0.3500000000000000D+00 * w(ijk,109) & + 0.3500000000000000D+00 * w(ijk,110) & + 0.4600000000000000D+00 * w(ijk,112) & + 0.3540000000000000D+00 * w(ijk,139) & + 0.3700000000000000D+00 * w(ijk,179) & + 0.2290000000000000D+00 * w(ijk,199) & + 0.4000000000000000D+00 * w(ijk,221) chem(ijk, 24) = - w(ijk, 68) & - w(ijk,100) & - w(ijk,108) & + 0.2500000000000000D+00 * w(ijk,141) & + 0.8000000000000000D-01 * w(ijk,179) & + 0.4000000000000000D+00 * w(ijk,181) & + 0.4000000000000000D+00 * w(ijk,201) & + 0.4000000000000000D+00 * w(ijk,223) chem(ijk, 25) = - w(ijk, 69) & - w(ijk,101) & - w(ijk,109) chem(ijk, 26) = - w(ijk, 59) & - w(ijk, 70) & - w(ijk,102) & - w(ijk,110) chem(ijk, 27) = - w(ijk, 71) & - w(ijk,103) & - w(ijk,111) chem(ijk, 28) = - w(ijk, 72) & - w(ijk,104) & - w(ijk,112) chem(ijk, 29) = - w(ijk, 73) chem(ijk, 30) = - w(ijk, 74) chem(ijk, 31) = - w(ijk, 75) & - w(ijk, 97) & + 0.1000000000000000D+00 * w(ijk,116) & + w(ijk,117) & + w(ijk,118) & + 0.2000000000000000D-01 * w(ijk,119) & + w(ijk,120) & + w(ijk,121) & + 0.2000000000000000D-01 * w(ijk,122) & + w(ijk,123) & + w(ijk,124) & + 0.2000000000000000D-01 * w(ijk,125) & + w(ijk,126) chem(ijk, 32) = - w(ijk, 10) & - w(ijk, 11) & + w(ijk, 13) & + 0.1300000000000000D+00 * w(ijk, 17) & + 0.4500000000000000D+00 * w(ijk, 18) & + w(ijk, 22) & + w(ijk, 23) & + 0.5000000000000000D-01 * w(ijk, 59) & + 0.1000000000000000D-01 * w(ijk, 63) & - w(ijk, 76) & + 0.8000000000000000D-01 * w(ijk, 82) & + 0.3500000000000000D+00 * w(ijk, 85) & + 0.3500000000000000D+00 * w(ijk, 87) & + w(ijk, 88) & + 0.4000000000000000D+00 * w(ijk, 89) & - w(ijk, 91) & + 0.4000000000000000D+00 * w(ijk,105) & + w(ijk,106) & + 0.6400000000000000D+00 * w(ijk,107) & + 0.2000000000000000D-01 * w(ijk,108) & + 0.9000000000000000D+00 * w(ijk,109) & + 0.9000000000000000D+00 * w(ijk,110) & + 0.4000000000000000D-01 * w(ijk,112) & + 0.4000000000000000D+00 * w(ijk,113) & + 0.7000000000000000D+00 * w(ijk,115) & + w(ijk,131) & + 0.4700000000000000D-01 * w(ijk,133) & + 0.2100000000000000D-01 * w(ijk,134) & + 0.1600000000000000D+01 * w(ijk,136) & + w(ijk,137) & + 0.6060000000000000D+00 * w(ijk,139) & + 0.2500000000000000D+00 * w(ijk,141) & + w(ijk,146) & + 0.2870000000000000D+00 * w(ijk,149) & + 0.1330000000000000D+01 * w(ijk,171) & + 0.7500000000000000D+00 * w(ijk,172) & + 0.8100000000000001D+00 * w(ijk,173) & + 0.8290000000000000D+00 * w(ijk,174) & + 0.7530000000000000D+00 * w(ijk,175) & + 0.1550000000000000D+01 * w(ijk,176) & + 0.1250000000000000D+01 * w(ijk,177) & + 0.7550000000000000D+00 * w(ijk,178) & + 0.1090000000000000D+01 * w(ijk,179) & + w(ijk,180) & + 0.1400000000000000D+01 * w(ijk,181) & + w(ijk,182) & + w(ijk,183) & + w(ijk,184) & + w(ijk,185) & + w(ijk,186) & + 0.2000000000000000D+01 * w(ijk,187) & + w(ijk,188) & + 0.7500000000000000D+00 * w(ijk,189) & + 0.7500000000000000D+00 * w(ijk,190) & + 0.9600000000000000D+00 * w(ijk,191) & + 0.9100000000000000D-01 * w(ijk,193) & + 0.7600000000000000D-01 * w(ijk,194) & + 0.8000000000000000D+00 * w(ijk,196) & + 0.5010000000000000D+00 * w(ijk,197) & + 0.3400000000000000D+00 * w(ijk,199) & + 0.4000000000000000D+00 * w(ijk,201) & + w(ijk,206) & + 0.2070000000000000D+00 * w(ijk,209) & + 0.2020000000000000D+00 * w(ijk,211) & + 0.5040000000000000D+00 * w(ijk,212) & + w(ijk,213) & + 0.4800000000000000D-01 * w(ijk,215) & + 0.2100000000000000D-01 * w(ijk,216) & + 0.1600000000000000D+01 * w(ijk,218) & + w(ijk,219) & + 0.6860000000000001D+00 * w(ijk,221) & + 0.4000000000000000D+00 * w(ijk,223) & + w(ijk,228) & + 0.2800000000000000D+00 * w(ijk,231) & + w(ijk,233) chem(ijk, 33) = - w(ijk, 12) & + w(ijk, 14) & + 0.2000000000000000D+00 * w(ijk, 21) & + w(ijk, 60) & + 0.3350000000000000D+00 * w(ijk, 63) & + 0.2500000000000000D-01 * w(ijk, 65) & - w(ijk, 77) & + 0.8800000000000000D+00 * w(ijk, 84) & + 0.8000000000000000D-01 * w(ijk, 86) & - w(ijk, 92) & + 0.2500000000000000D+00 * w(ijk, 96) & + 0.4400000000000000D+00 * w(ijk,107) & + 0.9900000000000000D+00 * w(ijk,108) & + 0.6500000000000000D+00 * w(ijk,111) & + 0.1600000000000000D+00 * w(ijk,114) & + w(ijk,132) & + 0.2330000000000000D+00 * w(ijk,133) & + 0.2110000000000000D+00 * w(ijk,134) & + 0.1500000000000000D+00 * w(ijk,135) & + 0.2000000000000000D+00 * w(ijk,136) & + 0.9399999999999999D+00 * w(ijk,137) & + 0.1710000000000000D+01 * w(ijk,138) & + 0.8000000000000000D+00 * w(ijk,140) & + 0.4600000000000000D+00 * w(ijk,147) & + 0.1240000000000000D+01 * w(ijk,149) & + 0.7500000000000000D+00 * w(ijk,172) & + 0.5800000000000000D+00 * w(ijk,173) & + 0.5230000000000000D+00 * w(ijk,174) & + 0.4110000000000000D+00 * w(ijk,175) & + 0.3500000000000000D+00 * w(ijk,176) & + 0.6690000000000000D+00 * w(ijk,177) & + 0.9320000000000001D+00 * w(ijk,178) & + w(ijk,180) & + 0.3000000000000000D+00 * w(ijk,189) & + 0.6400000000000000D+00 * w(ijk,191) & + w(ijk,192) & + 0.7240000000000000D+00 * w(ijk,193) & + 0.6770000000000000D+00 * w(ijk,194) & + 0.4970000000000000D+00 * w(ijk,195) & + 0.6000000000000000D+00 * w(ijk,196) & + 0.8590000000000000D+00 * w(ijk,197) & + 0.9409999999999999D+00 * w(ijk,198) & + w(ijk,200) & + 0.3500000000000000D+00 * w(ijk,207) & + 0.6500000000000000D+00 * w(ijk,209) & + 0.6400000000000000D+00 * w(ijk,211) & + 0.1210000000000000D+01 * w(ijk,212) & + w(ijk,214) & + 0.2430000000000000D+00 * w(ijk,215) & + 0.2390000000000000D+00 * w(ijk,216) & + 0.1870000000000000D+00 * w(ijk,217) & + 0.2000000000000000D+00 * w(ijk,218) & + 0.9399999999999999D+00 * w(ijk,219) & + 0.1710000000000000D+01 * w(ijk,220) & + w(ijk,222) & + 0.4600000000000000D+00 * w(ijk,229) & + 0.1240000000000000D+01 * w(ijk,231) chem(ijk, 34) = - w(ijk, 16) & + 0.8000000000000000D+00 * w(ijk, 21) & + 0.2500000000000000D+00 * w(ijk, 64) & - w(ijk, 78) & + 0.1200000000000000D+00 * w(ijk, 84) & + 0.4100000000000000D+00 * w(ijk, 86) & + 0.3000000000000000D-01 * w(ijk, 96) & + 0.3000000000000000D-01 * w(ijk,107) & + 0.1600000000000000D+00 * w(ijk,108) & + 0.5300000000000000D+00 * w(ijk,111) & + 0.6230000000000000D+00 * w(ijk,133) & + 0.7220000000000000D+00 * w(ijk,134) & + 0.6420000000000000D+00 * w(ijk,135) & + 0.6000000000000000D-01 * w(ijk,137) & + 0.2900000000000000D+00 * w(ijk,138) & + 0.8000000000000000D+00 * w(ijk,140) & + 0.4640000000000000D+00 * w(ijk,149) & + 0.1800000000000000D-01 * w(ijk,173) & + 0.2400000000000000D+00 * w(ijk,174) & + 0.4190000000000000D+00 * w(ijk,175) & + 0.8100000000000000D-01 * w(ijk,177) & + 0.3130000000000000D+00 * w(ijk,178) & + w(ijk,180) & + 0.1490000000000000D+00 * w(ijk,191) & + 0.1270000000000000D+00 * w(ijk,193) & + 0.3300000000000000D+00 * w(ijk,194) & + 0.5810000000000000D+00 * w(ijk,195) & + 0.1410000000000000D+00 * w(ijk,197) & + 0.5690000000000000D+00 * w(ijk,198) & + w(ijk,200) & + 0.1100000000000000D+00 * w(ijk,207) & + 0.1670000000000000D+00 * w(ijk,209) & + 0.1490000000000000D+00 * w(ijk,211) & + 0.2850000000000000D+00 * w(ijk,212) & + 0.6700000000000000D+00 * w(ijk,215) & + 0.8280000000000000D+00 * w(ijk,216) & + 0.8800000000000000D+00 * w(ijk,217) & + 0.6000000000000000D-01 * w(ijk,219) & + 0.2900000000000000D+00 * w(ijk,220) & + w(ijk,222) & + 0.4690000000000000D+00 * w(ijk,231) chem(ijk, 35) = - w(ijk, 17) & - w(ijk, 18) & + 0.3600000000000000D-01 * w(ijk, 63) & - w(ijk, 80) & + 0.1500000000000000D+00 * w(ijk, 83) & - w(ijk, 93) & + 0.2500000000000000D+00 * w(ijk, 96) & + 0.5000000000000000D+00 * w(ijk,114) & + 0.6300000000000000D-01 * w(ijk,133) & + 0.1200000000000000D+01 * w(ijk,142) & + 0.3500000000000000D+00 * w(ijk,143) & + w(ijk,144) & + 0.1190000000000000D+00 * w(ijk,173) & + 0.6500000000000000D+00 * w(ijk,182) & + 0.3700000000000000D+00 * w(ijk,183) & + w(ijk,184) & + 0.1000000000000000D+00 * w(ijk,193) & + 0.6500000000000000D+00 * w(ijk,202) & + 0.3700000000000000D+00 * w(ijk,203) & + w(ijk,204) & + 0.6300000000000000D-01 * w(ijk,215) & + 0.1300000000000000D+01 * w(ijk,224) & + 0.7400000000000000D+00 * w(ijk,225) & + w(ijk,226) chem(ijk, 36) = - w(ijk, 19) & + w(ijk, 79) & - w(ijk, 81) & + 0.8000000000000000D-01 * w(ijk, 82) & + 0.1500000000000000D+00 * w(ijk, 83) & - w(ijk, 94) & + 0.2500000000000000D+00 * w(ijk, 96) & + 0.6000000000000000D+00 * w(ijk,113) & + 0.6200000000000000D+00 * w(ijk,114) & + 0.6500000000000000D+00 * w(ijk,142) & + 0.6000000000000000D+00 * w(ijk,143) & + w(ijk,144) & + 0.5400000000000000D+00 * w(ijk,147) & + 0.5000000000000000D-02 * w(ijk,173) & + 0.3500000000000000D+00 * w(ijk,182) & + 0.6300000000000000D+00 * w(ijk,183) & + w(ijk,184) & + 0.4000000000000000D+00 * w(ijk,189) & + 0.4000000000000000D-02 * w(ijk,193) & + 0.3500000000000000D+00 * w(ijk,202) & + 0.6300000000000000D+00 * w(ijk,203) & + w(ijk,204) & + 0.5400000000000000D+00 * w(ijk,207) & + 0.7000000000000000D+00 * w(ijk,224) & + 0.1260000000000000D+01 * w(ijk,225) & + w(ijk,226) & + 0.5400000000000000D+00 * w(ijk,229) chem(ijk, 37) = - w(ijk, 20) & + 0.1300000000000000D+00 * w(ijk, 59) & - w(ijk, 83) & - w(ijk, 96) & - w(ijk,114) & + 0.5000000000000000D+00 * w(ijk,142) & + 0.9500000000000000D+00 * w(ijk,143) & + w(ijk,182) & + w(ijk,183) & + w(ijk,202) & + w(ijk,203) & + 0.5000000000000000D+00 * w(ijk,224) & + w(ijk,225) chem(ijk, 38) = - w(ijk, 22) & - w(ijk, 60) & - w(ijk, 82) & - w(ijk, 95) & + 0.9000000000000000D+00 * w(ijk,101) & + 0.9000000000000000D+00 * w(ijk,102) & + 0.3900000000000000D+00 * w(ijk,109) & + 0.3900000000000000D+00 * w(ijk,110) & + 0.7900000000000000D+00 * w(ijk,112) & - w(ijk,113) & + 0.4460000000000000D+00 * w(ijk,139) & + 0.4000000000000000D+00 * w(ijk,141) & + 0.5500000000000000D+00 * w(ijk,179) & + 0.6000000000000000D+00 * w(ijk,181) & + 0.7710000000000000D+00 * w(ijk,199) & + 0.6000000000000000D+00 * w(ijk,201) & + 0.6000000000000000D+00 * w(ijk,221) & + 0.6000000000000000D+00 * w(ijk,223) chem(ijk, 39) = + 0.3500000000000000D+00 * w(ijk, 83) & - w(ijk, 84) chem(ijk, 40) = - w(ijk, 23) & + 0.2400000000000000D-01 * w(ijk, 65) & - w(ijk, 79) & + 0.4100000000000000D+00 * w(ijk, 82) & + 0.6000000000000000D+00 * w(ijk, 89) & + 0.3000000000000000D+00 * w(ijk,189) chem(ijk, 41) = - w(ijk, 21) & - w(ijk, 90) & + 0.6000000000000000D+00 * w(ijk,105) & + w(ijk,116) & + 0.5900000000000000D-01 * w(ijk,133) & + 0.1240000000000000D+00 * w(ijk,134) & + 0.2610000000000000D+00 * w(ijk,135) & + 0.1530000000000000D+00 * w(ijk,139) & + 0.2000000000000000D+00 * w(ijk,140) & + 0.3500000000000000D+00 * w(ijk,141) & + 0.5000000000000000D-01 * w(ijk,142) & + 0.5000000000000000D-01 * w(ijk,143) & + w(ijk,148) & + w(ijk,169) & + w(ijk,170) & + w(ijk,190) & + 0.5000000000000000D+00 * w(ijk,191) & + w(ijk,208) & + 0.4840000000000000D+00 * w(ijk,209) & + 0.2000000000000000D+01 * w(ijk,210) & + 0.1500000000000000D+01 * w(ijk,211) & + w(ijk,212) & + w(ijk,230) chem(ijk, 42) = - w(ijk, 88) & + 0.4000000000000000D+00 * w(ijk, 89) & + 0.4000000000000000D+00 * w(ijk,105) & + 0.3000000000000000D+00 * w(ijk,115) & + w(ijk,127) & - w(ijk,128) chem(ijk, 43) = - w(ijk, 89) & - w(ijk,105) & - w(ijk,115) & + w(ijk,129) & - w(ijk,130) chem(ijk, 44) = - w(ijk, 13) & - w(ijk, 85) & + w(ijk,150) chem(ijk, 45) = - w(ijk, 14) & - w(ijk, 86) & + 0.1300000000000000D+00 * w(ijk,113) & + w(ijk,151) & + w(ijk,152) & + w(ijk,153) & + w(ijk,154) & + w(ijk,155) & + w(ijk,156) & + w(ijk,157) & + w(ijk,158) & + w(ijk,159) & + w(ijk,160) & + w(ijk,161) & + w(ijk,162) & + w(ijk,163) & + w(ijk,166) & + w(ijk,168) & + w(ijk,232) chem(ijk, 46) = - w(ijk, 15) & - w(ijk, 87) & + 0.1100000000000000D+00 * w(ijk,114) & + w(ijk,164) chem(ijk, 47) = + 0.3600000000000000D-01 * w(ijk, 63) & + 0.3700000000000000D+00 * w(ijk,106) & + 0.1400000000000000D+00 * w(ijk,107) & + 0.1500000000000000D+00 * w(ijk,109) & + 0.1500000000000000D+00 * w(ijk,110) & + 0.1000000000000000D-01 * w(ijk,112) & + 0.2200000000000000D+00 * w(ijk,113) & + 0.1100000000000000D+00 * w(ijk,114) & + 0.1100000000000000D+00 * w(ijk,115) chem(ijk, 48) = + 0.1000000000000000D+00 * w(ijk,107) & + 0.1400000000000000D+00 * w(ijk,108) & + 0.7000000000000001D-01 * w(ijk,112) & + 0.1300000000000000D+00 * w(ijk,113) & + 0.2100000000000000D+00 * w(ijk,114) & + w(ijk,165) & + w(ijk,167) & + w(ijk,186) & + w(ijk,188) & + 0.5000000000000000D+00 * w(ijk,192) & + 0.4990000000000000D+00 * w(ijk,193) & + 0.4950000000000000D+00 * w(ijk,194) & + 0.4950000000000000D+00 * w(ijk,195) & + 0.5000000000000000D+00 * w(ijk,196) & + 0.4990000000000000D+00 * w(ijk,197) & + 0.4900000000000000D+00 * w(ijk,198) & + 0.4940000000000000D+00 * w(ijk,199) & + 0.5000000000000000D+00 * w(ijk,207) & + 0.5000000000000000D+00 * w(ijk,208) & + 0.4840000000000000D+00 * w(ijk,209) chem(ijk, 49) = + w(ijk, 12) & + w(ijk, 15) & + w(ijk, 61) & + 0.6500000000000000D+00 * w(ijk, 85) & + 0.1900000000000000D+00 * w(ijk,107) & + 0.2300000000000000D+00 * w(ijk,108) & + 0.3000000000000000D-01 * w(ijk,109) & + 0.3000000000000000D-01 * w(ijk,110) & - w(ijk,131) & + 0.1500000000000000D+00 * w(ijk,133) & + 0.3100000000000000D-01 * w(ijk,134) & + w(ijk,145) & - w(ijk,150) & - 0.2000000000000000D+01 * w(ijk,171) & - w(ijk,172) & - 0.9930000000000000D+00 * w(ijk,173) & - 0.9510000000000000D+00 * w(ijk,174) & - w(ijk,175) & - w(ijk,176) & - w(ijk,177) & - w(ijk,178) & - w(ijk,179) & - w(ijk,180) & - w(ijk,181) & - w(ijk,182) & - w(ijk,183) & - w(ijk,184) & - w(ijk,186) & - w(ijk,187) & - w(ijk,188) & - w(ijk,189) & - w(ijk,190) & - w(ijk,191) & + 0.5000000000000000D+00 * w(ijk,192) & + 0.5080000000000000D+00 * w(ijk,193) & + 0.5540000000000000D+00 * w(ijk,194) & + 0.5070000000000000D+00 * w(ijk,195) & + 0.5000000000000000D+00 * w(ijk,196) & + 0.5010000000000000D+00 * w(ijk,197) & + 0.5100000000000000D+00 * w(ijk,198) & + 0.5060000000000000D+00 * w(ijk,199) & + w(ijk,200) & + w(ijk,201) & + w(ijk,202) & + w(ijk,203) & + w(ijk,204) & + 0.2000000000000000D+01 * w(ijk,205) & + w(ijk,206) & + 0.5000000000000000D+00 * w(ijk,207) & + 0.5000000000000000D+00 * w(ijk,208) & + 0.5160000000000000D+00 * w(ijk,209) & - w(ijk,213) & + 0.1550000000000000D+00 * w(ijk,215) & + 0.4000000000000000D-01 * w(ijk,216) & + w(ijk,227) & - w(ijk,233) & + w(ijk,234) chem(ijk, 50) = + w(ijk, 16) & + w(ijk, 62) & + 0.1000000000000000D+00 * w(ijk,107) & + 0.1800000000000000D+00 * w(ijk,108) & + 0.2000000000000000D+00 * w(ijk,111) & + 0.1600000000000000D+00 * w(ijk,112) & - w(ijk,132) & + 0.4800000000000000D-01 * w(ijk,133) & + 0.2450000000000000D+00 * w(ijk,134) & + 0.1330000000000000D+00 * w(ijk,135) & - w(ijk,151) & - w(ijk,172) & + 0.1400000000000000D-01 * w(ijk,174) & + 0.1300000000000000D-01 * w(ijk,175) & - w(ijk,192) & + 0.6000000000000000D-02 * w(ijk,193) & + 0.1800000000000000D-01 * w(ijk,194) & + 0.1500000000000000D-01 * w(ijk,195) & - w(ijk,214) & + 0.5300000000000000D-01 * w(ijk,215) & + 0.2620000000000000D+00 * w(ijk,216) & + 0.1550000000000000D+00 * w(ijk,217) chem(ijk, 51) = + 0.5830000000000000D+00 * w(ijk, 63) & + 0.4400000000000000D+00 * w(ijk, 86) & + w(ijk, 90) & - w(ijk,133) & - w(ijk,152) & - w(ijk,173) & - w(ijk,193) & - w(ijk,215) chem(ijk, 52) = + 0.7500000000000000D+00 * w(ijk, 64) & - w(ijk,134) & - w(ijk,153) & - w(ijk,174) & - w(ijk,194) & - w(ijk,216) chem(ijk, 53) = + 0.9510000000000000D+00 * w(ijk, 65) & - w(ijk,135) & - w(ijk,154) & - w(ijk,175) & - w(ijk,195) & - w(ijk,217) chem(ijk, 54) = + w(ijk, 66) & - w(ijk,136) & - w(ijk,155) & - w(ijk,176) & - w(ijk,196) & - w(ijk,218) chem(ijk, 55) = + w(ijk, 67) & - w(ijk,137) & - w(ijk,156) & - w(ijk,177) & - w(ijk,197) & - w(ijk,219) chem(ijk, 56) = + w(ijk, 68) & - w(ijk,138) & - w(ijk,157) & - w(ijk,178) & - w(ijk,198) & - w(ijk,220) chem(ijk, 57) = + w(ijk, 69) & + w(ijk, 70) & - w(ijk,139) & - w(ijk,158) & - w(ijk,179) & - w(ijk,199) & - w(ijk,221) chem(ijk, 58) = + w(ijk, 71) & - w(ijk,140) & - w(ijk,159) & - w(ijk,180) & - w(ijk,200) & - w(ijk,222) chem(ijk, 59) = + w(ijk, 72) & - w(ijk,141) & - w(ijk,160) & - w(ijk,181) & - w(ijk,201) & - w(ijk,223) chem(ijk, 60) = + 0.1000000000000000D+00 * w(ijk, 75) & + w(ijk, 97) & - w(ijk,116) & - w(ijk,117) chem(ijk, 61) = + 0.9000000000000000D+00 * w(ijk, 73) & - w(ijk,118) & - w(ijk,119) & - w(ijk,120) chem(ijk, 62) = + 0.9000000000000000D+00 * w(ijk, 74) & - w(ijk,121) & - w(ijk,122) & - w(ijk,123) chem(ijk, 63) = + 0.8500000000000000D+00 * w(ijk, 75) & - w(ijk,124) & - w(ijk,125) & - w(ijk,126) chem(ijk, 64) = + 0.9800000000000000D+00 * w(ijk,119) & - w(ijk,142) & - w(ijk,161) & - w(ijk,182) & - w(ijk,202) & - w(ijk,224) chem(ijk, 65) = + 0.9800000000000000D+00 * w(ijk,122) & - w(ijk,143) & - w(ijk,162) & - w(ijk,183) & - w(ijk,203) & - w(ijk,225) chem(ijk, 66) = + 0.9800000000000000D+00 * w(ijk,125) & - w(ijk,144) & - w(ijk,163) & - w(ijk,184) & - w(ijk,204) & - w(ijk,226) chem(ijk, 67) = + w(ijk, 16) & + w(ijk, 19) & + w(ijk, 22) & + w(ijk, 23) & + w(ijk, 77) & + w(ijk, 81) & + 0.6500000000000000D+00 * w(ijk, 87) & + w(ijk, 92) & + w(ijk, 94) & + 0.1500000000000000D+00 * w(ijk,109) & + 0.1500000000000000D+00 * w(ijk,110) & + 0.1300000000000000D+00 * w(ijk,113) & + 0.2800000000000000D+00 * w(ijk,114) & + 0.7000000000000000D+00 * w(ijk,115) & - w(ijk,127) & + w(ijk,128) & - w(ijk,145) & + w(ijk,146) & + 0.2300000000000000D+00 * w(ijk,147) & - w(ijk,164) & - w(ijk,165) & - w(ijk,185) & - w(ijk,186) & + w(ijk,187) & + 0.1200000000000000D+00 * w(ijk,189) & - w(ijk,192) & - w(ijk,193) & - w(ijk,194) & - w(ijk,195) & - w(ijk,196) & - w(ijk,197) & - w(ijk,198) & - w(ijk,199) & - w(ijk,200) & - w(ijk,201) & - w(ijk,202) & - w(ijk,203) & - w(ijk,204) & - 0.2000000000000000D+01 * w(ijk,205) & - 0.8800000000000000D+00 * w(ijk,207) & - w(ijk,208) & - w(ijk,209) & - w(ijk,227) & + w(ijk,228) & + 0.2300000000000000D+00 * w(ijk,229) & - w(ijk,234) chem(ijk, 68) = + w(ijk, 20) & + 0.5100000000000000D+00 * w(ijk, 82) & + 0.5000000000000000D+00 * w(ijk, 83) & + 0.2000000000000000D+00 * w(ijk, 95) & + 0.5000000000000000D+00 * w(ijk, 96) & - w(ijk,129) & + w(ijk,130) & - w(ijk,146) & - w(ijk,166) & - w(ijk,167) & - w(ijk,187) & - w(ijk,188) & - w(ijk,206) & - w(ijk,228) chem(ijk, 69) = + w(ijk, 78) & + 0.3000000000000000D-01 * w(ijk,107) & + 0.1200000000000000D+00 * w(ijk,108) & + 0.2000000000000000D-01 * w(ijk,109) & + 0.2000000000000000D-01 * w(ijk,110) & + 0.4200000000000000D+00 * w(ijk,111) & + 0.4200000000000000D+00 * w(ijk,112) & - w(ijk,147) & - w(ijk,168) & - w(ijk,189) & - w(ijk,207) & - w(ijk,229) chem(ijk, 70) = + 0.8000000000000000D+00 * w(ijk, 95) & + 0.8000000000000000D+00 * w(ijk, 98) & + 0.4300000000000000D+00 * w(ijk, 99) & + 0.1100000000000000D+00 * w(ijk,100) & + 0.9000000000000000D+00 * w(ijk,101) & + 0.9000000000000000D+00 * w(ijk,102) & + 0.1000000000000000D+00 * w(ijk,103) & + 0.1300000000000000D+00 * w(ijk,104) & - w(ijk,148) & - w(ijk,169) & - w(ijk,190) & - w(ijk,208) & - 0.2000000000000000D+01 * w(ijk,210) & - w(ijk,211) & - w(ijk,230) chem(ijk, 71) = + 0.2000000000000000D+00 * w(ijk, 98) & + 0.5700000000000000D+00 * w(ijk, 99) & + 0.8900000000000000D+00 * w(ijk,100) & + 0.1000000000000000D+00 * w(ijk,101) & + 0.1000000000000000D+00 * w(ijk,102) & + 0.9000000000000000D+00 * w(ijk,103) & + 0.8700000000000000D+00 * w(ijk,104) & - w(ijk,149) & - w(ijk,170) & - w(ijk,191) & - w(ijk,209) & - w(ijk,211) & - 0.2000000000000000D+01 * w(ijk,212) & - w(ijk,231) chem(ijk, 72) = + 0.1500000000000000D+00 * w(ijk, 59) & + 0.1000000000000000D+00 * w(ijk, 73) & + 0.1000000000000000D+00 * w(ijk, 74) & + 0.5000000000000000D-01 * w(ijk, 75) & + 0.4900000000000000D+00 * w(ijk, 82) & + 0.5000000000000000D+00 * w(ijk, 83) & + 0.7000000000000001D-01 * w(ijk, 86) & + 0.3500000000000000D+00 * w(ijk, 87) & + w(ijk, 88) & + w(ijk, 89) & + 0.5000000000000000D+00 * w(ijk, 96) & + w(ijk,105) & + 0.1300000000000000D+00 * w(ijk,109) & + 0.1300000000000000D+00 * w(ijk,110) & + 0.4800000000000000D-01 * w(ijk,133) & + 0.3340000000000000D+00 * w(ijk,134) & + 0.4160000000000000D+00 * w(ijk,135) & + 0.1600000000000000D+00 * w(ijk,147) & + 0.8500000000000001D-01 * w(ijk,173) & + 0.2450000000000000D+00 * w(ijk,174) & + 0.3220000000000000D+00 * w(ijk,175) & + 0.8000000000000000D-01 * w(ijk,189) & + 0.7099999999999999D-01 * w(ijk,193) & + 0.2370000000000000D+00 * w(ijk,194) & + 0.3180000000000000D+00 * w(ijk,195) & + 0.8000000000000000D-01 * w(ijk,207) & + 0.5100000000000000D-01 * w(ijk,215) & + 0.3910000000000000D+00 * w(ijk,216) & + 0.5870000000000000D+00 * w(ijk,217) & + 0.1600000000000000D+00 * w(ijk,229) & - w(ijk,232) & - w(ijk,233) & - w(ijk,234) & - 0.2000000000000000D+01 * w(ijk,235) & - w(ijk,236) & - w(ijk,237) END DO ! Volumic source terms. END SUBROUTINE fexchem END MODULE mod_chem_spack_fexchem