module mo_indprd use shr_kind_mod, only : r8 => shr_kind_r8 private public :: indprd contains subroutine indprd( class, prod, nprod, y, extfrc, rxt, ncol ) use chem_mods, only : gas_pcnst, extcnt, rxntot use ppgrid, only : pver implicit none !-------------------------------------------------------------------- ! ... dummy arguments !-------------------------------------------------------------------- integer, intent(in) :: class integer, intent(in) :: ncol integer, intent(in) :: nprod real(r8), intent(in) :: y(ncol,pver,gas_pcnst) real(r8), intent(in) :: rxt(ncol,pver,rxntot) real(r8), intent(in) :: extfrc(ncol,pver,extcnt) real(r8), intent(inout) :: prod(ncol,pver,nprod) !-------------------------------------------------------------------- ! ... "independent" production for Explicit species !-------------------------------------------------------------------- if( class == 1 ) then prod(:,:,1) =.080_r8*rxt(:,:,256)*y(:,:,56)*y(:,:,1) prod(:,:,2) =rxt(:,:,135)*y(:,:,7)*y(:,:,5) prod(:,:,3) = 0._r8 prod(:,:,4) = 0._r8 prod(:,:,5) = 0._r8 prod(:,:,6) = 0._r8 prod(:,:,7) = 0._r8 prod(:,:,8) = 0._r8 prod(:,:,9) = 0._r8 prod(:,:,10) = 0._r8 prod(:,:,11) = 0._r8 prod(:,:,12) = 0._r8 prod(:,:,13) = 0._r8 prod(:,:,14) = 0._r8 prod(:,:,15) = 0._r8 prod(:,:,16) = 0._r8 prod(:,:,17) = 0._r8 prod(:,:,18) = 0._r8 prod(:,:,19) = 0._r8 prod(:,:,20) = 0._r8 prod(:,:,21) = (rxt(:,:,211)*y(:,:,17) +rxt(:,:,212)*y(:,:,17) + & rxt(:,:,223)*y(:,:,119) +rxt(:,:,236)*y(:,:,46) + & .500_r8*rxt(:,:,249)*y(:,:,51) +.800_r8*rxt(:,:,250)*y(:,:,49) + & rxt(:,:,251)*y(:,:,50) +.500_r8*rxt(:,:,300)*y(:,:,78))*y(:,:,20) & + (rxt(:,:,244)*y(:,:,6) +.900_r8*rxt(:,:,247)*y(:,:,13) + & 2.000_r8*rxt(:,:,248)*y(:,:,45) +2.000_r8*rxt(:,:,296)*y(:,:,73) + & rxt(:,:,324)*y(:,:,88))*y(:,:,45) + (rxt(:,:,295)*y(:,:,13) + & 2.000_r8*rxt(:,:,297)*y(:,:,73))*y(:,:,73) +rxt(:,:,58)*y(:,:,51) & +.400_r8*rxt(:,:,59)*y(:,:,55) prod(:,:,22) = 0._r8 prod(:,:,23) = 0._r8 prod(:,:,24) = 0._r8 prod(:,:,25) = 0._r8 !-------------------------------------------------------------------- ! ... "independent" production for Implicit species !-------------------------------------------------------------------- else if( class == 4 ) then prod(:,:,128) = 0._r8 prod(:,:,133) =.180_r8*rxt(:,:,55)*y(:,:,12) +rxt(:,:,53)*y(:,:,117) & +2.000_r8*rxt(:,:,1) prod(:,:,131) =rxt(:,:,4)*y(:,:,4) prod(:,:,106) =1.440_r8*rxt(:,:,55)*y(:,:,12) prod(:,:,107) =.380_r8*rxt(:,:,55)*y(:,:,12) +rxt(:,:,53)*y(:,:,117) & + extfrc(:,:,3) prod(:,:,49) = 0._r8 prod(:,:,125) = + extfrc(:,:,1) prod(:,:,135) = + extfrc(:,:,2) prod(:,:,130) =.660_r8*rxt(:,:,55)*y(:,:,12) prod(:,:,127) = 0._r8 prod(:,:,98) = 0._r8 prod(:,:,68) = 0._r8 prod(:,:,55) = 0._r8 prod(:,:,121) =rxt(:,:,54)*y(:,:,12) +rxt(:,:,35)*y(:,:,99) +rxt(:,:,46) & *y(:,:,100) prod(:,:,69) = 0._r8 prod(:,:,50) = 0._r8 prod(:,:,38) = 0._r8 prod(:,:,123) =.180_r8*rxt(:,:,55)*y(:,:,12) prod(:,:,120) =rxt(:,:,54)*y(:,:,12) prod(:,:,132) = 0._r8 prod(:,:,86) = 0._r8 prod(:,:,119) =.050_r8*rxt(:,:,55)*y(:,:,12) prod(:,:,124) =rxt(:,:,35)*y(:,:,99) +3.000_r8*rxt(:,:,38)*y(:,:,101) & +2.000_r8*rxt(:,:,39)*y(:,:,102) +3.000_r8*rxt(:,:,40)*y(:,:,103) & +rxt(:,:,41)*y(:,:,104) +4.000_r8*rxt(:,:,36)*y(:,:,105) & +3.000_r8*rxt(:,:,37)*y(:,:,106) +rxt(:,:,50)*y(:,:,108) & +2.000_r8*rxt(:,:,44)*y(:,:,109) +rxt(:,:,45)*y(:,:,110) & +2.000_r8*rxt(:,:,42)*y(:,:,111) +rxt(:,:,43)*y(:,:,112) prod(:,:,47) = 0._r8 prod(:,:,126) = 0._r8 prod(:,:,45) = 0._r8 prod(:,:,39) = 0._r8 prod(:,:,129) = 0._r8 prod(:,:,100) = 0._r8 prod(:,:,101) = 0._r8 prod(:,:,52) = 0._r8 prod(:,:,122) =rxt(:,:,46)*y(:,:,100) +rxt(:,:,47)*y(:,:,107) +rxt(:,:,50) & *y(:,:,108) +2.000_r8*rxt(:,:,48)*y(:,:,113) +2.000_r8*rxt(:,:,49) & *y(:,:,114) +3.000_r8*rxt(:,:,51)*y(:,:,115) +2.000_r8*rxt(:,:,52) & *y(:,:,116) prod(:,:,134) = 0._r8 prod(:,:,96) = 0._r8 prod(:,:,94) = 0._r8 prod(:,:,80) = 0._r8 prod(:,:,104) = 0._r8 prod(:,:,99) = 0._r8 prod(:,:,92) = 0._r8 prod(:,:,110) = 0._r8 prod(:,:,85) = 0._r8 prod(:,:,81) = 0._r8 prod(:,:,118) = 0._r8 prod(:,:,75) = 0._r8 prod(:,:,74) = 0._r8 prod(:,:,66) = 0._r8 prod(:,:,56) = 0._r8 prod(:,:,62) = 0._r8 prod(:,:,35) = 0._r8 prod(:,:,82) = 0._r8 prod(:,:,36) = 0._r8 prod(:,:,61) = 0._r8 prod(:,:,90) = 0._r8 prod(:,:,88) = 0._r8 prod(:,:,71) = 0._r8 prod(:,:,89) = 0._r8 prod(:,:,57) = 0._r8 prod(:,:,42) = 0._r8 prod(:,:,43) = 0._r8 prod(:,:,77) = 0._r8 prod(:,:,70) = 0._r8 prod(:,:,48) = 0._r8 prod(:,:,103) = 0._r8 prod(:,:,72) = 0._r8 prod(:,:,78) = 0._r8 prod(:,:,87) = 0._r8 prod(:,:,116) = 0._r8 prod(:,:,114) = 0._r8 prod(:,:,111) = 0._r8 prod(:,:,115) = 0._r8 prod(:,:,58) = 0._r8 prod(:,:,117) = 0._r8 prod(:,:,93) = 0._r8 prod(:,:,59) = 0._r8 prod(:,:,84) = 0._r8 prod(:,:,40) = 0._r8 prod(:,:,97) = 0._r8 prod(:,:,64) = 0._r8 prod(:,:,91) = 0._r8 prod(:,:,65) = 0._r8 prod(:,:,76) = 0._r8 prod(:,:,53) = 0._r8 prod(:,:,95) = 0._r8 prod(:,:,109) = 0._r8 prod(:,:,83) = 0._r8 prod(:,:,54) = 0._r8 prod(:,:,44) = 0._r8 prod(:,:,63) = 0._r8 prod(:,:,108) = 0._r8 prod(:,:,112) = 0._r8 prod(:,:,105) = 0._r8 prod(:,:,102) = 0._r8 prod(:,:,113) = 0._r8 prod(:,:,60) = 0._r8 prod(:,:,79) = 0._r8 prod(:,:,46) = 0._r8 prod(:,:,67) = 0._r8 prod(:,:,73) = 0._r8 prod(:,:,41) = + extfrc(:,:,4) prod(:,:,51) = 0._r8 prod(:,:,37) = 0._r8 prod(:,:,1) = 0._r8 prod(:,:,2) = 0._r8 prod(:,:,3) = 0._r8 prod(:,:,4) = + extfrc(:,:,5) prod(:,:,5) = 0._r8 prod(:,:,6) = 0._r8 prod(:,:,7) = 0._r8 prod(:,:,8) = 0._r8 prod(:,:,9) = 0._r8 prod(:,:,10) = + extfrc(:,:,9) prod(:,:,11) = + extfrc(:,:,6) prod(:,:,12) = 0._r8 prod(:,:,13) = 0._r8 prod(:,:,14) = 0._r8 prod(:,:,15) = + extfrc(:,:,10) prod(:,:,16) = + extfrc(:,:,7) prod(:,:,17) = + extfrc(:,:,8) prod(:,:,18) = + extfrc(:,:,11) prod(:,:,19) = 0._r8 prod(:,:,20) = 0._r8 prod(:,:,21) = 0._r8 prod(:,:,22) = 0._r8 prod(:,:,23) = 0._r8 prod(:,:,24) = 0._r8 prod(:,:,25) = 0._r8 prod(:,:,26) = 0._r8 prod(:,:,27) = 0._r8 prod(:,:,28) = 0._r8 prod(:,:,29) = 0._r8 prod(:,:,30) = 0._r8 prod(:,:,31) = 0._r8 prod(:,:,32) = 0._r8 prod(:,:,33) = 0._r8 prod(:,:,34) = 0._r8 end if end subroutine indprd end module mo_indprd