module module_mosaic_sumpm 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 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 enddo enddo 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 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 enddo enddo 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 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)) 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)) enddo enddo enddo endif enddo enddo enddo 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 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)) 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)) enddo enddo enddo endif enddo enddo enddo 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. 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 enddo if( aero_diag_opt > 0 ) then do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) 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 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 enddo 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_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 enddo 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_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 enddo 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_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 enddo endif 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, & 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. 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 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 if( aero_diag_opt > 0 ) then iphase = 1 do itype=1,ntype_aer do n = 1, nsize_aer(itype) 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 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 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 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 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 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 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 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 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 endif end subroutine sum_vbs2 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. 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. iphase = 2 do itype=1,ntype_aer do n = 1, nsize_aer(itype) 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 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 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 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 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 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 enddo 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 enddo 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 enddo 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 enddo end subroutine sum_aq_vbs2 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. 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 enddo if( aero_diag_opt > 0 ) then do iphase=1,nphase_aer do itype=1,ntype_aer do n = 1, nsize_aer(itype) 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 enddo 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 enddo 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 enddo 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 enddo endif 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. 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) 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 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 enddo 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 enddo 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 enddo 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 enddo end subroutine sum_vbs9 end module module_mosaic_sumpm