module mo_adjrxt private public :: adjrxt contains subroutine adjrxt( rate, inv, m, ncol ) use ppgrid, only : pver use shr_kind_mod, only : r8 => shr_kind_r8 use chem_mods, only : nfs, rxntot implicit none !-------------------------------------------------------------------- ! ... dummy arguments !-------------------------------------------------------------------- integer, intent(in) :: ncol real(r8), intent(in) :: inv(ncol,pver,nfs) real(r8), intent(in) :: m(ncol,pver) real(r8), intent(inout) :: rate(ncol,pver,rxntot) !-------------------------------------------------------------------- ! ... local variables !-------------------------------------------------------------------- real(r8) :: im(ncol,pver) rate(:,:,115) = rate(:,:,115) * inv(:,:, 1) rate(:,:,117) = rate(:,:,117) * inv(:,:, 1) rate(:,:,120) = rate(:,:,120) * inv(:,:, 2) rate(:,:,126) = rate(:,:,126) * inv(:,:, 2) rate(:,:,128) = rate(:,:,128) * inv(:,:, 2) rate(:,:,160) = rate(:,:,160) * inv(:,:, 1) rate(:,:,169) = rate(:,:,169) * inv(:,:, 1) rate(:,:,177) = rate(:,:,177) * inv(:,:, 1) rate(:,:,187) = rate(:,:,187) * inv(:,:, 1) rate(:,:,191) = rate(:,:,191) * inv(:,:, 1) rate(:,:,193) = rate(:,:,193) * inv(:,:, 1) rate(:,:,194) = rate(:,:,194) * inv(:,:, 1) rate(:,:,195) = rate(:,:,195) * inv(:,:, 1) rate(:,:,201) = rate(:,:,201) * inv(:,:, 1) rate(:,:,203) = rate(:,:,203) * inv(:,:, 1) rate(:,:,217) = rate(:,:,217) * inv(:,:, 1) rate(:,:,221) = rate(:,:,221) * inv(:,:, 1) rate(:,:,222) = rate(:,:,222) * inv(:,:, 1) rate(:,:,238) = rate(:,:,238) * inv(:,:, 1) rate(:,:,265) = rate(:,:,265) * inv(:,:, 1) rate(:,:,280) = rate(:,:,280) * inv(:,:, 1) rate(:,:,281) = rate(:,:,281) * inv(:,:, 1) rate(:,:,283) = rate(:,:,283) * inv(:,:, 1) rate(:,:,285) = rate(:,:,285) * inv(:,:, 1) rate(:,:,300) = rate(:,:,300) * inv(:,:, 1) rate(:,:,308) = rate(:,:,308) * inv(:,:, 1) rate(:,:,310) = rate(:,:,310) * inv(:,:, 1) rate(:,:,353) = rate(:,:,353) * inv(:,:, 1) rate(:,:,354) = rate(:,:,354) * inv(:,:, 1) rate(:,:,355) = rate(:,:,355) * inv(:,:, 1) rate(:,:,450) = rate(:,:,450) * inv(:,:, 2) rate(:,:,460) = rate(:,:,460) * inv(:,:, 2) rate(:,:,115) = rate(:,:,115) * m(:,:) rate(:,:,116) = rate(:,:,116) * m(:,:) rate(:,:,117) = rate(:,:,117) * m(:,:) rate(:,:,118) = rate(:,:,118) * m(:,:) rate(:,:,119) = rate(:,:,119) * m(:,:) rate(:,:,121) = rate(:,:,121) * m(:,:) rate(:,:,122) = rate(:,:,122) * m(:,:) rate(:,:,124) = rate(:,:,124) * m(:,:) rate(:,:,125) = rate(:,:,125) * m(:,:) rate(:,:,129) = rate(:,:,129) * m(:,:) rate(:,:,130) = rate(:,:,130) * m(:,:) rate(:,:,131) = rate(:,:,131) * m(:,:) rate(:,:,132) = rate(:,:,132) * m(:,:) rate(:,:,133) = rate(:,:,133) * m(:,:) rate(:,:,134) = rate(:,:,134) * m(:,:) rate(:,:,135) = rate(:,:,135) * m(:,:) rate(:,:,136) = rate(:,:,136) * m(:,:) rate(:,:,137) = rate(:,:,137) * m(:,:) rate(:,:,138) = rate(:,:,138) * m(:,:) rate(:,:,139) = rate(:,:,139) * m(:,:) rate(:,:,140) = rate(:,:,140) * m(:,:) rate(:,:,141) = rate(:,:,141) * m(:,:) rate(:,:,142) = rate(:,:,142) * m(:,:) rate(:,:,143) = rate(:,:,143) * m(:,:) rate(:,:,144) = rate(:,:,144) * m(:,:) rate(:,:,145) = rate(:,:,145) * m(:,:) rate(:,:,146) = rate(:,:,146) * m(:,:) rate(:,:,147) = rate(:,:,147) * m(:,:) rate(:,:,148) = rate(:,:,148) * m(:,:) rate(:,:,149) = rate(:,:,149) * m(:,:) rate(:,:,150) = rate(:,:,150) * m(:,:) rate(:,:,151) = rate(:,:,151) * m(:,:) rate(:,:,152) = rate(:,:,152) * m(:,:) rate(:,:,153) = rate(:,:,153) * m(:,:) rate(:,:,154) = rate(:,:,154) * m(:,:) rate(:,:,155) = rate(:,:,155) * m(:,:) rate(:,:,156) = rate(:,:,156) * m(:,:) rate(:,:,157) = rate(:,:,157) * m(:,:) rate(:,:,158) = rate(:,:,158) * m(:,:) rate(:,:,159) = rate(:,:,159) * m(:,:) rate(:,:,160) = rate(:,:,160) * m(:,:) rate(:,:,161) = rate(:,:,161) * m(:,:) rate(:,:,162) = rate(:,:,162) * m(:,:) rate(:,:,163) = rate(:,:,163) * m(:,:) rate(:,:,164) = rate(:,:,164) * m(:,:) rate(:,:,165) = rate(:,:,165) * m(:,:) rate(:,:,166) = rate(:,:,166) * m(:,:) rate(:,:,167) = rate(:,:,167) * m(:,:) rate(:,:,168) = rate(:,:,168) * m(:,:) rate(:,:,169) = rate(:,:,169) * m(:,:) rate(:,:,170) = rate(:,:,170) * m(:,:) rate(:,:,171) = rate(:,:,171) * m(:,:) rate(:,:,172) = rate(:,:,172) * m(:,:) rate(:,:,173) = rate(:,:,173) * m(:,:) rate(:,:,174) = rate(:,:,174) * m(:,:) rate(:,:,175) = rate(:,:,175) * m(:,:) rate(:,:,176) = rate(:,:,176) * m(:,:) rate(:,:,177) = rate(:,:,177) * m(:,:) rate(:,:,178) = rate(:,:,178) * m(:,:) rate(:,:,179) = rate(:,:,179) * m(:,:) rate(:,:,180) = rate(:,:,180) * m(:,:) rate(:,:,181) = rate(:,:,181) * m(:,:) rate(:,:,182) = rate(:,:,182) * m(:,:) rate(:,:,183) = rate(:,:,183) * m(:,:) rate(:,:,184) = rate(:,:,184) * m(:,:) rate(:,:,185) = rate(:,:,185) * m(:,:) rate(:,:,186) = rate(:,:,186) * m(:,:) rate(:,:,187) = rate(:,:,187) * m(:,:) rate(:,:,188) = rate(:,:,188) * m(:,:) rate(:,:,189) = rate(:,:,189) * m(:,:) rate(:,:,190) = rate(:,:,190) * m(:,:) rate(:,:,191) = rate(:,:,191) * m(:,:) rate(:,:,192) = rate(:,:,192) * m(:,:) rate(:,:,193) = rate(:,:,193) * m(:,:) rate(:,:,195) = rate(:,:,195) * m(:,:) rate(:,:,196) = rate(:,:,196) * m(:,:) rate(:,:,197) = rate(:,:,197) * m(:,:) rate(:,:,198) = rate(:,:,198) * m(:,:) rate(:,:,199) = rate(:,:,199) * m(:,:) rate(:,:,200) = rate(:,:,200) * m(:,:) rate(:,:,201) = rate(:,:,201) * m(:,:) rate(:,:,202) = rate(:,:,202) * m(:,:) rate(:,:,204) = rate(:,:,204) * m(:,:) rate(:,:,205) = rate(:,:,205) * m(:,:) rate(:,:,206) = rate(:,:,206) * m(:,:) rate(:,:,207) = rate(:,:,207) * m(:,:) rate(:,:,208) = rate(:,:,208) * m(:,:) rate(:,:,209) = rate(:,:,209) * m(:,:) rate(:,:,210) = rate(:,:,210) * m(:,:) rate(:,:,211) = rate(:,:,211) * m(:,:) rate(:,:,212) = rate(:,:,212) * m(:,:) rate(:,:,213) = rate(:,:,213) * m(:,:) rate(:,:,214) = rate(:,:,214) * m(:,:) rate(:,:,215) = rate(:,:,215) * m(:,:) rate(:,:,216) = rate(:,:,216) * m(:,:) rate(:,:,217) = rate(:,:,217) * m(:,:) rate(:,:,218) = rate(:,:,218) * m(:,:) rate(:,:,219) = rate(:,:,219) * m(:,:) rate(:,:,220) = rate(:,:,220) * m(:,:) rate(:,:,221) = rate(:,:,221) * m(:,:) rate(:,:,223) = rate(:,:,223) * m(:,:) rate(:,:,224) = rate(:,:,224) * m(:,:) rate(:,:,225) = rate(:,:,225) * m(:,:) rate(:,:,226) = rate(:,:,226) * m(:,:) rate(:,:,227) = rate(:,:,227) * m(:,:) rate(:,:,228) = rate(:,:,228) * m(:,:) rate(:,:,229) = rate(:,:,229) * m(:,:) rate(:,:,230) = rate(:,:,230) * m(:,:) rate(:,:,231) = rate(:,:,231) * m(:,:) rate(:,:,232) = rate(:,:,232) * m(:,:) rate(:,:,233) = rate(:,:,233) * m(:,:) rate(:,:,234) = rate(:,:,234) * m(:,:) rate(:,:,235) = rate(:,:,235) * m(:,:) rate(:,:,236) = rate(:,:,236) * m(:,:) rate(:,:,237) = rate(:,:,237) * m(:,:) rate(:,:,238) = rate(:,:,238) * m(:,:) rate(:,:,239) = rate(:,:,239) * m(:,:) rate(:,:,240) = rate(:,:,240) * m(:,:) rate(:,:,241) = rate(:,:,241) * m(:,:) rate(:,:,242) = rate(:,:,242) * m(:,:) rate(:,:,243) = rate(:,:,243) * m(:,:) rate(:,:,244) = rate(:,:,244) * m(:,:) rate(:,:,245) = rate(:,:,245) * m(:,:) rate(:,:,246) = rate(:,:,246) * m(:,:) rate(:,:,247) = rate(:,:,247) * m(:,:) rate(:,:,248) = rate(:,:,248) * m(:,:) rate(:,:,249) = rate(:,:,249) * m(:,:) rate(:,:,250) = rate(:,:,250) * m(:,:) rate(:,:,251) = rate(:,:,251) * m(:,:) rate(:,:,252) = rate(:,:,252) * m(:,:) rate(:,:,253) = rate(:,:,253) * m(:,:) rate(:,:,254) = rate(:,:,254) * m(:,:) rate(:,:,255) = rate(:,:,255) * m(:,:) rate(:,:,256) = rate(:,:,256) * m(:,:) rate(:,:,257) = rate(:,:,257) * m(:,:) rate(:,:,258) = rate(:,:,258) * m(:,:) rate(:,:,259) = rate(:,:,259) * m(:,:) rate(:,:,260) = rate(:,:,260) * m(:,:) rate(:,:,261) = rate(:,:,261) * m(:,:) rate(:,:,262) = rate(:,:,262) * m(:,:) rate(:,:,263) = rate(:,:,263) * m(:,:) rate(:,:,264) = rate(:,:,264) * m(:,:) rate(:,:,265) = rate(:,:,265) * m(:,:) rate(:,:,266) = rate(:,:,266) * m(:,:) rate(:,:,267) = rate(:,:,267) * m(:,:) rate(:,:,268) = rate(:,:,268) * m(:,:) rate(:,:,269) = rate(:,:,269) * m(:,:) rate(:,:,270) = rate(:,:,270) * m(:,:) rate(:,:,271) = rate(:,:,271) * m(:,:) rate(:,:,272) = rate(:,:,272) * m(:,:) rate(:,:,273) = rate(:,:,273) * m(:,:) rate(:,:,274) = rate(:,:,274) * m(:,:) rate(:,:,275) = rate(:,:,275) * m(:,:) rate(:,:,276) = rate(:,:,276) * m(:,:) rate(:,:,278) = rate(:,:,278) * m(:,:) rate(:,:,279) = rate(:,:,279) * m(:,:) rate(:,:,280) = rate(:,:,280) * m(:,:) rate(:,:,281) = rate(:,:,281) * m(:,:) rate(:,:,282) = rate(:,:,282) * m(:,:) rate(:,:,283) = rate(:,:,283) * m(:,:) rate(:,:,284) = rate(:,:,284) * m(:,:) rate(:,:,285) = rate(:,:,285) * m(:,:) rate(:,:,286) = rate(:,:,286) * m(:,:) rate(:,:,287) = rate(:,:,287) * m(:,:) rate(:,:,288) = rate(:,:,288) * m(:,:) rate(:,:,290) = rate(:,:,290) * m(:,:) rate(:,:,291) = rate(:,:,291) * m(:,:) rate(:,:,292) = rate(:,:,292) * m(:,:) rate(:,:,293) = rate(:,:,293) * m(:,:) rate(:,:,294) = rate(:,:,294) * m(:,:) rate(:,:,295) = rate(:,:,295) * m(:,:) rate(:,:,296) = rate(:,:,296) * m(:,:) rate(:,:,297) = rate(:,:,297) * m(:,:) rate(:,:,298) = rate(:,:,298) * m(:,:) rate(:,:,299) = rate(:,:,299) * m(:,:) rate(:,:,300) = rate(:,:,300) * m(:,:) rate(:,:,301) = rate(:,:,301) * m(:,:) rate(:,:,302) = rate(:,:,302) * m(:,:) rate(:,:,303) = rate(:,:,303) * m(:,:) rate(:,:,304) = rate(:,:,304) * m(:,:) rate(:,:,305) = rate(:,:,305) * m(:,:) rate(:,:,306) = rate(:,:,306) * m(:,:) rate(:,:,307) = rate(:,:,307) * m(:,:) rate(:,:,309) = rate(:,:,309) * m(:,:) rate(:,:,310) = rate(:,:,310) * m(:,:) rate(:,:,311) = rate(:,:,311) * m(:,:) rate(:,:,312) = rate(:,:,312) * m(:,:) rate(:,:,313) = rate(:,:,313) * m(:,:) rate(:,:,314) = rate(:,:,314) * m(:,:) rate(:,:,315) = rate(:,:,315) * m(:,:) rate(:,:,316) = rate(:,:,316) * m(:,:) rate(:,:,317) = rate(:,:,317) * m(:,:) rate(:,:,318) = rate(:,:,318) * m(:,:) rate(:,:,319) = rate(:,:,319) * m(:,:) rate(:,:,320) = rate(:,:,320) * m(:,:) rate(:,:,321) = rate(:,:,321) * m(:,:) rate(:,:,322) = rate(:,:,322) * m(:,:) rate(:,:,323) = rate(:,:,323) * m(:,:) rate(:,:,324) = rate(:,:,324) * m(:,:) rate(:,:,325) = rate(:,:,325) * m(:,:) rate(:,:,326) = rate(:,:,326) * m(:,:) rate(:,:,327) = rate(:,:,327) * m(:,:) rate(:,:,328) = rate(:,:,328) * m(:,:) rate(:,:,329) = rate(:,:,329) * m(:,:) rate(:,:,330) = rate(:,:,330) * m(:,:) rate(:,:,331) = rate(:,:,331) * m(:,:) rate(:,:,332) = rate(:,:,332) * m(:,:) rate(:,:,333) = rate(:,:,333) * m(:,:) rate(:,:,334) = rate(:,:,334) * m(:,:) rate(:,:,335) = rate(:,:,335) * m(:,:) rate(:,:,336) = rate(:,:,336) * m(:,:) rate(:,:,337) = rate(:,:,337) * m(:,:) rate(:,:,338) = rate(:,:,338) * m(:,:) rate(:,:,339) = rate(:,:,339) * m(:,:) rate(:,:,340) = rate(:,:,340) * m(:,:) rate(:,:,341) = rate(:,:,341) * m(:,:) rate(:,:,342) = rate(:,:,342) * m(:,:) rate(:,:,343) = rate(:,:,343) * m(:,:) rate(:,:,344) = rate(:,:,344) * m(:,:) rate(:,:,345) = rate(:,:,345) * m(:,:) rate(:,:,346) = rate(:,:,346) * m(:,:) rate(:,:,347) = rate(:,:,347) * m(:,:) rate(:,:,348) = rate(:,:,348) * m(:,:) rate(:,:,349) = rate(:,:,349) * m(:,:) rate(:,:,350) = rate(:,:,350) * m(:,:) rate(:,:,351) = rate(:,:,351) * m(:,:) rate(:,:,352) = rate(:,:,352) * m(:,:) rate(:,:,353) = rate(:,:,353) * m(:,:) rate(:,:,355) = rate(:,:,355) * m(:,:) rate(:,:,356) = rate(:,:,356) * m(:,:) rate(:,:,357) = rate(:,:,357) * m(:,:) rate(:,:,358) = rate(:,:,358) * m(:,:) rate(:,:,359) = rate(:,:,359) * m(:,:) rate(:,:,360) = rate(:,:,360) * m(:,:) rate(:,:,361) = rate(:,:,361) * m(:,:) rate(:,:,362) = rate(:,:,362) * m(:,:) rate(:,:,363) = rate(:,:,363) * m(:,:) rate(:,:,364) = rate(:,:,364) * m(:,:) rate(:,:,365) = rate(:,:,365) * m(:,:) rate(:,:,366) = rate(:,:,366) * m(:,:) rate(:,:,367) = rate(:,:,367) * m(:,:) rate(:,:,368) = rate(:,:,368) * m(:,:) rate(:,:,369) = rate(:,:,369) * m(:,:) rate(:,:,370) = rate(:,:,370) * m(:,:) rate(:,:,371) = rate(:,:,371) * m(:,:) rate(:,:,372) = rate(:,:,372) * m(:,:) rate(:,:,373) = rate(:,:,373) * m(:,:) rate(:,:,374) = rate(:,:,374) * m(:,:) rate(:,:,375) = rate(:,:,375) * m(:,:) rate(:,:,376) = rate(:,:,376) * m(:,:) rate(:,:,377) = rate(:,:,377) * m(:,:) rate(:,:,378) = rate(:,:,378) * m(:,:) rate(:,:,379) = rate(:,:,379) * m(:,:) rate(:,:,380) = rate(:,:,380) * m(:,:) rate(:,:,381) = rate(:,:,381) * m(:,:) rate(:,:,382) = rate(:,:,382) * m(:,:) rate(:,:,383) = rate(:,:,383) * m(:,:) rate(:,:,384) = rate(:,:,384) * m(:,:) rate(:,:,385) = rate(:,:,385) * m(:,:) rate(:,:,386) = rate(:,:,386) * m(:,:) rate(:,:,387) = rate(:,:,387) * m(:,:) rate(:,:,388) = rate(:,:,388) * m(:,:) rate(:,:,389) = rate(:,:,389) * m(:,:) rate(:,:,390) = rate(:,:,390) * m(:,:) rate(:,:,391) = rate(:,:,391) * m(:,:) rate(:,:,392) = rate(:,:,392) * m(:,:) rate(:,:,393) = rate(:,:,393) * m(:,:) rate(:,:,394) = rate(:,:,394) * m(:,:) rate(:,:,395) = rate(:,:,395) * m(:,:) rate(:,:,396) = rate(:,:,396) * m(:,:) rate(:,:,397) = rate(:,:,397) * m(:,:) rate(:,:,398) = rate(:,:,398) * m(:,:) rate(:,:,399) = rate(:,:,399) * m(:,:) rate(:,:,404) = rate(:,:,404) * m(:,:) rate(:,:,405) = rate(:,:,405) * m(:,:) rate(:,:,406) = rate(:,:,406) * m(:,:) rate(:,:,407) = rate(:,:,407) * m(:,:) rate(:,:,408) = rate(:,:,408) * m(:,:) rate(:,:,435) = rate(:,:,435) * m(:,:) rate(:,:,436) = rate(:,:,436) * m(:,:) rate(:,:,437) = rate(:,:,437) * m(:,:) rate(:,:,440) = rate(:,:,440) * m(:,:) rate(:,:,441) = rate(:,:,441) * m(:,:) rate(:,:,446) = rate(:,:,446) * m(:,:) rate(:,:,447) = rate(:,:,447) * m(:,:) rate(:,:,448) = rate(:,:,448) * m(:,:) rate(:,:,449) = rate(:,:,449) * m(:,:) rate(:,:,451) = rate(:,:,451) * m(:,:) rate(:,:,452) = rate(:,:,452) * m(:,:) rate(:,:,453) = rate(:,:,453) * m(:,:) rate(:,:,454) = rate(:,:,454) * m(:,:) rate(:,:,455) = rate(:,:,455) * m(:,:) rate(:,:,456) = rate(:,:,456) * m(:,:) rate(:,:,457) = rate(:,:,457) * m(:,:) rate(:,:,458) = rate(:,:,458) * m(:,:) rate(:,:,459) = rate(:,:,459) * m(:,:) rate(:,:,461) = rate(:,:,461) * m(:,:) rate(:,:,462) = rate(:,:,462) * m(:,:) rate(:,:,463) = rate(:,:,463) * m(:,:) end subroutine adjrxt end module mo_adjrxt