! hereinafter the Contractor, under Contract No. DE-AC05-76RL0 1830 with ! the Department of Energy (DOE). NEITHER THE GOVERNMENT NOR THE ! CONTRACTOR MAKES ANY WARRANTY, EXPRESS OR IMPLIED, OR ASSUMES ANY ! LIABILITY FOR THE USE OF THIS SOFTWARE. ! ! Aerosol Option: MOSAIC (Model for Simulating Aerosol Interactions & ! Chemistry) ! * Primary investigator: Rahul A. Zaveri ! * Co-investigator: Richard C. Easter, William I. Gustafson Jr. ! Last update: February 2009 ! ! Contacts: ! Rahul A. Zaveri, PhD Jerome D. Fast, PhD ! Senior Research Scientist Staff Scientist ! Pacific Northwest National Laboratory Pacific Northwest National Laboratory ! P.O. Box 999, MSIN K9-30 P.O. Box 999, MSIN K9-30 ! Richland, WA 99352 Richland, WA, 99352 ! Phone: (509) 372-6159 Phone: (509) 372-6116 ! Email: Rahul.Zaveri@pnl.gov Email: Jerome.Fast@pnl.gov ! ! Please report any bugs or problems to Rahul Zaveri, the primary author ! of the code, or Jerome Fast, the WRF-chem implementation team leader ! for PNNL. ! ! Terms of Use: ! 1) Users are requested to consult the primary author prior to ! modifying the MOSAIC code or incorporating it or its submodules in ! another code. This is meant to ensure that the any linkages and/or ! assumptions will not adversely affect the operation of MOSAIC. ! 2) The MOSAIC source code is intended for research and educational ! purposes. Users are requested to contact the primary author ! regarding the use of MOSAIC code for any commercial application. ! 3) Users preparing publications resulting from the usage of MOSAIC are ! requested to cite one or more of the references below (depending on ! the application) for proper acknowledgement. ! ! References for MOSAIC (The first one is the primary MOSAIC reference): ! * Zaveri R.A., R.C. Easter, J.D. Fast, and L.K. Peters (2008), Model ! for Simulating Aerosol Interactions and Chemistry (MOSAIC), J. ! Geophys. Res., 113, D13204, doi:10.1029/2007JD008782. ! * Zaveri R.A., R.C. Easter, and A.S. Wexler (2005), A new method for ! multi-component activity coefficients of electrolytes in aqueous ! atmospheric aerosols, J. Geophys. Res., 110, D02201, ! doi:10.1029/2004JD004681. ! * Zaveri R.A., R.C. Easter, and L.K. Peters (2005), A computationally ! efficient multicomponent equilibrium solver for aerosols (MESA), J. ! Geophys. Res., 110, D24203, doi:24203, doi:10.1029/2004JD005618. ! ! Reference for volatility basis set SOA code: ! * Shrivastava, M., J. Fast, R. Easter, W. I. Gustafson Jr., R. A. ! Zaveri, J. L. Jimenez, P. Saide, and A. Hodzic, 2010: Modeling ! organic aerosols in a megacity: comparison of simple and complex ! representations of the volatility basis set approach. Atmos. Chem. ! Phys., in review. ! ! Reference for the alternative SOA treatment that is less expensive to run ! and is used in the package mozart_mosaic_4bin_vbs0 ! SOA is treated using a simplified approach for anthropogenic and biomass burning species based on Hodzic and Jimenez, GMD, 2011 ! Hodzic, A. and Jimenez, J. L.: Modeling anthropogenically controlled secondary organic aerosols in a megacity: ! a simplified framework for global and climate models, Geosci. Model Dev., 4, 901-917, doi:10.5194/gmd-4-901-2011, 2011. ! ! References for implementation of the aerosol optical properties and ! direct effect in WRF-Chem: ! * Fast, J.D., W.I. Gustafson Jr., R.C. Easter, R.A. Zaveri, J.C. ! Barnard, E.G. Chapman, G.A. Grell, and S.E. Peckham (2005), Evolution ! of ozone, particulates, and aerosol direct radiative forcing in the ! vicinity of Houston using a fully-coupled meteorology-chemistry- ! aerosol model, J. Geophys. Res., 111, D21305, ! doi:10.1029/2005JD006721. ! * Barnard, J.C., J.D. Fast, G. Paredes-Miranda, W.P. Arnott, and ! A. Laskin (2010), Technical note: evaluation of the WRF-Chem ! "aerosol chemical to aerosol optical properties" module using data ! from the MILAGRO campaign, Atmos. Chem. Phys., 10, 7325-7340, ! doi:10.5194/acp-10-7325-2010. ! ! References for implementation of the aerosol-cloud interactions and ! indirect effects in WRF-Chem: ! * Gustafson, W. I., E. G. Chapman, S. J. Ghan, R. C. Easter, and J. D. ! Fast (2007), Impact on modeled cloud characteristics due to ! simplified treatment of uniform cloud condensation nuclei during ! NEAQS 2004, Geophys. Res. Lett., 34, L19809, L19809, ! doi:10.1029/2007GL0300321. ! * Chapman, E.G., W. I. Gustafson Jr., R. C. Easter, J. C. Barnard, ! S. J. Ghan, M. S. Pekour, and J. D. Fast (2009): Coupling aerosol- ! cloud-radiative processes in the WRF-Chem model: Investigating the ! radiative impact of elevated point sources, Atmos. Chem. Phys., 9, ! 945-964, www.atmos-chem-phys.net/9/945/2009/. ! ! Contact Jerome Fast for updates on the status of manuscripts under ! review. ! ! Additional information: ! * www.pnl.gov/atmospheric/research/wrf-chem ! ! Support: ! Funding for developing and evaluating MOSAIC was provided by the U.S. ! Department of Energy under the auspices of Atmospheric Science Program ! of the Office of Biological and Environmental Research, the NASA Earth ! Sciences Enterprise under grant NAGW 3367, and PNNL Laboratory Directed ! Research and Development program. !************************************************************************ module module_mosaic_sumpm ! ! *** NOTE - when the cw species are NOT in the registry, then ! then the p_xxx_cwnn variables are not in module_state_description, ! and the following cpp directive should be commented out ! #define cw_species_are_in_registry contains !----------------------------------------------------------------------- subroutine sum_pm_mosaic ( & alt, chem, & pm2_5_dry, pm2_5_water, pm2_5_dry_ec, pm10, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: pm2_5_dry,pm2_5_water,pm2_5_dry_ec,pm10 REAL :: mass INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte ! ! Sum over bins with center diameter < 2.5e-4 cm for pm2_5_dry, ! pm2_5_dry_ec, and pm2_5_water. All bins go into pm10 ! pm2_5_dry(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_water(its:imax,kts:kmax,jts:jmax) = 0. pm10(its:imax,kts:kmax,jts:jmax) = 0. do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) if (dcen_sect(n,itype) .le. 2.5e-4) then do j=jts,jmax do k=kts,kmax do i=its,imax mass = chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(n,itype,iphase)) pm2_5_dry(i,k,j) = pm2_5_dry(i,k,j) + mass pm2_5_dry_ec(i,k,j) = pm2_5_dry_ec(i,k,j) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) pm2_5_water(i,k,j) = pm2_5_water(i,k,j) & + chem(i,k,j,waterptr_aer(n,itype)) pm10(i,k,j) = pm10(i,k,j) + mass enddo enddo enddo else do j=jts,jmax do k=kts,kmax do i=its,imax pm10(i,k,j) = pm10(i,k,j) & + chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(n,itype,iphase)) enddo enddo enddo endif enddo ! size enddo ! type enddo ! phase !Convert the units from mixing ratio to concentration (ug m^-3) pm2_5_dry(its:imax,kts:kmax,jts:jmax) = pm2_5_dry(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_water(its:imax,kts:kmax,jts:jmax) = pm2_5_water(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm10(its:imax,kts:kmax,jts:jmax) = pm10(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) end subroutine sum_pm_mosaic !----------------------------------------------------------------------- !----------------------------------------------------------------------- subroutine sum_pm_mosaic_vbs2 ( & alt, chem, & pm2_5_dry, pm2_5_water, pm2_5_dry_ec, pm10, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: pm2_5_dry,pm2_5_water,pm2_5_dry_ec,pm10 REAL :: mass INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte ! ! Sum over bins with center diameter < 2.5e-4 cm for pm2_5_dry, ! pm2_5_dry_ec, and pm2_5_water. All bins go into pm10 ! pm2_5_dry(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_water(its:imax,kts:kmax,jts:jmax) = 0. pm10(its:imax,kts:kmax,jts:jmax) = 0. do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) if (dcen_sect(n,itype) .le. 2.5e-4) then do j=jts,jmax do k=kts,kmax do i=its,imax mass = chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) pm2_5_dry(i,k,j) = pm2_5_dry(i,k,j) + mass pm2_5_dry_ec(i,k,j) = pm2_5_dry_ec(i,k,j) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) pm2_5_water(i,k,j) = pm2_5_water(i,k,j) & + chem(i,k,j,waterptr_aer(n,itype)) pm10(i,k,j) = pm10(i,k,j) + mass enddo enddo enddo else do j=jts,jmax do k=kts,kmax do i=its,imax pm10(i,k,j) = pm10(i,k,j) & + chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) enddo enddo enddo endif enddo ! size enddo ! type enddo ! phase !Convert the units from mixing ratio to concentration (ug m^-3) pm2_5_dry(its:imax,kts:kmax,jts:jmax) = pm2_5_dry(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_water(its:imax,kts:kmax,jts:jmax) = pm2_5_water(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm10(its:imax,kts:kmax,jts:jmax) = pm10(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) end subroutine sum_pm_mosaic_vbs2 !----------------------------------------------------------------------- !----------------------------------------------------------------------- subroutine sum_pm_mosaic_vbs0 ( & alt, chem, & pm2_5_dry, pm2_5_water, pm2_5_dry_ec, pm10, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: pm2_5_dry,pm2_5_water,pm2_5_dry_ec,pm10 REAL :: mass INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte ! ! ! Sum over bins with center diameter < 2.5e-4 cm for pm2_5_dry, ! pm2_5_dry_ec, and pm2_5_water. All bins go into pm10 ! pm2_5_dry(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_water(its:imax,kts:kmax,jts:jmax) = 0. pm10(its:imax,kts:kmax,jts:jmax) = 0. do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) if (dcen_sect(n,itype) .le. 2.5e-4) then do j=jts,jmax do k=kts,kmax do i=its,imax mass = chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_smpa_aer(n,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) ! SMPA and SMPBB do not participate to pm2.5 mass pm2_5_dry(i,k,j) = pm2_5_dry(i,k,j) + mass pm2_5_dry_ec(i,k,j) = pm2_5_dry_ec(i,k,j) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) pm2_5_water(i,k,j) = pm2_5_water(i,k,j) & + chem(i,k,j,waterptr_aer(n,itype)) pm10(i,k,j) = pm10(i,k,j) + mass enddo enddo enddo else do j=jts,jmax do k=kts,kmax do i=its,imax pm10(i,k,j) = pm10(i,k,j) & + chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_smpa_aer(n,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) ! SMPA and SMPBB do not participate to pm10 mass enddo enddo enddo endif enddo ! size enddo ! type enddo ! phase !Convert the units from mixing ratio to concentration (ug m^-3) pm2_5_dry(its:imax,kts:kmax,jts:jmax) = pm2_5_dry(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_water(its:imax,kts:kmax,jts:jmax) = pm2_5_water(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm10(its:imax,kts:kmax,jts:jmax) = pm10(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) end subroutine sum_pm_mosaic_vbs0 !----------------------------------------------------------------------- !----------------------------------------------------------------------- subroutine sum_pm_mosaic_vbs4 ( & alt, chem, & pm2_5_dry, pm2_5_water, pm2_5_dry_ec, pm10, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: pm2_5_dry,pm2_5_water,pm2_5_dry_ec,pm10 REAL :: mass INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte ! ! ! Sum over bins with center diameter < 2.5e-4 cm for pm2_5_dry, ! pm2_5_dry_ec, and pm2_5_water. All bins go into pm10 ! pm2_5_dry(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = 0. pm2_5_water(its:imax,kts:kmax,jts:jmax) = 0. pm10(its:imax,kts:kmax,jts:jmax) = 0. do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) if (dcen_sect(n,itype) .le. 2.5e-4) then do j=jts,jmax do k=kts,kmax do i=its,imax mass = chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoaX_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(n,itype,iphase)) ! SMPA and SMPBB do not participate to pm2.5 mass pm2_5_dry(i,k,j) = pm2_5_dry(i,k,j) + mass pm2_5_dry_ec(i,k,j) = pm2_5_dry_ec(i,k,j) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) pm2_5_water(i,k,j) = pm2_5_water(i,k,j) & + chem(i,k,j,waterptr_aer(n,itype)) pm10(i,k,j) = pm10(i,k,j) + mass enddo enddo enddo else do j=jts,jmax do k=kts,kmax do i=its,imax pm10(i,k,j) = pm10(i,k,j) & + chem(i,k,j,lptr_so4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_no3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_cl_aer(n,itype,iphase)) & + chem(i,k,j,lptr_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_na_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oin_aer(n,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(n,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoaX_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(n,itype,iphase)) ! SMPA and SMPBB do not participate to pm10 mass enddo enddo enddo endif enddo ! size enddo ! type enddo ! phase !Convert the units from mixing ratio to concentration (ug m^-3) pm2_5_dry(its:imax,kts:kmax,jts:jmax) = pm2_5_dry(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) = pm2_5_dry_ec(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm2_5_water(its:imax,kts:kmax,jts:jmax) = pm2_5_water(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) pm10(its:imax,kts:kmax,jts:jmax) = pm10(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) end subroutine sum_pm_mosaic_vbs4 !----------------------------------------------------------------------- subroutine sum_vbs0 ( aero_diag_opt, & alt, chem, & hoa_a01,hoa_a02,hoa_a03,hoa_a04, & bboa_a01,bboa_a02,bboa_a03,bboa_a04, & soa_a01,soa_a02,soa_a03,soa_a04, & bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04, & biog_a01,biog_a02,biog_a03,biog_a04, & asmpsoa_a01,asmpsoa_a02,asmpsoa_a03,asmpsoa_a04, & arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4, & smpa_v1,smpbb_v1, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: aero_diag_opt INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: hoa_a01,hoa_a02,hoa_a03,hoa_a04, & bboa_a01,bboa_a02,bboa_a03,bboa_a04, & soa_a01,soa_a02,soa_a03,soa_a04, & bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04, & biog_a01,biog_a02,biog_a03,biog_a04, & asmpsoa_a01,asmpsoa_a02,asmpsoa_a03,asmpsoa_a04, & arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4,smpa_v1,smpbb_v1 INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte totoa_a01(its:imax,kts:kmax,jts:jmax) = 0. totoa_a02(its:imax,kts:kmax,jts:jmax) = 0. totoa_a03(its:imax,kts:kmax,jts:jmax) = 0. totoa_a04(its:imax,kts:kmax,jts:jmax) = 0. if( aero_diag_opt > 0 ) then hoa_a01(its:imax,kts:kmax,jts:jmax) = 0. soa_a01(its:imax,kts:kmax,jts:jmax) = 0. bboa_a01(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a01(its:imax,kts:kmax,jts:jmax) = 0. biog_a01(its:imax,kts:kmax,jts:jmax) = 0. asmpsoa_a01(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a01(its:imax,kts:kmax,jts:jmax) = 0. hoa_a02(its:imax,kts:kmax,jts:jmax) = 0. soa_a02(its:imax,kts:kmax,jts:jmax) = 0. bboa_a02(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a02(its:imax,kts:kmax,jts:jmax) = 0. asmpsoa_a02(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a02(its:imax,kts:kmax,jts:jmax) = 0. biog_a02(its:imax,kts:kmax,jts:jmax) = 0. hoa_a03(its:imax,kts:kmax,jts:jmax) = 0. soa_a03(its:imax,kts:kmax,jts:jmax) = 0. bboa_a03(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a03(its:imax,kts:kmax,jts:jmax) = 0. asmpsoa_a03(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a03(its:imax,kts:kmax,jts:jmax) = 0. biog_a03(its:imax,kts:kmax,jts:jmax) = 0. hoa_a04(its:imax,kts:kmax,jts:jmax) = 0. soa_a04(its:imax,kts:kmax,jts:jmax) = 0. bboa_a04(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a04(its:imax,kts:kmax,jts:jmax) = 0. asmpsoa_a04(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a04(its:imax,kts:kmax,jts:jmax) = 0. biog_a04(its:imax,kts:kmax,jts:jmax) = 0. ! Species to calculate O:C ratios biog_v1(its:imax,kts:kmax,jts:jmax) = 0. biog_v2(its:imax,kts:kmax,jts:jmax) = 0. biog_v3(its:imax,kts:kmax,jts:jmax) = 0. biog_v4(its:imax,kts:kmax,jts:jmax) = 0. smpa_v1(its:imax,kts:kmax,jts:jmax) = 0. smpbb_v1(its:imax,kts:kmax,jts:jmax) = 0. ant_v1(its:imax,kts:kmax,jts:jmax) = 0. ant_v2(its:imax,kts:kmax,jts:jmax) = 0. ant_v3(its:imax,kts:kmax,jts:jmax) = 0. ant_v4(its:imax,kts:kmax,jts:jmax) = 0. endif do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax totoa_a01(i,k,j)= totoa_a01(i,k,j) & + chem(i,k,j,lptr_smpa_aer(1,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(1,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(1,itype,iphase)) totoa_a02(i,k,j)= totoa_a02(i,k,j) & + chem(i,k,j,lptr_oc_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(2,itype,iphase)) & + chem(i,k,j,lptr_smpa_aer(2,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) totoa_a03(i,k,j)= totoa_a03(i,k,j) & + chem(i,k,j,lptr_smpa_aer(3,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(3,itype,iphase)) totoa_a04(i,k,j)= totoa_a04(i,k,j) & + chem(i,k,j,lptr_smpa_aer(4,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase if( aero_diag_opt > 0 ) then do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) !The 4th bin is 2.5-10um and outside the AMS measurements do j=jts,jmax do k=kts,kmax do i=its,imax biog_v1(i,k,j)= biog_v1(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) ant_v1(i,k,j)= ant_v1(i,k,j) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) smpa_v1(i,k,j)= smpa_v1(i,k,j) & + chem(i,k,j,lptr_smpa_aer(n,itype,iphase)) smpbb_v1(i,k,j)= smpbb_v1(i,k,j) & + chem(i,k,j,lptr_smpbb_aer(n,itype,iphase)) enddo enddo enddo enddo enddo enddo biog_v1(its:imax,kts:kmax,jts:jmax) = biog_v1(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) ant_v1(its:imax,kts:kmax,jts:jmax) = ant_v1(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) smpa_v1(its:imax,kts:kmax,jts:jmax) = smpa_v1(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) smpbb_v1(its:imax,kts:kmax,jts:jmax) = smpbb_v1(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax ! missing summation! (hoa = hoa + x ...) ! hoa_a01(i,k,j)= chem(i,k,j,lptr_oc_aer(1,itype,iphase)) ! ! ! soa_a01(i,k,j)= (chem(i,k,j,lptr_smpa_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase))) ! ! ! bbsoa_a01(i,k,j)= chem(i,k,j,lptr_smpbb_aer(1,itype,iphase)) ! ! ! biog_a01(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase))) ! ! asmpsoa_a01(i,k,j)= chem(i,k,j,lptr_smpa_aer(1,itype,iphase)) ! ! totoa_a01(i,k,j)= ( & ! + chem(i,k,j,lptr_smpa_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) & ! + chem(i,k,j,lptr_oc_aer(1,itype,iphase))) hoa_a01(i,k,j)= hoa_a01(i,k,j) & + chem(i,k,j,lptr_oc_aer(1,itype,iphase)) soa_a01(i,k,j)= soa_a01(i,k,j) & + chem(i,k,j,lptr_smpa_aer(1,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) bbsoa_a01(i,k,j)= bbsoa_a01(i,k,j) & + chem(i,k,j,lptr_smpbb_aer(1,itype,iphase)) biog_a01(i,k,j)= biog_a01(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) asmpsoa_a01(i,k,j)= asmpsoa_a01(i,k,j) & + chem(i,k,j,lptr_smpa_aer(1,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax ! missing summation! (hoa = hoa + x ...) ! hoa_a02(i,k,j)= (chem(i,k,j,lptr_oc_aer(2,itype,iphase))) ! ! ! soa_a02(i,k,j)= (chem(i,k,j,lptr_smpa_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase))) ! ! ! bbsoa_a02(i,k,j)= (chem(i,k,j,lptr_smpbb_aer(2,itype,iphase))) ! ! ! biog_a02(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase))) ! ! ! asmpsoa_a02(i,k,j)= chem(i,k,j,lptr_smpa_aer(2,itype,iphase)) ! ! totoa_a02(i,k,j)= ( chem(i,k,j,lptr_oc_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_smpa_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase))) hoa_a02(i,k,j)= hoa_a02(i,k,j) & + chem(i,k,j,lptr_oc_aer(2,itype,iphase)) soa_a02(i,k,j)= soa_a02(i,k,j) & + chem(i,k,j,lptr_smpa_aer(2,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) bbsoa_a02(i,k,j)= bbsoa_a02(i,k,j) & + chem(i,k,j,lptr_smpbb_aer(2,itype,iphase)) biog_a02(i,k,j)= biog_a02(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) asmpsoa_a02(i,k,j)= asmpsoa_a02(i,k,j) & + chem(i,k,j,lptr_smpa_aer(2,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax ! missing summation! (hoa = hoa + x ...) ! hoa_a03(i,k,j)= (chem(i,k,j,lptr_oc_aer(3,itype,iphase))) ! ! soa_a03(i,k,j)= (chem(i,k,j,lptr_smpa_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase))) ! ! ! bbsoa_a03(i,k,j)= (chem(i,k,j,lptr_smpbb_aer(3,itype,iphase))) ! ! biog_a03(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase))) ! ! asmpsoa_a03(i,k,j)= chem(i,k,j,lptr_smpa_aer(3,itype,iphase)) ! ! ! totoa_a03(i,k,j)= ( & ! + chem(i,k,j,lptr_smpa_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) & ! + chem(i,k,j,lptr_oc_aer(3,itype,iphase))) hoa_a03(i,k,j)= hoa_a03(i,k,j) & + (chem(i,k,j,lptr_oc_aer(3,itype,iphase))) soa_a03(i,k,j)= soa_a03(i,k,j) & + chem(i,k,j,lptr_smpa_aer(3,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) bbsoa_a03(i,k,j)= bbsoa_a03(i,k,j) & + chem(i,k,j,lptr_smpbb_aer(3,itype,iphase)) biog_a03(i,k,j)= biog_a03(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) asmpsoa_a03(i,k,j)= asmpsoa_a03(i,k,j) & + chem(i,k,j,lptr_smpa_aer(3,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax ! missing summation! (hoa = hoa + x ...) ! hoa_a04(i,k,j)= (chem(i,k,j,lptr_oc_aer(4,itype,iphase))) ! ! ! soa_a04(i,k,j)= (chem(i,k,j,lptr_smpa_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase))) ! ! ! bbsoa_a04(i,k,j)= (chem(i,k,j,lptr_smpbb_aer(4,itype,iphase))) ! ! ! biog_a04(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase))) ! ! asmpsoa_a04(i,k,j)= chem(i,k,j,lptr_smpa_aer(4,itype,iphase)) ! ! ! ! totoa_a04(i,k,j)= ( & ! + chem(i,k,j,lptr_smpa_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_smpbb_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) & ! + chem(i,k,j,lptr_oc_aer(4,itype,iphase))) hoa_a04(i,k,j)= hoa_a04(i,k,j) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase)) soa_a04(i,k,j)= soa_a04(i,k,j) & + chem(i,k,j,lptr_smpa_aer(4,itype,iphase)) & + chem(i,k,j,lptr_smpbb_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) bbsoa_a04(i,k,j)= bbsoa_a04(i,k,j) & + chem(i,k,j,lptr_smpbb_aer(4,itype,iphase)) biog_a04(i,k,j)= biog_a04(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) asmpsoa_a04(i,k,j)= asmpsoa_a04(i,k,j) & + chem(i,k,j,lptr_smpa_aer(4,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase endif !Factor of 1.4 used below to convert OC to OA totoa_a01(its:imax,kts:kmax,jts:jmax) =totoa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) totoa_a02(its:imax,kts:kmax,jts:jmax) =totoa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) totoa_a03(its:imax,kts:kmax,jts:jmax) =totoa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) totoa_a04(its:imax,kts:kmax,jts:jmax) =totoa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) if( aero_diag_opt > 0 ) then hoa_a01(its:imax,kts:kmax,jts:jmax) =hoa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a01(its:imax,kts:kmax,jts:jmax) =soa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bboa_a01(its:imax,kts:kmax,jts:jmax) =bboa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bbsoa_a01(its:imax,kts:kmax,jts:jmax) =bbsoa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a01(its:imax,kts:kmax,jts:jmax) =biog_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) asmpsoa_a01(its:imax,kts:kmax,jts:jmax) =asmpsoa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) hoa_a02(its:imax,kts:kmax,jts:jmax) =hoa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a02(its:imax,kts:kmax,jts:jmax) =soa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bboa_a02(its:imax,kts:kmax,jts:jmax) =bboa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bbsoa_a02(its:imax,kts:kmax,jts:jmax) =bbsoa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a02(its:imax,kts:kmax,jts:jmax) =biog_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) asmpsoa_a02(its:imax,kts:kmax,jts:jmax) =asmpsoa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) hoa_a03(its:imax,kts:kmax,jts:jmax) =hoa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a03(its:imax,kts:kmax,jts:jmax) =soa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bboa_a03(its:imax,kts:kmax,jts:jmax) =bboa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bbsoa_a03(its:imax,kts:kmax,jts:jmax) =bbsoa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a03(its:imax,kts:kmax,jts:jmax) =biog_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) asmpsoa_a03(its:imax,kts:kmax,jts:jmax) =asmpsoa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) hoa_a04(its:imax,kts:kmax,jts:jmax) =hoa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a04(its:imax,kts:kmax,jts:jmax) =soa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bboa_a04(its:imax,kts:kmax,jts:jmax) =bboa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) bbsoa_a04(its:imax,kts:kmax,jts:jmax) =bbsoa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a04(its:imax,kts:kmax,jts:jmax) =biog_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) asmpsoa_a04(its:imax,kts:kmax,jts:jmax) =asmpsoa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) endif end subroutine sum_vbs0 !----------------------------------------------------------------------- subroutine sum_vbs2 ( aero_diag_opt, & alt, chem, & hoa_a01,hoa_a02,hoa_a03,hoa_a04, & hoa_a05,hoa_a06,hoa_a07,hoa_a08, & !BSINGH(12/04/2013): Added 4 more bins(5 to 8) for all apecies bboa_a01,bboa_a02,bboa_a03,bboa_a04, & bboa_a05,bboa_a06,bboa_a07,bboa_a08, & soa_a01,soa_a02,soa_a03,soa_a04, & soa_a05,soa_a06,soa_a07,soa_a08, & bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04, & bbsoa_a05,bbsoa_a06,bbsoa_a07,bbsoa_a08, & hsoa_a01,hsoa_a02,hsoa_a03,hsoa_a04, & hsoa_a05,hsoa_a06,hsoa_a07,hsoa_a08, & biog_a01,biog_a02,biog_a03,biog_a04, & biog_a05,biog_a06,biog_a07,biog_a08, & arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04, & arosoa_a05,arosoa_a06,arosoa_a07,arosoa_a08, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & totoa_a05,totoa_a06,totoa_a07,totoa_a08, & hsoa_c,hsoa_o,bbsoa_c,bbsoa_o, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: aero_diag_opt INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: hoa_a01,hoa_a02,hoa_a03,hoa_a04, & hoa_a05,hoa_a06,hoa_a07,hoa_a08, & bboa_a01,bboa_a02,bboa_a03,bboa_a04, & bboa_a05,bboa_a06,bboa_a07,bboa_a08, & soa_a01,soa_a02,soa_a03,soa_a04, & soa_a05,soa_a06,soa_a07,soa_a08, & bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04, & bbsoa_a05,bbsoa_a06,bbsoa_a07,bbsoa_a08, & biog_a01,biog_a02,biog_a03,biog_a04, & biog_a05,biog_a06,biog_a07,biog_a08, & hsoa_a01,hsoa_a02,hsoa_a03,hsoa_a04, & hsoa_a05,hsoa_a06,hsoa_a07,hsoa_a08, & arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04, & arosoa_a05,arosoa_a06,arosoa_a07,arosoa_a08, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & totoa_a05,totoa_a06,totoa_a07,totoa_a08, & hsoa_c,hsoa_o,bbsoa_c,bbsoa_o, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4 INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte totoa_a01(its:imax,kts:kmax,jts:jmax) = 0. totoa_a02(its:imax,kts:kmax,jts:jmax) = 0. totoa_a03(its:imax,kts:kmax,jts:jmax) = 0. totoa_a04(its:imax,kts:kmax,jts:jmax) = 0. totoa_a05(its:imax,kts:kmax,jts:jmax) = 0. totoa_a06(its:imax,kts:kmax,jts:jmax) = 0. totoa_a07(its:imax,kts:kmax,jts:jmax) = 0. totoa_a08(its:imax,kts:kmax,jts:jmax) = 0. if( aero_diag_opt > 0 ) then hoa_a01(its:imax,kts:kmax,jts:jmax) = 0. soa_a01(its:imax,kts:kmax,jts:jmax) = 0. bboa_a01(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a01(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a01(its:imax,kts:kmax,jts:jmax) = 0. biog_a01(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a01(its:imax,kts:kmax,jts:jmax) = 0. hoa_a02(its:imax,kts:kmax,jts:jmax) = 0. soa_a02(its:imax,kts:kmax,jts:jmax) = 0. bboa_a02(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a02(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a02(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a02(its:imax,kts:kmax,jts:jmax) = 0. biog_a02(its:imax,kts:kmax,jts:jmax) = 0. hoa_a03(its:imax,kts:kmax,jts:jmax) = 0. soa_a03(its:imax,kts:kmax,jts:jmax) = 0. bboa_a03(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a03(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a03(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a03(its:imax,kts:kmax,jts:jmax) = 0. biog_a03(its:imax,kts:kmax,jts:jmax) = 0. hoa_a04(its:imax,kts:kmax,jts:jmax) = 0. soa_a04(its:imax,kts:kmax,jts:jmax) = 0. bboa_a04(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a04(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a04(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a04(its:imax,kts:kmax,jts:jmax) = 0. biog_a04(its:imax,kts:kmax,jts:jmax) = 0. hoa_a05(its:imax,kts:kmax,jts:jmax) = 0. soa_a05(its:imax,kts:kmax,jts:jmax) = 0. bboa_a05(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a05(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a05(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a05(its:imax,kts:kmax,jts:jmax) = 0. biog_a05(its:imax,kts:kmax,jts:jmax) = 0. hoa_a06(its:imax,kts:kmax,jts:jmax) = 0. soa_a06(its:imax,kts:kmax,jts:jmax) = 0. bboa_a06(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a06(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a06(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a06(its:imax,kts:kmax,jts:jmax) = 0. biog_a06(its:imax,kts:kmax,jts:jmax) = 0. hoa_a07(its:imax,kts:kmax,jts:jmax) = 0. soa_a07(its:imax,kts:kmax,jts:jmax) = 0. bboa_a07(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a07(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a07(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a07(its:imax,kts:kmax,jts:jmax) = 0. biog_a07(its:imax,kts:kmax,jts:jmax) = 0. hoa_a08(its:imax,kts:kmax,jts:jmax) = 0. soa_a08(its:imax,kts:kmax,jts:jmax) = 0. bboa_a08(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a08(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a08(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a08(its:imax,kts:kmax,jts:jmax) = 0. biog_a08(its:imax,kts:kmax,jts:jmax) = 0. ! Species to calculate O:C ratios hsoa_c(its:imax,kts:kmax,jts:jmax) = 0. hsoa_o(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_c(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_o(its:imax,kts:kmax,jts:jmax) = 0. biog_v1(its:imax,kts:kmax,jts:jmax) = 0. biog_v2(its:imax,kts:kmax,jts:jmax) = 0. biog_v3(its:imax,kts:kmax,jts:jmax) = 0. biog_v4(its:imax,kts:kmax,jts:jmax) = 0. ant_v1(its:imax,kts:kmax,jts:jmax) = 0. ant_v2(its:imax,kts:kmax,jts:jmax) = 0. ant_v3(its:imax,kts:kmax,jts:jmax) = 0. ant_v4(its:imax,kts:kmax,jts:jmax) = 0. endif ! NOTE - summation seems to be wrong for code below, but not my code so did not fix... iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax totoa_a01(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(1,itype,iphase))) totoa_a02(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(2,itype,iphase))) totoa_a03(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(3,itype,iphase))) totoa_a04(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase))) totoa_a05(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(5,itype,iphase))) totoa_a06(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(6,itype,iphase))) totoa_a07(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(7,itype,iphase))) totoa_a08(i,k,j) = (chem(i,k,j,lptr_pcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(8,itype,iphase))) enddo enddo enddo enddo ! type if( aero_diag_opt > 0 ) then !do iphase=1,nphase_aer!BSINGH - Commented out as we need to add only phase 1 (interstitial) values species iphase = 1 do itype=1,ntype_aer do n = 1, nsize_aer(itype) !The 4th bin is 2.5-10um and outside the AMS measurements do j=jts,jmax do k=kts,kmax do i=its,imax hsoa_c(i,k,j)=hsoa_c(i,k,j) & + (chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)))*180/211 hsoa_o(i,k,j)= hsoa_o(i,k,j) & + (chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase))) bbsoa_c(i,k,j)= bbsoa_c(i,k,j) & + (chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)))*180/211 bbsoa_o(i,k,j)=bbsoa_o(i,k,j) & +(chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase))) biog_v1(i,k,j)= biog_v1(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) ant_v1(i,k,j)= ant_v1(i,k,j) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) enddo enddo enddo enddo enddo ! NOTE - summation also wrong for code below, but not my code so did not fix... !do iphase=1,nphase_aer !BSINGH - Commented out as we need to add only phase 1 (interstitial) values species iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a01(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(1,itype,iphase))) bboa_a01(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(1,itype,iphase))) soa_a01(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase))) arosoa_a01(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) bbsoa_a01(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase))) hsoa_a01(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase))) biog_a01(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=1,nphase_aer!BSINGH - Commented out as we need to add only phase 1 (interstitial) values species iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a02(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(2,itype,iphase))) bboa_a02(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(2,itype,iphase))) soa_a02(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase))) arosoa_a02(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) bbsoa_a02(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase))) hsoa_a02(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase))) biog_a02(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=1,nphase_aer!BSINGH - Commented out as we need to add only phase 1 (interstitial) values species iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a03(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(3,itype,iphase))) bboa_a03(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(3,itype,iphase))) soa_a03(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase))) arosoa_a03(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) bbsoa_a03(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase))) hsoa_a03(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase))) biog_a03(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=1,nphase_aer!BSINGH - Commented out as we need to add only phase 1 (interstitial) values species iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a04(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(4,itype,iphase))) bboa_a04(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(4,itype,iphase))) soa_a04(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase))) arosoa_a04(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) bbsoa_a04(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase))) hsoa_a04(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase))) biog_a04(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase))) enddo enddo enddo enddo ! type iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a05(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(5,itype,iphase))) bboa_a05(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(5,itype,iphase))) soa_a05(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(5,itype,iphase))) arosoa_a05(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(5,itype,iphase)) bbsoa_a05(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(5,itype,iphase))) hsoa_a05(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(5,itype,iphase))) biog_a05(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(5,itype,iphase))) enddo enddo enddo enddo ! type iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a06(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(6,itype,iphase))) bboa_a06(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(6,itype,iphase))) soa_a06(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(6,itype,iphase))) arosoa_a06(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(6,itype,iphase)) bbsoa_a06(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(6,itype,iphase))) hsoa_a06(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(6,itype,iphase))) biog_a06(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(6,itype,iphase))) enddo enddo enddo enddo ! type iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a07(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(7,itype,iphase))) bboa_a07(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(7,itype,iphase))) soa_a07(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(7,itype,iphase))) arosoa_a07(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(7,itype,iphase)) bbsoa_a07(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(7,itype,iphase))) hsoa_a07(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(7,itype,iphase))) biog_a07(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(7,itype,iphase))) enddo enddo enddo enddo ! type iphase = 1 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a08(i,k,j)= (chem(i,k,j,lptr_pcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(8,itype,iphase))) bboa_a08(i,k,j)= (chem(i,k,j,lptr_pcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(8,itype,iphase))) soa_a08(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(8,itype,iphase))) arosoa_a08(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(8,itype,iphase)) bbsoa_a08(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(8,itype,iphase))) hsoa_a08(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(8,itype,iphase))) biog_a08(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(8,itype,iphase))) enddo enddo enddo enddo ! type endif end subroutine sum_vbs2 !BSINGH(12/12/2013): Added following function for SAPRC 8 bin !---------------------------------------------------------------------- subroutine sum_aq_vbs2 ( & alt, chem, & hoa_cw01,hoa_cw02,hoa_cw03,hoa_cw04,hoa_cw05,hoa_cw06,hoa_cw07,hoa_cw08, & bboa_cw01,bboa_cw02,bboa_cw03,bboa_cw04,bboa_cw05,bboa_cw06,bboa_cw07,bboa_cw08, & soa_cw01,soa_cw02,soa_cw03,soa_cw04,soa_cw05,soa_cw06,soa_cw07,soa_cw08, & bbsoa_cw01,bbsoa_cw02,bbsoa_cw03,bbsoa_cw04,bbsoa_cw05,bbsoa_cw06,bbsoa_cw07,bbsoa_cw08, & hsoa_cw01,hsoa_cw02,hsoa_cw03,hsoa_cw04,hsoa_cw05,hsoa_cw06,hsoa_cw07,hsoa_cw08, & biog_cw01,biog_cw02,biog_cw03,biog_cw04,biog_cw05,biog_cw06,biog_cw07,biog_cw08, & arosoa_cw01,arosoa_cw02,arosoa_cw03,arosoa_cw04,arosoa_cw05,arosoa_cw06,arosoa_cw07,arosoa_cw08, & totoa_cw01,totoa_cw02,totoa_cw03,totoa_cw04,totoa_cw05,totoa_cw06,totoa_cw07,totoa_cw08, & hsoa_cw_c,hsoa_cw_o,bbsoa_cw_c,bbsoa_cw_o, & biog_cw_v1, & ant_cw_v1, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: hoa_cw01,hoa_cw02,hoa_cw03,hoa_cw04, & hoa_cw05,hoa_cw06,hoa_cw07,hoa_cw08, & bboa_cw01,bboa_cw02,bboa_cw03,bboa_cw04, & bboa_cw05,bboa_cw06,bboa_cw07,bboa_cw08, & soa_cw01,soa_cw02,soa_cw03,soa_cw04, & soa_cw05,soa_cw06,soa_cw07,soa_cw08, & bbsoa_cw01,bbsoa_cw02,bbsoa_cw03,bbsoa_cw04, & bbsoa_cw05,bbsoa_cw06,bbsoa_cw07,bbsoa_cw08, & biog_cw01,biog_cw02,biog_cw03,biog_cw04, & biog_cw05,biog_cw06,biog_cw07,biog_cw08, & hsoa_cw01,hsoa_cw02,hsoa_cw03,hsoa_cw04, & hsoa_cw05,hsoa_cw06,hsoa_cw07,hsoa_cw08, & arosoa_cw01,arosoa_cw02,arosoa_cw03,arosoa_cw04, & arosoa_cw05,arosoa_cw06,arosoa_cw07,arosoa_cw08, & totoa_cw01,totoa_cw02,totoa_cw03,totoa_cw04, & totoa_cw05,totoa_cw06,totoa_cw07,totoa_cw08, & hsoa_cw_c,hsoa_cw_o,bbsoa_cw_c,bbsoa_cw_o, & biog_cw_v1, & ant_cw_v1 INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte hoa_cw01(its:imax,kts:kmax,jts:jmax) = 0. soa_cw01(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw01(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw01(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw01(its:imax,kts:kmax,jts:jmax) = 0. biog_cw01(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw01(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw01(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw02(its:imax,kts:kmax,jts:jmax) = 0. soa_cw02(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw02(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw02(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw02(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw02(its:imax,kts:kmax,jts:jmax) = 0. biog_cw02(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw02(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw03(its:imax,kts:kmax,jts:jmax) = 0. soa_cw03(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw03(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw03(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw03(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw03(its:imax,kts:kmax,jts:jmax) = 0. biog_cw03(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw03(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw04(its:imax,kts:kmax,jts:jmax) = 0. soa_cw04(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw04(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw04(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw04(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw04(its:imax,kts:kmax,jts:jmax) = 0. biog_cw04(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw04(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw05(its:imax,kts:kmax,jts:jmax) = 0. soa_cw05(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw05(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw05(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw05(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw05(its:imax,kts:kmax,jts:jmax) = 0. biog_cw05(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw05(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw06(its:imax,kts:kmax,jts:jmax) = 0. soa_cw06(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw06(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw06(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw06(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw06(its:imax,kts:kmax,jts:jmax) = 0. biog_cw06(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw06(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw07(its:imax,kts:kmax,jts:jmax) = 0. soa_cw07(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw07(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw07(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw07(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw07(its:imax,kts:kmax,jts:jmax) = 0. biog_cw07(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw07(its:imax,kts:kmax,jts:jmax) = 0. hoa_cw08(its:imax,kts:kmax,jts:jmax) = 0. soa_cw08(its:imax,kts:kmax,jts:jmax) = 0. bboa_cw08(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw08(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw08(its:imax,kts:kmax,jts:jmax) = 0. arosoa_cw08(its:imax,kts:kmax,jts:jmax) = 0. biog_cw08(its:imax,kts:kmax,jts:jmax) = 0. totoa_cw08(its:imax,kts:kmax,jts:jmax) = 0. ! Species to calculate O:C ratios hsoa_cw_c(its:imax,kts:kmax,jts:jmax) = 0. hsoa_cw_o(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw_c(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_cw_o(its:imax,kts:kmax,jts:jmax) = 0. biog_cw_v1(its:imax,kts:kmax,jts:jmax) = 0. ant_cw_v1(its:imax,kts:kmax,jts:jmax) = 0. !do iphase=2,2 !set nphase_aer=2 for cloud-borne aerosols !BSINGH - Commented out iphase = 2 do itype=1,ntype_aer do n = 1, nsize_aer(itype) !The 4th bin is 2.5-10um and outside the AMS measurements do j=jts,jmax do k=kts,kmax do i=its,imax hsoa_cw_c(i,k,j)=hsoa_cw_c(i,k,j) & + (chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)))*180/211 hsoa_cw_o(i,k,j)= hsoa_cw_o(i,k,j) & + (chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase))) bbsoa_cw_c(i,k,j)= bbsoa_cw_c(i,k,j) & + (chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)))*180/211 bbsoa_cw_o(i,k,j)=bbsoa_cw_o(i,k,j) & +(chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase))) biog_cw_v1(i,k,j)= biog_cw_v1(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) ant_cw_v1(i,k,j)= ant_cw_v1(i,k,j) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) enddo enddo enddo enddo enddo !do iphase=2,2 !set nphase_aer=2 for cloud-borne aerosols !BSINGH - Commented out iphase = 2 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw01(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(1,itype,iphase)) bboa_cw01(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(1,itype,iphase)) soa_cw01(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) arosoa_cw01(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) bbsoa_cw01(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase))) hsoa_cw01(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase))) biog_cw01(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) totoa_cw01(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(1,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=2,2 !set nphase_aer=2 for cloud-borne aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw02(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(2,itype,iphase)) bboa_cw02(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(2,itype,iphase)) soa_cw02(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) arosoa_cw02(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) bbsoa_cw02(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase))) hsoa_cw02(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase))) biog_cw02(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) totoa_cw02(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(2,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=2,2 !set nphase_aer=2 for cloud-borne aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw03(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(3,itype,iphase)) bboa_cw03(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(3,itype,iphase)) soa_cw03(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) arosoa_cw03(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) bbsoa_cw03(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase))) hsoa_cw03(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase))) biog_cw03(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) totoa_cw03(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(3,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=2,2 !set nphase_aer=2 for cloud-borne aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw04(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(4,itype,iphase)) bboa_cw04(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(4,itype,iphase)) soa_cw04(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) arosoa_cw04(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) bbsoa_cw04(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase))) hsoa_cw04(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase))) biog_cw04(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) totoa_cw04(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase))) enddo enddo enddo enddo ! type !do iphase=2,2 !set nphase_aer=1 for interstitial aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer if(nsize_aer(itype).ge.5) then do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw05(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(5,itype,iphase)) bboa_cw05(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(5,itype,iphase)) soa_cw05(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(5,itype,iphase)) arosoa_cw05(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(5,itype,iphase)) bbsoa_cw05(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(5,itype,iphase))) hsoa_cw05(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(5,itype,iphase))) biog_cw05(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(5,itype,iphase)) totoa_cw05(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(5,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(5,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(5,itype,iphase))) enddo enddo enddo endif ! check for number of size bins enddo ! type !do iphase=2,2 !set nphase_aer=1 for interstitial aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer if(nsize_aer(itype).ge.6) then do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw06(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(6,itype,iphase)) bboa_cw06(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(6,itype,iphase)) soa_cw06(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(6,itype,iphase)) arosoa_cw06(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(6,itype,iphase)) bbsoa_cw06(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(6,itype,iphase))) hsoa_cw06(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(6,itype,iphase))) biog_cw06(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(6,itype,iphase)) totoa_cw06(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(6,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(6,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(6,itype,iphase))) enddo enddo enddo endif ! size bins enddo ! type !do iphase=2,2 !set nphase_aer=1 for interstitial aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer if(nsize_aer(itype).ge.7) then do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw07(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(7,itype,iphase)) bboa_cw07(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(7,itype,iphase)) soa_cw07(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(7,itype,iphase)) arosoa_cw07(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(7,itype,iphase)) bbsoa_cw07(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(7,itype,iphase))) hsoa_cw07(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(7,itype,iphase))) biog_cw07(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(7,itype,iphase)) totoa_cw07(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(7,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(7,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(7,itype,iphase))) enddo enddo enddo endif ! size bins enddo ! type !do iphase=2,2 !set nphase_aer=1 for interstitial aerosols!BSINGH - Commented out iphase = 2 do itype=1,ntype_aer if(nsize_aer(itype).ge.8) then do j=jts,jmax do k=kts,kmax do i=its,imax hoa_cw08(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(8,itype,iphase)) bboa_cw08(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(8,itype,iphase)) soa_cw08(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(8,itype,iphase)) arosoa_cw08(i,k,j)= chem(i,k,j,lptr_ant1_c_aer(8,itype,iphase)) bbsoa_cw08(i,k,j)= (chem(i,k,j,lptr_opcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(8,itype,iphase))) hsoa_cw08(i,k,j)= ( chem(i,k,j,lptr_opcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(8,itype,iphase))) biog_cw08(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(8,itype,iphase)) totoa_cw08(i,k,j)= ( chem(i,k,j,lptr_pcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(8,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(8,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(8,itype,iphase))) enddo enddo enddo endif ! size bins enddo ! type end subroutine sum_aq_vbs2 !BSINGH -ENDS !---------------------------------------------------------------------- !----------------------------------------------------------------------- subroutine sum_vbs4 ( aero_diag_opt, & alt, chem, & hoa_a01,hoa_a02,hoa_a03,hoa_a04, & soa_a01,soa_a02,soa_a03,soa_a04, & biog_a01,biog_a02,biog_a03,biog_a04, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: aero_diag_opt INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: hoa_a01,hoa_a02,hoa_a03,hoa_a04, & soa_a01,soa_a02,soa_a03,soa_a04, & biog_a01,biog_a02,biog_a03,biog_a04, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4 INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte totoa_a01(its:imax,kts:kmax,jts:jmax) = 0. totoa_a02(its:imax,kts:kmax,jts:jmax) = 0. totoa_a03(its:imax,kts:kmax,jts:jmax) = 0. totoa_a04(its:imax,kts:kmax,jts:jmax) = 0. if( aero_diag_opt > 0 ) then hoa_a01(its:imax,kts:kmax,jts:jmax) = 0. soa_a01(its:imax,kts:kmax,jts:jmax) = 0. biog_a01(its:imax,kts:kmax,jts:jmax) = 0. hoa_a02(its:imax,kts:kmax,jts:jmax) = 0. soa_a02(its:imax,kts:kmax,jts:jmax) = 0. biog_a02(its:imax,kts:kmax,jts:jmax) = 0. hoa_a03(its:imax,kts:kmax,jts:jmax) = 0. soa_a03(its:imax,kts:kmax,jts:jmax) = 0. biog_a03(its:imax,kts:kmax,jts:jmax) = 0. hoa_a04(its:imax,kts:kmax,jts:jmax) = 0. soa_a04(its:imax,kts:kmax,jts:jmax) = 0. biog_a04(its:imax,kts:kmax,jts:jmax) = 0. ! Species to calculate O:C ratios biog_v1(its:imax,kts:kmax,jts:jmax) = 0. biog_v2(its:imax,kts:kmax,jts:jmax) = 0. biog_v3(its:imax,kts:kmax,jts:jmax) = 0. biog_v4(its:imax,kts:kmax,jts:jmax) = 0. ant_v1(its:imax,kts:kmax,jts:jmax) = 0. ant_v2(its:imax,kts:kmax,jts:jmax) = 0. ant_v3(its:imax,kts:kmax,jts:jmax) = 0. ant_v4(its:imax,kts:kmax,jts:jmax) = 0. endif do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax totoa_a01(i,k,j)= totoa_a01(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(1,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(1,itype,iphase)) totoa_a02(i,k,j)= totoa_a02(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(2,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(2,itype,iphase)) totoa_a03(i,k,j)= totoa_a03(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(3,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(3,itype,iphase)) totoa_a04(i,k,j)= totoa_a04(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(4,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase if( aero_diag_opt > 0 ) then do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) !The 4th bin is 2.5-10um and outside the AMS measurements do j=jts,jmax do k=kts,kmax do i=its,imax biog_v1(i,k,j)= biog_v1(i,k,j) & + chem(i,k,j,lptr_bsoaX_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(n,itype,iphase)) ant_v1(i,k,j)= ant_v1(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(n,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(n,itype,iphase)) enddo enddo enddo enddo enddo enddo biog_v1(its:imax,kts:kmax,jts:jmax) = biog_v1(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) ant_v1(its:imax,kts:kmax,jts:jmax) = ant_v1(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a01(i,k,j)= hoa_a01(i,k,j) & +chem(i,k,j,lptr_oc_aer(1,itype,iphase)) soa_a01(i,k,j)= soa_a01(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(1,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(1,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(1,itype,iphase)) biog_a01(i,k,j)= biog_a01(i,k,j) & + chem(i,k,j,lptr_bsoa1_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(1,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(1,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a02(i,k,j)= hoa_a02(i,k,j) & +chem(i,k,j,lptr_oc_aer(2,itype,iphase)) soa_a02(i,k,j)= soa_a02(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(2,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(2,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(2,itype,iphase)) biog_a02(i,k,j)= biog_a02(i,k,j) & + chem(i,k,j,lptr_bsoa1_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(2,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(2,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a03(i,k,j)= hoa_a03(i,k,j) & + chem(i,k,j,lptr_oc_aer(3,itype,iphase)) soa_a03(i,k,j)= soa_a03(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(3,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(3,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(3,itype,iphase)) biog_a03(i,k,j)= biog_a03(i,k,j) & + chem(i,k,j,lptr_bsoa1_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(3,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(3,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a04(i,k,j)= hoa_a04(i,k,j) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase)) soa_a04(i,k,j)= soa_a04(i,k,j) & + chem(i,k,j,lptr_asoaX_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa3_aer(4,itype,iphase)) & + chem(i,k,j,lptr_asoa4_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoaX_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_r2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_sfc_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_oh_aer(4,itype,iphase)) & + chem(i,k,j,lptr_glysoa_nh4_aer(4,itype,iphase)) biog_a04(i,k,j)= biog_a04(i,k,j) & + chem(i,k,j,lptr_bsoa1_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa2_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa3_aer(4,itype,iphase)) & + chem(i,k,j,lptr_bsoa4_aer(4,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase endif !Factor of 1.4 used below to convert OC to OA totoa_a01(its:imax,kts:kmax,jts:jmax) =totoa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) totoa_a02(its:imax,kts:kmax,jts:jmax) =totoa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) totoa_a03(its:imax,kts:kmax,jts:jmax) =totoa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) totoa_a04(its:imax,kts:kmax,jts:jmax) =totoa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) if( aero_diag_opt > 0 ) then hoa_a01(its:imax,kts:kmax,jts:jmax) =hoa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a01(its:imax,kts:kmax,jts:jmax) =soa_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a01(its:imax,kts:kmax,jts:jmax) =biog_a01(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) hoa_a02(its:imax,kts:kmax,jts:jmax) =hoa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a02(its:imax,kts:kmax,jts:jmax) =soa_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a02(its:imax,kts:kmax,jts:jmax) =biog_a02(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) hoa_a03(its:imax,kts:kmax,jts:jmax) =hoa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a03(its:imax,kts:kmax,jts:jmax) =soa_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a03(its:imax,kts:kmax,jts:jmax) =biog_a03(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) hoa_a04(its:imax,kts:kmax,jts:jmax) =hoa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) soa_a04(its:imax,kts:kmax,jts:jmax) =soa_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) biog_a04(its:imax,kts:kmax,jts:jmax) =biog_a04(its:imax,kts:kmax,jts:jmax) & / alt(its:imax,kts:kmax,jts:jmax) endif end subroutine sum_vbs4 !---------------------------------------------------------------------- subroutine sum_vbs9 ( & alt, chem, & hoa_a01,hoa_a02,hoa_a03,hoa_a04, & bboa_a01,bboa_a02,bboa_a03,bboa_a04, & soa_a01,soa_a02,soa_a03,soa_a04, & bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04, & hsoa_a01,hsoa_a02,hsoa_a03,hsoa_a04, & biog_a01,biog_a02,biog_a03,biog_a04, & arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & hsoa_c,hsoa_o,bbsoa_c,bbsoa_o, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4, & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte ) USE module_state_description, only: num_chem USE module_data_mosaic_asect IMPLICIT NONE INTEGER, INTENT(IN ) :: & ids,ide, jds,jde, kds,kde, & ims,ime, jms,jme, kms,kme, & its,ite, jts,jte, kts,kte REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(IN) :: alt REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_chem ), & INTENT(IN ) :: chem REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), & INTENT(OUT) :: hoa_a01,hoa_a02,hoa_a03,hoa_a04, & bboa_a01,bboa_a02,bboa_a03,bboa_a04, & soa_a01,soa_a02,soa_a03,soa_a04, & bbsoa_a01,bbsoa_a02,bbsoa_a03,bbsoa_a04, & biog_a01,biog_a02,biog_a03,biog_a04, & hsoa_a01,hsoa_a02,hsoa_a03,hsoa_a04, & arosoa_a01,arosoa_a02,arosoa_a03,arosoa_a04, & totoa_a01,totoa_a02,totoa_a03,totoa_a04, & hsoa_c,hsoa_o,bbsoa_c,bbsoa_o, & biog_v1,biog_v2,biog_v3,biog_v4, & ant_v1,ant_v2,ant_v3,ant_v4 INTEGER :: i,imax,j,jmax,k,kmax,n,itype,iphase imax = min(ite,ide-1) jmax = min(jte,jde-1) kmax = kte hoa_a01(its:imax,kts:kmax,jts:jmax) = 0. soa_a01(its:imax,kts:kmax,jts:jmax) = 0. bboa_a01(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a01(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a01(its:imax,kts:kmax,jts:jmax) = 0. biog_a01(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a01(its:imax,kts:kmax,jts:jmax) = 0. totoa_a01(its:imax,kts:kmax,jts:jmax) = 0. hoa_a02(its:imax,kts:kmax,jts:jmax) = 0. soa_a02(its:imax,kts:kmax,jts:jmax) = 0. bboa_a02(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a02(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a02(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a02(its:imax,kts:kmax,jts:jmax) = 0. biog_a02(its:imax,kts:kmax,jts:jmax) = 0. totoa_a02(its:imax,kts:kmax,jts:jmax) = 0. hoa_a03(its:imax,kts:kmax,jts:jmax) = 0. soa_a03(its:imax,kts:kmax,jts:jmax) = 0. bboa_a03(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a03(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a03(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a03(its:imax,kts:kmax,jts:jmax) = 0. biog_a03(its:imax,kts:kmax,jts:jmax) = 0. totoa_a03(its:imax,kts:kmax,jts:jmax) = 0. hoa_a04(its:imax,kts:kmax,jts:jmax) = 0. soa_a04(its:imax,kts:kmax,jts:jmax) = 0. bboa_a04(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_a04(its:imax,kts:kmax,jts:jmax) = 0. hsoa_a04(its:imax,kts:kmax,jts:jmax) = 0. arosoa_a04(its:imax,kts:kmax,jts:jmax) = 0. biog_a04(its:imax,kts:kmax,jts:jmax) = 0. totoa_a04(its:imax,kts:kmax,jts:jmax) = 0. ! Species to calculate O:C ratios hsoa_c(its:imax,kts:kmax,jts:jmax) = 0. hsoa_o(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_c(its:imax,kts:kmax,jts:jmax) = 0. bbsoa_o(its:imax,kts:kmax,jts:jmax) = 0. biog_v1(its:imax,kts:kmax,jts:jmax) = 0. biog_v2(its:imax,kts:kmax,jts:jmax) = 0. biog_v3(its:imax,kts:kmax,jts:jmax) = 0. biog_v4(its:imax,kts:kmax,jts:jmax) = 0. ant_v1(its:imax,kts:kmax,jts:jmax) = 0. ant_v2(its:imax,kts:kmax,jts:jmax) = 0. ant_v3(its:imax,kts:kmax,jts:jmax) = 0. ant_v4(its:imax,kts:kmax,jts:jmax) = 0. do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) !The 4th bin is 2.5-10um and outside the AMS measurements do j=jts,jmax do k=kts,kmax do i=its,imax hsoa_c(i,k,j)=hsoa_c(i,k,j) & + (chem(i,k,j,lptr_pcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(n,itype,iphase)))*180/211 hsoa_o(i,k,j)= hsoa_o(i,k,j) & + chem(i,k,j,lptr_pcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(n,itype,iphase)) bbsoa_c(i,k,j)= bbsoa_c(i,k,j) & + (chem(i,k,j,lptr_pcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(n,itype,iphase)))*180/211 bbsoa_o(i,k,j)=bbsoa_o(i,k,j) & + chem(i,k,j,lptr_pcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(n,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(n,itype,iphase)) biog_v1(i,k,j)= biog_v1(i,k,j) & + chem(i,k,j,lptr_biog1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(n,itype,iphase)) biog_v2(i,k,j)= biog_v2(i,k,j) & + chem(i,k,j,lptr_biog2_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(n,itype,iphase)) biog_v3(i,k,j)= biog_v3(i,k,j) & + chem(i,k,j,lptr_biog3_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(n,itype,iphase)) biog_v4(i,k,j)= biog_v4(i,k,j) & +chem(i,k,j,lptr_biog4_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(n,itype,iphase)) ant_v1(i,k,j)= ant_v1(i,k,j) & + chem(i,k,j,lptr_ant1_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(n,itype,iphase)) ant_v2(i,k,j)= ant_v2(i,k,j) & + chem(i,k,j,lptr_ant2_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(n,itype,iphase)) ant_v3(i,k,j)= ant_v3(i,k,j) & + chem(i,k,j,lptr_ant3_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(n,itype,iphase)) ant_v4(i,k,j)= ant_v4(i,k,j) & + chem(i,k,j,lptr_ant4_c_aer(n,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(n,itype,iphase)) enddo enddo enddo enddo enddo enddo ! NOTE - summation also wrong for code below, but not my code so did not fix... do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a01(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(1,itype,iphase)) bboa_a01(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(1,itype,iphase)) soa_a01(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(1,itype,iphase)) & + (chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(1,itype,iphase))) & + chem(i,k,j,lptr_ant1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(1,itype,iphase)) & + (chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(1,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(1,itype,iphase)) arosoa_a01(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(1,itype,iphase)) bbsoa_a01(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(1,itype,iphase)) hsoa_a01(i,k,j)= chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(1,itype,iphase)) biog_a01(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(1,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(1,itype,iphase)) totoa_a01(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(1,itype,iphase)) & + (chem(i,k,j,lptr_ant1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(1,itype,iphase))) & + chem(i,k,j,lptr_ant1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(1,itype,iphase)) & + (chem(i,k,j,lptr_biog1_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(1,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(1,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(1,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a02(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(2,itype,iphase)) bboa_a02(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(2,itype,iphase)) soa_a02(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(2,itype,iphase)) & + (chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(2,itype,iphase))) & + chem(i,k,j,lptr_ant1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(2,itype,iphase)) & + (chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(2,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(2,itype,iphase)) arosoa_a02(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(2,itype,iphase)) bbsoa_a02(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(2,itype,iphase)) hsoa_a02(i,k,j)= chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(2,itype,iphase)) biog_a02(i,k,j)= chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(2,itype,iphase)) totoa_a02(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(2,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(2,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a03(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(3,itype,iphase)) bboa_a03(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(3,itype,iphase)) soa_a03(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(3,itype,iphase)) arosoa_a03(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(3,itype,iphase)) bbsoa_a03(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(3,itype,iphase)) hsoa_a03(i,k,j)= chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(3,itype,iphase)) biog_a03(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(3,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(3,itype,iphase)) totoa_a03(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(3,itype,iphase)) & + (chem(i,k,j,lptr_ant1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(3,itype,iphase))) & + chem(i,k,j,lptr_ant1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(3,itype,iphase)) & + (chem(i,k,j,lptr_biog1_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(3,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(3,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(3,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase do iphase=1,nphase_aer do itype=1,ntype_aer do j=jts,jmax do k=kts,kmax do i=its,imax hoa_a04(i,k,j)= chem(i,k,j,lptr_pcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(4,itype,iphase)) bboa_a04(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(4,itype,iphase)) soa_a04(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(4,itype,iphase)) & + (chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(4,itype,iphase))) & + chem(i,k,j,lptr_ant1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(4,itype,iphase)) & + (chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(4,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(4,itype,iphase)) arosoa_a04(i,k,j)= chem(i,k,j,lptr_ant1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(4,itype,iphase)) bbsoa_a04(i,k,j)= chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(4,itype,iphase)) hsoa_a04(i,k,j)= chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(4,itype,iphase)) biog_a04(i,k,j)= (chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(4,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(4,itype,iphase)) totoa_a04(i,k,j)= chem(i,k,j,lptr_pcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_b_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg2_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg3_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg4_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg5_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg6_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_pcg7_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg1_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg2_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg3_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg4_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg5_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg6_f_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_opcg7_f_o_aer(4,itype,iphase)) & + (chem(i,k,j,lptr_ant1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant2_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant3_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant4_c_aer(4,itype,iphase))) & + chem(i,k,j,lptr_ant1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant2_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant3_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_ant4_o_aer(4,itype,iphase)) & + (chem(i,k,j,lptr_biog1_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog2_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog3_c_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog4_c_aer(4,itype,iphase))) & + chem(i,k,j,lptr_biog1_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog2_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog3_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_biog4_o_aer(4,itype,iphase)) & + chem(i,k,j,lptr_oc_aer(4,itype,iphase)) enddo enddo enddo enddo ! type enddo ! phase end subroutine sum_vbs9 !----------------------------------------------------------------------- end module module_mosaic_sumpm